Modelling of the interstellar hydrogen distribution in the heliosphere

The detailed knowledge of the distribution of neutral interstellar hydrogen in the interplanetary space is necessary for a reliable interpretation of optical and H⁺ pickup ions observations. In the paper, we review the status of the modelling efforts with the emphasis on recent improvements in that field. We discuss in particular the role of the nonstationary, solar cycle-related effects and the consequences of hydrogen filtration through the heliospheric interface region for its distribution in the inner Solar System. We demonstrate also that the use of the simple cold model, neglecting the thermal character of the hydrogen gas (T 8000 K), is generally incorrect for the whole region of the inner heliosphere (R < 5 AU) since it leads to a substantial underestimation of the local hydrogen density and thus influences the derivation of the H properties in the outer heliosphere/LISM. Referring to recent Ulysses measurements, we point out also the need to consider in the modelling the effects of the latitudinal asymmetry of the ionization rate.     

Relations between ISM inside and outside the heliosphere

Thanks to remarkable new tools, such as the Goddard High Resolution Spectrograph (GHRS) on board the HST and the EUVE spectrometer on the interstellar side, and Ulysses particle detectors on the heliospheric side, it is possible now to begin to compare abundances and physical properties of the interstellar matter outside the heliosphere (from absorption features in the stellar spectra), and inside the heliosphere (from in situ or remote detection of the interstellar neutrals or their derivatives, the pick-up ions or the Anomalous Cosmic Rays detected by the two Voyager spacecraft).Ground-based and UV spectra of nearby stars show that the Sun is located between two volumes of gas of different heliocentric velocities V and temperatures T (see also Linsky et al, this issue). One of these clouds has the same velocity (V= 25.6 km s^–1 from = 255 and =8) and temperature (6700 K) as the heliospheric helium of interstellar origin probed by Ulysses, and is certainly surrounding our star (and then the Local Interstellar Cloud or LIC). This Identification allows comparisons between interstellar constituents on both sides of the heliospheric interface.Ly-alpha background data (absorption cell and recent HST-GHRS spectra) suggest that the heliospheric neutral H velocity is smaller by 5–6 km s^–1 than the local cloud velocity, and therefore that H is decelerated at its entrance into the heliosphere, in agreement with interaction models between the heliosphere and the ISM which include the coupling with the plasma. This is in favor of a non negligible electron density (at least 0.05 cm³). There are other indications of a rather large ionization of the ambient ISM, such as the ionization equilibrium of interstellar magnesium and of sodium. However the resulting range for the plasma density is still broad.The heliospheric neutral hydrogen number density (0.08–0.16 cm^–3) is now less precisely determined than the helium density (0.013–0.017 cm^–3, see Gloeckler, Witte et al, Mobius, this issue). The comparison between the neutral hydrogen to neutral helium ratios in the ISM (recent EUVE findings) and in the heliosphere, suggests that 15 to 70% of H does not enter the heliosphere. The comparison between the interstellar oxygen relative abundance (with respect to H and He) in the ISM and the heliospheric abundance deduced from pick-up ions is also in favor of some filtration, and thus of a non-negligible ionization.For a significant ISM plasma density, one expects a Hydrogen wall to be present as an intermediate state of the interstellar H around the interface between inside and outside. Since 1993, the two UVS instruments on board Voyager 1 and 2 indeed reveal clearly the existence of an additional Ly-alpha emission, probably due to a combination of light from the compressed H wall, and from a galactic source. On the other hand, the decelerated and heated neutral hydrogen of this H wall has recently been detected in absorption in the spectra of nearby stars (see Linsky, this issue).     

Physical and chemical characteristics of the ISM inside and outside the heliosphere

This paper summarizes some of the discussions of working group 8–9 during the ISSI Conference on The Heliosphere in the Local Interstellar Medium. Because the subject of these working groups has become significantly broader during the last ten years, we have selected three topics for which recent observations have modified and improved our knowledge of the heliosphere and the surrounding interstellar medium. These topics are the number densities and ISM ionization states of hydrogen and helium, the newly discovered hot gas from the H wall seen in absorption, and the comparison between ISM and heliospheric minor element abundances. Papers from this volume in which more details on these topics can be found are quoted throughout the report.     

Possible shock wave in the local interstellar plasma, very close to the heliosphere

