Simulations of the behaviour of the polar thermosphere and ionosphere for northward IMF

The dynamics and structure of the polar thermosphere and ionosphere within the polar regions are strongly influenced by the magnetospheric electric field. The convection of ionospheric plasma imposed by this electric field generates a large-scale thermospheric circulation which tends to follow the pattern of the ionospheric circulation itself. The magnetospheric electric field pattern is strongly influenced by the magnitude and direction of the interplanetary magnetic field (IMF), and by the dynamic pressure of the solar wind. Previous numerical simulations of the thermospheric response to magnetospheric activity have used available models of auroral precipitation and magnetospheric electric fields appropriate for a southward-directed IMF. In this study, the UCL/Sheffield coupled thermosphere/ionosphere model has been used, including convection electric field models for a northward IMF configuration. During periods of persistent strong northward IMF B_z, regions of sunward thermospheric winds (up to 200 m s⁻¹) may occur deep within the polar cap, reversing the generally anti-sunward polar cap winds driven by low-latitude solar EUV heating and enhanced by geomagnetic forcing under all conditions of southward IMF B_z. The development of sunward polar cap winds requires persistent northward IMF and enhanced solar wind dynamic pressure for at least 2–4 h, and the magnitude of the northward IMF component should exceed approximately 5 nT. Sunward winds will occur preferentially on the dawn (dusk) side of the polar cap for IMF B_y negative (positive) in the northern hemisphere (reverse in the southern hemisphere). The magnitude of sunward polar cap winds will be significantly modulated by UT and season, reflecting E-and F-region plasma densities. For example, in northern mid-winter, sunward polar cap winds will tend to be a factor of two stronger around 1800 UT, when the geomagnetic polar cusp is sunlit, then at 0600 UT, when the entire polar cap is in darkness.     

Equatorial penetration of magnetic disturbance effects in the thermosphere and ionosphere

The neutral dynamic and electrodynamic coupling between high and low latitudes, and the mutual interactions between these two processes, are investigated. For 22 March 1979, when a sudden increase in magnetic activity occurred, we have analyzed the following experimental data: (a) neutral densities and cross-track neutral winds as a function of latitude (0°–80°) near 200 km from a satellite-borne accelerometer; (b) hourly mean H-component magnetic data from the Huancayo Observatory (0.72°S, 4.78°E; dipole geomagnetic coordinates) magnetometer; and (c) hourly mean foF2 measurements from the ionosonde at Huancayo. Comparisons are also made with a self-consistent thermosphere-ionosphere general circulation model and with observationally-based empirical models of winds and density.In concert with the increase in magnetic activity to K_p levels of 5–7, a nighttime (2230 LT) westward intensification of the neutral wind approaching 400 ± 100 ms⁻¹ occurred near the magnetic equator on 22 March 1979, accompanied by a 35% increase in neutral mass density. About 2 h after each of two substorm commencements associated with periods of southward IMF, ∼100γ and ∼200γ reductions in the daytime Huancayo H-component (corrected for ring current effects) are interpreted in terms of ∼0.5 and ∼1.0 mVm⁻¹ westward perturbation electric fields, respectively. An intervening 2-hour period of northward IMF preceded a positive equatorial magnetic perturbation of about 200γ. Time scales for field variations are a few hours, suggesting that processes other than Alfven shielding are involved. Variations in f₀F2 (∼ ± 1.0 MHz) over Huancayo are consistent with the inferred electric fields and magnetic variations. Similar equatorial perturbations are found through examination of other magnetic disturbances during 1979.     

A survey of simultaneous observations of the high-latitude ionosphere and interplanetary magnetic field with EISCAT and AMPTE-UKS

