Modeling equatorial ionospheric electric fields

This paper gives a brief overview of the processes responsible for the equatorial electric field, and reviews relevant modeling work of these processes, with emphases on basic aspects and recent progress. Modeling studies have been able to explain most of the observed features of equatorial electric fields, although some uncertainties remain. The strong anisotropy of the conductivity and the presence of an east-west electric field lead to a strong vertical polarization electric field in the lower ionosphere at the magnetic equator, whose magnitude can be limited by plasma irregularities. Local winds influence the structure of the equatorial polarization field in both the E and F regions. The evening pre-reversal enhancement of the eastward electric field has been modeled by considering a combination of effects due to the presence of a strong eastward wind in the F region and to east-west gradients of the conductivity, current, and wind. Models of coupled thermosphere-ionosphere dynamics and electrodynamics have demonstrated the importance of mutual-coupling effects. The low-latitude east-west electric field arises mainly from the global ionospheric wind dynamo and from the magnetospheric dynamo, but models of these dynamos and of their coupling have not yet attained accurate predictive capability.     

Plasmaspheric zonal electric fields and coupling fluxes from the Dunedin VLF Doppler experiment,for 180°E,L ≈ 2.3, at solstice and equinox

Group delays and Doppler shifts from ducted whistler-mode signals are measured using the VLF Doppler experiment at Dunedin, New Zealand (45.8°S, 170.5°E). Equatorial zonal electric field and plasmasphere-ionosphere coupling fluxes are determined for L ≈ 2.3 at June solstice and equinox during magnetically quiet periods. The general features of the electric field measured at Dunedin agree with those predicted from ionospheric dynamo theory with a (1,−2) tidal component. Some seasonal variations are observed, with the electric field measured during equinox being smaller and predominantly westward during the night. The electric field at June solstice is also westward during the evening and for part of the night, but turns sharply eastward during the pre-dawn and dawn period at the duct entry site. The June electric field appears to follow a diurnal variation whereas the equinox electric field shows a possible 4-hourly periodic variation. Seasonal variations in the neutral wind pattern, altering the configuration of the ionospheric dynamo field, are the probable cause of the seasonal differences in the electric field. The seasonal variation of the coupling fluxes can be explained by the alteration of the E x B drift pattern, caused by the changes in the electric field.     

Dynamo action in sporadic-E formation

A theoretical study of the effects of background winds on wind shear-produced sporadic-E layers requires an account of the dynamo electric fields which result from the plasma motion produced by these winds. When a sporadic-E layer is carried downward by a descending wind shear the final height of the layer may vary by some 10 km depending on the background wind. Due to a loading effect on the dynamo, for a given background wind, the final height may vary by about 6 km depending on the degree of ionization in the E-layer in the magnetically conjugate hemisphere. The time scale for dispersal of a sporadic-E layer by drifts in wind-induced electric fields may be as little as 2 h.     

Low latitude electrodynamic plasma drifts : a review

Radar and radio measurements have provided detailed information on the dependence of F-region electrodynamic drifts on height, season, solar cycle and magnetic activity. Recently, satellite ion drift and electric field probes have determined the variation of low latitude ionospheric drifts over a large range of altitudes and latitudes. The general characteristics of the quiet time plasma can be explained as resulting from E- and F-region dynamo and interhemispheric coupling processes. The low latitude and equatorial zonal and upward/poleward components of the plasma drift respond differently to geomagnetic activity. Disturbance dynamo effects are responsible for the drift perturbations following periods of enhanced magnetic activity. The prompt penetration of high latitude electric fields to lower latitudes produces large perturbations on the upward/poleward drifts, but has no significant effect on the low latitude and equatorial zonal drifts. A number of processes such as ‘overshielding’, ‘fossil wind’ and magnetic reconfiguration were suggested as being responsible for the direct penetration of high latitude electric fields to lower latitudes. Detailed low latitude and global numerical models were used to study the characteristics of low latitude and equatorial plasma drifts and their response to changes in the polar cap potential drop or in the high latitude field-aligned currents. These models can reproduce the latitudinal variation of the perturbation electric fields and their diurnal variations, but are still unable to account for several aspects of the experimental data as a result of the complexity of the high latitude and magnetospheric processes involved.     