We show, using the HST — GHRS data on velocity and temperature in the nearby interstellar medium, that the observed 3 – 4 km s^–1 relative velocity between the Local Interstellar Cloud (LIC) and the so-called G-cloud located in the Galactic Center hemisphere can be quite naturally explained assuming that the two clouds do interact with each other. In the proposed interpretation the two media are separated by a (quasiperpendicular) MHD shock front propagating from the LIC into the G-cloud. The LIC plasma is then nothing else but the shocked (compression 1.3 – 1.4) gas of the G-cloud. A 1-D single-fluid solution of the Rankine — Hugoniot equations can fit the most probable observed values of the relative velocity (3.75 km/s), LIC (6700 K) and G-cloud (5400 K) kinetic temperatures, if the plasma-beta of the LIC plasma is in the range 1.3 – 1.5 (Table 1). This corresponds to a super — fast magnetosonic motion of the heliosphere through the LIC, independently of LIC density. The LIC magnetic field strength is 1.9 (3.1) G for the LIC electron density n_e = 0.04 (0.10) cm^–3. In this case the shock is less than 30 000 AU away and moves at about 10 km s^–1 relative to the LIC plasma. The Sun is chasing the shock and should catch up with it in about 10⁴ years. If the heliospheric VLP emissions cutoff at 1.8 kHz is indicative of n_e (LIC) = 0.04 cm^–3 (Gurnett et al., 1993), the (pure plasma) bowshock ahead of the heliopause could be the source of quasi-continuous heliospheric 2-kHz emission band. We believe that with the expected increase in the performance of modern spectroscopic instrumentation the proposed method of magnetic field evaluation may in the future find wider application in the studies of the interstellar medium.     

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.     

Interstellar gas in the heliosphere and the solar wind anisotropies

A critical review of the interstellar hydrogen in the heliosphere will be presented. Recent Sun-interstellar matter interaction model improvements, a non-stationary flow and a flexible latitude dependence, will be discussed. We also consider the influence of heliospheric interface on neutral flow and the remaining refinements, which could help to better interpret the results of the SWAN experiment on board SOHO.     

Relative ionizations in the nearest interstellar gas

We compare CLOUDY predictions for the equilibrium ionization in the interstellar cloud surrounding the solar system with pick-up ion data. The incident radiation field includes contributions from hot stars, the emission from the conductive cloud boundary and the diffuse FUV back-ground. To within the observational uncertainties, CLOUDY predictions for the ratios n(He)/n(O), n(N)/n(O), n(Ne)/n(O), and n(He)/n(Ne) are consistent with pick-up ion data, provided that O and N are filtered by 50% in the heliopause region and the outer heliosphere as predicted by others. Thus, the steady-state ionization model and assumed radiation field appear approximately valid. However, the youth and low intervening column density towards the Vela pulsar leave open the possibility that the parent supernova explosion 10,500 years ago, and 200 pc distant, may also have affected LISM ionization, although the mechanism is uncertain. Support for this last possibility is provided by the apparent signature of the Vela explosion in the terrestrial geological record.Abbreviations ISM Interstellar Medium - FUV Far Ultraviolet - EUV Extreme Ultraviolet - SNR SN remnant - SXRB SXR Background - LISM Local Interstellar Matter     

Modelling the heliosphere

An overview of our present efforts at the Bartol Research Institute in modelling the largescale interaction of the solar wind with the local interstellar medium is presented. Particular stress is placed on the self-consistent inclusion of neutral hydrogen in the models and both 2D and 3D structure is discussed. Observational implications are noted.     

Pickup protons in the heliosphere

We calculate the conditions of pickup protons inside the termination shock. Outside 50 AU the partial pressure of pickup protons is greater than the magnetic pressure by a factor of > 10, and greater than the partial pressure of solar wind protons by a factor of > 100. Thus, pickup protons have a significant dynamical influence on the structures of the solar wind in the outer heliosphere.     

MHD processes in the outer heliosphere

Conclusion Much has been learned about the structure and dynamics of the outer heliosphere during the last decade as a result observations from the Voyager and Pioneer spacecraft. The large scale of the observations forces one to consider the heliosphere from a new perspective, to think of new dynamical processes, and to introduce new concepts. The early studies of isolated gas dynamic flows must be replaced by MHD dynamics of interacting flows and flow systems. The simple deterministic models that have been dominant in early studies of the solar wind are now seen to have limited applicability, and statistical approaches are being developed. New concepts that have been introduced, such as inverse cascades, filtering, entrainment, etc., must be further explored and clarified, to make them more precise and quantitative. MHD turbulence is probably very important in solar wind dynamics, but the subject is poorly developed from a theoretical point of view. The statistical analysis of solar wind parameters has scarcely begun, but it is clearly necessary for an understanding of complex, large-scale flows. The multitude of possible interactions among shocks and flows of various types needs to be explored systematically with observations, models and analytical theory. Voyagers 1 and 2 and Pioneers 10 and 11 are continuing to move through the outer heliosphere and gather data. The lengthy data reduction procedures require even more care in dealing with the low field strengths, densities and temperatures at large heliocentric distances, and the analysis of the complex flows and fields in the outer heliosphere becomes increasingly difficult. Thus one can expect continued growth of our knowledge of the heliosphere, but comprehensive understanding of the data will take some time. If this review stimulates the specialists in solar wind physics to think critically about the results presented and to remedy the deficiencies of current knowledge of the heliosphere, then it will have served its purpose. It is also hoped that this review will serve to encourage specialists in other fields to bring their talents to bear on heliospheric problems and to transfer results of heliospheric physics to their fields.     