This paper surveys the results of simultaneous observations by the EISCAT incoherent scatter radar and the AMPTE-UKS satellite, made during three periods in September and October 1984, when AMPTE-UKS was in the solar wind on the dayside of the Earth and the UK-POLAR EISCAT experiment was measuring ionospheric parameters at invariant latitudes 70.8–75.0°. A total of 42 h of EISCAT convection velocity data, with 2.5 min resolution, were obtained, together with 28 h of simultaneous 5 s resolution AMPTE-UKS observations of the solar wind and interplanetary magnetic field (IMF). The general features of the AMPTE-UKS data are described in Section 2 and those of the EISCAT data are described in Sections 3 and 4. The main subjects discussed are the form of the plasma convection patterns and their dependence on all three components of the IMF (Section 5), the ionospheric response to abrupt changes in the IMF (Section 6), in particular a sharp ‘southward turning’ of the IMF on 27 October 1984, and a crossing of an IMF sector boundary. Section 7 describes ‘short lived rapid flow burst’, which are believed to be related to flux transfer events at the magnetopause.     

Control of diurnal and seasonal variation of particle precipitation by regular changes of the interplanetary magnetic field

Systematic changes of the position of the dipole axis of the Earth's magnetic field with respect to the solar axis induce distinct daily and seasonal variations of the vertical B_z-component in the solarmagnetospheric coordinate system (B_ZSM). Depending on the direction of the interplanetary magnetic field (IMF), negative B_ZSM- values are produced in spring by T polarity and in autumn by A polarity, whereas in the diurnal variation lowest B_ZSM-values have been calculated to occur near 23 UT for T, and near 11 UT for A polarity, respectively. In different ionospheric and geomagnetic parameters measured at high and midlatitudes increased precipitation of high energetic particles into the lower thermosphere and upper mesosphere has been detected during periods with negative B_ZSM-components. The seasonal variation of the parameters investigated, with maximum values near the equinoxes, as well as a part of their diurnal variations, can thus be explained by particle precipitation being markedly controlled by the IMF sector structure.     

Interplanetary magnetic field structure and the variations of the ELF and VLF emissions in the topside ionosphere

The Intercosmos-13 data obtained when measuring ELF and VLF emission amplitudes during the vernal equinox of 1975 at auroral latitudes and over the polar caps are compared with certain IMF parameters [the polarity of the sector structure, the signs and magnitudes of the B_y, and B_z, components of the IMF as defined in the Solar Ecliptic coordinate system (Nishida, 1978)]. The comparison shows that:

• 1.(i) the positive polarity of the IMF sector structure (when IMF vector is directed toward the Earth) involves an enhanced probability of the detection of larger emission field intensities (>25–30dB);
• 2.(ii) the emission median intensity is ~20dB higher at B_y > 0 compared with B_y < 0;
• 3.(iii) the 0.72 kHz emission median intensity in the polar caps and at night-side auroral latitudes is lower when B_z > 0 as compared with B_z < 0;
• 4.(iv) at vernal equinox there is no north-south asymmetry in the dependence of ELF and VLF emission intensity on the IMF parameters.

     

Theoretical models of polar-cap convection under the influence of a northward interplanetary magnetic field

The unexpected patterns of high-latitude auroral luminosity and ionospheric convection that are observed when the interplanetary magnetic field (IMF) has a northward orientation have inspired a variety of theoretical interpretations. The existing models, all referring to steady-state conditions, can be classified according to the topology of the polar magnetic field lines and of the polar-cap convection streamlines. The classes of model include: (1) a closed magnetosphere model, (2) a conventional open model with a distorted, but topologically unchanged, polar-cap boundary, (3) a conventional open model with distorted, but topologically unchanged, polar-cap convection cells, (4) a modified open model with ‘lobe convection cells’ contained wholly on open magnetic-field lines, and (5) a modified open model with a bifurcated polar cap. The third and fourth types require significant regions of sunward flow on open polar-cap field lines, a concept that presents serious theoretical difficulties. The other three types appear equally viable from a theoretical point of view, and the comparison against observations is an ongoing enterprise. Outstanding theoretical questions include (a) how do observed structures in the polar ionosphere map along magnetic field lines into the magnetosphere?, (b) what is the mechanism that drives the observed sunward convection at highest latitudes on the day side?, and (c) what role does time dependence play in the observed phenomena?     

Dynamical coupling of the auroral F-region ionosphere and thermosphere: case studies

