Effects on the plasmasphere of a time-varying convection electric field

The time evolution of the plasmasphere has been investigated theoretically, using simple computational models. The magnetic field is assumed to be dipolar and time-independent, but the convection electric field is allowed to vary in time. For purposes of comparison, various spatial distributions of the magnetospheric electric field are considered. Plasmasphere flux tubes are assumed to be filled by diffusion of plasma upwards from the dayside ionosphere. Following a reduction in the convection field, the bulge of the original plasmasphere develops into a long tail that gradually wraps itself around the main plasmasphere. Periodic gusts in a spatially uniform convection field produce extremely complicated fine structure that depends strongly on both local time and universal time. Each large gust produces a distinct tail of cold plasma that stretches from the main body of the plasmasphere to the magnetopause, and causes a peak in density, outside the main plasmapause; similar features have been observed by OGO satellites. The calculations indicate that a periodic gusty field has a major effect on the size of the plasmasphere if the field has large Fourier components close to the drift period of cold plasma near the plasmapause. Gusts occurring randomly, at an average rate of several a day, can also cause substantial reduction in the size of the plasmasphere. The assumption that the convection field is spatially uniform, but gusty, leads to better agreement with the observed average shape of the plasmasphere than the assumption of a constant, uniform electric field. The theory indicates that the thickness of the plasmasphere boundary should be inversely correlated with magnetic activity, in general agreement with OGO 5 observations.     

Effects on the plasmasphere of irregular electric fields

A conservative convection electric field model developed by Volland (1973) to describe the solar wind induced plasma flow within the inner magnetosphere is modified to include a noisy spatial component. Under steady state conditions such a random component will result in spatial irregularities in the thermal plasma density distribution in the vicinity of the plasmapause—particularly near dusk. Spatial irregularities in the convection can produce longitudinally restricted perturbations near the plasmapause some of which are detached from the main body of the plasmasphere. Temporal variations in the midnight to noon flow intensity are shown to produce elongated extensions of the plasmasphere known as plasmatails but even short period variations of the overall magnitude of the convection cannot produce longitudinally localized perturbations in the thermal plasma distribution. Convection models based on the 3 hr magnetic index K_p yield plasmasphere structures which are qualitatively similar to those based on shorter period variations, but the exact location at any given time of the plasmapause is dependent upon the characteristic time scale employed.     

Thermal proton flow in the plasmasphere: The morning sector

Vertical profiles of electron density obtained in the vicinity of the plasmapause using the Alouette-II topside sounder have been analyzed to assess the presence of H⁺ flow in the topside ionosphere. The observations in the midnight sector show clearly the presence of the plasmapause; i.e. there is a sharp boundary separating the poleward regions of polar wind H⁺ flow and the more gentle conditions of the plasmasphere where light ions are present in abundance. In contrast, in the sunlit morning sector upwards H⁺ flow is deduced to be present to invariant latitudes as low as 48° (L = 2·2) in the regions normally known to be well inside the plasmasphere. The upwards H⁺ flux is sufficiently large (3 × 10⁸ ions cm^?2 sec^?1) that the plasmapause cannot be seen in the latitudinal electron density contours of the topside ionosphere. The cause for this flow remains unknown but it may be a result of a diurnal refilling process.     

Effects of convection electric field on the thermal plasma flow between the ionosphere and the protonosphere

Dynamic behavior of the coupled ionosphere-protonosphere system in the magnetospheric convection electric field has been theoretically studied for two plasmasphere models. In the first model, it is assumed that the whole plasmasphere is in equilibrium with the underlying ionosphere in a diurnal average sense. The result for this model shows that the plasma flow between the ionosphere and the protonosphere is strongly affected by the convection electric field as a result of changes in the volume of magnetic flux tubes associated with the convective cross-L motion. Since the convection electric field is assumed to be directed from dawn to dusk, magnetic flux tubes expand on the dusk side and contract on the dawn side when rotating around the earth. The expansion of magnetic flux tubes on the dusk side causes the enhancement of the upward H⁺ flow, whereas the contraction on the dawn side causes the enhancement of the downward H⁺ flow. Consequently, the H⁺ density decreases on the dusk side and increases on the dawn side. It is also found that significant latitudinal variations in the ionospheric structures result from the L-dependency of these effects. In particular, the H⁺ density at 1000 km level becomes very low in the region of the plasmasphere bulge on the dusk side. In the second model, it is assumed that the outer portion of the plasmasphere is in the recovery state after depletions during geomagnetically disturbed periods. The result for this model shows that the upward H⁺ flux increases with latitude and consequently the H⁺ density decreases with latitude in the region of the outer plasmasphere. In summary, the present theoretical study provides a basis for comparison between the equatorial plasmapause and the trough features in the topside ionosphere.     

