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.     

A new theoretical approach to the modelling of toroidal Pc5 pulsations

A model is developed to represent a toroidal mode of Pc5 geomagnetic pulsations. It is shown that this model is consistent in its predictions, such as the latitude profiles of amplitude and phase and their dependence on the height integrated Pedersen conductivity, Σ_p, with those of Walker's (1980) theory. It is also shown that this theory is relatively easily capable of accommodating (i) a variety of field line plasma mass density distributions, (ii) a variety of external excitation schemes, (iii) unequal Σ_p's at each end of the field lines and (iv) non-dipolar geomagnetic fields. The theory yields the transient as well as the steady state response, an important feature permitting application to short-lived events or to those for which the generator is amplitude modulated. It is shown, for instance, that the amplitude-latitude profile varies during the transient. It is also shown that the steady state latitude profiles of amplitude and phase are the dual of those observed as a function of frequency when the excitation frequency is scanned through a resonance. A more realistic steady state energy flow from a generator along the field lines to the ionosphere is inherent in this theory compared with that from the mode to the ionosphere which is inherent in Walker's theory.     

Infrared brightness temperature of Saturn's rings based on the inhomogeneous multilayer assumption

Models of Saturn's rings based on the classical multilayer assumption have been studied in the infrared. Thermal energy balance of Saturn's rings is treated rigorously by solving the infrared radiative transfer equations. It was found that a homogeneous multilayer model is incompatible with the observed infrared brightness variation of the A and B rings, although it can fit that of the C ring. The alternative inhomogeneous multilayer model with dark particles within a bright haze of small icy particles is presented in order to satisfy the available infrared data of the A, B, and C rings. The results based on the inhomogeneous multilayer model may be summarized as follows: The observed infrared brightness data of the three rings are explained in terms of the different optical thickness without having significant differences in the ring-particle properties, such as albedo, spin rate, and sizes. But each ring contains a different amount of bright haze particles and their concentration within the rings depends on whether or not dark particles emit radiation mostly from one hemisphere (slow rotator and/or low thermal inertia). If a dark particle is an isothermal radiator, the possible ranges of A₁ and A₂ for all three rings are given by 0.9 ? A₁ ? 0.95 and 0.0 ? A₂ ? 0.15, where A₁ and A₂ are the bolometric bond albedos of a bright haze and a dark particle, respectively. The possible ranges of the optical thickness ratio X of the dark particle layer to the total ring layer for the rings A, B, and C are given by 0.65 ? X ? 0.75, 0.8 ? X ? 0.9, and 0.8 ? X ? 1.0, respectively. If a dark particle is a slow rotator, we obtain 0.9 ? A₁ ? 0.95 and 0.0 ? A₂ ? 0.4 for all three rings. The ranges of X for the rings A, B, and C are given by 0.35 ? X ? 0.7, 0.65 ? X ? 0.9, and 0.35 ? X ? 1.0, respectively. In this paper, for the first time, a consistent model is presented which is applicable to all three rings from the multilayer point of view.     

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     

Saturn's rings: II. Condensations of light and optical thickness of Cassini's division

“Condensations” of light have been observed when Saturn's rings are seen almost edge on, and the Sun and the Earth are on opposite sides of the ring plane. These condensations are associated with ring C and Cassini's division. If the relative brightness between the two condensations and the optical thickness of ring C are known, we can calculate the optical thickness of Cassini's division, τ_CASS. Using Barnard's and Sekiguchi's measurements, we have obtained 0.01 ? τ_CASS ? 0.05. A brightness profile of the condensations which agrees well with visual observations is also presented.We are able to set an upper limit of about 0.01 for the optical thickness of any hypothetical outer ring. This rules out a ring observed by C. Cragg in 1954, but does not eliminate the D′ ring observed by Feibelman in 1967.It is known that the outer edge of ring B is almost at the position of the 1/2 resonance with Mimas. Franklin, Colombo, and Cook explained this fact in 1971, postulating a total mass of ring B of 10^?6M_SATURN. We have derived a formula for the mass of the rings, which is a linear function of the mean particle size. We find that 10^?6M_SATURN implies large particles (～70m). If the particles are small (～10cm), as currently believed, the total mass of ring B is not enough to shift the outer edge. We conclude that the above explanation and current size estimates are inconsistent.     

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.     

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).     

A mathematical magnetospheric field model with independent physical parameters

