Search

A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere- ionosphere coupling loop. First, in addition to the RC ions, we solve an electron kinetic equation in our model. Second, using the relation of Galand and Richmond [2001], we calculate the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our code. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has has been comprehensively studied, and compare our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and, their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. This changes the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere-ionosphere coupling loop. First, in addition to the RC ions, we have solved an electron kinetic equation in our model. Second, using the relation of Galand and Richmond, we have calculated the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our model. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has been comprehensively studied by Fok et al., and have compared our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. These effects change the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

The initial results of a model of inner magnetospheric convection are presented. The model employs the Rice convection model with a magnetic field computed with the constraint of magnetostatic equilibrium. The approach computes equilibria from a magnetofriction code which is a modified version of the Hesse-Birn equilibrium code adopted for use in the inner magnetosphere. The code uses the pressure distribution computed from the Rice convection model to update the magnetic field. The algorithm used to compute the inner magnetospheric equilibria is outlined, and the coupling of the equilibrium code with the convection model is described.

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

The Rice Convection Model (RCM) is used to illustrate theoretical possibilities for generating region 1 Birkeland currents by pressure gradients on closed field lines in the Earth's magnetosphere. Inertial effects and viscous forces are neglected. The RCM is applied to idealized cases, to emphasize the basic physical ideas rather than realistic representation of the actual magnetosphere. Ionospheric conductance is taken to be uniform, and the simplest possible representations of the magnetospheric plasma are used. Three basic cases are considered: (1) the case of pure northward Interplanetary Magnetic Field (IMF), with cusp merging assumed to create new closed field lines near the nose of the magnetosphere, following the suggestion by Song and Russell (1992); (2) the case where Dungey-type reconnection occurs at the nose, but magnetosheath plasma somehow enters closed field lines on the dawnside and duskside of the merging region, causing a pressure-driven low-latitude boundary layer; and (3) the case where Dungey-type reconnection occurs at the nose, but region 1 currents flow on sunward drifting plasma sheet field lines. In case 1, currents of region 1 sense are generated by pressure gradients, but those currents do not supply the power for ionospheric convection. Results for case 2 suggest that pressure gradients at the inner edge of the low-latitude boundary layer might generate a large fraction of the region 1 Birkeland currents that drive magnetospheric convection. Results for case 3 indicate that pressure gradients in the plasma sheet could provide part of the region 1 current.

The Rice convection model (RCM) has been modified to a form suitable for Jupiter (RCM-J) to study plasma interchange motion in and near the Io plasma torus. The net result of the interchange is that flux tubes, heavily loaded with torus plasma, are transported outward, to be replaced by tubes containing little low-energy (less than 1 keV) plasma. The process is numerically simulated in terms of time evolution from an initial torus that is longitudinally asymmetric and with gradually decreasing density outward from Io's orbit. In the simulations, the nonlinear stage of the instability characteristically exhibits outreaching fingers of heavily-loaded flux tubes that lengthen at an accelerating rate. The principal finding is that the primary geometrical form of outward transport of torus plasma in Jupiter's magnetosphere is through long, outward-moving fingers of plasma. In the simulations, the fingers mainly form in the active sector of the Io torus (the heavier side of the asymmetric torus), and they are spaced longitudinally roughly 20 deg apart.

This paper describes the physics of the Harang discontinuity and its application to forecasting of the magnetospheric environment. Gradient/curvature depletion of energetic ions from earthward drifting flux tubes results in upward Birkeland currents. Ionospheric closure rotates the dawn-dusk convection electric field resulting in the Harang reversal. This physical mechanism has been incorporated into the Rice Convection Model (RCM) and can serve as a fundamental part of the forecasting of the magnetospheric system. The inclusion of this physics within the RCM results in a physical picture of the Harang discontinuity which did not previously exist. Model results are close to observations and represent a significant improvement in predicted electric field and Birkeland current patterns.