Thermal protons in the morning magnetosphere: Filling and heating near the equatorial plasmapause

Thermal H⁺ distributions have been measured as the European Space Agency GEOS-1 satellite passed through the late morning equatorial magnetosphere, plasmapause and plasmasphere. The unique capabilities of the on-board Supralhermal Plasma Analysers (SPA) have been used to overcome the retarding floating potential of the satellite and measure the velocity distribution of the cold protons. In the magnetosphere an enhanced source cone of such ions with a temperature of ~ 0.5 eV is a signature of the filling process occurring outside the plasmapause where flux tubes are relatively empty. In the plasmasphere the thermal H⁺ is essentially isotropic with a temperature less than 0.5 eV but the motion of the satellite introduces apparent drift.These measurements of cold proton velocity distribution now permit a reappraisal of the definition of the “plasmapause”. It becomes inappropriate to use an arbitrary empirical density, e.g. the conventional 10 cm ^?3, in order to establish a boundary. It is now possible to identify a plasmapause interaction region where the two cold proton populations co-exist. This region generally lies Earthward of the 10 cm ^?3 density level, has a width which is strongly dependent on magnetic activity and the temperature is typically between 0.5 and 1.5 eV. The change from “filled” to “unfilled” flux tubes relates to the physical processes which are occurring and the controlling electric field configuration; in particular, the last closed equipotential. Throughout this region, in going from the plasmasphere to the magnetosphere, the plasma drift motion is expected to change from corotation to a convection which is controlled by E ×B, and is predominantly Sunward due to the dawn-dusk electric field. Crossing the plasmapause on the morning side, little change in drift direction should occur but subtle variations in the ionic velocity distribution do reflect the change in the degree of flux tube density equilibrium.Our first direct measurement of the magnetospheric E × B drift has been reported previously but here measurements from a selected six day period show how the plasma in the plasmapause region responds to changing magnetospheric activity. The drift velocities cannot he derived with high accuracy but the analysis shows that the technique can provide a valid mapping of the magnelospheric electric field. In addition, since the magnetospheric cold plasma distribution is observed after it has come from the ionosphere, a distance of many Earth radii, the scattering and accelerating mechanisms along the flux tube can be studied. For this particular data-set taken in the late morning, the maximum potential drops along the flux tubes were less than a volt. The ionospheric proton source cone is observed to be broad, pitch angle scattering persists up to 40 or even 70°.Although these results throw new light on the plasmaspheric filling process one must recognise that, however the plasmapause is defined, it is not a simple matter to map this boundary from the equatorial plane down to low altitudes and the mid-latitude trough.     

Ionospheric and magnetospheric “plasmapauses”

During August 1972, Explorer 45 orbiting near the equatorial plane with an apogee of ～5.2 R_e traversed magnetic field lines in close proximity to those simultaneously traversed by the topside ionospheric satellite ISIS 2 near dusk in the L range 2.0–5.4. The locations of the Explorer 45 plasmapause crossings (determined by the saturation of the d.c. electric field double probe) during this month were compared to the latitudinal decreases of the H⁺ density observed on ISIS 2 (by the magnetic ion mass spectrometer) near the same magnetic field lines. The equatorially determined plasmapause field lines typically passed through or poleward of the minimum of the ionospheric light ion trough, with coincident satellite passes occurring for which the L separation between the plasmapause and trough field lines was between 1 and 2. Hence, the abruptly decreasing H⁺ density on the low latitude side of the ionospheric trough is not a near earth signature of the equatorial plasmapause. Vertical flows of the H⁺ ions in the light ion trough as detected by the magnetic ion mass spectrometer on ISIS were directed upward with velocities between 1 and 2 km s^?1 near dusk on these passes. These velocities decreased to lower values on the low latitude side of the H⁺ trough but did not show any noticeable change across the field lines corresponding to the magnetospheric plasmapause. The existence of upward accelerated H⁺ flows to possibly supersonic speeds during the refilling of magnetic flux tubes in the outer plasmasphere could produce an equatorial plasmapause whose field lines map into the ionosphere at latitudes which are poleward of the H⁺ density decrease.     

Comparison of solutions to the thirteen-moment and standard transport equations for low speed thermal proton flows