A numerical model of the ionospheric dynamo—I. Formulation and numerical technique

A new method of numerically solving a suitably formulated ionospheric wind dynamo equation for electrostatic potential and field is developed. Unlike in many other dynamo models, the upper boundary does not exist and the formulation asymptotically approaches the equatorial boundary condition. Therefore, it naturally incorporates the symmetric, asymmetric E- and F-region dynamo actions in any given ionosphere and any given global or local wind field. It also enables the equation to be posed as an initial value problem and solved numerically using an efficient, accurate, stable and fast integration method of ordinary differential equations. The numerical technique can be extended to compute three dimensional dynamo-generated electric currents in the ionosphere.     

A numerical model of the ionospheric dynamo—III electric current at equatorial and low latitudes

Using the general dynamo model and its special cases derived in a previous paper, the distributions of three dimensional electric current density in a magnetic meridional plane in the equatorial and low latitude ionosphere are computed. The winds generating the ionospheric dynamo are tide-like and locally periodic, similar to those in an internal gravity wave. Very large (several μA m⁻²) field-aligned current density is obtained in the equatorial region at places of sharp vertical gradients of the wind velocity. The currents generated by locally periodic winds of latitudinal wavelength less than several hundred kilometers do not significantly affect the normal equatorial electrojet.     

Electric field sources in the quiet plasmasphere from whistler observations

Existing evidence for the ionospheric dynamo being the source of quiet time electric fields in the plasmasphere is reviewed. Part of a 24 h set of whistler data recorded continuously at Sanae, Antarctica (L = 4), during quiet magnetic (average K_p = 1) is analysed to obtain westward electric fields in the equatorial plane. These electric fields are examined as a function of L-value in order to infer their source. It is found that for periods of outward flow of plasma during the noon-midnight local time period, the electric fields are consistent with the dominant source being the ionospheric dynamo. There is some evidence that during the evening period of inward flow the electric fields are magnetospheric in origin, although this could also be consistent with a refined dynamo model. The observed whistler duct convection patterns do not fit either of two theoretical models, which invoke a magnetospheric field but not a dynamo field.     

On the global atmospheric electrical circuit

A model of global atmospheric electric circuit has been developed, which incorporates the pollution due to aerosol particles of anthropogenic and volcanic origins and the ionization caused by the coronal discharges due to intense electric field beneath thunderclouds in the atmospheric layers close to the earth's surface. Also the effects of solar activity and of Stratospheric Aerosol Particles (SAP) on the parameters of the circuit have been studied. The global distribution of Aerosol Particles (AP) of man-made origin has been assumed on the basis of world population density and that of volcanic origin on the basis of global distribution of volcanic activities. The thunderstorms are assumed to be the current generators of the global circuit. The latitudinal, longitudinal and height variations of atmospheric conductivity have been calculated from the known variations of ionization caused by cosmic rays, the radio-active emanations and intense electric fields beneath thunderclouds (over continents), along with the assumed variation of AP concentration.On increasing the AP concentration over the northern hemisphere by about 5 times that of the southern hemisphere, the ionospheric potential increases to about 6% and electric field increases in non-mountainous regions more than that in the high mountain regions. A large increase in SAP serves to increase the global resistance while both global current and ionospheric potential decrease. However, for low SAP concentration, the ionospheric potential increases. The SAP affect the electrical structure of the stratosphere without much influencing the troposphere except in the volcanically active regions, where conductivity is low due to high AP concentration. The major influence on the global atmospheric electric circuit occurs in response to solar activity. In the absence of local effects, the calculated variation of ionospheric potential is 14%. However, in real atmosphere (where both local and global processes act simultaneously), the ionospheric potential is found to vary by only about 3%. This has little effect on the ground electrical properties where more than 30% of the variations have been found to be caused by the local effect.     

Midlatitude measurements of the ionospheric electric field during the ALADDIN programme