The Rice Convection Model (RCM) is extended and merged with empirical models so as to cover the entire high-latitude ionosphere with the aim of providing precipitation and electric field inputs for ionosphere and thermosphere modelers and producing a model in which the boundaries of the precipitation and electric field patterns maintain physically consistent relationships to each other. The computed auroral electron energy flux, plotted as a function of latitude, exhibited an exaggerated two-peak structure. When no floor was placed under the precipitation rate, the minimum between the two peaks was much too deep to be consistent with typical observations. The regions of excessively weak precipitation map to equatorial distances of 15-35 RE and thus to the regions of the plasma sheet that were not included in previous self-consistent convection calculations.

It is found that the Harang discontinuity results from a dawn/dusk asymmetry in total energetic plasma content across the nightside of the magnetosphere. This effect, referred to as the dawnside depletion effect, is discussed in detail. The basic physical assumptions usually regarded as appropriate for the description of the closed field line plasma sheet are outlined and a quantitative, bounce-averaged description of the particle drifts is presented. A simplified, illustrative model is given of the magnetospheric geometry and convection electric field in order to demonstrate the dawnside depletion effect and examine its ionospheric and magnetospheric implications. For a quantitative examination of the implications of the dawnside depletion effect and the resulting ionospheric electric field structure, the Rice convection model Z(RCM) is used. A brief description is provided of the RCM, the specific run setups, and the RCM results. A general discussion of the basic physical connection between cross-tail gradients in plasma contents and the Harang discontinuity is presented.

F-region incoherent scatter radar drift observations from Millstone Hill and Jicamarca, h-prime F observations from Huancayo, and high latitude ground-magnetometer measurements taken during the Sundial 1986 campaign are used to study the relationship between plasmaspheric electric field perturbations and high latitude currents during disturbed periods. The observations are in good agreement with numerical results from a Rice Covection Model run that involved a sharp increase in the polar cap potential drop followed by a subsequent decrease. The zonal disturbance electric field pattern is latitude independent, and the corresponding amplitudes change approximately as L exp n (where n is about 1.5). The meridional electric field patterns and amplitudes have larger latitudinal variations. The mid-, low, and equatorial electric fields from the Rice Convection Model are in good agreement with previous results from the semianalytic, Senior-Blanc (1987) model. Also discussed are three physical mechanisms (over-shielding, fossil winds, and magnetic reconfiguration) that contribute to the long lasting (1-2 h) equatorial zonal electric field perturbations associated with a sudden northward turning of the IMF. It is predicted that the penetration of high latitude electric fields to low latitudes should, in general, be closely related to the rate of motion of the shielding layer and the equatorward edge of the diffuse aurora.

Electric-field-penetration events have been identified using F-region vertical-drift measurements obtained in the October 6-13, 1984 period by the Jicamarcan incoherent-backscatter radar and corresponding h-prime F measurements from ionosondes at Fortaleza, Cachoeira Paulista, and Dakar. Predictions made using the Rice Convection Model for the pattern, strength, and duration of the low-latitude electric field occurring in response to an increasing high-latitude convection agree with observations. The observed 1-2 h duration of the low-latitude response to decreased convection can be explained by the fossil-wind theory of Richmond (1983).

The traditional theoretical interpretations of the observed plasmapause are compared, namely, the plasmapause as: (1) The boundary between closed flux tubes that have been in the inner magnetosphere for several days and those that have recently drifted in from the magnetotail or (2) the last closed electric equipotential. Although the two interpretations become equivalent in the case where the electric-field pattern is steady for several days, interpretation 1 seems theoretically more secure for typical magnetospheric conditions, due to the essentially time-dependent nature of the magnetospheric electric field. The results of old theoretical studies of the effects of time variations in the electric-field pattern on the shape of the plasmapause are reviewed briefly. The formulation of the present version of the Rice Convection Model is also reviewed. Preliminary results of recent computations of quiet-time electric fields, carried out with this model, are presented and discussed. Quiet-time thermospheric winds are found to have only minor effects on drift paths of magnetospheric particles.

Data recorded by Defense Meteorological Satellite Program, Tiros and P-78-1 satellites for the CDAW 6 event on March 22, 1979, have been compared with a statistical model of precipitating electron fluxes. Comparisons have been made on both an orbit-by-orbit basis and on a global basis by sorting and binning the data by AE index, invariant latitude, and magnetic local time in a manner similar to which the model was generated. It is concluded that the model flux agrees with the data to within a factor of two, although small features and the exact locations of features are not consistently reproduced. In addition, the latitude of highest electron precipitation usually occurs about 3 deg more poleward in the model than in the data. This discrepancy is attributed to ring current inflation of the storm time magnetosphere.