A quantitative magnetospheric magnetic field model has been calculated in three dimensions. The model is based on an analytical solution of the Chapman-Ferraro problem. For this solution, the magnetopause was assumed to be an infinitesimally thin discontinuity with given geometry. The shape of the dayside magnetopause is in agreement with measurements derived from spacecraft boundary crossings.The magnetic field of the magnetopause currents can be derived from scalar potentials. The scalar potentials result from solutions of Laplace's equation with Neumann's boundary conditions. The boundary values and the magnetic flux through the magnetopause are determined by all magnetic sources which are located inside and outside the magnetospheric cavity. They include the Earth's dipole field, the fields of the equatorial ring current and tail current systems, and the homogeneous interplanetary magnetic field. In addition, the flux through the magnetopause depends on two constants of interconnection which provide the possibility of calculating static interconnection between magnetospheric and interplanetary field lines. Realistic numerical values for both constants have been derived empirically from observed displacements of the polar cusps which are due to changes in the orientation of the interplanetary field. The transition from a closed to an open magnetosphere and vice versa can be computed in terms of a change of the magnetic boundary conditions on the magnetopause. The magnetic field configuration of the closed magnetosphere is independent of the amount and orientation of the interplanetary field. In contrast, the configuration of the open magnetosphere confirms the observational finding that field line interconnection occurs primarily in the polar cusp and high latitude tail regions.The tail current system reflects explicitly the effect of dayside magnetospheric compression which is caused by the solar wind. In addition, the position of the plasma sheet relative to the ecliptic plane depends explicitly on the tilt angle of the Earth's dipole. Near the tail axis, the tail field is approximately in a self-consistent equilibrium with the tail currents and the isotropic thermal plasma.The models for the equatorial ring current depend on the D_st-parameter. They are self-consistent with respect to measured energy distributions of ring current protons and the axially symmetric part of the magnetospheric field.     

Fluid description of a collisionless plasma and electric current systems in cosmic plasmas

The large-scale structure of a collisonless, two-component plasma with a typical Larmor radius of ions ? and scale-lengthL is discussed using Maxwell transport equations. Special attention is paid to the situations in which the usual one-fluid model of plasma based on the expansion of the transport equations in the powers of the ratio ?/L is not a satisfactory approximation. The one-fluid model fails if the magnetic-field-aligned component of the mass velocity or the magnetic-field-aligned component of the typical random velocity of particles is much larger than the other components of the mass and random velocities. The model also fails if the component of the typical random velocity of particles, which is perpendicular to the field lines, substantially exceeds the mass velocity of particles across the field lines. A quasi-static plasma is discussed as an example of plasmas on which the expansion in the powers of ?/L is not applicable. The relation between the electric current flowing in a quasi-static plasma (or in a hot plasma streaming along the field lines) and the topology of the magnetic lines of force is analysed. There are two distinguishable currents of different origin in such a plasma. Magnetic field generated by the currents acquires a geometry in which one current flows in the surfaces perpendicular to the binormals to the field lines while the other current flows along the binormals.     

Saturn's rings: Particle size distributions for thin layer models

The small physical thickness of Saturn's rings requires that radio occultation observations be interpreted using scattering models with limited amounts of multiple scatter. A new model in which the possible order of near-forward scatter is strictly limited allows for the small physical thickness, and can be used to relate Voyager 1 observations of 3.6-and 13-cm wavelength microwave scatter from Saturn's rings to the ring particle size distribution function n(a), for particles with radius 0.001 ≤ a ≤ 20 m. This limited-scatter model yields solutions for particle size distribution functions for eight regions in Saturn's rings, which exhibit approximately inverse-cubic power-law behavior, with large-size cutoffs in particle radius ranging from about 5 m in ring C to about 10 m in parts of ring A. The power-law index is about 3.1 in ring C, about 2.8 in the Cassini division, and increases systematically with radial location in ring A from 2.7 at 2.10R_s to slightly more than 3.0 at 2.24R_s. Corresponding mass densities are 32–43 kg/m² in ring C, 188 kg/m² in the Cassini division, and 244–344 kg/m² in ring A, under the assumption that the material density of the particles is 0.9 g/cm³. These values are a factor of 1 to 2 lower than first-order mass loading estimates derived from resonance phenomena. In view of the uncertainties in the measurements and in the linear density wave model, and the strong arguments for icy particles with specific gravity not greater than about 1, we interpret this discrepancy as being indicative of possible differences in the regions studied, or systematic errors in the interpretation of the scattering results, the density wave phenomena, or some combination of the above.     

Precipitation of charged particles by a parallel electric field

