Voyager energetic particle observations at interplanetary shocks and upstream of planetary bow shocks: 1977–1990

The Voyager 1 and 2 spacecraft include instrumentation that makes comprehensive ion (E 28 keV) and electron (E 22 keV) measurements in several energy channels with good temporal, energy, and compositional resolution. Data collected over the past decade (1977–1988), including observations upstream and downstream of four planetary bow shocks (Earth, Jupiter, Saturn, Uranus) and numerous interplanetary shocks to 30 AU, are reviewed and analyzed in the context of the Fermi and shock drift acceleration (SDA) models. Principal findings upstream of planetary bow shocks include the simultaneous presence of ions and electrons, detection of tracer ions characteristic of the parent magnetosphere (O, S, O⁺), power-law energy spectra extending to 5 MeV, and large (up to 100:1) anisotropies. Results from interplanetary shocks include observation of acceleration to the highest energies ever seen in a shock ( 22 MeV for protons, 220 MeV for oxygen), the saturation in energy gain to 300 keV at quasi-parallel shocks, the observation of shock-accelerated relativistic electrons, and separation of high-energy (upstream) from low-energy (downstream) populations to within 1 particle gyroradius in a near-perpendicular shock. The overall results suggest that ions and electrons observed upstream of planetary bow shocks have their source inside the parent magnetosphere, with first order Fermi acceleration playing a secondary role at best. Further, that quasi-perpendicular interplanetary shocks accelerate ions and electrons most efficiently to high energies through the shock-drift process. These findings suggest that great care must be exercised in the application of concepts developed for heliosphere shocks to cosmic ray acceleration through shocks at supernova remnants.     

Initial ISEE magnetometer results: Shock observation

ISEE-1 and -2 magnetic field profiles across 6 terrestrial bow shock and one interplanetary shock are examined. The interplanetary shock illustrates the behavior of a low Mach number shock. It had an upstream whistler wave precursor with an apparent wavelength of 180 km. The shock thickness was about 90 km for the thickness of the final field jump or 270 km for the exponential growth of the precursor wave packet. The ion inertial length was 50 km, upstream of the shock.Three examples of low or moderate , high Mach number, quasiperpendicular shocks are examined. These did not have upstream waves, but rather had waves growing in the field gradient. The growth length for these waves and the shock profile was of the order of the ion inertial length.Two examples of high shocks showed little coherence in field variation even though the two vehicles were only a few hundred kilometers apart. Thus, we cannot gauge their velocity and turn the time profiles into distances. The final crossing examined shows clearly the effect of changing the orientation of the interplanetary magnetic field. Initially the upstream magnetic field made an angle of about 80° to the shock normal and the shock position remained fairly steady. Then the field rotated to 45° to the normal and the field profiles became very irregular and the shock position very unstable. Discrete wave packets appeared.Finally, we present the joint behavior of wave, particle and field data across some of these shocks to show some of the myriad of shock features whose behavior we are now beginning to investigate.     

Pulsars as gamma ray sources: Nebular shocks and magnetospheric gaps

I summarize the results of recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is primarily electrons and positrons with an admixture of heavy ions. Shocks which contain heavy ions that are a minority constituent by number but which carry most of the energy density in the upstream medium put 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) E ^-2, where N(E)dE is the number of particles with energy between E and E+dE. Synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front, provides the mechanism of thermalization and non-thermal particle acceleration. The maximum energy achievable by the pairs is _± m _± c ² = m _i c ² ₁/Z _i, where ₁ is the Lorentz factor of the upstream flow and Z _i is the atomic number of the ions. The shock's spatial structure contains a series of overshoots in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. These overshoots provide a new interpretation of the structure of the inner regions of the Crab Nebula, in particular of the wisps, surface brightness enhancements near the pulsar. The wisps appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure ₁ 4 × 10⁶, the upstream magnetic field B ₁ to be 3 × 10^-5 Gauss, as well as to show that the total ion flow is 3 × 10³⁴ elementary charges/sec, in good agreement with the total current flow predicted by the early Goldreich and Julian (1969) model. The total pair outflow is shown to be about 5 × 10³⁷ pairs per second, in good agreement with the particle flux required to explain the nebular X—ray source.The energetics of particle acceleration within the magnetospheres of rotation powered pulsars and the consequences for pulsed gamma ray emission are also briefly discussed. The gamma ray luminosity above 100 MeV is shown to scale in proportion to _R ^1/2 , as is in accord with some of the simplest ideas about polar cap models. Models based on acceleration in the outer magnetosphere are also briefly discussed.     

