The morphology of dayside Birkeland currents

Intense (10⁵ A) electric currents flow into and from the Earth's two polar ionospheres near magnetic noon. These currents, called Birkeland or magnetic field-aligned currents, are the agent by which momentum couples from the flowing solar wind plasma to drive plasma motions in the high latitude ionosphere. Coupling is strongest when the interplanetary magnetic field (IMF) has a southward component and when this occurs there exist two principal regions of Birkeland current near magnetic noon called the region 1 and the cusp systems. We present a simple model bringing theoretical order to the many patterns proposed previously for the morphology of these dayside Birkeland currents as observed by orbiting satellites in the topside polar ionosphere. Specifically we show that the cusp Birkeland current system is not a latitudinally separate region but is instead the extension in longitude of the region 1 Birkeland current from either dawn or dusk; which particular one depends on the sign of the east-west (Y) component of the IMF. The presence of an IMF Y-component therefore leads to two region 1 current systems near magnetic noon, with the poleward one being that previously called the ‘cusp’ system.     

The influence of ionization and recombination processes on the dynamics of ionospheric plasma clouds

The dynamics of a one-dimensional ionospheric irregularity interacting with the magnetosphere is studied by numerical simulation. The polarization electric field produced by charge separation within the irregularity propagates along magnetic field lines with the Alfvén velocity V_A and drives polarizational and field-aligned currents in the magnetosphere. Their values and localization are controlled by motion and deformation of the irregularity resulting from its electrostatic coupling to the background ionosphere. The pattern of the field-aligned currents varies with time and depends primarily on gradients of the polarization electric field. The latter is controlled by the ambient electric field, diffusion, recombination process, intensity of the initial perturbation, etc. Feedback effect of the magnetospheric conductance on the development of the irregularity is examined.     

Modelling the penetration of thundercloud electric fields into the ionosphere

In a previous paper, we considered the penetration of DC thundercloud electric fields E into the ionosphere and also into the region between the ionosphere and the ground (Velinov and Tonev, 1994). In the present paper, we extend the analysis by making a more precise approximation of the electric conductivity profiles by 5–10 piecewise exponential functions of altitude instead of the two functions used up to now. This allows a much more realistic representation of the atmospheric conductivity profile. Besides, Maxwell's equations are solved for more general boundary conditions, taking into account that the electrosphere is not a perfect conductor. This leads to the appearance not only of the transverse E_r (as had been assumed until now), but also of the geomagnetic field-aligned E_z component of the penetrating thundercloud electric fields. The computations show that both E_r and E_z cause significant variations of the electron density profiles N(z) in the ionosphere.     

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.     

A unified view on convection and field-aligned current patterns in the polar cap

During the last two decades measurements of polar cap ionospheric electric fields and currents, field-aligned currents, and global auroral forms have been made from ground-based and space-based platforms. An attempt is made to unify these observations into a large-scale view of polar phenomena. In this view, plasma convection patterns and the corresponding electrodynamics in the polar region can consistently be ordered by the orientation of the interplanetary magnetic field (IMF). The different patterns of the electric potential and of field-aligned currents depend on where the main interaction between the terrestrial and interplanetary fields occurs, on the morning or evening side of the central polar cap, or on the dayside portion of the ‘closed’ cusp region, or on the nightside portion of the ‘open’ cusp region. One of the essential elements of this unified view is that it is possible to account for various convection patterns ranging from the four-cell pattern (during periods of strong northward IMF and B_y ~ 0), to the three-cell pattern (B_z > 0 and |B_y| 2> 0), to the conventional two-cell pattern (B_z < 0) with its possible deformation into a convection throat near the dayside cusp (during southward IMF). We also discuss the way in which the complicated field-aligned current systems can consistently be accounted for in terms of these convection patterns.     

Electric field generation at the magnetospheric boundary for northward IMF

We discuss three different processes which generate electric fields at the magnetopause during northward interplanetary magnetic field (IMF) conditions. These are (1) Petschek-type magnetic field reconnection, (2) magnetic field diffusion, and (3) viscous-like interaction resulting from the Kelvin-Helmholtz instability. For northward IMF all three processes lead to the formation of a boundary layer on closed magnetic field lines adjacent to the magnetospheric boundary. The thickness of the boundary layer depend on Petschek's parameter in the first case, the magnetic Reynolds number in the second case, and an effective Reynolds number in the third case. In each case coupling between the boundary layer and the ionosphere occurs via field-aligned currents. These field-aligned currents result from the penetration into the polar ionosphere of the electric field generated at the magnetospheric boundary. These currents are closed by a transverse current in the boundary layer and the associated Lorentz force causes a decrease of the kinetic energy of the solar wind plasma inside the boundary layer. As a result of this velocity decrease the thickness of the boundary layer increases on both flanks of the magnetosphere near the equatorial plane. The convergence of the boundary layer on the dawn and dusk sides leads to antisunward plasma flow in the magnetospheric tail.     