A time-dependent model of the effect of a parallel electric field on particle precipitation from a closed field-line has been constructed and the results are presented. A pattern of field-aligned pitch-angle distributions and energy peaks develops rapidly and then persists unchanged in shape while the intensity decreases for a time of the order of the bounce period of the energetic particles. It is shown that the structures in velocity space are created by the juxtaposition of particles from different source populations. Four sources are found to be sufficient to reproduce the principal features observed frequently by rockets and satellites. They are, a trapped plasma sheet distribution, a loss-cone partially filled by pitch-angle diffusion at the equator, cold ionospheric plasma which has flowed outward along the field line and particles backscattered from the precipitation into the atmosphere.The model develops density gradients and discontinuities far sharper than any observed, so that any parallel electric field actually occurring in an aurora must be accompanied by strong wave-particle interactions either as part of the accelerating mechanism or as a result of the density gradients produced by it.     

Criticism of reconnection models of the magnetosphere

Reconnection involves singular lines called X-lines on the day and night sides of the magnetosphere, and the reconnection rate is proportional to the component of the electric field along the X-line. Although there is some indirect support for this model, nevertheless direct support is totally lacking. However, there are two distinct pieces of clearly contradictory observational evidence on the dayside. First is the failure to account for the implied energy dissipation by the magnetopause current, over 10¹¹ W, which should be easily observable as heating or enhanced flow of the plasma near the magnetopause. In marked contrast to this prediction, HEOS-2 satellite data reveal a plasma with decreased energy density and reduced flow. Second, the boundary of closed magnetic field lines is in the wrong location. In the reconnection process the plasma outflow would cut across open field lines toward higher latitudes; there should be a band of open field lines equatorward of the cleft. Observations of trapped energetic particles indicate closed field lines within the entry layer and cleft. Either one of these pieces of evidence is sufficient by itself to require drastic revision, even rejection, of the reconnection model. There is also contradictory evidence on the night side. The last closed field line capable of trapping energetic particles is poleward of auroral arcs. The implication is that the X-line is at the distant magnetopause, and not in the plasma sheet. Consequently, even if the reconnection process were operative at the nightside X-line, it would be isolated from steady state plasma sheet and auroral processes. On the other hand, substorm phenomena, in which stored magnetic energy is converted into particle kinetic energy, necessarily involve an induced electric field; that is excluded in theories of the reconnection process in which it is assumed that curl E = 0. Nevertheless, the observed easy access of energetic solar flare particles to the polar caps, and especially the preservation of interplanetary anisotropies as differences between the two polar caps, argues strongly for an open magnetosphere, with interconnection between geomagnetic and inter-planetary magnetic field lines. It is suggested that the resolution of this apparent paradox involves electric fields parallel to the magnetic field lines somewhere on the dawn and dusk sides of the magnetosphere, with an equipotential dayside magnetopause.     

Some properties of geomagnetic field pulsations in the 0.1–1.0-Hz frequency range: A quantitative description and comparison with the satellite and ground-based observations

By using an image-dipole magnetic field model for a variety of plasma density profiles we have studied the latitude effect of the 0.1–1.0-Hz hydromagnetic wave propagation in the Earth's magnetosphere. On comparing the results of signal group delay time calculations for dipole and model magnetic fields with ground and satellite observations we obtain some propagation characteristics of Pc1s and localize the regions of their generation. Our results show that most high-latitude Pc1 events are generated in the outer magnetosphere in accordance with ground and satellite observations and theoretical considerations. The non-dipole geometry of the geomagnetic field in the outer magnetosphere (at geomagnetic latitudes φ₀ > 66°, L > 6) has a significant effect on the hydromagnetic wave propagation.     

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.     

Jupiter's radiation belts

A model for the production and loss of energetic electrons in Jupiter's radiation belt is presented. It is postulated that the electrons originate in the solar wind and are diffused in toward the planet by perturbations which violate the particles' third adiabatic invariant. At large distances, magnetic perturbations, electric fields associated with magnotospheric convection, or interchange instabilities driven by thermal plasma gradients may drive the diffusion. Inside about 10 R_J the diffusion is probably driven by electric fields associated with the upper atmosphere dynamo which is driven by neutral winds in the ionosphere. The diurnal component of the dynamo wind fields produces a dawn-dusk asymmetry in the decimetric radiation from the electrons in the belts, and the lack of obvious measured asymmetries in the decimetric radiation measurements provides estimates of upper limits for these Jovian ionospheric neutral winds. The average diurnal winds are less than or comparable to those on earth, but only modest fluctuating winds are required to drive the energetic electron diffusion referred to above.The winds required to diffuse the energetic particles across the orbit of the satellite lo in a time equal to their drift period are also estimated. If Io is non-conducting, modest winds are required, but if Io is conducting, only small winds are needed. It is concluded that both protons and electrons are diffused in from the solar wind to small distances without serious losses occurring due to the particles being swept up by the satellites.Consideration of proton and electron diffusion in energy shows that once the electrons become relativistic, the ratio of proton to electron energy increases. Thus, if protons and electrons have the same energy in the solar wind, when the electrons reach nMeV, the protons will be nMeV if n ? 1 or n² MeV if n ? 1. If the proton-to-electron energy ratio is initially, e.g., 5, then these figures are 5n and 5n², respectively.     