It is pointed out that the maximum electrostatic potential difference across the polar cap, Phi, is a fundamental measure of the coupling between the solar wind and the earth's magnetosphere/ionosphere sytem. During the Coordinated Data Analysis Workshop (CDAW) 6 intervals, no suitably instrumented spacecraft was in an appropriate orbit to determine the polar-cap potential drop directly. However, two recently developed independent techniques make it possible to estimate the polar-cap potential drop for times when direct spacecraft data are not available. The present investigation is concerned with a comparison of cross-polar-cap potential drop estimates calculated for the two CDAW 6 intervals on the basis of these two techniques. In the case of one interval, the agreement between the potential drops and Joule heating rates is relatively good. In the second interval, however, the agreement is not very good. Explanations for this discrepancy are discussed.

The Rice Convection Model deals with large-scale processes in the earth's inner and middle magnetosphere, including coupling to the ionosphere. Starting from appropriate initial and boundary conditions, the model computes the following physical parameters: ionospheric electric fields and currents; magnetospheric particle distributions, electric fields, and electric currents; and magnetic-field-aligned (Birkeland) currents connecting the two regions. This paper evaluates work on the model, with emphasis on the assumptions made, the basic equations, and the numerical methods. The theoretical basis of the model is compared and contrasted with standard magnetohydrodynamics. The limitations imposed by the major assumptions are discussed. Model inputs and boundary conditions are listed, and the methods of specifying them discussed. Some physical conclusions and insights that have been gained from the model are listed and described very briefly. References are given to published discussions of the major points of physics.

The following international Magnetospheric Study quantitative models of observed ionosphere-magnetosphere events are reviewed: (1) a theoretical model of convection; (2) algorithms for deducing ionospheric current and electric-field patterns from sets of ground magnetograms and ionospheric conductivity information; and (3) empirical models of ionospheric conductances and polar cap potential drop. Research into magnetic-field-aligned electric fields is reviewed, particularly magnetic-mirror effects and double layers.

On the basis of Birkeland current observations, analytic stability arguments demonstrate that there may be a sector within the plasma sheet in the night side magnetosphere where plasma pressures are reduced relative to the surrounding regions. Physical arguments, in addition to this analysis, imply that a depleted region near the center of the magnetotail would be stable against the interchange instability, thereby generating currents of the same sense as the region-1 currents. The depleted region is presently simulated by means of computer experiments based on the Rice Convection Model. Results indicate that flux tube content gradients across the tailward boundary can cause region-1 currents to flow along sunward-convecting flux tubes, in general agreement with satellite data.

During the past ten years, substantial progress has been made in the development of quantitative models of convection in the magnetosphere and of the electrodynamic processes that couple that magnetosphere and ionosphere. Using a computational scheme first proposed by Vasyliunas, the convection models under consideration separate the three-dimensional problem of convection in the inner magnetosphere/ionosphere into a pair of two-dimensional problems coupled by Birkeland currents flowing between the two regions. The logic, development, and major results of the inner magnetosphere convection model are reviewed with emphasis on ionospheric and magnetospheric currents. A major theoretical result of the models has been the clarification of the relationship between the region 1/region 2 picture of field-aligned currents and the older partial ring current/tail current interruption picture of substorm dynamics.

It is pointed out that the energy which is released during a magnetospheric substorm has large and obvious effects on the earth's auroral ionosphere. The present investigation is concerned with the question of the effect of the ionosphere on substorm phenomena which occur near the magnetospheric equatorial plane. Plasma transfer between ionosphere and magnetosphere are discussed along with general theoretical considerations regarding magnetospheric currents, plasma sheet dynamics, and ring current injection. Attention is given to questions concerning the occurrence of substorms, ion circulation in the magnetospheric convection cycle, time-dependent plasma ejection from plasma-sheet flux tubes, and a schematic of the magnetospheric-convection current circuit.