The importance of anisotropic interstellar turbulence and molecular-cloud magnetic mirrors for galactic cosmic-ray propagation

Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale fluctuations satisfy the inequality k k , where k and k are, respectively, the components of a fluctuations wave vector and to the background magnetic field. Such anisotropic fluctuations are very inefficient at scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in confining and isotropizing cosmic rays when scattering is weak is also discussed.     

LISM structure — Fragmented superbubble shell?

Small scale structure in local interstellar matter (LISM) is considered. Overall morphology of the local cloud complex is inferred from Ca II absorption lines and observations of H I in white dwarf stars. Clouds with column densities ranging from 2–100 × 10¹⁷ cm^–2 are found within 20 pc of the Sun. Cold (50 K) dense (10⁵ cm^–3) small (5–10 au) clouds could be embedded and currently undetected in the upwind gas. The Sun appears to be embedded in a filament of gas with thickness 0.7 pc, and cross-wise column density 2 × 10¹⁷ cm^–2. The local magnetic field direction is parallel to the filament, suggesting that the physical process causing the filamentation is MHD related. Enhanced abundances of refractory elements and LISM kinematics indicate outflowing gas from the Scorpius-Centaurus Association. The local flow vector and Sco data are consistent with a 4,000,000 year old superbubble shell at –22 km s^–1, which is a shock front passing through preshock gas at –12 km s^–1, and yielding cooled postshock gas at –26 km s^–1in the upwind direction. A preshock magnetic field strength of 1.6 G, and postshock field strength of 5.2 G embedded in the superbubble shell, are consistent with the data.Abbreviations LISM Local ISM - SIC Surrounding Interstellar Cloud - LIC Local Interstellar Cloud     

Initial observations of low energy charged particles near the Earth's bow shock on ISEE-1

We report initial measurements from the ULECA sensor of the Max-Planck-Institut/University of Maryland experiment on ISEE-1. ULECA is an electrostatic deflection — total energy sensor consisting of a collimator, deflection analyzer and an array of solid state detectors. The position of a given detector, which determines the energy per charge of an incident particle, together with the measured energy determine the particle's charge state. We find that a rich variety of phenomena are operative in the transthermal energy regime (10 keV/Q to 100 keV/Q) covered by ULECA. Specifically, we present observations of locally accelerated protons, alpha particles, and heavier ions in the magnetosheath and upstream of the Earth's bow shock. Preliminary analysis indicates that the behavior of these locally accelerated particles is most similar at the same energy per charge.     

Temperature Fluctuations,H II Region Abundances,and Primordial Helium

There is evidence for temperature fluctuations in Planetary Nebulae and in some Galactic H II regions. If such fluctuations occur in the low metallicity, extragalactic H II regions used to probe the primordial helium abundance, the derived 4He mass fraction, YP, could be systematically different from the true primordial value. Although this effect could be large, there are no data which allow us to estimate the size of the temperature fluctuations for the extragalactic H II regions. Therefore, we have explored this effect via Monte Carlo simulations of the data in which the abundances derived from a fiducial data set are modified by T chosen from a distribution with 0 T Tmax where Tmax is varied from 500 K to 4000 K.     

Dual spacecraft observations of energetic particles in the vicinity of the magnetopause,bow shock,and the interplanetary medium

This article presents some of the new and important particle features that have been detected in the energy range 1 keV to 290 keV by the ISEE-1 and -2 spacecraft near the magnetopause, bow shock, and the interplanetary space. Only examples of data from the first few orbits, when the spacecraft were on the front side, are shown.Paper presented at 13th ESLAB Symposium, Innsbruck, Austria (June 5, 1978).     

