Recent Progress in Physics-Based Models of the Plasmasphere

We describe recent progress in physics-based models of the plasmasphere using the fluid and the kinetic approaches. Global modeling of the dynamics and influence of the plasmasphere is presented. Results from global plasmasphere simulations are used to understand and quantify (i) the electric potential pattern and evolution during geomagnetic storms, and (ii) the influence of the plasmasphere on the excitation of electromagnetic ion cyclotron (EMIC) waves and precipitation of energetic ions in the inner magnetosphere. The interactions of the plasmasphere with the ionosphere and the other regions of the magnetosphere are pointed out. We show the results of simulations for the formation of the plasmapause and discuss the influence of plasmaspheric wind and of ultra low frequency (ULF) waves for transport of plasmaspheric material. Theoretical models used to describe the electric field and plasma distribution in the plasmasphere are presented. Model predictions are compared to recent Cluster and Image observations, but also to results of earlier models and satellite observations.     

An Overview of Results From rpi on Image

The Radio Plasma Imager (RPI) on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) spacecraft was designed as a long-range magnetospheric radio sounder, relaxation sounder, and a passive plasma wave instrument. The RPI is a highly flexible instrument that can be programmed to perform these types of measurements at times when IMAGE is located in key regions of the magnetosphere. RPI is the first radio sounder ever flown to large radial distances into the magnetosphere. The long-range sounder echoes from RPI allow remote sensing of a variety of plasmas structures and boundaries in the magnetosphere. A profile inversion technique for RPI echo traces has been developed and provides a method for determining the density distribution of the plasma from either direct or field-aligned echoes. This technique has enabled the determination of the evolving density structure of the polar cap and the plasmasphere under a variety of geomagnetic conditions. New results from RPI show that the plasmasphere refills in slightly greater than a day at L values of 2.8 and that ion heating is probably playing a major role in the overall density distribution along the field-line. In addition, RPI's plasma resonance observations at large radial distances over the polar cap provided in situ measurements of the plasma density with an accuracy of a few percent. For the first time in the magnetosphere, RPI has also observed the plasma D resonances. RPI's long antennas and its very low noise receivers provide excellent observations in the passive receive-only mode when the instrument measures the thermal plasma noise as well as natural emissions such as the continuum radiation and auroral kilometric radiation (AKR). Recent passive measurements from RPI have been compared extensively with images from the Extreme Ultraviolet (EUV) imager on IMAGE resulting in a number of new discoveries. For instance, these combined observations show that kilometric continuum can be generated at the plasmapause from sources in or very near the magnetic equator, within a bite-out region of the plasmasphere. The process by which plasmaspheric bite-out structures are produced is not completely understood at this time. Finally, RPI has been used to successfully test the feasibility of magnetospheric tomography. During perigee passages of the Wind spacecraft, RPI radio transmissions at one and two frequencies have been observed by the Waves instrument. The received electric field vector was observed to rotate with time due to the changing density of plasma, and thus Faraday rotation was measured. Many future multi-spacecraft missions propose to use Faraday rotation to obtain global density pictures of the magnetosphere.     

Quantifying Global Plasmaspheric Images With in situ Observations

Simultaneous IMAGE EUV plasmaspheric images and Magnetospheric Plasma Analyzer (MPA) data from the Los Alamos National Laboratory's geosynchronous satellites are combined to understand plasmaspheric behavior and to quantify the global images. A brief review of the understanding of the plasmasphere as learned from in situ observations prior to the launch of IMAGE is given to place the results presented here into context.     

Plasmasphere Response: Tutorial and Review of Recent Imaging Results

The plasmasphere is the cold, dense innermost region of the magnetosphere that is populated by upflow of ionospheric plasma along geomagnetic field lines. Driven directly by dayside magnetopause reconnection, enhanced sunward convection erodes the outer layers of the plasmasphere. Erosion causes the plasmasphere outer boundary, the plasmapause, to move inward on the nightside and outward on the dayside to form plumes of dense plasma extending sunward into the outer magnetosphere. Coupling between the inner magnetosphere and ionosphere can significantly modify the convection field, either enhancing sunward flows near dusk or shielding them on the night side. The plasmaspheric configuration plays a crucial role in the inner magnetosphere; wave-particle interactions inside the plasmasphere can cause scattering and loss of warmer space plasmas such as the ring current and radiation belts.     

