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.     

Geomagnetic field variations at low latitudes and ionospheric electric fields

The solar cycle, seasonal and daily variations of the geomagnetic H field at an equatorial station, Kodaikanal, and at a tropical latitude station, Alibag, are compared with corresponding variations of the E-region ionization densities. The solar cycle variation of the daily range of H at either of the stations is shown to be primarily contributed to by the corresponding variation of the electron density in the E-region of the ionosphere. The seasonal variation of the ΔH at equatorial stations, with maxima during equinoxes, is attributed primarily to the corresponding variation of the index of horizontal electric field in the E-region. The solar daily variation of ΔH at the equatorial station is attributed to the combined effects of the electron density with the maximum very close to noon and the index of electric field with the maximum around 1030 LT, the resulting current being maximum at about 1110 LT. These results are consistent with the ionosphere E-region electron horizontal velocity measurements at the equatorial electrojet station, Thumba in India.     

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.     

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

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.     

Middle atmosphere and lower thermosphere processes at high latitudes studied with the EISCAT radars

The middle and upper atmosphere and the ionosphere at high latitudes are studied with the EISCAT incoherent scatter radars in northern Scandinavia. We describe here the investigations of the lower thermosphere and the E-region, and the mesosphere and the D-region. In the auroral zone both these altitude regions are influenced by magnetospheric processes, such as charged particle precipitation and electric fields, which are measured with the incoherent scatter technique. Electron density, neutral density, temperature and composition are determined from the EISCAT data. By measuring the ion drifts, electric fields, mean winds, tides and gravity waves are deduced. Sporadic E-layers and their relation to gravity waves, electric fields and sudden sodium layers are also investigated with EISCAT. In the mesosphere coherent scatter occurs from unique ionization irregularities. This scatter causes the polar mesosphere summer echoes (PMSE), which are examined in detail with the EISCAT radars. We describe the dynamics of the PMSE, as well as the combination with aeronomical processes, which could give rise to the irregularities. We finally outline the future direction which is to construct the EISCAT Svalbard Radar for studying the ionosphere and the upper, middle and lower atmosphere in the polar cap region.     

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.     

Magnetometer and incoherent scatter observations of an intense Ps 6 pulsation event

In the morning sector of 21 April 1985, during the recovery phase of a geomagnetic storm, a Ps 6 pulsation event was recorded by the EISCAT magnetometer cross in northern Scandinavia. Simultaneously, the EISCAT incoherent scatter radar measured E- and F-region plasma parameters with a latitudinal scanning program. Electric fields and height-integrated Hall and Pedersen conductivities are derived. Two-dimensional patterns of these quantities are constructed for one Ps 6 period. The conductance patterns closely resemble the typical auroral forms of eastward drifting Ω bands with low and high conductances at the northern and southern edges of the scanned area, respectively. From the equatorward region a tongue of high ionization extends poleward into the dark area. The location of the maximum southward current is slightly displaced towards the west from the centre of the conductance tongue. The east-west disturbance electric field points towards the tongue; the north-south fields are enhanced outside and reduced inside the high conductance region. As has been previously suggested, the observations can be explained with a model which superposes currents caused by conductance variations and electric fields. Both effects need to be taken into account for this event. The current structures move within a few degrees in the direction of the background E×B drift, but their speed is about 15% lower than the average F-region plasma drift.     

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.     

Altitude and latitude dependence of the equatorial electrojet

The effects of day-to-day or seasonal variation of altitude and latitude profiles of the Elayer plasma density in the equatorial ionosphere on equatorial electrojet (EEJ) structure are examined numerically using a self-consistent and high resolution dynamo model. It is found that variations in the E-layer peak altitude and amplitude and its gradient below significantly affect EEJ structure. For any realistic shape, the EEJ peak appears at or below the E-layer peak altitude. Distinct double peaks appear in the EEJ structure, such as revealed by rocket measurements, if the E-layer peak is above 105 km or the gradient is large, as when sporadic-E is present. The influence of the latitudinal variation of ionospheric field line integrated conductivities upon the amplitude and altitude of the EEJ peak is demonstrated.     

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.     

Solar activity effects on equatorial plasma bubble zonal velocity and its latitude gradient as measured by airglow scanning photometers

Using a new mode of scanning 630-nm photometer operation the zonal velocities of ionospheric plasma depletions were measured over Cachoeira Paulista in Brasil in two east-west planes tilted 30°N and 30° S with respect to zenith. The measurements cover a time period of approximately 2 years, from January 1988 to January 1990, a period marked by significant increase in solar activity of the ongoing cycle. The results have permitted a rather detailed evaluation of the local time and latitude variations in the zonal plasma bubble velocity as a function of solar activity. Although the mean trend in the velocity local time variation is a decrease from early evening to post-midnight hours, a strong tendency for velocity peaks is observed near 21 LT and midnight. The velocities as well as their height (latitude) gradients show perceivable increases with solar activity represented as sunspot numbers. The present results are compared with the ambient plasma velocities measured using the Jicamarca radar by Fejer el al. (1985), J. Geophys. Res. 90, 12249, with that measured on board the DE 2 satellite on the equatorial latitudes by Coley and Heelis (1989), J. geophys. Res. 94, 6751, and with various theoretical calculations, in an attempt to bring out the salient features of the plasma dynamics of the equatorial ionosphere.     