The ionospheric plasma flow and current patterns of travelling convection vortices: a case study

Following a short duration density enhancement in the solar wind, observed by the AMPTE/IRM spacecraft, transient disturbances appeared in the polar ionosphere in the prenoon local time sector which were identified as Travelling Convection Vortices (TCV). This event has been studied intensively by combining radar and magnetometer observations. EISCAT radar was operated in the special programme U.K.-POLAR which provides F-region plasma parameters from invariant latitudes around 72° at a rate of one sample per 15 s. The combined data set provides a detailed picture of the drift pattern of the plasma and the three-dimensional current distribution. There are two Hall current eddies drifting westward at a speed of 0.15° s⁻¹. The leading one circulating clockwise is associated with a downward field-aligned current and the oppositely circulating eddy with an upward current. The ionospheric conductivity seems to be enhanced in the leading vortex compared to the trailing, although the latter is connected to an upward field-aligned current. Still unexplained is the mechanism generating the electric field which drives the vortices. The direction of the electric field observed in the ionosphere is opposite to that expected if the source were a compression of the magnetosphere.     

A quasi-DC model of electric fields in the ionosphere-ground region due to electrified clouds

The penetration of the quasi-DC electric fields, E due to electrified clouds, into the middle- and high-latitude ionosphere is theoretically studied during the initial stage of cloud charge separation. The electrification process is characterized by a source function S(t), whose variations are on the time scale of the relaxation process in the cloud. A first-order approximation solution for the time variation of E is obtained as an explicit function of the electric field in the steady-state case. Some features of the time variations of E at different altitudes (in the ionosphere and also in the atmosphere) are investigated, depending on the source function and the relaxation time constants. This result can also be applied for slowly electrifying clouds which do not produce lightning.     

Magnetopause plasma transients: signatures in the auroral ionosphere

The magnetopause and adjacent boundary layers of the Earth's magnetosphere play important roles in transferring momentum and energy from the solar wind to the magnetosphere-ionosphere system. The details of the different boundary processes, their ionospheric signatures and relative importance are not well known at present. Particle precipitation, field-aligned current, auroral emission, ionospheric ion drift and ground magnetic perturbations are among the low-altitude parameters that show signatures of various plasma processes in the LLBL and the magnetopause current layer. Magnetic merging events, Kelvin-Helmholtz waves, and pressure pulses excited by the variable solar wind/magnetosheath plasma are examples of boundary phenomena that may be coupled to the ionosphere via field-aligned currents. In this paper, attention is focussed on a specific category of auroral activity occurring in the cusp/cleft region predominantly during the southward directed interplanetary magnetic field (IMF). Co-ordinated observations from the ground and satellites in polar orbit have been used to study the temporal/spatial development of the events in relation to the background patterns of particle precipitation and ionospheric convection as well as the field-aligned current and ion drift characteristics of the individual events. The auroral phenomenon is characterized by a sequence of elongated forms moving laterally into the polar cap. Spatial scales of major events repeating every 5–10 min are ∼200 km (N-S) times 300–1000 km (E-W). Smaller scale auroral structures with more irregular occurrence rates are observed at times. The preliminary evidence suggests that the motion pattern is regulated by the IMF orientation, that is, the direction of longitudinal motion along the polar cap boundary is determined by the IMF B_Y polarity. The examples reported here occurred within 1000–1400 MLT, near the zero point potential line separating the morning and post-noon convection cells. During nonzero IMF B_Y the auroral structures are associated with channels of enhanced zonal ionospheric ion flow and Birkeland current sheets of opposite polarity, imbedded within the larger scale IMF B_Y-related cusp-mantle current system. These characteristics are discussed in relation to model predictions of ionospheric signatures of magnetopause plasma transients, with particular emphasis placed on impulsive magnetic merging events.     

Effect of the large-scale convection electric field structure on the formation of thin ionization layers at high latitudes

A three-dimensional simulation of the high-latitude ionosphere was applied to investigate the geographical distribution of E-region thin ionization layers which may be formed by the action of the convection electric field. The simulation model computes the ion densities (O⁺, O⁺₂, N⁺, N⁺₂, NO⁺, Fe⁺), and temperatures as a function of altitude, latitude, and longitude. The stationary state momentum and continuity equations are solved for each ion species, then the energy equation is solved for electrons, neutrals, and a generic ion having the mean ion mass and velocity. The various electric field patterns of the Heppner and Maynard [(1987) J. geophys. Res.92, 4467–4489] convection electric field model were applied and the ionization density pattern was examined after a time sufficient for the formation of thin layers (≈2000 s). It was found that large areas of thin ionization layers were formed for each of the electric field patterns examined. Southward IMF B_z conditions resulted in thin layers forming in the pre-midnight sector in the latitude range north of about 70° to about 80°, and after midnight between 60 and 70°. For northward B_z conditions, the layers were mainly in the pre-midnight sector and covered a latitude range from about 60 to 80°.     