Data from the low energy electron (LEE) experiments on the Atmosphere Explorer C and D satellites have been used to determine the average global distribution of the energy flux of precipitating auroral electrons and their average energy for different levels of geomagnetic activity. Measurements from the Atmosphere Explorer unified abstract file (15-s resolution) have been binned according to invariant latitude (in the range 50-90 deg), magnetic local time, and geomagnetic activity as measured by the Kp and auroral electrojet (AE) indices, separately. Bin-averaged values of precipitating energy flux and average energy have been calculated, and a smoothing algorithm used to reduce stochastic variations in the raw data. The results indicate that, for the parameters studied, the AE inces does a superior job of ordering the data with regard to geomagnetic activity. The global distribution of the auroral enhancement porition of the Pedersen and Hall conductances were inferred from the data by means of an empirical fit to detailed energy deposition calculations.

Mathematical expressions have been constructed that allow the large-scale global convection characteristics of the high-latitude ionosphere to be reproduced. The model contains no discontinuities in the ion convection velocity and as such should be useful in F region chemical models. The number of variables in the model allow such features as the dayside throat and the Harang discontinuity to be modeled. The applicability of the model to magnetospheric physics is limited by the exclusion of large-magnitude small-scale flow features associated with discrete arcs and by the inability of the model to produce separate flow cells at the same local time.

The Rice University convection model is applied to the early main phase of the July 29, 1977 magnetic storm through a computer implementation that self-consistently calculates electric fields, currents, and plasma distributions and velocities in the inner magnetosphere/ionosphere system. On the basis of solar wind parameters and AL index as inputs, the model predicts the injection of plasma sheet plasma to form a substantial storm time ring current whose total predicted strength agrees with the observed Dst index. The possibility that the magnetic field may be sufficiently inflated to make 60 deg field lines extend to the outer magnetosphere is examined. In the model, distortion of the plasma sheet inner edge by magnetospheric compression associated with the sudden commencement temporarily disrupts the normal Birkeland current pattern. Normal Birkeland currents and shielding reassert themselves in about an hour.

The currents and electric fields in the area covered by the poleward set of field-aligned currents are computed by means of a quantitative model, in which the high-latitude band considered carries most of the poleward electrojet, and lies poleward of the region covering the inner magnetosphere and corresponding ionosphere. The Birkeland current and aurorally enhanced conductivity are assumed to be uniformly distributed across the band. The time-dependent conductivity model used is based on electron fluxes and mean energies measured from the S3-2 satellite. Joule heating of the upper atmosphere was found to be about 2 x 10 to the 11th W during the substorm period. Model values for the strength of the electrojet and the amount of Joule heating agree to within about 20% with results based on a simple Cowling conductivity band.

The convection potential drop across the polar cap is computed from data obtained on high-inclination low-altitude satellites. Potential measurements are correlated with various combinations of parameters measured simultaneously in the upstream solar wind. Most of the potential drop is successfully predicted by merging theory, although a significant background potential drop of 35 kV does not depend on IMF parameters and is attributed to a process other than merging. Results indicate that small values of the IMF are amplified by a factor of 5-10 at the dayside magnetopause, which, when taken into account, improves correlations between IMF parameters and polar cap potential drop. Potential drop is better correlated with IMF parameters than with geomagnetic indices, due to nonlinear response of the magnetosphere affecting geomagnetic activity indices.

Results of a substorm simulation are used to investigate the penetration of substorm-associated electric fields into the plasmasphere. Near 4 earth radii in the equatorial plane, the time-dependent electric field model is characterized by eastward components in the dusk-midnight local time sector and westward components after midnight. With the exception of a small region just before dusk, the model predicts eastward electric field components throughout the daytime sector. The characteristic radial component is directed inward at all local times with the exception of a small region just after dawn. It is noted that these results compare favorably with available whistler and incoherent-scatter radar measurements obtained during magnetically disturbed periods. By assuming an initial plasmapause shape and by following the computed E x B drift trajectories of plasma flux tubes from that initial boundary, the short-term evolution of the plasmapause during the substorm-like event of September 19, 1976, is examined.