The effect of rotation on the hydrodynamics of stellar collapse

Summary Using values of d, _min, and _max that Van Riper (1978) has found most promising for a hydrodynamic envelope ejection, we have shown that even a small amount of rotation in the initial core can stop its collapse before nuclear densities are reached. We expected _i > 0.02 to produce significant deviations from a spherically symmetric collapse, but have found that _i as much as ten times smaller than this will not allow the core to reach densities as high as in the spherical collapse. In no case, however, does the core flatten very much, nor does the value of become very large. Low final 's preclude the formation of an axisymmetric torus. They also indicate that deformation of an iron core into a triaxial configuration or fragmentation of the core during its collapse is an extremely unlikely event. (Note: Classically, must exceed 0.27 before a dynamic instability to non-axisymmetric perturbations is encountered.)The small degree of flattening of the core also suggests that the reduced moment of inertia I of the core will always be relatively small in magnitude and hence that the third time derivative of I, which is proportional to the energy emitted in gravity wave radiation, will not be very significant. Numerically calculated estimates of I- during some of these model evolutions supports this suspicion. If the _min and used here are found to be realistic values after the detailed physics of the core collapse is well understood, it is clear that gravitational radiation from a core collapse will be difficult to measure.Finally, we should point out that it is the relatively large values of Y_min (near 4/3) combined with values of d near unity that (a) prevented the core from flattening significantly in these models and (b) prevented the core from reaching high configurations. If realistic values of either one (or both) of these parameters are found to be much smaller in more complete models of the core collapse, then the core will have to become flatter (and denser) before pressure gradients will support it along the rotation axis. All of the conclusions drawn here would be modified accordingly under those circumstances. It should also be noted that in general relativistic models, the critical for spherical collapse is somewhat larger than 4/3 (Van Riper, 1979). Therefore, we predict that when fully general relativistic core collapses are performed including rotation, a given choice of _min and _i will produce a slightly flatter and slightly denser core than the corresponding model that has been presented here.     

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.     

Unusually Distant Bow Shock Encounters at Mars: Analysis of March 24, 1989 Event

Detailed analysis of disturbances observed on 24 March, 1989 far upstream of the usual Martian bow shock position was completed with the use of the planetary obstacle and bow shock models relevant for the period of Phobos 2 observations and for low Mach numbers, respectively. It is proven that the system of discontinuities observed in the solar wind between 18:42 and 19:36 UT was the consequence of unusually distant planetary bow shock excursions. The cause was unusually small ρV ² and M _a values in the solar wind flow.     

Fourier parameters of heliospheric current sheet and their significance

If the path of the neutral line on the coronal source surface is expressible as a singlevalued function (colatitude vs longitude ), then Fourier analysis of ctn with respect to leads to a simple algorithm for realistically mapping the neutral line outward to model the heliospheric current sheet (HCS) at distancesr1 AU. To be compatible with MHD, the source surface used for this mapping should be prolate (aligned with dipole axis) rather than spherical. Orientation of the Sun's magnetic-dipole moment is indicated by them=1 Fourier amplitude (a ₁ sin +b ₁ cos ) of ctn on the source surface. Physical features (including the neutral line) on a prolate source surface intrinsically map to lower dipole latitudes atr1 AU in the heliosphere, and Ulysses observations of a unipolar field at latitudes beyond 30°S (when the neutral line on the source surface still reached 39°S) confirm the expected geometry.     

Emission mechanism for low-frequency radiation in the outer heliosphere