What determines the latitudinal extent of the auroral acceleration region?

Characteristic scales associated with auroral precipitation are investigated on the basis of quasistatic magnetotail models, resistive MHD simulations of magnetotail dynamics, and a general relation between parallel electric fields and velocity shear. Since the inverted-V precipitation region of discrete auroras (on the dusk side) is associated with upward flowing, region 1, currents, we investigate the distribution of these currents first. The overall distribution of region 1 type field-aligned currents and their dynamic changes can be explained by characteristic scales in the magnetotail and their mapping to the ionosphere. The quiet time region 1 currents are associated with the decrease of tail flaring. Their overall extent in the north-south direction is closely related to the scale height of the cross-tail current. Dynamic region 1 currents are related to the velocity shear of earthward flow, which can be generated by a tail instability. In that case the peaks of the enhanced region 1 currents are found to map closer to midnight and to lower latitudes than the quiet region 1 currents, consistent with average observations [Iijima and Potemra (1976a), J. geophys. Res.81, 2165]. On the basis of a general relation between parallel electric fields and ‘slippage’ in the plasma transport [Hesse and Schindler (1988), J. geophys. Res.93, 5559; Schindleret al. (1991), Astrophys. J.380, 293], we make estimates of the spatial extent of nonideal regions, where parallel electric fields may exist. For a plausible model of substorm reconfiguration, we find a latitudinal extent of about 7 km for a time scale of 1 min and a integrated parallel electric field of 5 kV. The length scale is proportional to the time scale. The sign of this parallel electric field is consistent with downward acceleration of electrons on the dusk side. The spatial extent of the parallel electric field region is independent of the microscopic generation mechanism if the time scale and the characteristic parallel potential difference (i.e. the integrated parallel electric field) are independent of this mechanism.     

Topside and intelhemispheric ion flows in the mid-latitude plasmasphere

A modelling study has been carried out of field-aligned ion flows in the topside ionospheres of conjugate hemispheres under solstice conditions at mid to low latitudes. In the model calculations coupled time-dependent O⁺, H⁺ and electron continuity, momentum and heat balance equations are solved along dipole magnetic field lines at L = 1.5 and 3.0 Sunspot medium and sunspot minimum atmospheric conditions are considered.It has been found that thermal coupling between conjugate hemispheres gives rise to strong flows of O⁺ in the topside ionosphere of the summer hemisphere that are directed upwards at conjugate sunrise and directed downwards at conjugate sunset. At conjugate sunrise in the winter hemisphere there is a small upward-directed signature in the O⁺ field-aligned flux; there is no observable signature in the O⁺ field-aligned flux in the winter hemisphere at conjugate sunset. There are strong upward and downward flows of O⁺ at local sunrise and local sunset, respectively, in both the summer and winter hemispheres.At both L = 1.5 and 3.0 the 24 h time-integrated interhemispheric H⁺ flux is in the direction summer hemisphere to winter hemisphere. At L = 1.5 its magnitude is in good agreement with the magnitude of the 24 h time-integrated plasma (O⁺ + H⁺) field-aligned flux at 1000 km altitude; there are no such agreements at L = 3.0.A study of the roles played by the individual terms of the O⁺ momentum equation has demonstrated the complex structure of momentum balance. Certain of the terms may be orders of magnitude greater than the combined total of the individual terms, i.e. the O⁺ field-aligned flux.     

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.     

Comparison of ionospheric electrical conductances inferred from coincident radar and spacecraft measurements and photoionization models

Height-integrated electrical conductivities (conductances) inferred from coincident Sondrestrom incoherent scatter radar and DMSP-F7 observations in the high-latitude ionosphere during solar minimum are compared with results from photoionization models. We use radar and spacecraft measurements in combination with atmospheric and ionospheric models to distinguish between the contributions of the two main sources of ionization of the thermosphere, namely, solar UV/EUV radiation and auroral electron precipitation. The model of Robinsonet al. (1987, J. geophys. Res.89, 3951) of Pedersen and Hall conductances resulting from electron precipitation appears to be in accordance with radar measurements. Published models of the conductances resulting from photoionization that use the solar zenith angle and the solar 10.7-cm radio flux as scaling parameters are, however, in discrepancy with radar observations. At solar zenith angles of less than 90°, the solar radiation components of the Pedersen and Hall conductances are systematically overestimated by most of these models. Geophysical conditions that have some bearing on the state of the high-latitude thermosphere (e.g. geomagnetic and substorm activity and a seasonal variation of the neutral gas distribution) seem to influence the conductivity distribution but are to our knowledge not yet sufficiently well modelled.     