We have compared solutions obtained from the general 13-moment system of transport equations with those obtained from the standard collision-dominated transport equations for conditions corresponding to low speed thermal proton flow in the topside ionosphere in the vicinity of the plasmapause. In general, the solutions obtained from the 13-moment system of equations, which allows for different species temperatures parallel and perpendicular to the geomagnetic field and non-classical heat flows, are different from those obtained from the standard transport equations, which account for isotropic temperatures and classical collision-dominated heat flows. Within the plasmasphere, where the electron density is high, the differences between the 13-moment and standard solutions are typically small. However, outside the plasmasphere where the electron density is lower and in the ionosphere above SAR-arcs, where substantial electron and proton heat flows occur, there can be significant differences between the 13-moment and standard solutions. Generally, the differences are much larger for the protons than for the electrons. Our 13-moment solutions indicate that the proton and electron distributions are anisotropic with the difference between parallel and perpendicular temperatures approaching 4000 K for the protons and 2500 K for the electrons in the ionosphere above SAR-arcs. Also, above SAR-arcs the 13-moment heat flow equations yield proton heat flows as much as a factor of 10 lower and electron heat flows as much as a factor of 2 lower than those predicted by the classical collision-dominated heat flow expressions for the same boundary conditions.     

Correlations between ground observations of Pi 2 geomagnetic pulsations and satellite plasma density observations

Ground observations of Pi 2 geomagnetic pulsations are correlated with satellite measurements of plasma density for three time intervals. The pulsations were recorded using the IGS network of magnetometer stations and the plasma density measurements were made on board GEOS-1 and ISEE-1. Using the technique of complex demodulation, the amplitude, phase and polarisation characteristics of the Pi 2 pulsations are observed along two meridional profiles; one from Eidar, Iceland (L = 6.7) to Cambridge, U.K. (L = 2.5) and the other from Tromso, Norway (tL = 6.2) to Nurmijarvi, Finland (L = 3.3). The observed characteristics of the Pi 2 pulsations are then compared with the plasma density measurements. Close relationships between the plasmapause position and the position of an ellipticity reversal and a variation in H component phase are observed. A small, secondary amplitude maximum is observed on the U.K./Iceland meridian well inside the position of the projection of the equatorial plasmapause. The primary maxima on the two meridians, in general occur close to the estimated position of the equatorward edge of a westward electrojet. Using the plasma density measurements, the periods of surface waves at the plasmapause for two intervals are estimated and found to be in good agreement with the dominant spectral peaks observed at the ground stations near the plasmapause latitude and within the plasmasphere. The polarisation reversal, together with phase characteristics, spectral evidence and the agreement between the theoretical and observed periods leads to the suggestion that on occasions a surface wave is excited on the plasmapause as an intermediate stage in the propagation of Pi 2 pulsations from the auroral zone to lower latitudes.     

Dynamical interpretation of observed plasmasphere deformations

Based on all of the OGO-5 light ion density measurements (covering the period from March, 1968 to May, 1969), a definition of “isolated plasma regions” was employed to locate the most prominent patches of enhanced light ion densities in the midst of the depleted region, outside of the main plasmasphere. On the dayside, the distribution of these isolated plasma in L.T. vs. L coordinates was quite similar to that of the “detached plasma regions” by Chappell (1974a). On the nightside, however, the new distribution revealed more frequent occurrence of these regions. Elongated thick plasmatails produced during periods of sudden enhancement of convection electric fields and subsequentially thinning and corotating of the plasmatails during quieting periods, in general, could account for the statistical distribution as well as the individual events, such as those between March 27 and April 2, 1968 and Oct. 21 and Oct. 24, 1968. As demonstrated by Kivelson (1976), wave-particle interactions could produce tremendously complicated structures observed in the near vicinity of the plasmapause and far away from the plasmasphere. Examination of H⁺ and He⁺ density measurements for period of Aug. 12–Aug. 20, 1968 indicated that the density reduction of the plasmasphere during a magnetic storm was on the same order of magnitude as that obtained from whistler techniques during a magnetospheric substorm.     

The role of ground arrays of magnetometers in the study of pulsation resonance regions

The last ten years have seen significant progress in the understanding of magnetic pulsations. Results from ground arrays of magnetometers were instrumental in the establishment of the concept of field line resonance, which has proved essential to the ordering of observed pulsation phenomena. Ground arrays have continued to be indispensable in the study of pulsation resonance regions. In this paper we will review the major contributions that ground arrays have made in the areas of field line resonance inside and outside the plasmasphere, the role of the plasmapause, azimuthal propagation, the effect of the ionosphere and correlation with satellite observations. Emphasis will be on recent developments and we will also make some suggestions for possible array work in the future.     