Radio emissions from the outer heliosphere

For nearly fifteen years the Voyager 1 and 2 spacecraft have been detecting an unusual radio emission in the outer heliosphere in the frequency range from about 2 to 3 kHz, Two major events have been observed, the first in 1983–84 and the second in 1992–93. In both cases the onset of the radio emission occurred about 400 days after a period of intense solar activity, the first in mid-July 1982, and the second in May–June 1991. These two periods of solar activity produced the two deepest cosmic ray Forbush decreases ever observed. Forbush decreases are indicative of a system of strong shocks and associated disturbances propagating outward through the heliosphere. The radio emission is believed to have been produced when this system of shocks and disturbances interacted with one of the outer boundaries of the heliosphere, most likely in the vicinity of the the heliopause. The emission is believed to be generated by the shock-driven Langmuir-wave mode conversion mechanism, which produces radiation at the plasma frequency (f _p ) and at twice the plasma frequency (2f _p ). From the 400-day travel time and the known speed of the shocks, the distance to the interaction region can be computed, and is estimated to be in the range from about 110 to 160 AU.Abbreviations PWS Plasma Wave Subsystem - AU Astronomical Unit - DSN Deep Space Network - NASA National Aeronautics and Space Administration - GMIR Global Merged Interaction Region - MHD Magnetohydrodynamic - CME coronal mass ejection - f _p plasma frequency - R radial distance - AGC automatic gain control     

Shock interactions in the outer heliosphere

Observations of plasma and magnetic fields by Pioneer 10 and 11 and Voyager 1 and 2 reveal that MHD shocks are an important component of the large-scale solar wind structures in the outer heliosphere. This review discusses recent progress in simulation studies of the nonlinear evolution of the solar wind structures, and in particular concentrates on the theoretical development and applications of the shock interactions model. Various stream propagation models, which do not use the Rankine-Hugoniot relations to calculate the jump conditions at shock crossings, have been used to simulate the essential evolution process of isolated streams and the formation and propagation of corotating and transient shocks. They produce fairly good results in the region up to a few AU. In 1984, the shock interactions model was introduced to study the evolution of large-scale solar wind structures in the region outside 1 AU up to several tens of AU. The model uses the exact Rankine-Hugoniot relations to calculate the shock speed and shock strength at all shock crossings. So that the model can more accurately calculate the shock speeds and the accumulated irreversible shock heating of plasma at several tens of AU. The applications of the shock interactions model are presented in three groups. (a) The first group covers the basic interaction of a shock with the ambient solar wind, the formation and propagation of shock pairs, and the collision and merging of shocks. (b) The second group covers the use of the shock interactions model to simulate the nonlinear evolution of large-scale solar wind structures in the outer heliosphere. These simulation results can provide the detailed evolution process for large-scale solar wind structures in the vast region not directly observed. Two selected studies are reported. (c) Finally, the shock interactions model is applied to studying the heating of the solar wind in the outer heliosphere. The model calculations support shocks being chiefly responsible for the heating of the solar wind plasma in the outer heliosphere at least up to 30 AU.     

Energetic particle transport coefficients in the heliosphere

It is the purpose of this review to summarize and discuss recent research done in the field of particle propagation in the heliosphere. Several lines of approach have been followed to treat this problem. As a starting point the different forms of the transport equation are discussed. Quasi-Linear Theory (QLT) relates the power contained in fluctuations of the Interplanetary Magnetic Field (IMF) to the transport coefficients of energetic particles, an outline of the basic results of this theory is presented followed by a discussion of subsequent corrections made to the original formulation with an emphasis in recent developments where the effects of wave polarization, its propagation respect to the solar wind and the dissipation of power at large frequencies have been taken into account. The numerical approach using test particle trajectory integrations to obtain transport coefficients based on in situ satellite measureents is also discussed. It is well known that the determination of the particles mean free path for solar particle events by alternative methods leads to conflicting results, corrections made to original QLT are attempts to bridge the gap. Determination of the transport parameters from different lines of approach in a comparative basis have been done recently by calculating power spectra of IMF measured at the time solar particles were detected on the same spaceprobe, and performing numerical simulations with equivalent IMF data. Some of the results of such studies point to the solution of the conflicting determinations of the mean free path which has existed for nearly 30 years. An assesment of the present situation in this respect is given. Numerical determinations of transport parameters in the outer heliosphere are also reviewed and its consequences for solar modulation of galactic cosmic rays discussed. Space Science Reviews 62: Printed in Belgium.     