During two 24 h periods of EISCAT observations in the summer of 1982, the F-region ion temperature and density responded differently before and after midnight to large ion convective flows. Such observations were recently reported at Chatanika (Alaska), however, the mechanism invoked to interpret these measurements (large day-to-night variation in electron density affecting the coupling between ions and neutrals) appears insufficient, for summer conditions, to account for the EISCAT observations. Hence, it is proposed, with the support of Fabry-Perot observations and numerical models, that in addition to the electron density asymmetry, the presence of a large southward neutral wind around midnight induces, through Coriolis coupling, a zonal neutral wind of an opposite direction to the convective flow. This enhances considerably the frictional energy and momentum transfer between ions and neutrals in the post-midnight sector.     

Inversion of oblique ionograms including the Earth's magnetic field

A technique for determining ionospheric electron distribution from oblique ionograms is presented, based on the inversion method of Reilly and Kolesar (Radio Sci. 24, 575, 1989). It makes use of an equivalent operating frequency and an additional term to account for magnetoionic effects associated with the Earth's magnetic field. The technique is demonstrated by application to synthetic oblique ionograms, and to an experimentally obtained ionogram.     

On the global and regional behaviour of the mid-latitude ionosphere

Data from a chain of seven ionosondes in the range of 56–38 N and 1–38° E geographic coordinates were analysed to illustrate the global and regional behaviour of the mid-latitude F-region for some selected geomagnetic storms that occurred during the solar cycle 21. It was found that there are different spatial scales in the response of the mid-latitude ionosphere to the disturbance in the magnetosphere-ionosphere thermosphere system. The physical mechanisms and processes are discussed in relation to the relevance of various theories in the understanding of the dynamics of ionospheric storms.     

Mechanisms and time-scales of the magnetospheric response to the interplanetary magnetic field changes during the 8 May 1986 substorm

On 8 May 1986, between 1113 and 1600 UT, an isolated magnetospheric substorm was observed, during which the AE-index exceeded 700 nT (CDAW 9E event). Three available sets of measurements (a) of the solar-wind parameters (IMP-8 satellite), (b) of the magnetotail energy flux (ISEE-1 spacecraft), and (c) of ground magnetic observatories, allowed us to make a detailed study of the overall magnetospheric response to changes of the interplanetary magnetic field (IMF) direction, during this event of weak solar-wind coupling.In order to study the mechanisms and time-delays of the magnetospheric response to the abrupt increase of the solar-wind energy input, we have evaluated the total magnetospheric energy output U_T following two different methods: (a) Akasofu's method, taking the ring current decay time τ_R constant, and (b) Vasyliunas' method where the values of u_t are independent of the solar-wind energy input as determined from the epsilon parameter. Both methods suggest that the driven system has been considerably developed during this substorm, while an unloading event has been superposed at the expansion onset.     

Seasonal and latitudinal variations of eddy diffusion coefficient in the mesosphere and lower thermosphere

Two approaches to turbulence parameter determination around the turbopause are discussed. It is shown that there is a contradiction in the results concerning seasonal variations of the eddy diffusion coefficient between direct measurements and estimations based on minor constituents. Taking account of the vertical mean transport in the continuity equations for the above constituents might eliminate that contradiction. Results of the re-examination of the authors' data on the Ar to N₂ ratio published earlier are presented. The results show stronger turbulence in winter than in summer, the amplitude of the effect increasing towards higher latitudes.     

Simulation of the global structure of stationary planetary waves in the mesosphere and lower thermosphere

When simulating the global structure of stationary planetary waves (SPW) the problem of obtaining the numerical solution in the equatorial region appears. It results from the presence of apparent singularities in the operator of the SPW latitudinal structure when the Coriolis parameter is small. The new method based on SPW latitudinal operator inversion is proposed. This method permits the difficulties arising from the simulation of stationary large scale disturbances at low latitudes to be avoided. The global structure of SPW with zonal wave number m = 1 at the mesosphere and lower thermosphere heights has been calculated for the background zonal wind distribution representing a climatic picture of the solstice conditions. In the region of the mean zonal westerlies the SPW penetration across the equator is obtained. The SPW at low latitudes are shown to appear most significantly in the zonal component of the wind velocity. The influence of planetary wave motions on the distribution of longlived species in the ionospheric D-region and at the heights of lower thermosphere are discussed.     