The day-to-day variability of Sq in geomagnetically conjugate areas

The correlation of daily ranges of Sq(H) has been found to be slightly but significantly greater for pairs of stations which are approximately at magnetically conjugate locations than for other, non-conjugate pairs of similar distance of separation. It is suggested that the effect is due to field-aligned currents.     

Small-scale ionospheric structures associated with mid-latitude spread-F

Slant-F traces on ionograms recorded by a modern ionosonde in a sunspot-minimum period have revealed the existence of field-aligned irregularities at times of spread-F occurrence. This appears to be the first investigation in a mid-latitude region around 36° (geomagnetic) to detect these irregularities at F₂-region heights using an ionosonde. Although such traces were observed frequently near sunspot minimum they were seldom recorded for periods close to sunspot maximum. Also, for a specific spread-F event in August 1989, both the ionograms from the modern ionosonde and scintillations of 150 MHz transmissions from a Transit satellite indicate the existence in the ionosphere of periodic structures (period around 11 min). The scintillation recording also included rapidly fading signals indicative of small-scale structures. The satellite had a path close to the magnetic meridian which passed through the recording station (Brisbane, Australia). Because of the enhanced signal fluctuations in the scintillation recording on this occasion it seems likely (with the support of other evidence on the ionograms) that the small-scale structures present were field-aligned.     

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.     

Boundary populations in the polar caps

The morphology of precipitating particles, measured at low altitude in the polar regions, varies systematically with the strength and direction of IMF B_z and with solar wind speed V_sw. We use particle data taken onboard the DMSP satellites to determine these variations. Both individual satellite passes during the storm/quieting period of 26 and 27 August 1990, and statistical maps compiled from a data base over 4.5 yr are presented. We focus attention on those magnetospheric populations that have magnetosheath characteristics, the boundary populations. We show that the precipitating ion boundary population, whose down-coming spectra can be fitted to streaming Maxwellians, expands from a region confined near the dayside cusp for southward IMF, to a thick, annular region, including the dayside cusp, for northward IMF. The expansion in local time is inhibited by increasing solar wind speed. Boundary electrons behave somewhat differently. They have easier access to the polar regions and their variations have shorter spatial/temporal scale lengths than the boundary ions. For strongly northward IMF, intense, agitated boundary electrons can be found over all or part of the polar cap. Broad regions (up to ~ 100 km) of strongly accelerated electrons (several keV) that produce visible arcs are embedded in this population. Two features of the ion boundary population help identify its source. (1) The spectra of the boundary ions expanding into the polar cap exhibit field-aligned streaming, which, downtail, is toward the Earth. (2) The region into which the boundary ions expand best maps magnetically to a dawn-dusk cut across the neutral sheet, rather than to the low-latitude boundary layer. Therefore, we conclude that the immediate source for boundary ions in the polar regions during northward IMF is the plasma sheet boundary layer. These ions reach tail lobe field lines by convection whose direction when mapped to the ionosphere is sunward. Significant change in the topology of the magnetospheric magnetic field, and, in particular, the closing of high-latitude field lines, is not required to explain the data.     

Atmospheric winds calculated from diurnal changes in the mid-latitude ionosphere

Diurnal variations in the electron content (N_t) and peak density (N_m) of the ionosphere are calculated using a full time-varying model which includes the effects of electric fields, interhemispheric fluxes and neutral winds. The calculation is iterated, adjusting the assumed hourly values of neutral wind until a good match is obtained with mean experimental values of N_t and N_m. Using accurate ionospheric data for quiet conditions at 35°S and 43°S, winds are derived for summer, equinox and winter conditions near solar maximum and solar minimum. Solar maximum results are also obtained at 35°N. Changes in the neutral wind are found to be the major cause of seasonal changes in the ionosphere, and of differences between the two hemispheres. Calculated winds show little variation with latitude, but the winds increase by about 30% at solar minimum (in equinox and winter). The HWM90 wind model gives daytime winds which are nearly twice too large near solar maximum. The theoretical VSH model agrees better with observed daytime variations, and both models fit the observed winds reasonably well at night. Results indicate that modelling of the quiet, mid-latitude ionosphere should be adequate for many purposes when improved wind models are available. Model values for the peak height of the ionosphere are also provided; these show that wind calculations using servo theory are unreliable from sunrise to noon and for several hours after sunset.