The interplanetary medium and its interaction with the Earth's magnetosphere

This paper reviews the principal results of direct measurements of the plasma and magnetic field by spacecraft close to the Earth (within the heliocentric distance range 0.7–1.5 AU). The paper gives an interpretation of the results for periods of decrease, minimum and increase of the solar activity. The following problems are discussed: the interplanetary plasma (chemical composition, density, solar wind flow speed, temperature, temporal and spatial variation of these parameters), the interplanetary magnetic field (intensity, direction, fluctuations and its origin), some derived parameters characterizing the physical condition of the interplanetary medium; the quasi-stationary sector structure and its connection with solar and terrestrial phenomena; the magnetohydrodynamic discontinuities in the interplanetary medium (tangential discontinuities and collisionless shock waves); the solar magnetoplasma interaction with the geomagnetic field (the collisionless bow shock wave, the magnetosheath, the magnetopause, the Earth's magnetic tail, the internal magnetosphere characteristics), the connection between the geomagnetic activity and the interplanetary medium and magnetosphere parameters; peculiarities in behaviour of the interplanetary medium and magnetosphere during geomagnetic storms; energetic aspects of the geomagnetic storms.     

Magnetohydrodynamic discontinuities in space plasmas: Interrelation between stability and structure

It is believed that shock waves and other discontinuous flows form the basis of a wide range of phenomena in space medium. We review the results concerning the particular property of MHD discontinuities, the interrelation between their stability and structure. Such an interrelation is associated primarily with the requirement of their evolutionarity. For a non-evolutionary discontinuity the amplitudes of reflected and refracted waves caused by a small amplitude incident wave are not determined unambiguously from the conservation laws at the discontinuity surface. Since the problem of the further time evolution of the initial small perturbation does not have a unique solution, such a discontinuity cannot exist in a real medium as a single steady-state configuration. Therefore it is unstable with respect to a disintegration into several discontinuities or to a transition to some more general unsteady flow. This is confirmed by the fact that, in the studied cases, the non-evolutionary shock transitions can, while the evolutionary ones cannot, be represented as a set of several discontinuities and waves of rarefaction. One more fundamental property of the non-evolutionary shocks that argues for their non-existence is that they do not have a unique structure for all values of the dissipative transport coefficients. At the same time, the possibility of their existence cannot be excluded when the shock has a unique structure or when the disintegration is forbidden for some reasons. Besides the non-evolutionarity, which is a direct reason for a disintegration, there is an indirect one. It can be shown that the hydrodynamic shocks without magnetic field that are corrugationally unstable also allow the shock transitions through more than one discontinuity. This suggests that the shocks unstable in the ordinary sense, for which the small perturbation grows with time, do not exist at all, but they disintegrate into stable ones. However, the physical mechanism that distinguishes between these two scenarios remains unclear. The interrelation between the stability and structure of MHD discontinuities thus governs their nonlinear evolution. This fact is essential when describing shock waves and other discontinuous flows in the space medium.     

Dayside Proton Aurora: Comparisons between Global MHD Simulations and IMAGE Observations

The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the IMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model, though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside ionosphere are enhanced. Overall the results of these initial comparisons between global MHD simulation results and IMAGE observations emphasize the interplay between reconnection and dynamic pressure processes at the dayside magnetopause, as well as the intricate connection between the bow shock and the auroral region.     

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.     

Theoretical studies of the solar wind phenomenon

This paper is a review of current theoretical topics concerning the solar wind. Broadly speaking the questions outstanding at the present time concern the loss of angular momentum to the sun, the origin of the fluctuations observed in the wind at the orbit of earth, conditions in the wind in regions yet unvisited by spacecraft (inside the orbit of Venus, beyond the orbit of Mars, and out of the plane of the ecliptic), conditions at the terminus of the wind, etc. The question of angular momentum loss is important in understanding the evolution of the sun to its present form with a slowly rotating surface. Evidence from both comet and spacecraft observations of the wind indicate that the rate at which angular momentum is being carried away by the solar wind is very large, of the order of 10³¹ dyne/cm in the gas flow and half as much by the interplanetary magnetic field. But theory cannot account for more than about 10³⁰ dynes/cm in the gas without special assumptions.The fluctuations presently observed in the wind at the orbit of earth have scales ranging upward from 10² km. Their presence is puzzling because fluctuations with scales less than about 10⁶ km are not expected to survive from the sun. Presumably, therefore, the fluctuations are generated by the velocity differences of more than 100 km/sec in the wind from different regions in the solar corona and by instabilities produced by the anisotropy of the electrons of the wind plasma.Conditions in the wind at places far removed from the orbit of earth can be inferred from the behavior of cosmic rays. The evidence is that the wind becomes relatively placid beyond about 5 AU, extending from there out to 30–300 AU without much small-scale turbulence. There are also some suggestions that the wind may perhaps be less turbulent toward the sun from 1 AU, and that the wind may be faster and more turbulent at higher solar latitudes. But the ambiguity of the situation does not permit a firm conclusion on this yet.     

Structure of reconnection layers in the magnetosphere

Magnetic reconnection can lead to the formation of observed boundary layers at the dayside magnetopause and in the nightside plasma sheet of the earth's magnetosphere. In this paper, the structure of these reconnection layers is studied by solving the one-dimensional Riemann problem for the evolution of a current sheet. Analytical method, resistive MHD simulations, and hybrid simulations are used. Based on the ideal MHD formulation, rotational discontinuities, slow shocks, slow expansion waves, and contact discontinuity are present in the dayside reconnection layer. Fast expansion waves are also present in the solution of the Riemann problem, but they quickly propagate out of the reconnection layer. Our study provides a coherent picture for the transition from the reconnection layer with two slow shocks in Petschek's model to the reconnection layer with a rotational discontinuity and a slow expansion wave in Levy et al's model. In the resistive MHD simulations, the rotational discontinuities are replaced by intermediate shocks or time-dependent intermediate shocks. In the hybrid simulations, the time-dependent intermediate shock quickly evolves to a steady rotational discontinuity, and the contact discontinuity does not exist. The magnetotail reconnection layer consists of two slow shocks. Hybrid simulations of slow shocks indicate that there exists a critical number,M _c, such that for slow shocks with an intermediate Mach numberM _IM _c, a large-amplitude rotational wavetrain is present in the downstream region. For slow shocks withM _I<M _c, the downstream wavetrain does not exist. Chaotic ion orbits in the downstream wave provide an efficient mechanism for ion heating and wave damping and explain the existence of the critical numberM _c in slow shocks.     

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.     

The acceleration of charged particles in interplanetary shock waves

The theory and observations of energetic ion acceleration in interplanetary shock waves is reviewed. The shock acceleration of the solar wind plasma and particle transport effects are discussed. Suggestions are offered for future research in shock acceleration physics.An invited paper presented at STIP Workshop on Shock Waves in the Solar Corona and Interplanetary Space, 15–19 June, 1980, Smolenice, Czechoslovakia.NAS/NRC Research Associate.     

Global Response of Martian Plasma Environment to an Interplanetary Structure: From Ena and Plasma Observations at Mars

As a part of the global plasma environment study of Mars and its response to the solar wind, we have analyzed a peculiar case of the subsolar energetic neutral atom (ENA) jet observed on June 7, 2004 by the Neutral Particle Detector (NPD) on board the Mars Express satellite. The “subsolar ENA jet” is generated by the interaction between the solar wind and the Martian exosphere, and is one of the most intense sources of ENA flux observed in the vicinity of Mars. On June 7, 2004 (orbit 485 of Mars Express), the NPD observed a very intense subsolar ENA jet, which then abruptly decreased within ∼10 sec followed by quasi-periodic (∼1 min) flux variations. Simultaneously, the plasma sensors detected a solar wind structure, which was most likely an interplanetary shock surface. The abrupt decrease of the ENA flux and the quasi-periodic flux variations can be understood in the framework of the global response of the Martian plasma obstacle to the interplanetary shock. The generation region of the subsolar ENA jet was pushed towards the planet by the interplanetary shock; and therefore, Mars Express went out of the ENA jet region. Associated global vibrations of the Martian plasma obstacle may have been the cause of the quasi-periodic flux variations of the ENA flux at the spacecraft location.     

3-D magnetic field and current system in the heliosphere

In this paper a global system of the magnetic field and current from the interaction of the solar wind plasma and the interstellar medium is modeled using a 3-D MHD simulation. The terminal shock, the heliopause and the outer shock are clearly determined in our simulation. In the heliosheath the toroidal magnetic field is found to increase with the distance from the sun. The magnetic field increases rapidly in the upstream region of the heliosheath and becomes maximum between the terminal shock and the heliopause. Hence a shell-type magnetic wall is found to be formed in the heliosheath. Because of this magnetic wall the radially expanding solar wind plasma changes its direction tailward in all latitudes except the equatorial region. Only the equatorial disk-like plasma flow is found to extend to the heliopause through the weak magnetic-field region around the equator. Two kinds of global current loops which sustain the toroidal magnetic field in the heliosphere are found in our simulation.The influence of the 11-year solar cycle variation of the magnetic polarity is also examined. It is found that the polarity of the toroidal magnetic field in the heliosheath switches at every solar cycle change. Hence the heliosheath is found to consist of such magnetized plasma bubbles. The neutral sheets are found to extend between such magnetized plasma bubbles in the 3-D heliosheath in a complicated form. The magnetic-pressure effect on the heliosheath plasma structure is also examined.     

MHD Waves and Heating in Coronal Holes

Coronal holes have been identified as source regions of the fast solar wind, and MHD wave activity has been detected in coronal holes by remote sensing, and in situ in fast solar wind streams. I review some of the most suggestive wave observations, and discuss the theoretical aspects of MHD wave heating and solar wind acceleration in coronal holes. I review the results of single fluid 2.5D MHD, as well as multi-fluid 2.5D MHD models of waves in coronal holes, the heating, and the acceleration of the solar wind be these waves.     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

The interplanetary magnetic field. Solar origin and terrestrial effects

Many observations related to the large-scale structure of the interplanetary magnetic field, its solar origin and terrestrial effects are discussed. During the period observed by spacecraft the interplanetary field was dominated by a sector structure corotating with the sun in which the field is predominantly away from the sun (on the average in the Archimedes spiral direction) for several days (as observed near the earth), and then toward the sun for several days, etc. The average sector appears to be a coherent entity with internal structure such that its preceding portion is more active than its following portion. Cosmic rays corotate with the interplanetary field, and there are differential flows associated with the sector pattern. Profound effects on geomagnetic activity and the radiation belts are produced as the sector pattern rotates past the earth. The solar origin of the sector pattern is discussed. The solar source may be associated with the large-scale weak background photospheric fields observed with the solar magnetograph. It is suggested that there may be a rather continual relation between this solar structure and terrestrial responses, of which the recurring M-Region geomagnetic storms are just the most prominent example.     

Ground-based solar radio observations of the August 1972 events

Ground-based observations of the variable solar radio emission ranging from few millimetres to decametres have been used here as a diagnostic tool to gain coherent phenomenological understanding of the great 2, 4 and 7 August, 1972 solar events in terms of dominant physical processes like generation and propagation of shock waves in the solar atmosphere, particle acceleration and trapping.The basic data used in this review have been collected by many workers throughout the world utilizing a variety of instruments such as fixed frequency radiometers, multi-element interferometers, dynamic spectrum analysers and polarimeters. Four major flares are selected for detailed analysis on the basis of their ability to produce energetic protons, shock waves, polar cap absorptions (PCA) and sudden commencement (SC) geomagnetic storms. A comparative study of their radio characteristics is made. Evidence is seen for the pulsations during microwave bursts by the mechanism similar to that proposed by McLean et al. (1971), to explain the pulsations in the metre wavelength continuum radiation. It is suggested that the multiple peaks observed in some microwave bursts may be attributable to individual flares occurring sequentially due to a single initiating flare. Attempts have been made to establish identification of Type II bursts with the interplanetary shock waves and SC geomagnetic storms. Furthermore, it is suggested that it is the mass behind the shock front which is the deciding factor for the detection of shock waves in the interplanetary space. It appears to us that more work is necessary in order to identify which of the three moving Type IV bursts (Wild and Smerd, 1972), namely, advancing shock front, expanding magnetic arch and ejected plasma blob serves as the piston-driver behind the interplanetary shocks. The existing criteria for proton flare prediction have been summarized and two new criteria have been proposed. Observational limitations of the current ground-based experimental techniques have been pointed out and a suggestion has been made to evolve appropriate observational facilities for solar work before the next Solar Maximum Year (SMY).     

Microscopic structure of the solar wind

As the solar wind flows out from the coronal base the coulomb collision frequencies rapidly become small and particle-particle collisions can no longer maintain local statistical equilibrium. At 1 AU the particle distribution functions have important non-Maxwellian characteristics and the firehose instability, a cyclotron resonance whistler-mode instability, and several heat flux current instabilities should be operative. Superthermal particle populations also provide large wave levels, and other forms of enhanced plasma turbulence develop at shock fronts and discontinuities. This report contains a review of the theoretical concepts and a progress report on the experimental study of interplanetary wave-particle interactions.Prepared for Space Science Reviews.     

Ion acceleration at shocks in interplanetary space: A brief review of recent observations

The mechanism by which ions are accelerated near the Earth's bow shock and near shocks propagating outward from the Sun in response to solar activity appears to be essentially the same. For both types of shock the solar wind thermal distribution acts as a seed population. Leaked magnetospheric ions and resident flare ions are additional seed populations for the bow shock and outward propagating shocks respectively. The acceleration of solar wind ions at these shocks begins with either the reflection of ions off the shock or leakage of shocked plasma back through the shock. Interaction with a disruption wave field self-generated by these backstreaming ions is responsible for the remainder of the acceleration at the bow shock. Both the disruption wave field and the ambient interplanetary wave field play important roles in accelerating ions at outward propagating shocks, but on different time scales. The geometry of the shock and the duration of field line connection to the shock play decisive roles in determining what is observed.