Plasmasphere response to the onset of quiet magnetic conditions: Plasma convection patterns

A whistler study has been made of plasma convection within the plasmasphere during a transition from steady moderate geomagnetic activity to quiet conditions. Continuous whistler data recorded at Sanae, Antarctica (L= 3.98) for the period 0400 UT, 10 July to 0400 UT, 11 July 1973 have been analyzed in 15 min intervals.This study has revealed two distinct bulges in the plasmasphere centred on 1700 and 0100 UT. The bulges appear to result from the outward flow of plasma rather than the addition of new plasma. We tentatively interpret the late bulge at 0100 UT as being the duskside bulge of earlier studies rotated into the midnight region. In this bulge, plasma above L = 3.8 appears to convect outwards to form the bulge whereas plasma at lower L-values is relatively undisturbed. For the early bulge (1700 UT) the plasma convection pattern is similar over all observable L-values and closely reflects the shape of the estimated plasmapause in that region. Comparison of the bulges, with those obtained by Carpenter (1966) suggests that the onset of quiet conditions results in a general displacement of the bulges in an eastward direction by about 3 hr.     

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.     

Horizontal electric fields in the middle latitude

Three dimensional electric fields were measured at the altitude of about 27 km in the stratosphere over the Pacific Ocean about 200–400 km away from the Sanriku coast of Honsyu Island (L = 1·4) on 16–17 October 1973, which was magnetically disturbed. The average horizontal electric field thus measured is about 10 mV/m, and the electric field vectors made clockwise semidiurnal rotations rather than diurnal. Daily variation of this electric field was compared with data at L = 2·7–3·5 published by Mozer (1973) and was found to be very similar. This suggests that these electric fields are of common origin in the plasmasphere. From their mean daily variation it is estimated that the plasmaspheric convection is decreased in the night side and is increased in the day side by 200–300 m/sec, and there is an outward flow in the first half of the afternoon and an inward flow in the plasma bulge region of about 500 m/sec.     

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.     

Spectral studies of geomagnetic pulsations with periods between 20 and 120 sec,and their relationship to the plasmapause region

Digital spectrograms have been computed for 18 days of geomagnetic pulsation activity at three UK Earth current stations (L = 2.6?3.6).Three main conclusions are drawn: (1) There are days when the period of the dominant spectral amplitude is ordered according to the observatory latitude. The most frequently observed large amplitude spectral peaks are centred on 80, 60 and 45s for South Uist (L = 3.6), Eskdalemuir (L = 3.1) and East Anglia (L = 2.6). respectively. (2) There are other days when the period of the dominant spectral amplitude is the same at all the observatories. (3) When Pc 3 and 4 period waves have been detected together, the latitude dependence of the amplitudes supports the theory that the shorter period pulsation is enhanced in the plasmatrough while the longer period wave is enhanced within the plasmasphere.     

On the diurnal period variation of midlatitude ULF pulsations

Magnetometer studies of the periods of mid-latitude ULF pulsations have produced conflicting results on the variation of the pulsation periods with both latitude and local time. Since the mid-latitude geomagnetic field is not expected to be significantly distorted by the solar wind, the observed diurnal period variations should be determined by changes in the ambient plasma density. We have applied a physically realistic plasmasphere model to the determination of pulsation eigenperiods over a 24-h interval at L=2.3 (appropriate to Wellington, New Zealand). The resulting model pulsation eigenperiods are largest during the day, with minimum and maximum values at 05.00 and 18.00 L.T. respectively. The model predicts a general increase in the eigenperiods during the replenishment of the protonosphere after a period of geomagnetic activity.     

Statistical studies of geomagnetic pulsations with periods between 10 and 70 sec and their relationship to the plasmapause region

Geomagnetic pulsations, in the period range 10–150 sec, have been analysed from five stations; Eskdalemuir (L = 3.1), Lerwick (L = 4.0), St. Anthony (L = 4.9), Sodankyla (L = 5.3) and Tromsø (L = 6.6). The results of 12 observatory years' worth of data are presented in the form of contour maps showing the frequency of occurrence of the pulsations as a function of K_p index and of local time. The maps show that a ground based observatory is more likely to record shorter period oscillations (pc 3) when the geomagnetic field line linking the station with the southern hemisphere passes through the plasmatrough than when the observatory field line links the plasmasphere. The peak occurrence of pc 3 for the observatories considered is at 08:45 hr ± 1 hr LT and is related to the observatory L value and the average night-time K_p index by the equation, L = 8.1 ? 1.2K_p. At Eskdalemuir, the spectrum is broader band than the other stations and tends to divide into two peaks; the pc 3 (20 sec) peak tends to occur when the plasmapause has moved in close to the observatory; while the pc 4 (60 sec) peak occurs when the K_p values have been lower and the plasmapause is further away at higher 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.     

Meridional characteristics of a Pc 4 micropulsation event in the plasmasphere