The Extreme Ultraviolet Imager Investigation for the IMAGE Mission

The Extreme Ultraviolet Imager (EUV) of the IMAGE Mission will study the distribution of He⁺ in Earth's plasmasphere by detecting its resonantly-scattered emission at 30.4 nm. It will record the structure and dynamics of the cold plasma in Earth's plasmasphere on a global scale. The 30.4-nm feature is relatively easy to measure because it is the brightest ion emission from the plasmasphere, it is spectrally isolated, and the background at that wavelength is negligible. Measurements are easy to interpret because the plasmaspheric He⁺ emission is optically thin, so its brightness is directly proportional to the He⁺ column abundance. Effective imaging of the plasmaspheric He⁺ requires global `snapshots in which the high apogee and the wide field of view of EUV provide in a single exposure a map of the entire plasmasphere. EUV consists of three identical sensor heads, each having a field of view 30° in diameter. These sensors are tilted relative to one another to cover a fan-shaped field of 84°×30°, which is swept across the plasmasphere by the spin of the satellite. EUVs spatial resolution is 0.6° or 0.1 R _E in the equatorial plane seen from apogee. The sensitivity is 1.9 count s^–1 Rayleigh^–1, sufficient to map the position of the plasmapause with a time resolution of 10 min.     

Extreme Ultraviolet Imager Observations of the Structure and Dynamics of the Plasmasphere

The IMAGE Extreme Ultraviolet Imager (EUV) provides our first global images of the plasmasphere by imaging the distribution of He⁺ in its 30.4-nm resonance line. The images reveal the details of a highly structured and dynamic entity. Comparing EUV images and selected in-situ observations has helped to validate the remote sensing measurements. The brightness in the EUV images is heavily weighted by the He⁺ density near the plane of the magnetic equator, but two lines of evidence emphasize that the features seen by EUV extend far from the equator, and in at least some cases reach the ionosphere. Certain features and behaviors, including shoulders, channels, notches, and plasma erosion events, appear frequently in the EUV images. These are keys to understanding the ways that electric fields in the inner magnetosphere affect the large and meso-scale distribution of plasma, and their study can elucidate the mechanisms by which the solar wind and interplanetary magnetic field couple to the inner magnetosphere.     

The First two Years of Image

The Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) is the first satellite mission that is dedicated to imaging the Earth's magnetosphere. Using advanced multispectral imaging techniques along with omnidirectional radio sounding, IMAGE has provided the first glimpses into the global structure and behavior of plasmas in the inner magnetosphere. Scientific results from the two-year prime mission include the confirmation of the theory of plasmaspheric tails and the discovery of several new and unpredicted features of the plasmasphere. Neutral-atom imaging has shown how the ring current develops during magnetic storms and how ionospheric ions are injected into the ring current during substorms. The first global imaging of proton auroras has allowed the identification of the ionospheric footprint of the polar cusp and its response to changes in the interplanetary magnetic field. Detached subauroral proton arcs have been found to appear in the afternoon sector following south-north and east-west rotations of the IMF. Low-energy neutral atom imaging has shown global-scale ionospheric outflow to be an immediate response to solar-wind pressure pulses. Such imaging has also provided the first measurements of solar wind and interstellar neutral atoms from inside the magnetosphere. Radio sounding has revealed the internal structure of the plasmasphere and identified plasma cavities as the source of kilometric continuum radiation. These and numerous other scientific results now set the stage for the extended mission of IMAGE in which the imaging perspective will change markedly owing to orbital evolution while the magnetospheric environment undergoes a transition from solar maximum toward solar minimum.     

IMAGE mission overview

The Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission is the first mission in NASA's MIDEX (Mid-size Explorer) program. It is the first satellite mission that is dedicated to imaging the Earth's magnetosphere. IMAGE will utilize the techniques of ultraviolet imaging, neutral atom imaging, and radio plasma imaging to map out global distributions of the electron and proton aurora; the helium ions of the plasmasphere; the ionospheric ion outflow; the medium-energy ions of the near-Earth plasma sheet, ring current, and polar cusp; the high-energy ions of the ring current and trapped radiation belts; and the total plasma density from the ionosphere out to the magnetopause. The imaging perspective is from an elliptical polar orbit with apogee at latitudes from 40° to 90° in the northern hemisphere. For ultraviolet and neutral atom imaging, the time resolution is set by the two-minute spin period of the IMAGE spacecraft, which will be sufficient to track the development of magnetospheric substorms. An important feature of the IMAGE mission is its completely open data set with no proprietary data or intervals. All data, along with software needed for plotting and analysis, will be available within 24 hours of acquisition.     

Investigations of irregularities in remote plasma regions by radio sounding: applications of the radio plasma imager on image

The Radio Plasma Imager (RPI) on the IMAGE mission operates like a radar by transmitting and receiving coherent electromagnetic pulses. The RPI is designed to receive mirror-like (specular) reflections and coherent scatter returns. Long-range echoes of electromagnetic sounder waves are reflected at remote plasma cutoffs. Thus, analyses of RPI observations will yield the plasma parameters and distances to the remote reflection points. The RPI will employ pulse compression and spectral integration techniques, perfected in ground-based ionospheric digital sounders, in order to enhance the signal-to-noise ratio in long-range magnetospheric sounding. When plasma irregularities exist in the remote magnetospheric plasmas being probed by the sounder waves, echo signatures may become complicated. Experience in ionospheric sounding under such conditions indicates that sounding echo strengths can actually be enhanced by the presence of irregularities, and ground-based sounding indicates that coherent detection techniques can still be employed. In this paper we investigate the conditions that will allow coherent signals to be detected by the RPI and the signatures of scattering to be expected in the presence of multi-scale irregularities. Sounding of irregular plasma structures in the plasmasphere, plasmapause and magnetopause are also discussed.     

Outer Plasmaspheric Plasma Properties: What We Know from Satellite Data

Our understanding of the composition, density, and temperature structure of the outer plasmasphere has undergone several major revisions since the discovery of the plasmasphere over 35 years ago. Each new era of understanding was brought about by the launching of satellites carrying suites of plasma diagnostic experiments each more sophisticated than its predecessor. Instead of each new generation of missions providing incremental advancement in our knowledge of the plasmasphere's plasma properties, they most often than not, initiated a complete re-thinking of our understanding. There are still a number of observational "puzzles" and large gaps in our knowledge of how plasmaspheric structure and plasma properties change with geomagnetic storm and substorm activity. As we develop new missions to attempt to globally image the outer plasmasphere in order to put simultaneous multi-spacecraft observations into context, we undoubtedly will uncover new tantalizing observations that will cause us to re-think our mental picture of the plasmasphere. What we have learned from satellite data about the chemical composition, density, temperature and pitch-angle distribution properties and their variability is outlined in this review paper.     

Radio plasma imager simulations and measurements

The Radio Plasma Imager (RPI) will be the first-of-its kind instrument designed to use radio wave sounding techniques to perform repetitive remote sensing measurements of electron number density (N _e) structures and the dynamics of the magnetosphere and plasmasphere. RPI will fly on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission to be launched early in the year 2000. The design of the RPI is based on recent advances in radio transmitter and receiver design and modern digital processing techniques perfected for ground-based ionospheric sounding over the last two decades. Free-space electromagnetic waves transmitted by the RPI located in the low-density magnetospheric cavity will be reflected at distant plasma cutoffs. The location and characteristics of the plasma at those remote reflection points can then be derived from measurements of the echo amplitude, phase, delay time, frequency, polarization, Doppler shift, and echo direction. The 500 m tip-to-tip X and Y (spin plane) antennas and 20 m Z axis antenna on RPI will be used to measures echoes coming from distances of several R _E. RPI will operate at frequencies between 3 kHz to 3 MHz and will provide quantitative N _e values from 10^–1 to 10⁵ cm^–3. Ray tracing calculations, combined with specific radio imager instrument characteristics, enables simulations of RPI measurements. These simulations have been performed throughout an IMAGE orbit and under different model magnetospheric conditions. They dramatically show that radio sounding can be used quite successfully to measure a wealth of magnetospheric phenomena such as magnetopause boundary motions and plasmapause dynamics. The radio imaging technique will provide a truly exciting opportunity to study global magnetospheric dynamics in a way that was never before possible.     

Extraction of ion distributions from magnetospheric ENA and EUV images

Energetic neutral atom (ENA) and extreme ultra-violet photon (EUV) imagers will soon be probing magnetospheric ion distributions from the NASA space missions IMAGE and TWINS. Although ENA and EUV images will differ greatly, the same basic mathematical approach can be applied to deducing the ion distributions: extracting the parameters of a model ion distribution in a model magnetic field (and, in the case of ENA, interacting with a model cold neutral population). The model ion distribution is highly non-linear in its many parameters (as many as 38 have been used) in order to describe the strong spatial gradients of ion intensities in the magnetosphere. We have developed several new computer algorithms to accomplish the extraction by minimizing the differences between a simulated (instrument-specific) image and an observed image (or set of images). Towards the goal of a truly automated `hands-off extraction algorithm, we have combined three algorithms into a Hierarchical Simplex Algorithm. At each step of the minimization, it first tries a sophisticated and efficient Adaptive Conjugate Gradient algorithm. Then, if the error function is not reduced, it defers to an intermediate Analytic Simplex algorithm, and (if this step also fails) it finally defaults to the robust but inefficient Downhill Simplex algorithm. Whenever a step is successful, the algorithm begins the next step at the top of the hierarchy. We also offer a completely different approach (without minimization) for the interpretation of sharp `edges in the images (e.g., the plasmapause in He ⁺ 30.4 nm EUV images of the plasmasphere). We demonstrate mathematically that the equatorial shape of the plasmapause can be constructed directly from the image using a simple graphical algorithm.     

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.     

Ring Currents and Internal Plasma Sources

The discovery of terrestrial O⁺ and other heavy ions in magnetospheric hot plasmas, combined with the association of energetic ionospheric outflows with geomagnetic activity, led to the conclusion that increasing geomagnetic activity is responsible for filling the magnetosphere with ionospheric plasma. Recently it has been discovered that a major source of ionospheric heavy ion plasma outflow is responsive to the earliest impact of coronal mass ejecta upon the dayside ionosphere. Thus a large increase in ionospheric outflows begins promptly during the initial phase of geomagnetic storms, and is already present during the main phase development of such storms. We hypothesize that enhancement of the internal source of plasma actually supports the transition from substorm enhancements of aurora to storm-time ring current development in the inner magnetosphere. Other planets known to have ring current-like plasmas also have substantial internal sources of plasma, notably Jupiter and Saturn. One planet having a small magnetosphere, but very little internal source of plasma, is Mercury. Observations suggest that Mercury has substorms, but are ambiguous with regard to the possibility of magnetic storms of the planet. The Messenger mission to Mercury should provide an interesting test of our hypothesis. Mercury should support at most a modest ring current if its internal plasma source is as small as is currently believed. If substantiated, this hypothesis would support a general conclusion that the magnetospheric inflationary response is a characteristic of magnetospheres with substantial internal plasma sources. We quantitatively define this hypothesis and pose it as a problem in comparative magnetospheres.     

Global Imaging of Proton and Electron Aurorae in the far Ultraviolet

The IMAGE spacecraft carries three FUV photon imagers, the Wideband Imaging Camera (WIC) and two channels, SI-12 and SI-13, of the Spectrographic Imager. These provide simultaneous global images, which can be interpreted in terms of the precipitating particle types (protons and electrons) and their energies. IMAGE FUV is the first space-borne global imager that can provide instantaneous global images of the proton precipitation. At times a bright auroral spot, rich in proton precipitation, is observed on the dayside, several degrees poleward of the auroral zone. The spot was identified as the footprint of the merging region of the cusp that is located on lobe field lines when IMF B_z was northward. This identification was based on compelling statistical evidence showing that the appearance and location of the spot is consistent with the IMF B_z and B_y directions. The intensity of the spot is well correlated with the solar wind dynamic pressure and it was found that the direct entry of solar wind particles could account for the intensity of the observed spot without the need for any additional acceleration. Another discovery was the observation of dayside sub-auroral proton arcs. These arcs were observed in the midday to afternoon MLT sector. Conjugate satellite observations showed that these arcs were generated by pure proton precipitation. Nightside auroras and their relationship to substorm phases were studied through single case studies and in a superimposed epoch analysis. It was found that generally there is substantial proton precipitation prior to substorms and the proton intensity only doubles at substorm onset while the electron auroral brightness increases on average by a factor of 5 and sometimes by as much as a factor of 10. Substorm onset occurs in the central region of the pre-existing proton precipitation. Assuming that nightside protons are precipitating from a quasi-stable ring current at its outer regions where the field lines are distorted by neutral sheet currents we can associate the onset location with this region of closed but distorted field lines relatively close to the earth. Our results also show that protons are present in the initial poleward substorm expansion however later they are over taken by the electrons. We also find that the intensity of the substorms as quantified by the intensity of the post onset electron precipitation is correlated with the intensity of the proton precipitation prior to the substorms, highlighting the role of the pre-existing near earth plasma in the production of the next substorm.     

The IMAGE science and mission operations center

The Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) will produce forefront science by quantifying the response of the magnetosphere to the time variable solar wind. It will acquire, for the first time, a variety of three-dimensional images of magnetospheric boundaries and plasma distributions extending from the magnetopause to the inner plasmasphere. The images will be produced on time scales needed to answer important questions about the interactions of the solar wind and the magnetosphere. The IMAGE team will provide open access to all IMAGE data. Thus there will be no proprietary rights or periods. All IMAGE data products will be archived and available to the scientific research community. The IMAGE mission will operate with a near 100% duty cycle with all instruments in their baseline operational modes. A Science and Mission Operations Control Center or SMOC has been developed at the NASA Goddard Space Flight Center (GSFC) to be the main data and command processing system for IMAGE. The IMAGE Level-0 data will be processed into Level 0.5 and Level-1 data and browse products within 24 hours after their receipt of raw data in the SMOC. These data products will be transferred to the NSSDC, for long-term archiving, and posted immediately on the world-wide-web for use by the international scientific community and the public.     

Polarity Reversals from Paleomagnetic Observations and Numerical Dynamo Simulations

Recent advances in the study of geomagnetic field reversals are reviewed. These include studies of the transitional field during the last geomagnetic reversal and the last geomagnetic excursion based on paleomagnetic observations, and analysis of reversals in self-consistent 3D numerical dynamo simulations. Field models inferred from observations estimate reversal duration in the range of 1–10 kyr (depending on site location). The transitional fields during both the Matuyama/Brunhes reversal and the Laschamp excursion are characterized by low-latitude reversed flux formation and subsequent poleward migration. During both events the dipole as well as the non-dipole field energies decrease. However, while the non-dipole energy dominates the dipole energy for a period of 2 kyr in the reversal, the non-dipole energy merely exceeds the dipole energy for a very brief period during the excursion. Numerical dynamo simulations show that stronger convection, slower rotation, and lower electrical conductivity provide more favorable conditions for reversals. A non-dimensional number that depends on the typical length scale of the flow and represents the relative importance of inertial effects, termed the local Rossby number, seems to determine whether a dynamo will reverse or not. Stable polarity periods in numerical dynamos may last about 1 Myr, whereas reversals may last about 10 kyr. Numerical dynamo reversals often involve prolonged dipole collapse followed by shorter directional instability of the dipole axis, with advective processes governing the field variation. Magnetic upwellings from the equatorial inner-core boundary that produce reversed flux patches at low-latitudes of the core-mantle boundary could be significant in triggering reversals. Inferences from the observational and modeling sides are compared. We summarize with an outlook on some open questions and future prospects.     

Interplanetary studies: Propagation of disturbances between the Sun and the magnetosphere

This review is concerned with the interplanetary ‘transmission line’ between the Sun and the Earth's magnetosphere. It starts with comments about coronal mass ejections (CMEs) that are associated with various forms of solar activities. It then continues with some of the current views about their continuation through the heliosphere to Earth and elsewhere. The evolution of energy, mass, and momentum transfer is of prime interest since the temporal/spatial/magnitude behavior of the interplanetary electric field and transient solar wind dynamic pressure is relevant to the magnetospheric response (the presence or absence of geomagnetic storms and substorms) at Earth. Energetec particle flux predictions are discussed in the context of solar activity (flares, prominence eruptions) at various positions on the solar disk relative to Earth's central meridian. A number of multi-dimensional magnetohydrodynamic (MHD) models, applied to the solar, near-Sun, and interplanetary portions of the ‘transmission line’, are discussed. These model simulations, necessary to advancing our understanding beyond the phenomenological or morphological stages, are directed to deceptively simple questions such as the following: can one-to-one associations be made between specific forms of solar activity and magnetosphere response?