Three measurements of ionospheric electric field were made during the 24 h ALADDIN rocket programme at Wallops Island (37°50′N, 75°29′W) on June 29–30, 1974. The first of these used a double probe instrument, flown at 1500 Local Solar Time, and the second and third measurements were made by barium cloud releases at evening and morning twilight. These three electric field vectors have been compared with the predictions of a number of models of electric field due to the dynamo effects of various atmospheric tides, and also of a possible magnetospheric origin. On the assumption that the measurements were made at a location equatorward of the afternoon convergence and poleward of the morning divergence in the electric field patterns related to the Sq current cystem, Stening's model of the diurnal variation of the electric field induced by the (1, −2) tidal model at the time of the Summer solstice correctly predicts the directions of the observed electric field. Forbes and Lindzen's model, incorporating the three major propagating tidal modes as well as the evanescent (1, −2) mode, also bears an acceptable relationship to the ALADDIN electric field directions. The ALADDIN E-field magnitudes are comparable with those obtained by ground-based observations (incoherent scatter) from Millstone Hill and from Saint Santin but are about half of Stening's model values, and three times those of Forbes and Lindzen.While the Millstone Hill E-field directions are compatible with the ALADDIN observations, Saint Santin E-field directions, at the same latitude but 75° difference in longitude, are distinctly different from ALADDIN, implying that longitudinal differences are significant.     

Responses of atmospheric electric field and air-earth current to variations of conductivity profiles

A global circuit model is constructed to study responses of air-earth current and electric field to a variation of atmospheric electrical conductivity profile. The model includes the orography and the global distribution of thunderstorm generators. The conductivity varies with latitude and exponentially with altitude. The thunderstorm cloud is assumed to be a current generator with a positive source at the top and a negative one at the bottom. The UT diurnal variations of the global current and the ionospheric potential are evaluated considering the local-time dependence of thunderstorm activity. The global distributions of the electric field and the air-earth current are affected by the orography and latitudinal effects. Assuming a variation of conductivity profile, responses of atmospheric electrical parameters are investigated. The non-uniform decrement of the conductivity with altitude increases both the electric field and the air-earth current. The result suggests a possibility that the increment of the electric field and the air-earth current after a solar flare may be caused by this scheme, due to Forbush decrease.     

A new role for the plasmaphere in “quiet” ionospheric electrodynamics

Global scale longitudinal gradients of pressure in the plasmasphere may be formed naturally by ionospheric processes, or caused by electrostatic fields of ionospheric dynamo origin. It is shown that plasmaspheric gradients of pressure, orthogonal both to the magnetic field (B) and to grad B, generate geophysically significant field-aligned currents. Considering the ionosphere and plasmasphere as a coupled electrodynamic system, these currents alter non-negligibly the self-consistent ionospheric electric field and current. Criteria are established for this coupling mechanism (a kind of plasmaspheric impedance) to be significant. This has implications for the relationships of ionospheric electric fields and currents, F-region drifts, and magnetic variations, due to upper atmosphere tides and winds.     

Magnetosheath,magnetopause and low latitude boundary layer research, 1987–1989

Important developments over the past two years (1987–1989) in studies of the magnetosheath, the magnetopause and the magnetopause boundary layers are reviewed. The pace of research shows no sign of abating, but progress has become focused on certain topics. These include flux transfer event reconnection and the relevance of solar wind dynamic pressure variations for magnetopause and polar cusp dynamics. A promising development in both these areas is the use of ground observations under the polar cusp ionosphere for monitoring solar wind-magnetosphere coupling processes at the magnetopause. In contrast to these progress areas, the magnetosheath and the low latitude boundary layer remain topics which are being unjustifiably neglected. The controversies and the key remaining questions are highlighted and suggestions are made on how to proceed further.     

Electrical structure in the high-latitude middle atmosphere

Middle atmosphere electrodynamics at high latitudes differs significantly from the normally assumed picture of a passive region through which electric fields of external origin couple. Large Vm ⁻¹ electric fields, both horizontal and vertical, have been observed within bounded regions of the upper stratosphere and lower mesosphere. They seem to occur only in regions where the electrical conductivity is a few times 10⁻¹⁰ S m⁻¹ or less and appear to be current limned. While low conductivity is necessary, it is not a sufficient condition for occurrence. The observed large horizontal electric fields were found to be anticorrelated with the local neutral wind. However, a generation mechanism of these electric fields is as yet unknown but must involve space charge separation rather than dynamo effects. Large variations in the conductivity were also observed to occur with fluctuations in magnetic activity, and these were found to be consistent with measured variations in energy deposition during auroral phenomena. Theoretical concepts of mapping of electric fields downward from the thermosphere along equipotential magnetic field lines were shown to hold qualitatively in the D-region at the mV m⁻¹ level. Perturbations affecting such models were determined to be small.     

