Determination of electric fields and neutral wind velocities in the ionospheric E- and F-regions from observations of artificial ion clouds

A new method of interpreting the behaviour of artificial ion clouds released in the Earth's ionosphere is presented. It is shown that values for the ionospheric electric field, neutral wind velocities and, in some circumstances, ion collision frequency, can be deduced from a study of the motion and deformation of the ion clouds, including those released in the E-region.     

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.     

Behaviour of ion velocity distributions for a simple collision model

For application to studies of the high latitude ionosphere, we have calculated ion velocity distributions for a weekly-ionized plasma subjected to crossed electric and magnetic fields. An exact solution to Boltzmann's equation has been obtained by replacing the Boltzmann collision integral with a simple relaxation model. At altitudes above about 150 km, where the ion collision frequency is much less than the ion cyclotron frequency, the ion distribution takes the shape of a torus in velocity space for electric fields greater than 40 mV m^?1. This shape persists for 1–2 hr after application of the electric field. At altitudes where the ion collision and cyclotron frequencies are approximately equal (about 120 km), the ion velocity distribution is shaped like a bean for large electric field strengths. This bean-shaped distribution persists throughout the lifetime of ionospheric electric fileds. These highly non-Maxwellian ion velocity distributions may have an appreciable affect on the interpretation of ion temperature measurements.     

The influence of convection electric fields on thermal proton outflow from the ionosphere

The continuity, momentum and energy hydrodynamic equations for an O⁺-H⁺ ionosphere have been solved self-consistently for steady state conditions when a perpendicular (convection) electric field is present. Comparison of the H⁺ temperature profiles obtained with and without the electric field show that the effect of the electric field is to enhance the H⁺ temperature at high altitudes from about 3600 to 6400 K. Due to ion heating by the electric field, there is a net reduction of O⁺ in the F2-region as compared with the case of a non-convecting ionosphere. When the reduction of O⁺ is neglected, the electric field acts to increase the H⁺ outward flux from 8.3 × 10⁷ to 2.7 × 10⁸ cm^?2 sec^?1 for average ionospheric conditions. However, when the reduction of O⁺ is included, there is a net reduction in the outward H⁺ flux. Nevertheless, the convection electric field still results in an increase in the rate of depletion of the F-re m^?1 electric field.     

On the divergence of the ionospheric electric field around the westward travelling surge

The pattern of the ionospheric electric field around the westward travelling surge (WTS) is theoretically studied. This is obtained by solving the current continuity equation at the ionospheric altitude for temporal and spatial development of the field-aligned current density modelled as the WTS phenomenon. The results show that the divergence of the ionospheric electric field is significantly changed depending on the dawn-to-dusk convection electric field E₀ because of non-uniformity in the ionospheric conductivity: the ionospheric electric field diverges in the upward current region (around the head of the WTS) when a westward electric field E₀ of 10 mV m⁻¹ is uniformly applied. On the other hand, the ionospheric electric field converges without E₀. From the observational inference that the ionospheric electric field converges around the head of the WTS, it is suggested that the WTS phenomenon may not be accounted for by the discharging process in the presence of the enhanced dawn-to-dusk convection electric field and non-uniform conductivity as was studied by previous authors.     

Comments on the proton whistler generation process

The problem of the propagation of an electromagnetic wave originating for instance in a lightning flash through the ionospheric medium is analysed in order to understand the formation at high ionospheric altitudes of the so-called proton whistler. It is shown that the accessibility of the hydrodynamic (or kinetic) proton resonance at the satellite altitude requires that a mode conversion process must take place slightly above the transition region separating the one ion (O⁺) from the two ion (O⁺ + H⁺) component plasmas. Moreover, the transformation conditions in the wave conversion region imply that the magnetic field should be (almost) perpendicular to the density gradient. Otherwise, the incident electromagnetic wave will never reach the satellite altitude in the frequency range of the proton whistler. However, some former proton whistler theories have postulated that the signal is the result of simple ionospheric propagation effects, in contradiction with the above results. These former proton whistler theories are reviewed and it is shown that the basic flaw in these theories lies in that the incident electromagnetic wave has been supposed from the beginning to have reached the high ionospheric altitudes where is located the satellite without being influenced by the lower ionospheric layers. Some various aspects, like the high variability of the wave electric to magnetic field ratio and the harmonics bands as observed by Injun are analysed in the light of the obtained results. Finally, numerical solutions of the wave dispersion relation for both the fast hydrodynamic mode (the extraordinary mode) and the slow ion kinetic mode are presented which shows that a coupling process between the two modes may take place at various frequencies between the O⁺ and the H⁺ gyrofrequencies.     

The source of midlatitude metallic ions at F-region altitudes