An analysis was made of a complex large amplitude Pc 4 micropulsation, of four hours duration around local noon, observed at five ground stations in the United Kingdom (2.4? L ?3.8). The final pulsation waveform was shown to be the results of the superposition of wave packets of different periods. The meridional variation of the amplitude of the different period wavepackets was consistent with their being fundamental “toroidal” field line resonances within the plasmasphere, rotated through 90° in their transmission through the ionosphere in accordance with recent theoretical predictions. Other predicted ionospheric effects, such as the loss of the sense-of-polarization reversal across the amplitude maximum, were apparent in the meridional variation of the polarization characteristics.     

Electron content modelling: The significance of protonospheric contents

The principal advance of the ATS-6 satellite beacon experiment was the ability to deduce continuously the electron content along the entire slant path from ground-based measurements of the signal group delay. This feature has been exploited in conjunction with the more usual Faraday rotation technique to separate the total electron content into ionospheric and protonospheric components. The physical validity of the deduced quantities is investigated using a mathematical model of the plasmasphere in which integration of the time-dependent continuity and momentum equations for oxygen and hydrogen ions along selected L shells yields the ion concentrations and field-aligned fluxes. The ion concentrations are then integrated along the propagation path to various ground stations from ATS-6 to give computed values for comparison with observations. The mathematical model is used with different sets of atmospheric parameters to investigate the significance of ionospheric and protonospheric contents as derived from beacon data.The calculated electron concentrations are able to reproduce mid-latitude equinoctial electron content observations. The shape parameters τ and F can also be simulated by day, but night-time values do not match the observations well, a greater protonospheric content being required. The calculations show that the quantity N_p, which is readily derived from ATS-6 observations, may be interpreted as the slant H⁺ content above some fixed height in the case of some stations (but not others) if the plasmasphere is reasonably full. The total slant content of H⁺ is approx. twice the value of N_p, though it appears that for the Lancaster raypath a closer relationship exists between N_p and the H⁺ tube content at L = 1.8. In general,N_p is most closely related to the tube content for an L value slightly greater than the minimum L intersected along the raypath.     

Ring current protons in the upper atmosphere within the plasmasphere

keV protons observed by the ESRO 1A satellite in the upper atmosphere equatorward of the main precipitation zone are described and discussed. The protons are highly anisotropic (empty loss cone) between the low-latitude boundary of the main precipitation zone and approximately L=4 during quiet and moderately disturbed conditions (K_p=0?4). Between L=4 and L=2.7 the proton flux is generally enhanced compared to that at L values somewhat above 4 and only moderately anisotropic. Substorms push the outer main precipitation zone equatorwards, but the boundaries of the innermost, moderately anisotropic region (at L=2.7 and L=4) move only when strong magnetic storms compress the plasmasphere to within L=4. It is suggested that the moderately anisotropic zone is caused by the ion-cyclotron instability for which the growth rate may have a broad maximum between about L=2.7 and L=4. For proton energies in the keV range the instability is excited only in regions with cold plasma densities above several hundred ions per cubic centimeter. It is finally concluded that the observations of low-latitude proton precipitation lend further support to the mechanism of ion-cyclotron instability as the cause of proton pitch angle diffusion, as proposed by Cornwall et al. (1970).     

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.     

Diurnal variation of thermal plasma in the plasmasphere

All of the OGO-5 light ion density measurements (covering the period from March 1968 to May 1969) obtained from the Lockheed Light Ion Mass Spectrometer were used to determine the average global topology of the equatorial plasmasphere density distribution. The variation of the light ion equatorial density at L?3.2 with local time was deduced by determining the average density observed within one hour of a specific local time and within 0.1 of a given L coordinate. The average H⁺ density showed a semidiurnal variation with peaks near noon and midnight. The He⁺ observations also revealed multiple peaks throughout the day but with smaller amplitudes than those of H⁺. At L>3.2 plasma trough conditions increase the scatter of densities. The average variation of the H⁺ density with L within the plasmasphere is found to be steepest near midnight and can be least squares fitted equally well to either an exponential variation exp (?bL) where b is between 0.85 and 1.5 or to a power law L^?a where a varies from 3.2 to 5.     

Outlying plasmasphere structure detected by whistlers