Results of a comprehensive computer simulation of the behavior of the earth's inner magnetosphere during a substorm-type event are reported. It is pointed out that the computer model self-consistently computes electric fields, currents, and plasma distributions and velocities in the inner-magnetosphere/ionosphere system; parallel electric fields and ionospheric neutral winds, however, are not included. The basic equations of the model are derived, and the inputs are described. An overview of the results is also given. The first appendix contains derivations of general, useful laws of bounce-averaged gradient, curvature, and E x B drifts in a plasma with isotropic pitch angle distribution. The second appendix gives a description of the numerical method used in the simulation.

Results of the computer simulation of the behavior of the inner magnetosphere during the substorm-type event of September 19, 1976, are discussed. The computed electric fields are found to compare satisfactorily with electric fields measured from S3-2, although there are detailed differences. The three general features on which the model and observations are in good agreement are (1) the magnitude and direction of the high-latitude electric field, (2) the degree to which the low-latitude ionosphere is shielded from the high-latitude convection electric field, and (3) the fact that the poleward electric field on the duskside is significantly larger, on the average, than the equatorward electric field on the dawnside. Simple formulas are presented that give rough estimates of global Joule heating rates from observable parameters.

Data from Atmosphere Explorer D for periods of strong northward interplanetary magnetic field are discussed. In the dayside magnetospheric cusp region energy time spectrograms of suprathermal positive ion fluxes exhibit a characteristic 'V' pattern as the spacecraft moves toward higher latitudes; that is, with the peak in the energy spectrum falling in energy and then rising again. Convection velocities follow this pattern closely with strong eastwest flows (with antisunward components) occurring in the equatorward half of the 'V' and significant sunward flows occurring in the poleward half of the 'V'. These patterns can be understood qualitatively in terms of a model of ionospheric electric potential produced by the known dependence of Birkeland current densities on magnetic activity.

The variation of proton average energy with latitude in the cusp has been suggested as an indicator of the means of particle entry. If magnetic merging is the principal means of particle entry, the proton average energy should fall with increasing latitude; if diffusion is the principal means, the average energy should first fall and then briefly rise as a function of latitude, showing a 'V' signature. In the present work, 60 selected cusp passes of the AE-D satellite for which IMP-J interplanetary or magnetosheath magnetic-field data were available were examined. About one-third of these passes showed clear or likely merging-type energy dispersions; a third showed clear or likely V-type energy dispersions, and a third showed unclear or no energy dispersions. The results are strongly correlated with the IMF - the merging signatures are associated with southward IMF and the V signatures with northward IMF. Unclear cases are associated with unsteady or weakly northward IMF.

Data from Atmosphere Explorer D for periods of strong northward interplanetary magnetic field show the following characteristic behavior in the dayside magnetospheric cusp region: energy-time spectrograms of suprathermal positive ion fluxes exhibit a characteristic 'V' pattern as the spacecraft moves toward higher latitudes; that is, with the peak in the energy spectrum falling in energy and then rising again. Convection velocities follow this pattern closely with strong east-west flows (with antisunward components) occurring in the equatorward half of the 'V' and significant sunward flows occurring in the poleward half of the 'V'. These patterns can be understood qualitatively in terms of a model of ionospheric electric potential produced by the known dependence of Birkeland current densities on magnetic activity.

Results are presented for an investigation of rapid subauroral ion drift features using data obtained over a nearly five-year period from the ion RPA/drift meter on Atmosphere Explorer-C. These latitudinally narrow features are found to be confined predominantly to the local time sector between 18:00 and 02:00 hr. They occur either singly or as multiple events, one of which almost always straddles the equatorward edge of the auroral zone. Their occurrence probability suggests a dependence of magnetic substorm activity.

Total ion concentration and ion drift velocity data from Atmospheric Explorer C are used to examine the convective flow of plasma in the vicinity of the midlatitude F region ionization trough. It is found that the low-latitude portion of the premidnight trough is a region of slow eastward plasma drift while the poleward portion is characterized by generally westward flow. The reversal from eastward to westward flow is described. It is proposed that the plasma flux tubes within the premidnight trough have corotated eastward into the evening local time sector from dusk before stagnating and flowing westward at higher latitudes in the nighttime sector. The role of F region recombination in trough formation is considered.