Simultaneous observation of the quasi-two-day variations in the lower and upper ionosphere over Europe

The large scale character of the observed quasi-two-day fluctuations in the whole ionosphere (from D- uptoF-region maximum) over Europe is shown. The study is based on the lower and upper ionospheric data obtained in Sofia (42.9°, 23.4°E), Ebre Observatory (40.9°N, 0.5°E) and El Arenosillo (37.1°N, 6.7°W) during two summer intervals: June–August 1980 and 1983. The obtained prevailing periods for the F-region fluctuations are 52–55 h and the mean amplitude is higher than 1 MHz. It was found that the fluctuations propagate westward with a mean phase velocity between 4.6 and 6° /h. The quasi-two-day variations in the F-region maximum are probably generated by flucutations in the mesospheric, neutral wind. During the time when well developed quasi-two-day fluctuations exist in the mesospheric neutral wind, similar variations are observed in the lower ionosphere also. Possible mechanisms for generating the D- andF-region electron density fluctuations from these oscillations in the neutral wind are proposed.     

Seasonal and tidal variations of neutral temperatures and densities in the high latitude lower thermosphere measured by EISCAT

Measurements of ion temperature, ion-neutral collision frequency and ion drift in the E-region from the period December 1984 to November 1985 are used to derive neutral temperatures, densities and meridional winds in the altitude intervals 92–120 km, 92–105 km and 92–120 km, respectively. Altitude profiles of temperature and density and their seasonal variations are compared with the CIRA 1972 and MSIS 1983 models and the effects of geomagnetic activity are demonstrated. Semi-diurnal tidal variations in all three parameters are derived and the comparison with lower latitude measurements is discussed.     

Interactions between whistler-mode waves and energetic electrons in the coupled system formed by the magnetosphere,ionosphere and atmosphere

The physical mechanism of a cyclotron resonance interaction between trapped energetic electrons and whistler-mode waves in the magnetosphere is discussed. Not only do the electrons have their pitch angles reduced in this interaction, so that they may be precipitated into the upper atmosphere, but also the waves can be amplified. Such a flux of precipitating electrons can, either by direct ionisation or via bremsstrahlung radiation, cause a pimple to be produced on the bottom of the ionosphere. That can significantly modify the amplitude and/or phase of very low frequency radio signals propagating in the Earth-ionosphere waveguide. Various experimental observations that demonstrate the reality of such effects are reviewed. The conditions necessary for a positive feedback situation are discussed, and some evidence for its existence assessed.     

Small-scale fluctuations of magnetic and electric components of the ELF and VLF wave fields in the sub-auroral topside ionosphere : stochastic characteristics of the wave field

When the Interkosmos-14 and Interkosmos-19 satellites crossed the region of spatially varying electron concentration in the topside ionosphere adjacent to the high-latitude boundary of the main ionospheric trough, it was discovered that there were simultaneous fluctuations of plasma density, temperature and the amplitudes (H_x and E_y) of the ELF and VLF radio/plasma emissions. The probability characteristics of the naturally perpendicular H_x and E_y fluctuations are analysed. The correlation coefficient R(_H, E_y) turned out to be less than 0.6 at frequencies of F ⩽ 4.65 kHz, while at higher frequencies R increases, up to 0.9 at 15 kHz. The following interpretations are proposed:

1. 1.1. While measuring noise emissions, as a rule a mixture of numerous elementary waves is recorded.
2. 2.2. At frequencies exceeding the local lower hybrid resonance frequency (in our case f_LHR ≈ 5 kHz), a mixture of electromagnetic waves experiencing the influence of the inhomogeneous electron concentration N_e is registered.
3. 3.3. At frequencies which are lower than the local value f_LHR the mixture mainly consists of ELF waves. The wave field has a complicated structure, and the dynamical coherence between electric and magnetic field components is not as simple as at VLF frequencies (f ≈ 15 kHz).
4. 4.4. It is shown that the wave components for a mixture of electromagnetic and electrostatic waves (for instance a mixture of VLF and lower hybrid frequency waves) have a lower correlation coefficient because the electrostatic waves are unrelated to the electromagnetic waves.
5. 5.5. The correlation analysis offers an opportunity to detect the presence of waves of various types in the wave mixture.

     

Application of a simplified theory of ELF propagation to a simplified worldwide model of the ionosphere

The approximate theory of ELF propagation in the Earth-ionosphere transmission line described by Booker (1980) is applied to a simplified worldwide model of the D- and E-regions, and of the Earth's magnetic field. At 1000 Hz by day, reflection is primarily from the gradient on the underside of the D-region. At 300 Hz by day, reflection is primarily from the D-region at low latitudes, but it is from the E-region at high latitudes. Below 100 Hz by day, reflection is primarily from the gradient on the underside of the E-region at all latitudes. By night, reflection from the gradient on the topside of the E-region is important. There is then a resonant frequency (~300 Hz) at which the optical thickness of the E-region for the whistler mode is half a wavelength. At the Schumann resonant frequency in the Earth-ionosphere cavity (~8Hz) the nocturnal E-region is almost completely transparent for the whistler mode and is semi-transparent for the Alfvén mode. Reflection then takes place from the F-region. ELF propagation in the Earth-ionosphere transmission line by night is quite dependent on the magnitude of the drop in ionization density between the E- and F-regions. Nocturnal propagation at ELF therefore depends significantly on an ionospheric feature whose magnitude and variability are not well understood. A comparison is made with results based on the computer program of the United States Naval Ocean Systems Center.     

Solar cycle and seasonal variations of the low latitude OI 630 nm nightglow

Regular zenith measurements of the OI 630 nm nightglow emission have been carried out at Cachoeira Paulista (22.7°S, 45.0°W; geomag. 11.9°S), Brazil, since 1975. The long series of observations during the period 1975–1982, including the ascending phase of the last solar cycle, permitted studies of solar cycle effects and seasonal variations. A large intensity increase, about seven times, from low solar activity to high solar activity has been observed. Also, the seasonal-nocturnal intensity variations show large changes between years of low and high solar activity. The characteristics of the variations observed are closely related to the equatorial electric field variations, since the observation site is under the southern equatorial ionospheric anomaly crest.     

Comparisons between EISCAT observations and model calculations of the high latitude ionosphere

Calculations using a numerical model of the convection dominated high latitude ionosphere are compared with observations made by EISCAT as part of the UK-POLAR Special Programme. The data used were for 24–25 October 1984, which was characterized by an unusually steady IMF, with Bz < 0 and By > 0; in the calculations it was assumed that a steady IMF implies steady convection conditions. Using the electric field models of Heppner and Maynard (1983) appropriate to By > 0 and precipitation data taken from Spiroet al. (1982), we calculated the velocities and electron densities appropriate to the EISCAT observations. Many of the general features of the velocity data were reproduced by the model. In particular, the phasing of the change from eastward to westward flow in the vicinity of the Harang discontinuity, flows near the dayside throat and a region of slow flow at higher latitudes near dusk were well reproduced. In the afternoon sector modelled velocity values were significantly less than those observed. Electron density calculations showed good agreement with EISCAT observations near the F-peak, but compared poorly with observations near 211 km. In both cases, the greatest disagreement occurred in the early part of the observations, where the convection pattern was poorly known and showed some evidence of long term temporal change. Possible causes for the disagreement between observations and calculations are discussed and shown to raise interesting and, as yet, unresolved questions concerning the interpretation of the data. For the data set used, the late afternoon dip in electron density observed near the F-peak and interpreted as the signature of the mid-latitude trough is well reproduced by the calculations. Calculations indicate that it does not arise from long residence times of plasma on the nightside, but is the signature of a gap between two major ionization sources, viz. photoionization and particle precipitation.