Whistlers recorded at Eights ($\text{L ? 4}$) and Byrd ($\text{f ? 7}$), Antarctica have been used to study large-scale structure in equatorial plasma density at geocentric distances $\text{?3–6}\text{R}{}_{\mathrm{E}}$. The observations were made during conditions of magnetic quieting following moderate disturbance. The structures were detected by a “scanning” process involving relative motion, at about one tenth of the Earth's angular velocity or greater, between the observed density features and the observing whistler station or stations. Three case studies are described, from 26 March 1965, 11 May 1965 and 29 August 1966. The cases support satellite results by showing outlying high density regions at $\text{?4–6}\text{R}{}_{\mathrm{E}}$ that are separated from the main plasmasphere by trough-like depressions ranging in width from $\text{?0.2}$ to 1 R_E. The structures evidently endured for periods of 12 hr or more. In the cases of deepest quieting their slow east-west motions with respect to the Earth are probably of dynamo origin. The cases observed during deep quieting (11 May 1965 and 29 August 1966) suggest the approximate rotation with the Earth of structure formed during previous moderate disturbance activity in the dusk sector. The third case, from 26 March 1965, may represent a structure formed near local midnight. The reported structures appear to be closely related to the bulge phenomenon. The present work supports other experimental and theoretical evidence that the dusk sector is one of major importance in the generation of outlying density structure. It is inferred that irregularities of the type reported here regularly develop near 4–5 R_E during moderate substorm activity. This research suggests that at least a major class of the density structures that develop near 4 R_E are tail-like in nature, joined to the main body of the plasmasphere. The apparent disagreement with Chappell's results from OGO 5, which are interpreted as showing regions of “detached” plasma beyond 5 R_E, may be related to the pronounced spatial structure of electric fields observed in high-latitude ionospheric regions that are conjugate to the magnetospheric regions in which the OGO-5 observations were made.     

Temporal variations in the dawn and dusk midlatitude trough and plasmapause position

The temporal development of the latitudinal position of a 600 km midlatitude electron density trough at dawn and dusk during the period 25–27 May 1967, which encompassed a large magnetic storm, was measured by the RF capacitive probe on the polar orbiting Ariel 3 satellite. The substorm-related changes in the L coordinate of the trough minimum and the point of most rapid change of density gradient on the low latitude side of the trough are similar. Oscillations of the trough position at dusk are in phase with substorm activity whereas movement of the trough at dawn is only apparent with the onset of the large storm. Detailed model calculations of the plasmasphere dynamics assuming a spatially invariant equatorial convection E-field which varies in step with the K_p index produces a plasmapause motion which parallels the observed trough behaviour, particularly at dusk, and shows that the outer plasmasphere and possibly the trough region are characterized by complex fine structured variations due to the past history of the magnetosphere convection.     

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     

IN-SITU observations of irregular ionospheric structure associated with the plasmapause

Additional studies of the ion composition results obtained from the OGO-6 satellite support earlier observations of irregularities in the distribution of H⁺ and He⁺ within the light ion trough near L = 4, which has been associated with the plasmapause. These irregularities are in the form of sub-troughs superimposed upon the major mid latitude decrease of the light ions. In the sub-troughs, ionization depletions and recoveries of as much as an order of magnitude are observed within a few degrees of latitude, usually exhibited in a pattern which changes significantly with longitude as the Earth rotates beneath the relatively fixed satellite orbit. The location and properties exhibited by these sub-troughs appear to be consistent with the concept of a plasmasphere distortion in the form of “plasmatails” resulting from the combined effects of magnetospheric convection plus corotation. Like the light ion trough, the “plasmatail” irregularity in H⁺ may be obscured on the day side by the dominant topside distribution of O⁺. Consequently, these light ion irregularities are seen as an important factor for studies of plasmapause-trough relationships.     

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.     

Quiet-time convection electric field properties derived from keV electron measurements at the inner edge of the plasma sheet by means of GEOS-2

From an analysis of the local time distribution of the electron upper energy limit reached by the geostationary satellite GEOS-2 in cutting through the innermost part of the electron plasma sheet during fairly quiet conditions the following results have been obtained, among others. An electric field model given by $\text{E}\text{= ?▽{AR}{}^4\text{sin}\text{(φ+}\text{π}\text{4}\text{)}}$, with the dusk singular point of the forbidden region boundary at 1500, instead of at 1800 M.L.T., is in quite good agreement with the observations. This means that effects due to the shielding by the hot plasma of the inner magnetosphere from the convection electric field are quite strong in situations of low disturbance level. The quiet-time convection electric field strength at 2100 M.L.T. in the geostationary orbit obtained from this analysis varies in the range 0.15–0.3 kV/R_e. Six hours earlier or later in the satellite orbit the convection field is four times stronger. Also when the convection field varies, some information about its magnitude can be obtained from the keV electron measurements.     

Giant pulsations in the plasmasphere

It has been generally accepted up to now that giant pulsations (Pg) are auroral zone phenomena but here we present observations of a sequence of three Pg events on successive days at three stations well within the plasmasphere. Field line resonance behaviour is exhibited with one of the events clearly resonating at L ? 2.8. From the resonant frequency (10.4 mHz) equatorial mass densities are calculated and from these, and the measured azimuthai polarization at resonance, the inference is drawn that Pgs are oscillations in the fundamental guided poloidal mode. We suggest that the drift wave instability of the compressional Alfvén wave may be the source mechanism for Pgs and speculate how conditions for the instability may have arisen.