A survey of metallic ions detected by the Bennett Ion Mass Spectrometers flown on the Atmosphere Explorer satellites, including both circular and eccentric orbital configurations, shows that patches of these ions of meteoric origin are frequently present during magnetically active periods on the bottomside of the F-layer at middle and high latitudes. In particular the F-region metals statistically tend to appear at night in the vicinity of the main ionospheric trough (in a band of invariant latitudes approx. 10 degrees wide) and on the day side of the polar cap. These distributions were previously associated with the expected dynamics of ions in the F-region above 140 km where meridional neutral wind drag and convection electric fields are the dominant ion transport mechanisms. However, the main meteor deposition layer—the presumed source region of the metals—is located below 100 km where these transport mechanisms do not prevail. It is demonstrated that the Pedersen ion drifts driven by intense electric fields such as those associated with sub-auroral ion drifts (SAID) are sufficient to transport the long-lived metallic ions upward from the main meteor layer to altitudes where the drag of equatorial directed neutral winds and electric field convection can support them against the downward pull of gravity and transport them to other locations. The spatial and temporal distribution of the middle and high latitude F-region metals are consistent with the known characteristics of the electric fields and with the expected F-region ion dynamics.     

Ground-based observations of an isolated substorm

An isolated substorm occurred in Northern Scandinavia on 1 March, 1977 around magnetic midnight. The ionospheric phenomena associated with this substorm were studied by ground magnetometers, the Scandinavian Twin Auroral Radar Experiment (STARE), riometers and an all-sky camera. The physical properties of the auroral electrojet are determined from the ground magnetic field and the ionospheric electric field data. Mid and low latitude magnetic field data show evidence of field-aligned current flow. It is shown that the enhancement of the electrojet's current density is essentially determined by an increase in the ionospheric conductivity. The current system derived from the data of this study corresponds to a model of Yasuhara et al. (1975a).     

Theoretical studies of storm effects in ionospheric total electron content

Investigations have been made of the effects produced by thermospheric winds, composition changes and magnetospheric electric fields on the ionospheric F-layer during disturbed conditions. The results of the computations suggest that a combination of realistic temperature and electric field changes would explain fairly satisfactorily the observed changes in total electron content.     

Ionosphere-plasmasheet field-aligned currents and parallel electric fields

Two kinetic models for the auroral topside ionosphere are compared. The collisionless plasma distributed along an auroral magnetic field line behaves like a non-Ohmic conducting medium with highly non-linear characteristic curves relating the parallel current density to the potential difference between the cold ionosphere and the hot plasmasheet region. The (zero-electric current) potential difference, required to balance the current carried by the precipitating plasmasheet particles and the current transported by the outflowing ionospheric particles, depends on the ratio n_ps.e/n_th.e and T_ps.e/T_th.e of the plasmasheet and ionospheric electron densities and temperatures. When in the E-region the magnetic field lines are interconnected by a high conductivity plasma the resulting field-aligned currents driven by the magnetospheric potential distribution are limited by the integrated Pedersen conductivity of the ionospheric layers. These currents are not related to the parallel electric field intensity as they would be in Ohmic materials. The parallel electric field intensity is necessarily determined by the local quasi-neutrality of the plasma.     

The calculated and observed diurnal variation of the ionosphere over Millstone Hill on 23–24 March 1970

A numerical model of current F-region theory is use to calculate the diurnal variation of the mid-latitude ionospheric F-region over Millstone Hill on 23–24 March 1970, during quiet geomagnetic conditions. From the solar EUV flux, the model calculates at each altitude and time step primary photoelectron spectra and ionization rates of various ion species. The photoelectron transport equation is solved for the secondary ionization rates, photoelectron spectra, and various airglow excitation rates. Five ion continuity equations that include the effects of transport by diffusion, magnetospheric-ionospheric plasma transport, electric fields, and neutral winds are solved for the ion composition and electron density. The electron and ion temperatures are also calculated using the heating rates determined from chemical reactions, photoelectron collisions, and magnetospheric-ionospheric energy transport. The calculations are performed for a diurnal cycle considering a stationary field tube co-rotating with the Earth; only the vertical plasma drift caused by electric fields perpendicular to the geomagnetic field line is allowed but not the horizontal drift. The boundary conditions used in the model are determined from the incoherent scatter radar measurements of T_e, T_i and O⁺ flux at 800km over Millstone Hill (Evans, 1971a). The component of the neutral thermospheric winds along the geomagnetic field has an important influence on the overall ionospheric structure. It is determined from a separate dynamic model of the neutral thermosphere, using incoherent scatter radar measurements.The calculated diurnal variation of the ionospheric structure agrees well with the values measured by the incoherent scatter radar when certain restrictions are placed on the solar EUV flux and model neutral atmospheric compositions. Namely, the solar EUV fluxes of Hinteregger (1970) are doubled and an atomic oxygen concentration of at least $\text{10}{}^{11}\text{cm}{}^3$ at 120 km is required for the neutral model atmosphere. Calculations also show that the topside thermal structure of the ionosphere is primarily maintained by a flow of heat from the magnetosphere and the night-time F2-region is maintained in part by neutral winds, diffusion, electric fields, and plasma flow from the magnetosphere. The problem of maintaining the calculated night-time ionosphere at the observed values is also discussed.     