Hydromagnetic wave coupling in the magnetosphere—Plasmapause effects on impulse-excited resonances

Recent analytical and numerical modelling has demonstrated the possibility that impulsively-stimulated compressional hydromagnetic cavity resonances can drive local field-line resonances in the magnetosphere. This paper extends the modelling to include axisymmetric plasmapause structures with realistic radial variation in the magnetospheric cavity. The results show that: (a) the plasmapause plays an important rôle in determining which cavity resonances are dominant; (b) when the wave fields are significantly non-axisymmetric, additional cavity resonances are evident which are at least partly trapped within the plasmasphere; (c) the position of the plasmapause determines where (and whether) cavity resonances couple significantly to field-line resonances; (d) for the small “azimuthal” wavenumber chosen, there is no evidence of a compressional surface wave on the plasmapause.     

A high resolution study of continuous pulsations in the European sector

Complex demodulation has been described in detail and applied to Pi2 pulsations in a previous paper by Beamish et al. (1979). The technique is now extended to demonstrate spatiotemporal variations in the fundamental characteristics of Pc3 and Pc4 pulsations along a meridional profile extending from the U.K. to Iceland. With the exception of a high latitude Pc4 coupled resonance the results are consistent with a ?90° Hughes rotation (introduced by the ionosphere) of magnetospheric toroidal line resonances. Furthermore, the ionosphere appears capable of smoothing away the polarisation reversal which would be expected across such amplitude maxima within the plasmasphere. However, a toroidal line resonance in the Pc3 period range about which a sense of polarisation reversal is clearly observed on the ground is suggested as occurring at the plasmapause. This is accounted for in terms of the width of the resonance structure.     

The ‘Roche-Limit’ of ionospheric plasma and the formation of the plasmapause

The plasmapause position is determined by the innermost equipotential surface which is tangent to the ‘Roche-Limit’ surface of the ionospheric plasma filling the magnetosphere. When the thermal particles corotate with the Earth's angular velocity, the ‘Roche-Limit’ equatorial distance is L_c=5.78 [R_E]. When the angular convection velocity is evaluated from the quiet time electric field distribution E3 of McIIwain (1972), L_c depends on the local time. Its minimum value is then L_C=4.5Near 2400 LT, and the plasmapause shape and position satisfactorily fit the observations. The diffusive equilibrium dnesity distribution appropriated inside the plasmasphere, becomes convectively unstable beyond L = L_c, where the collisions type of model satisfactorily represents the observations. In the intermediate region between the plasmapause and the last closed magnetic field line, contimues ionization fluxes are expected to flow out of the midlatitude ionosphere     

Monitoring the Dynamics of the Ionosphere–Plasmasphere System by Ground-Based ULF Wave Observations

Cross-spectral analysis of ULF wave measurements recorded at ground magnetometer stations closely spaced in latitude allows accurate determinations of magnetospheric field line resonance (FLR) frequencies. This is a useful tool for remote sensing temporal and spatial variations of the magnetospheric plasma mass density. The spatial configuration of the South European GeoMagnetic Array (SEGMA, 1.56 <  L <  1.89) offers the possibility to perform such studies at low latitudes allowing to monitor the dynamical coupling between the ionosphere and the inner plasmasphere. As an example of this capability we present the results of a cross-correlation analysis between FLR frequencies and solar EUV irradiance (as monitored by the 10.7-cm solar radio flux F10.7) suggesting that changes in the inner plasmasphere density follow the short-term (27-day) variations of the solar irradiance with a time delay of 1–2 days. As an additional example we present the results of a comparative analysis of FLR measurements, ionospheric vertical soundings and vertical TEC measurements during the development of a geomagnetic storm.     

Magnetospheric plasma discontinuity by the geostationary version of the differential phase method during storms

It is shown that the dynamics of the plasmapause, the plasmasphere plasma tails, the plasma sheet and the magnetosheath boundaries of the geomagnetosphere may be investigated by means of the geostationary version of the differential phase method, by which a signal transmitted from a sounding station (a geostationary satellite) and received by a response station on the Earth may be transformed, allowing the sign of the frequency shift and of the phase lag to be changed. Information on the location, the motion of the magnetospheric plasma discontinuities and the concentration drop at their boundaries may be obtained from measurements carried out on board the geostationary satellite of the phase difference of the sounding and response signals ΔΦ, the time of its increase Δt and the phase difference change rate (fast beating frequency Δƒ = ΔΦ/2π Δt). The establishment of communication between appropriately spaced ground stations and a satellite with a quasi-polar orbit allows the midlatitude plasmapause dynamics, and those of the ionosphere trough, polar cusp boundaries and of polar cap inhomogeneities to be studied. Equipment with a stability of 10⁻¹¹–10⁻¹² is needed for the most dynamical events (for ΔΦ= 10⁻⁴ tens of rad. and for Δƒ= 10⁻⁵ tens of Hz) occurring in the radio path during storms.     