The question of how low-frequency radio emissions in the outer heliosphere might be generated is considered. It is argued that the free energy contained in an electron beam distribution is first transformed into electrostatic Langmuir waves. The nonlinear interactions of these waves which can produce electromagnetic waves are then treated in the semi-classical formalism. Comparison of the results of the discussed model with electromagnetic radiation coming from upstream of the Earth's bow shock shows that the model adequately explains the generation of plasma waves at planetary shocks. By analogy, this model can provide a quantitative explanation of intensity of radio emissions at 2 to 3 kHz detected by the Voyager plasma wave instrument in the outer heliosphere provided that the electron beams generating Langmuir waves exist also in the postshock plasma due to secondary shocks in the compressed solar wind beyond the termination shock. The field strength of Langmuir waves required to generate the second harmonic emissions are approximately of 100–200 V m^–1 for the primary and 50–100 V m^–1 for the secondary foreshocks. However, only in the secondary foreshock the expected density is consistent with the observed frequency.     

Finite Larmor radius effects in the magnetosphere

This article reviews theories and observations related to effects produced by finite (and large) Larmor radii of charged particles in the magnetosphere. The FLR effects depend on =r _H /L, wherer _H is the Larmor radius andL is the spatial scale for field/plasma inhomogeneity. The parameter is a basic expansion parameter for most equations describing plasma dynamics in the magnetosphere. The FLR effects enter naturally the drift approximation for particle motion and represent also non-ideal MHD terms in the fluid formalism. The linear and higher order terms in lead to charge separation, energization of particles, and produce viscosity without collisions. The FLR effects introduce also important corrections to the dispersion relations for MHD waves and drift instabilities. Expansion of plasma into magnetic field leads to filamentation of the plasma boundary and to creation of structures with thickness less than an ion gyroradius. Large Larmor radius effects (1) in curved magnetic field geometry lead to stochastic behaviour of particle trajectories and to deterministic chaos. The tiny scale of the electron and ion gyroradii does not necessarily mean that FLR/LLR phenomena have negligible effect on the macroscopic dynamics and energetics of the whole magnetosphere. On the contrary, the small scale gyro-effects may provide the physical mechanism for gyroviscous coupling between the solar wind and the magnetosphere, the mechanism for triggering disruption of the magnetotail current layer, and the mechanism for parallel electric field that accelerate auroral particles.     

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.     

Turbulent transition in solar surges

It has been suggested that a surge can be modelled as a jet travelling in a sheared magnetic field, and that the transition to turbulence of this MHD tearing jet can explain several key observed features. In this paper we present our preliminary results of the transition to turbulencevia secondary instabilities of the MHD tearing jet. Our results confirm that turbulent transition can decelerate the surge, with decay times which compare well with surge data. Furthermore, we find that the turbulent MHD tearing jet forms magnetic field-aligned velocity filaments similar to those often observed in the surge flow field.     

High temporal resolution observations of electron heating at the bow shock

High temporal resolution measurements of solar wind electrons at the Earth's bow shock on the dawn side have been made with the LASL/MPI fast plasma experiments on ISEE-1 and 2. One dimensional, 1-d, temperatures, T _e , and densities, N _e , are obtained every 0.3 s and 2-d values are obtained every 3 s. Profiles of T _e and N _e at the shock usually are found to be similar to one another and also to the profile of the magnetic field magnitude. The time scale of electron thermalization varies from about 0.5 s to greater than 1 min, depending importantly on the shock motion and the orientation of the magnetic field. Typical thermalization times from 05:00–12:00 LT are 10 s, considerably shorter than proton thermalization times at the shock. This time scale corresponds to a distance of 100 km, comparable to but somewhat larger than the typical ion inertial length. The electron thermalization times are significantly longer than some of the values frequently cited in the past. At the end of the electron thermalization there typically is an overshoot in electron thermal pressure followed by an undershoot which give the pressure profile of the shock the appearance of a damped wave. Ion thermalization is essentially completed by the time the electron pressure wave is damped. The most probable value of the electron temperature ratio across the shock is 1.7, and this value is relatively independent of the Sun-Earth-satellite angle, _ss , for _ss between 25° and 100°.The Los Alamos Scientific Laboratory requests that the publisher identify this article as work performed under the auspices of the Department of Energy.By acceptance of this article, the publisher recognizes that the U.S. Government retains a non-exclusive, royalty-free license to publish or reproduce the published form of this contribution, or to allow others to do so, for U.S. Government purposes.     