Ionospheric electrodynamic model with an eccentric dipole geomagnetic field

A computer model of ionospheric electrodynamic processes using an eccentric dipole (ED) for the geomagnetic field has been developed. This is a development from existing models which are based on the centred dipole (CD) coaxial with the geographic axis. The ED dynamo model introduces or modifies the effects of hemispherical asymmetry and longitudinal variation in the dynamo processes through two explicit parameters—the geomagnetic field intensity and the length of the field lines. These parameters of the ED field have been quantified and displayed. An additional contribution to the above effects comes implicitly from the ionospheric parameters—plasma density and atmospheric tidal winds—which become asymmetric relative to the ED dip equator. The integrated effect of the geomagnetic and ionospheric parameters produces significant variation in the field line integrated ionospheric conductivity. The ED dynamo model shows that the peak height of the equatorial electrojet (EEJ) moves by over 2 km and height profiles of the EEJ display strong hemispherical asymmetry.     

The F-region dynamo

This paper reviews the theory of the F-region dynamo which drives about 10–15% of the total mid-latitude ionospheric current by day, and the major part at night (Section 2). Polarization fields associated with the dynamo cause marked effects in the night-time F-region, notably the mean eastward wind (Section 3). The paper also discusses the equipotentiality of geomagnetic field lines (Section 4 and Appendix) and the question of location of Sq and L current systems (Section 5).     

Global dynamics of the magnetosphere for northward IMF conditions

Ground-based and spacecraft observations of polar cap geophysical phenomena during periods of northward interplanetary magnetic field (IMF) show specific patterns of electric fields, field-aligned currents, aurora and particle precipitation. These are basically different from those when the IMF is southward. The total combination of observational data for northward IMF indicates rather a closed magnetosphere. This topology has led to the formation of a specific convection pattern in the distant plasma sheet. As different theoretical studies show, the connection of the IMF to geomagnetic flux tubes poleward of the cusp region may serve as the driving mechanism for plasma sheet convection and as the dynamo of current systems. Unfortunately, the direct observations of processes in the distant magnetosphere are too scarce either to accept or reject the concept of a closed magnetosphere. There are also some experimental data that are inconsistent with the closed magnetosphere topology. Definitive open or closed models must await future measurements.     

Multiple pathways across past landscapes: circuit theory as a complementary geospatial method to least cost path for modeling past movement

The incorporation of geospatial technologies in archaeology has resulted in productive advances in the analysis of past behavioral processes. Archaeologists have relied on the Geographic Information Systems (GIS) application of cost surface analysis and the computation of least cost paths (LCP) to study movement, one key social process. Recent research has identified limitations with LCP for modeling nonhuman species movement, notably the inability of LCP to accommodate multiple pathways. Archaeology must consider the implications of these critiques for models of past human movement. In this paper, I apply a different approach, circuit theory modeling, enacted through the program Circuitscape, to an archaeological case study previously analyzed with LCP modeling, travel to a regionally significant ceremonial earthwork center during Late Prehistory (ca. AD 1200–1600) in the Northern Great Lakes. Through this comparative analysis, I evaluate the advantages and disadvantages of these methods for modeling past movement. The results suggest LCP modeling and circuit modeling offer archaeology complementary geospatial methods for conceptualizing past mobility. Combining circuit theory and LCP allows archaeologists to produce richer models of past movement by appreciating scenarios where multiple pathways are important as well as scenarios where optimum single travel routes have priority.     

Electromagnetic scattering from single electrified raindrops

On the basis of recent experimental and theoretical works carried out to study the shape of electrified raindrops, their radiophysical and polar diagram characteristics are mathematically modelled. Both a numerical implementation and analyses are made using models for liquid drops, and considering factors determining their shape such as electric fields, surface charge, surface tension, hydrostatic pressure, dynamic pressure and sizes, as well as the angle of the drop's axis of symmetry relative to the vertical. Analytical approximations for the generators of the shapes of such drops are obtained.