Some possible effects of Jupiter's rings on the Jovian inner plasmasphere

The ionospheric plasma density on magnetic field lines threading the Jovian rings which are located inside ～1.8 R_J on the jovigraphic equatorial plane, is calculated by using a rotating ion exosphere model. It is found that the bulk of the ionospheric particles on these field lines are on ballistic trajectories. On field lines approximately symmetric with respect to the jovigraphic equator, the ring, which to a first approximation would absorb the population of trapped particles, consequently has little effect. On field lines which are made asymmetric by the higher-order multipoles of Jupiter's field and the tilt of the dipole axis, the rings may have a significant effect. It is suggested that better definition of the rings' atmospheric and ionospheric properties is required to model these localized effects. If the rings are found to be an important plasma source for the inner magnetosphere, the present exospheric model will have to be revised.     

Three-dimensional current system in different phases of a substorm

It is assumed that the three-dimensional current system of a substorm passes three successive stages. (1) When a dawn-to-dusk magnetospheric electric field appears, a current system with field-aligned currents at the poleward boundary of the auroral zone arises. An equivalent ionospheric current system calculated, taking into account a day-night asymmetry of ionospheric conductivity, looks like the well-known DP-2 system including an eastward low-latitude current and a greater magnitude of the dusk vortex in comparison with the dawn one. (2) An electric drift of plasma towards the Earth leads to the appearance of a westward partial ring current increasing in time. This current is closed by field-aligned currents at the equatorward boundary of the auroral zone. The calculated equivalent current system is similar to the well-known one of the precursory phase. (3) An increase of the auroral ionospheric conductivity during the expansive phase produces an increase of all currents and a turning of field-aligned currents at the equatorward boundary of the auroral zone relative to those at the poleward one. The calculated equivalent current system is similar to the DP-1 system.     

On the origin of the hot ions in the disturbed dayside magnetosphere

On the basis of observations in the dayside magnetosphere of the O⁺ and H⁺ ion densities as function of radial distance under fairly undisturbed and under storm conditions it is argued that acceleration of the hot magnetospheric ions of ionospheric origin cannot be limited to the outer parts of the field tubes. The extraction process seems to work below 1000 km altitude in storm conditions and to have a fairly small extension in altitude. The acceleration mechanism(s) do(es) not affect only one ion species. Variation in the altitude of the extraction of ionospheric ions is the most likely reason for the observed variations in the n(O⁺)/n(H⁺) ratio. Extraction of ionospheric ions into the magnetosphere does not seem to be a main cause of the storm time density decrease of the ionosphere.     

Ionospheric plasma acceleration at Mars: ASPERA-3 results

The Analyzer of Space Plasma and Energetic Atoms (ASPERA) on-board the Mars Express spacecraft (MEX) measured penetrating solar wind plasma and escaping/accelerated ionospheric plasma at very low altitudes (250 km) in the dayside subsolar region. This implies a direct exposure of the martian topside atmosphere to solar wind plasma forcing leading to energization of ionospheric plasma. The ion and electron energization and the ion outflow from Mars is surprisingly similar to that over the magnetized Earth. Narrow “monoenergetic” cold ion beams, ion beams with broad energy distributions, sharply peaked electron energy spectra, and bidirectional streaming electrons are particle features also observed near Mars. Energized martian ionospheric ions (O⁺, O⁺₂, CO⁺₂, etc.) flow in essentially the same direction as the external sheath flow. This suggests that the planetary ion energization couples directly to processes in the magnetosheath/solar wind. On the other hand, the beam-like distribution of the energized plasma implies more indirect energization processes like those near the Earth, i.e., energization in a magnetized environment by waves and/or parallel (to B) electric fields. The general conditions for martian plasma energization are, however, different from those in the Earth's magnetosphere. Mars has a weak intrinsic magnetic field and solar wind plasma may therefore penetrate deep into the dense ionospheric plasma. Local crustal magnetization, discovered by Acuña et al. [Acuña, M.J., Connerey, J., Ness, N., Lin, R., Mitchell, D., Carlsson, C., McFadden, J., Anderson, K., Rème, H., Mazelle, C., Vignes, D., Wasilewski, P., Cloutier, P., 1999. Science 284, 790-793], provide some dayside shielding against the solar wind. On the other hand, multiple magnetic anomalies may also lead to “hot spots” facilitating ionospheric plasma energization. We discuss the ASPERA-3 findings of martian ionospheric ion energization and present evidences for two types of plasma energization processes responsible for the low- and mid-altitude plasma energization near Mars: magnetic field-aligned acceleration by parallel electric fields and plasma energization by low frequency waves.     