Data comparison of EISCAT and the energetic particle instrument on GEOS 2

Flux variations of high energy (E > 16 keV) electrons, measured by the energetic particle spectrometer onboard the geosynchronous satellite GEOS 2, have been compared to simultaneously observed electron density structures in the lower ionosphere, measured with the incoherent scatter facility EISCAT. The very localized E-region electron density enhancements caused by the precipitation allow a preliminary estimate of the location of the footpoint of the magnetic field line through GEOS 2. Various other results of the intercomparison are discussed.     

Applicability of the two-layer model for simulation of non-linear evolution of a large scale plasma cloud in the ionospheric F-region

The adequacy of the two-layer model of Lloyd and Haerendel for describing the behaviour of an ionospheric irregularity is verified by numerical simulation of large plasma cloud dynamics. The background ionosphere is approximated by a set of conductive layers with ion mobilities and concentrations corresponding to the real ionospheric conditions. Polarization electric field produces positive and negative image clouds, i.e. plasma density enhancements and depletions in each layer. Their intensity, form and orientation turn out to change with height depending on the local conditions. However, the drift and deformation of the released cloud slightly differ from the case when the ionosphere is characterized by constant, height averaged parameters, at least if altitude dependent neutral wind and photochemical processes are ignored.     

Incoherent scatter observations of thin ionization layers at Sondrestrom

High resolution incoherent-scatter observations of E-region thin (1–3 km) metallic ion layers are presented. Data were collected during three different periods from August 1990 to August 1991, in three different experimental modes. First, the antenna was directed vertically and the entire duty cycle was devoted to Barker coded multi-pulse [Zamlutti (1980) J. atmos. terr. Phys.42, 975–982] measurements to determine the densities and temperatures in the E-region with 300 m resolution. The second experiment measured the F-region electric field as well as the high resolution E-region densities. For the third experiment the antenna was scanned magnetic north-south while only the E-region densities were measured. The experiments were carried out on 16 different nights for a period of 4 h each night at a time near magnetic midnight. Thin ionization layers were observed on 12 of the 16 nights. The first experiment demonstrated that the thin layers are composed of a significant fraction of heavy metallic ions; assuming the layers are composed of a mixture of Fe⁺ and Mg⁺ a composition estimate of 63% Fe⁺ was obtained in one example. The second experiment investigated the relationship between the direction of the electric field and the presence of the thin layers. In these observations thin layers were only present when the electric field was pointed in the magnetic north-west or south-west quadrants, most frequently when the field was near magnetic west. Correlation between layer altitude and field direction was also observed, layers occurring at higher altitudes for fields directed in the north-west, and lower altitudes for fields directed to the south-west. The observations are compatible with the electric field mechanism for thin ionization layer formation. The scanning experiment showed that the layers were of a limited latitudinal extent, typically about 100 km up to a maximum of about 200 km.     

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.     

Simultaneous observations of E and F region electric field fluctuations at the magnetic equator

Simultaneous daytime observations of E region horizontal irregularity drift velocities in the equatorial electrojet and F region vertical plasma drifts were made on a few magnetically quiet days at the magnetic equatorial station of Trivandrum (dip 0.5°N). Measurements of the electrojet irregularity velocities by VHF backscatter radar and the F region vertical plasma drifts by HF Doppier radar are used to deduce the daytime East-West electric fields in the E and F regions, respectively. The fluctuating components of the electric fields are separated and subjected to power spectral analysis. The E and F region electric field fluctuations are found to be well correlated; the estimated correlation coefficient is in the range of 0.52–0.8. The fluctuation amplitudes are of the order of 15% over the background for the E region and 25% for the F region. The spectral analysis reveals dominant components in the range of 30–90 min with F region components stronger than those of the E region by a factor of about 1.5 on the average. The F region electric fields during daytime being coupled from the low latitude E region, the good correlation observed between the E and F region perturbations suggests that the electric fields in the E region at low and equatorial latitudes are coherent for the temporal scales of the order of few tens of minutes. The spectral characteristics are such that the commonly occurring medium scale gravity waves could possibly be the source for the observed fluctuations in the E and F region electric fields.     

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.     

EISCAT measurements of ion temperatures which indicate non-isotropic ion velocity distributions

Substantial increases of the ion temperature can be observed at high latitudes as a consequence of strong convection electric fields. We have measured, with EISCAT, three independent components of the ion velocity vector and temperature in the same scattering volume, at about 300 km. During periods of strong variations in ion velocity (consequently of the E-field), the ion temperatures derived at the 3 sites are different. This difference, which appears to be systematic for the two experiments studied, can be interpreted in terms of different ion temperature perpendicular and parallel to the magnetic field, i.e. T_i⊥ greater than T_i∥. Assuming that a bi-Maxwellian distribution is present for convection electric field strengths as large as 50 mV m⁻¹, one obtains an anisotropy factor of approximately 1.5. It also appears that resonant charge exchange is the dominant collision process. During the evening sector events studied, the electron density was decreasing, whereas the electron temperature was generally increasing. Such events are strongly related to variations in the magnetic H component detected on the ground.     

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.