The formation of saturn's satellites and rings,as influenced by saturn's contraction history

We have used Pollack et al.'s 1976 calculations of the quasi-equilibrium contraction of Saturn to study the influence of the planet's early high luminosity on the formation of its satellites and rings. Assuming that the condensation of ices ceased at the same time within Jupiter's and Saturn's primordial nebulae, and using limits for the time of cessation derived for Jupiter's system by Pollack and Reynolds (1974) and Cameron and Pollack (1975), we arrive at the following tentative conclusions. Titan is the innermost satellite at whose position a methane-containing ice could condense, a result consistent with the presence of methane in this satellite's atmosphere. Water ice may have been able to condense at the position of all the satellites, a result consistent with the occurrence of low-density satellites close to Saturn. The systematic decrease in the mass of Saturn's regular satellites with decreasing distance from Saturn may have been caused partially by the larger time intervals for the closer satellites between the start of contraction and the first condensation of ices at their positions and between the start of contraction and the time at which Saturn's radius became less than a satellite's orbital radius. Ammonia ices, principally NH₄SH, were able to condense at the positions of all but the innermost satellites.Water ice may bave been able to condense in the region of the rings close to the end of the condensation period. We speculate that the rings are unique to Saturn because on the one hand, temperatures within Jupiter's Roche limit never became cool enough for ice particles to form before the end of the condensation period and on the other hand, ice particles formed only very early within Uranus' and Neptune's Roche limits, and were eliminated by gas drag effects that caused them to spiral into the planet before the gas of these planets' nebula was eliminated. Gas drag would also have eliminated any rocky particles initially present inside the Roche limit.We also derive an independent estimate of several million years for the time between the start of the quasi-equilibrium contraction of Saturn and the cessation of condensation. This estimate is based on the density and mass characteristics of Saturn's satellites. Using this value rather than the one found for Jupiter's satellites, we find that the above conclusions about the rings and the condensation of methane-and ammonia-containing ices remain valid.     

Size distribution of particles in planetary rings

Harris (Icarus24, 190–192) has suggested that the maximum size of particles in a planetary ring is controlled by collisional fragmentation rather than by tidal stress. While this conclusion is probably true, estimated radius limits must be revised upward from Harris' values of a few kilometers by at least an order of magnitude. Accretion of particles within Roche's limit is also possible. These considerations affect theories concerning the evolution of Saturn's rings, of the Moon, and of possible former satellites of Mercury and Venus. In the case of Saturn's rings, comparison of various theoretical scenarios with available observational evidence suggests that the rings formed from the breakup of larger particles rather than from original condensation as small particles. This process implies a distribution of particle sizes in Saturn's rings possibly ranging up to ～100 km but with most cross-section in cm-scale particles.     

Local plasma processes and enhanced electron densities in the lower ionosphere in magnetic cusp regions on Mars

Both the MARSIS ionospheric sounder and the charged particle instrument package ASPERA-3 are experiments on board the Mars Express spacecraft. Joint observations have shown that events of intense ionospheric electron density enhancements occur in the lower ionosphere of magnetic cusp regions, and that these enhancements are not associated with precipitation of charged particles above a few hundred electron volts (<300 eV). To account for the enhancement by particle precipitation, electron fluxes are required with mean energy between 1 and 10 keV. No ionizing radiation, neither energetic particles nor X-rays, could be identified, which could produce the observed density enhancement only in the spatially limited cusp regions. Actually, no increase in ionizing radiation, localized or not, was observed during these events. It is argued that the process causing the increase in density is controlled mainly by convection of ionosphere plasma driven by the interaction between the solar wind and crustal magnetic field lines leading to excitation of two-stream plasma waves in the cusp ionosphere. The result is to heat the plasma, reduce the electron–ion recombination coefficient and thereby increase the equilibrium electron density.     

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.     

Results from a swept-frequency ionospheric impedance probe

Experimental results are presented for the impedance of a rocket-borne dipole antenna immersed in the ionospheric plasma. The dependence of several interesting impedance artifacts upon the antenna position relative to the Earth's magnetic field and rocket motion through the ionospheric plasma are shown. Possible evidence for plasma compressibility is provided by an impedance discontinuity occurring consistently at approximately twice the electron cyclotron frequency and a frequency shifted cyclotron-resonance impedance minimum.