Energetic particle transport coefficients in the heliosphere

It is the purpose of this review to summarize and discuss recent research done in the field of particle propagation in the heliosphere. Several lines of approach have been followed to treat this problem. As a starting point the different forms of the transport equation are discussed. Quasi-Linear Theory (QLT) relates the power contained in fluctuations of the Interplanetary Magnetic Field (IMF) to the transport coefficients of energetic particles, an outline of the basic results of this theory is presented followed by a discussion of subsequent corrections made to the original formulation with an emphasis in recent developments where the effects of wave polarization, its propagation respect to the solar wind and the dissipation of power at large frequencies have been taken into account. The numerical approach using test particle trajectory integrations to obtain transport coefficients based on in situ satellite measureents is also discussed. It is well known that the determination of the particles mean free path for solar particle events by alternative methods leads to conflicting results, corrections made to original QLT are attempts to bridge the gap. Determination of the transport parameters from different lines of approach in a comparative basis have been done recently by calculating power spectra of IMF measured at the time solar particles were detected on the same spaceprobe, and performing numerical simulations with equivalent IMF data. Some of the results of such studies point to the solution of the conflicting determinations of the mean free path which has existed for nearly 30 years. An assesment of the present situation in this respect is given. Numerical determinations of transport parameters in the outer heliosphere are also reviewed and its consequences for solar modulation of galactic cosmic rays discussed. Space Science Reviews 62: Printed in Belgium.     

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.     

Yohkoh results in the context of the high-latitude heliosphere

Designed primarily to study solar activity, Yohkoh includes an X-ray telescope that obtains full-sun coronal images which show a range of features. Coronal X-ray emission-exclusive of flares, is notable for its variability even in the largest structures. A mass ejection event is related to magnetic field reconnection. Such events exhibit both accelerated and decelerated behaviour. Coronal hole temperatures are estimated from the filter ratio method. A plasma component at around 2.10⁶ K is identified. X-ray emission is detected from the South polar coronal hole. A preliminary comparison of Spartan coronagraph images with Yohkoh data suggests that polar plumes or rays are not connected to bright points.     

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.     

Cosmic-ray intensity variations in the 3-dimensional heliosphere

The study of cosmic-ray intensity variations have been carried out with data registered by ground-based and balloon-borne equipment for the past 50 years or more. The International Geophysical Year (IGY) from July 1957 to December 1958 gave an impetus to global collaborations. A world-wide network of concerted measurements became available with the advent of the space age.In situ measurements by satellite-borne detectors led to deep-space exploration. The spacecraft Pioneers and Voyagers, during the past 15 years, traversing farther out into the heliosphere at increasing radial distances from the sun have changed the study of time variations into one of time and spatial variations.Furthermore, with the Voyager 1, proceeding asymptotically towards heliolatitudes of 35° north since its encounter with Saturn and the anticipated direction of Voyager 2 after its encounter with Neptune in late-1989 towards 48° south heliolatitude, is converting the study into a truly three-dimensional exploration of the heliosphere. Thus, the investigation of galactic cosmic-ray intensity variations fromin situ measurements deep in the heliosphere in distance, latitude, and over solar cycles is indeed a remarkable achievement.The various cosmic-ray intensity variations over different time-scales, the modulation of the intensity by the evolving solar activity and the role of the electromagnetic state of the interplanetary medium (otherwise called heliosphere) can now be investigated as never before; these studies contribute immensely to our knowledge of the solar neighbourhood. This article essentially deals with the studies of time and spatial variations of cosmic-ray intensity that have been conducted especially over the past two decades.     

Cosmic ray modulation in the 3-D heliosphere

The basic physical processes that lead to the long-term modulation of cosmic rays by the solar wind have been known for many years. However our knowledge of the structure of the heliosphere, which determines which processes are most important for the modulation, and of the variation of this structure with time and solar activity level is still incomplete. Study of the modulation provides a tool for probing the scale and structure of the heliosphere. While the Pioneer and Voyager spacecraft are surveying the radial structure and extent of the heliosphere at modest heliographic latitudes, theUlysses mission is the first to undertake a nearly complete scan of the latitudinal structure of the modulated cosmic ray intensity in the inner heliosphere (R<5.4 AU).Ulysses will reach latitudes of 80°S in September 1994 and 80°N in July 1995 during the approach to minimum activity in the 11 year solar cycle. We present a first report of measurements extending to latitudes of 52°S, which show surprisingly little latitudinal effect in the modulated intensities and suggest that at this time modulation in the inner heliosphere may be much more spherically symmetric than had generally been believed based upon models and previous observations.