Modulation of type Pi waves by temporal variations in ionospheric conductivity in a three-dimensional magnetosphere-ionosphere current system

It is suggested that Pi ULF waves are generated from magnetosphere-ionosphere current systems. This current system is modeled by an R (resistance), C (capacitance) and L (inductance) circuit loop in which R = R(t). We studied three cases of modification of ULF waves by variations in ionospheric conductivity: (1) ω ? ω', (2) ω ≈ ω' and (3) ω ? ω', where ω and ω' are the frequency of the driving electric field and ionospheric conductivity variations, respectively, assuming that both variations are sinusoidal. The characteristics of the modification are very different in these three cases. In case 1, the envelope of the ULF wave intensity correlates with the variations in ionospheric conductivity. In case 2, the wave form of the ULF waves is slightly modified from a sinusoidal wave. In case 3, high frequency components are generated in the ULF wave form due to rapid oscillations in ionospheric resistance. We present observational evidence for the existence of the three types of modifications.     

A numerical study of polar ionospheric currents

Numerical calculations for the electric current in the polar ionosphere have been made by assuming some realistic distributions of the electric field and conductivity. Two dynamo actions are taken into account; one of which is induced by ionospheric winds and the other by the solar wind. For the solar wind dynamo action, it is found that the secondary polarization field caused by non-uniform distribution of ionospheric conductivity is much larger than the primary field induced by the solar wind, suggesting its important effect on charged particles in the magnetosphere, and that the irrotational current having a source and sink is of the same order of magnitude as the solenoidal current closing its circuit in the ionosphere. It is also found that the solar wind is, in general, more effective than the ionospheric winds in producing polar current systems such as DP 1 and 2, but in some cases the ionospheric winds have a significant effect on the current distribution.     

Modeling the midlatitude F-region ionospheric storm using east-west drift and a meridional wind

Under magnetically quiet conditions, ionospheric plasma in the midlatitude F-region corotates with the Earth and relative east-west drifts are small compared to the corotation velocity. During magnetic storms, however, the enhanced dawn-to-dusk magnetospheric convection electric field often penetrates into the midlatitude region, where it maps into the ionosphere as a poleward electric field in the 18:00 LT sector, producing a strong westward plasma drift. To evaluate the ionospheric response to this east-west drift, the time-dependent O⁺ continuity equation is solved numerically, including the effects of production by photoionization, loss by charge exchange and transport by diffusion, neutral wind and $\text{E}\text{×}\text{B}$ drift. In this investigation only the neutral wind's meridional component and east-west $\text{E}\text{×}\text{B}$ drift are included. It is found that an enhanced equatorward wind coupled with westward drift produces an enhancement in the peak electron density (NMAX(F2)) and in the electron content (up to 1000 km) in the afternoon sector and a subsequent greater-than-normal decay in ionization after 18:00 LT. These results agree in general with midlatitude F-region ionospheric storm observations of NMAX(F2) and electron content which show an afternoon enhancement over quiet-time values followed by an abrupt transition to lower-than-normal values. Westward drift appears to be a sufficient mechanism in bringing about this sharp transition.     

Mapping electrostatic potentials from the ionosphere to the magnetosphere

With the advent of long duration incoherent scatter radar experiments measuring ionospheric plasma convection over a wide range of latitudes and at all local times, the mapping and study of the spatial and temporal distribution of electric fields within the magnetosphere becomes possible. We consider the problems of mapping the ionospheric electrostatic potential distribution into the magnetosphere under the assumption that magnetic field lines are electrostatic equipotentials. We address the practical problem of developing a mapping technique which can adequately project ionospheric observations, acquired at different geographic longitudes, into the magnetosphere. The mapping must include the effects of a magnetotail when considering auroral latitudes and is a function of diurnal and seasonal effects and is ionospheric longitude dependent. Mapping observed ionospheric potential distributions into the magnetosphere yields parameters such as: the distribution of electric fields in the magnetosphere, the cross tail potential, the location of plasmapause, the local time of the ‘stagnation’ bulge region and the gradients of the electric field in the vicinity of the plasmapause.