Dynamical behavior of thermal protons in the mid-latitude ionosphere and magnetosphere

Satellite and other observations have shown that H⁺ densities in the mid-latitude topside ionosphere are greatly reduced during magnetic storms when the plasmapause and magnetic field convection move to relatively low L-values. In the recovery phase of the magnetic storm the convection region moves to higher L-values and replenishment of H⁺ in the empty magnetospheric field tubes begins. The upwards flow of H⁺, which arises from O⁺—H charge exchange, is initially supersonic. However, as the field tubes fill with plasma, a shock front moves downwards towards the ionosphere, eventually converting the upwards flow to subsonic speeds. The duration of this supersonic recovery depends strongly on the volume of the field tube; for example calculations indicate that for L = 5 the time is approximately 22 hours. The subsonic flow continues until diffusive equilibrium is reached or a new magnetic storm begins. Calculations of the density and flux profiles expected during the subsonic phase of the recovery show that diffusive equilibrium is still not reached after an elapsed time of 10 days and correspondingly there is still a net loss of plasma from the ionosphere to the magnetosphere at that time. This slow recovery of the H⁺ density and flux patterns, following magnetic storms, indicates that the mid-latitude topside ionosphere may be in a continual dynamic state if the storms occur sufficiently often.     

An ionosphere-plasmasphere global electron density model

Global Positioning System (GPS) satellites now provide a major source of data used in the study of the ionized regions that surround the Earth. With the greater use of GPS satellites in this role, understanding the effect the plasmasphere has on GPS radio signals has become increasingly important. The plasmasphere, with similar properties to the underlying ionosphere, can potentially contribute a significant fraction of the overall ionospheric error in radio signals transmitted from a GPS satellite to ground or space based receivers. A global plasmasphere-ionosphere electron density model has been developed to investigate this contribution. A variety of output examples that the model can produce are shown.     

Coincident observations of ionospheric troughs and the equatorial plasmapause

Explorer 45 traversed the plasmapause (determined approximately via the saturation of the d.c. electric field experiment) at near-equatorial latitudes on field lines which were crossed by Ariel 4 (~600km altitude) near dusk in May 1972 and on field lines which were crossed by Isis II (~1400km altitude) near midnight in December 1971 and January 1972. Many examples were found in which the field line through the near-equatorial plasmapause was traversed by Explorer 45 within one hour local time and one hour universal time of Ariel and Isis crossings of the same L coordinate. For the coincident passes near dusk, the RF electron density probe on Ariel detected electron density depletions near the plasmapause L coordinates when Ariel was in darkness. When the Ariel passes were in sunlight, however, electron depletions were not discernable near the plasmapause field line. On the selected near-midnight passes of Isis II, electron density depressions were typically detected (via the topside sounder) near the plasmapause L coordinate. The dusk Ariel electron density profiles are observed to reflect O⁺ density variations. Even at the high altitude of Isis near midnight, O⁺ is found to be the dominant ion in the trough region whereas H⁺ is dominant at lower latitudes as is evident from the measured electron density scale heights. In neither local time sector was it possible to single out a distinctive topside ionosphere feature as an indicator of the plasmapause field line as identified near the equator. At both local times the equator-determined plasmapause L coordinate showed a tendency to lay equatorward of the trough minimum.     

Standing slow magnetosonic waves in a dipole-like plasmasphere

A problem of the structure and spectrum of standing slow magnetosonic waves in a dipole plasmasphere is solved. Both an analytical (in WKB approximation) and numerical solutions are found to the problem, for a distribution of the plasma parameters typical of the Earth's plasmasphere. The solutions allow us to treat the total electronic content oscillations registered above Japan as oscillations of one of the first harmonics of standing slow magnetosonic waves. Near the ionosphere the main components of the field of registered standing SMS waves are the plasma oscillations along magnetic field lines, plasma concentration oscillation and the related oscillations of the gas-kinetic pressure. The velocity of the plasma oscillations increases dramatically near the ionospheric conductive layer, which should result in precipitation of the background plasma particles. This may be accompanied by ionospheric F2 region airglows modulated with the periods of standing slow magnetosonic waves.