Interplanetary shock waves generated by solar flares

Recent observational and theoretical studies of interplanetary shock waves associated with solar flares are reviewed. An attempt is made to outline the framework for the genesis, life and demise of these shocks. Thus, suggestions are made regarding their birth within the flare generation process, MHD wave propagation through the chromosphere and inner corona, and maturity to fully-developed coronal shock waves. Their subsequent propagation into the ambient interplanetary medium and disturbing effects within the solar wind are discussed within the context of theoretical and phenomenological models. The latter — based essentially on observations — are useful for a limited interpretation of shock geometric and kinematic characteristics. The former — upon which ultimate physical understanding depends — are used for clarification and classification of the shocks and their consequences within the solar wind. Classification of limiting cases of blast-produced shocks (as in an explosion) or longer lasting ejecta (or piston-driven shocks) will hopefully be combined with the study of the flare process itself.The theoretical approach, in spite of its contribution to clarification of various concepts, contains some fundamental limitations and requires further study. Numerical simulations, for example, depend upon a non-unique set of multi-parameter initial conditions at or near the Sun. Additionally, the subtle but important influence of magnetic fields upon energy transport processes within the solar wind has not been considered in the numerical simulation approach. Similarity solutions are limited to geometrical symmetries and have not exploited their potential beyond the special cases of the blast and the constant-velocity, piston-driven shock waves. These continuum fluid studies will probably require augmentation or even replacement by plasma kinetic theory in special situations when observations indicate the presence of anomalous transport processes. Presently, for example, efforts are directed toward identification of detailed shock structures (as in the case of Earth's bow shock) and of the disturbed solar wind (such as the piston).Further progress is expected with extensive in situ and remote monitoring of the solar wind over a wide range of heliographic radii, longitudes and latitudes.This paper is a revised and updated version of an invited review originally presented at the IUGG XV General Assembly, Moscow, U.S.S.R., 2–14 August 1971.     

Large electric fields in the magnetosphere

In several regions of the magnetosphere, perpendicular and/or parallel electric fields are found to be orders-of-magnitude larger than expected from simple considerations. Problems associated with these large fields that may be amenable to study through computer simulations are discussed. Regions in which large electric fields are observed include: a) The auroral ionosphere, where Langmuir soliton-like structures have been measured to contain plasma frequency oscillations as large as 500 mV/m, the envelopes of which have parallel electric fields of 100 mV/m lasting for fractions of a millisecond; b) The auroral acceleration region, where electrostatic shocks have been observed to contain perpendicular fields as large as 1000 mV/m and parallel fields as large as 100 mV/m, and where double layers having parallel fields up to 10 mV/m have been observed; c) The high latitude boundary of the plasma sheet, where turbulent electric fields as large as 100 mV/m have been seen along with quasi-static fields of 5–10 mV/m; d) Inside the plasma sheet, where fields of 5–10 mV/m have frequently been observed; e) The bow shock, where turbulent fields as large as 100 mV/m and d.c. fields of 5 mV/m normal to the shock have been seen.also Physics Department     

Anomalous Cosmic Rays and Solar Modulation

We review aspects of anomalous cosmic rays (ACRs) that bear on the solar modulation of energetic particles in the heliosphere. We show that the latitudinal and radial gradients of these particles exhibit a 22-year periodicity in concert with the reversal of the Sun's magnetic field. The power-law index of the low energy portion of the energy spectrum of ACRs at the shock in 1996 appears to be -1.3, suggesting that the strength of the solar wind termination shock at the helioequatorial plane is relatively weak, with s 2.8. The rigidity dependence of the perpendicular interplanetary mean free path in the outer heliosphere for particles with rigidities between 0.2 and 0.7 GV varies approximately as R2, where R is particle rigidity. There is evidence that ACR oxygen is primarily multiply charged above 20 MeV/nuc and primarily singly-charged below 16 MeV/nuc. The location of the termination shock was at 65 AU in 1987 and 85 AU in 1994.