A numerical study of polar ionospheric currents

Numerical calculations for the electric current in the polar ionosphere have been made by assuming some realistic distributions of the electric field and conductivity. Two dynamo actions are taken into account; one of which is induced by ionospheric winds and the other by the solar wind. For the solar wind dynamo action, it is found that the secondary polarization field caused by non-uniform distribution of ionospheric conductivity is much larger than the primary field induced by the solar wind, suggesting its important effect on charged particles in the magnetosphere, and that the irrotational current having a source and sink is of the same order of magnitude as the solenoidal current closing its circuit in the ionosphere. It is also found that the solar wind is, in general, more effective than the ionospheric winds in producing polar current systems such as DP 1 and 2, but in some cases the ionospheric winds have a significant effect on the current distribution.     

Equivalent ionospheric current systems representing IMF sector effects on the polar geomagnetic field

Equivalent ionospheric current systems representing IMF sector effects on the geomagnetic field in high latitudes are examined for each of the twelve calendar months by spherical harmonic analyses of geomagnetic hourly data at 13 northern polar stations for seven years. The main feature of obtained equivalent current systems includes circular currents at about 80° invariant latitude mostly in the daytime in summer and reversed circular currents at about 70° invariant latitude mainly at night in winter. Field-aligned current distributions responsible for equivalent currents, as well as vector distributions of electric fields and ionospheric currents, are approximated numerically from current functions of equivalent current systems by taking assumed distributions of the ionospheric conductivity. Two sets of upward and downward field-aligned current pairs in the auroral region, and also a field-aligned current region near the pole show seasonal variations. Also, ionospheric electric-field propagation along geomagnetic field lines from the summer hemisphere to the winter hemisphere with auroral Hall-conductivity effects may provide an explanation for the winter reversal of sector effects.     

Indication of anomalous conductivity at the nigerian dip equator

The geomagnetic daily variations at the Nigerian dip equator have been analyzed with the methodology introduced in a previous paper. It has been found that the height integrated current presents a notoriously higher amplification in Nigeria than in Peru. It has also been found that there exists a strong and inhomogeneous anomaly in the Earth's conductivity in Nigeria. And contrary to what is usually accepted, it is shown that its latitudinal distribution can not be precisely determined until the distribution and magnitude of the ionospheric currents at F-region heights is more accurately known.     

中低纬度变化磁场的分析(英文)

地球变化磁场呈现复杂时空特点,这是由引起该磁场的磁层一电离层电流以及地球内部感应电流的特性决定的。为了研究变化磁场的物理成因及其在日地物理事件中的特性。首先必须将组成变化磁场的各种成分分离开来,然后逐一加以研究。 我们采用自然正交分量法对我国八个地磁台站的时均值序列进行了分析,这些台站展布在27°12′48″到49°36′的中低纬度带内,正是Sq电流体系焦点所在的纬度带。分析结果表明,由发电机过程产生的Sq电流体系是这一纬度带主要的电流体系,与磁暴环电流有关的扰动电流体系也是十分重要的电流体系,在冬季月份,它往往超过Sq程度。此外与UT有关的磁扰变化也被明显地分离出来,它的成因可能与地球磁场的偏心结构有关。这些成份的相对大小随季节变化,而且有确定的纬度分布。 我们提出了一套单台分析和多台分析的方法。考虑到自然正交分量法收效快,稳定性好,所需资料列序列短的特点,这种方法可以推广到台站使用。自然正交分量法可以从成因上分离不同成因,使它在理论研究中具有优于一般付氏分析、时序迭加等方法,可为中低纬电流系成因研究提供有用的结果。     

Turbulent transport and heating in the auroral plasma of the topside ionosphere

Using plasma parameters from a typical stormtime ionospheric energy balance model, we have investigated the effects of plasma turbulence on the auroral magnetoplasma. The turbulence is assumed to be comprised of electrostatic ion cyclotron waves. These waves have been driven to a nonthermal level by a geomagnetic field-aligned, current-driven instability. The evolution of this instability is shown to proceed in two stages and indicates an anomalous increase in field-aligned electrical resistivity and cross-field ion thermal conductivity as well as a decrease in electron thermal conductivity along the geomagnetic field. In addition, this turbulence heats ions perpendicular to the geomagnetic field and hence leads to a significant ion temperature anisotropy.     

Numerical analysis of equatorial enhancement of geomagnetic sudden commencement

Equatorial behaviour of a polar-originating ionospheric current is examined by solving numerically the continuity equation on a two-dimensional spherical shell with appropriate assumptions for the ionospheric conductivity and the field-aligned source currents. The results show a clear daytime equatorial enhancement of the ionospheric currents in spite of much reduced electric field due to shielding effects of the enhanced Cowling conductivity there. The results are used for interpretation of the preliminary impulse of the geomagnetic sudden commencement.     

Description of the coastal effect at equatorial latitudes with applications to the Peruvian and Nigerian zones

The problem of the electromagnetic induction produced by a localized and an extended ionospheric current near an ocean coast, over a mantle of infinite conductivity, has been reduced to the solution of an integral equation where the induced current density appears in an implicit form. This formalism is applied to calculate the field induced by the geomagnetic daily variation due to the presence of the ocean at the Peruvian and Nigerian equatorial zones.     

The effect of realistic conductivities on the high-latitude neutral thermospheric circulation

The dynamics of the high latitude thermosphere are dominated by the ion circulation pattern driven by magnetospheric convection. The reaction of the neutral thermosphere is influenced by both the magnitude of the ion convection velocity and by the conductivity of the thermosphere. Using a threedimensional, time-dependent, thermospheric, neutral model together with different ionospheric models, the effect of changes in conductivity can be assessed. The ion density is described by two models: the first is the empirical model of Chiu (1975) appropriate for very quiet geomagnetic conditions, and the second is a modified version of the theoretical model of Quegan et al. (1982). The differences in the neutral circulation resulting from the use of these two ionospheric models emphasizes the need for realistic high latitude conductivities when attempting to model average or disturbed geomagnetic conditions, and a requirement that models should couple realistically the ionosphere and the neutral thermosphere. An attempt is made to qualitatively interpret some of the features of the neutral circulation produced at high latitudes by magnetospheric processes.     

The electrostatic field in the ionosphere

Theoretical calculations of the electrostatic field in the ionosphere are presented for different seasons and longitude zones. The corresponding current systems have been shown earlier to give good agreement with the geomagnetic Sq variations. The question of whether the electrostatic field is generated by winds in the ‘dynamo region’ or by other processes higher in the magnetosphere is evaluated in the light of recent observations. Several details of the electrostatic field variation, such as an increase near sunset, are noted.     

Extension of a polar ionospheric current to the nightside equator

A detailed analysis of rapid-run magnetograms from Guam (geomagnetic latitude = 4.2°) revealed that there are two kinds of geomagnetic sudden commencement (SC) observed in nighttime. One is the ordinary SC consisting of a main impulse only which has a smooth rise of the H-component. The other is a superposition by a small positive impulse on the very beginning part of the smooth rise of the main impulse and consequently the SC starts with a small stepwise increase of the H-component. The latter type of SC occurs between 20 and 08 h L.T. and its occurrence rate takes the maximum value of about 50% around 03 h L.T. Corresponding magnetograms from a dayside equatorial station (Huancayo, geomagnetic latitude = ?0.7°) were examined and a good correlation was found between the stepwise SC at the nightside (Guam) and SC¹ with a preliminary reverse impulse (PRI) at the dayside (Huancayo). Since PRI observed at the dayside equator may be interpreted as an extension of an ionospheric current due to an dusk-to-dawn electric field impressed on the polar ionosphere, our results show that a polar originating ionospheric current can extend to the nightside equator and produce a small but observable magnetic effect in spite of much reduced nighttime ionospheric conductivity.     

Optimal dynamos in the cores of terrestrial exoplanets: Magnetic field generation and detectability

A potentially promising way to gain knowledge about the internal dynamics of extrasolar planets is by remote measurement of an intrinsic magnetic field. Strong planetary magnetic fields, maintained by internal dynamo action in an electrically conducting fluid layer, are helpful for shielding the upper atmosphere from stellar wind induced mass loss and retaining water over long (Gyr) time scales. Here we present a whole planet dynamo model that consists of three main components: an internal structure model with composition and layers similar to the Earth, an optimal mantle convection model that is designed to maximize the heat flow available to drive convective dynamo action in the core, and a scaling law to estimate the magnetic field intensity at the surface of a terrestrial exoplanet. We find that the magnetic field intensity at the core surface can be up to twice the present-day geomagnetic field intensity, while the magnetic moment varies by a factor of 20 over the models considered. Assuming electron cyclotron emission is produced from the interaction between the stellar wind and the exoplanet magnetic field we estimate the cyclotron frequencies around the ionospheric cutoff at 10 MHz with emission fluxes in the range 10⁻⁴-10⁻⁷ Jy, below the current detection threshold of radio telescopes. However, we propose that anomalous boosts and modulations to the magnetic field intensity and cyclotron emission may allow for their detection in the future.     

Jupiter's outer atmosphere

The current state of the theory of Jupiter's outer atmosphere is briefly reviewed. The similarities and dissimilarities between the terrestrial and Jovian upper atmospheres are discussed, including the interaction of the solar wind with the planetary magnetic fields. Estimates of Jovian parameters are given, including magnetosphere and auroral zone sizes, ionospheric conductivity, energy inputs, and solar wind parameters at Jupiter. The influence of the large centrifugal force on the cold plasma distribution is considered. The Jovian Van Alien belt is attributed to solar wind particles diffused in towards the planet by dynamo electric fields from ionospheric neutral winds and consequences of this theory are given.     

Accurate approximate formulae for toroidal standing hydromagnetic oscillations in a dipolar geomagnetic field

Calculations of the toroidal eigenmodes of oscillation of the magnetospheric plasma have been important in explaining the nature of Pc3, Pc4 and Pc5 geomagnetic pulsations. In this paper perturbation solutions of the governing equations are presented. These are much more accurate than the WKB approximation which has often been used, and much simpler to compute than the numerical solutions which have been used. A method of including the finite ionospheric conductivity is also presented.     

Precession of the lunar core

Goldreich (Goldreich, P. [1967]. J. Geophys. Res. 72, 3135) showed that a lunar core of low viscosity would not precess with the mantle. We show that this is also the case for much of lunar history. But when the Moon was close to the Earth, the Moon’s core was forced to follow closely the precessing mantle, in that the rotation axis of the core remained nearly aligned with the symmetry axis of the mantle. The transition from locked to unlocked core precession occurred between 26.0 and 29.0 Earth radii, thus it is likely that the lunar core did not follow the mantle during the Cassini transition. Dwyer and Stevenson (Dwyer, C.A., Stevenson, D.J. [2005]. An Early Nutation-Driven Lunar Dynamo. AGU Fall Meeting Abstracts GP42A-06) suggested that the lunar dynamo needs mechanical stirring to power it. The stirring is caused by the lack of locked precession of the lunar core. So, we do not expect a lunar dynamo powered by mechanical stirring when the Moon was closer to the Earth than 26.0-29.0 Earth radii. A lunar dynamo powered by mechanical stirring might have been strongest near the Cassini transition.     

Possible relations between the rotational axis of the Earth's inner core and the magnetic dipole axis

The geomagnetic field is maintained by amagnetohydrodynamic dynamo process within the liquid outer core. The distribution of the associated electric currents is modified if the outer core is bounded by electrically conducting material. Then, eddy currents and the related magnetic fields are generated within these regions. In particular, the relative rigid rotation of the inner core produces a secondary magnetic field, which is superimposed on the dynamo field. The angle between the dipole axis of the total field and the rotational axis of the inner core is an important quantity needed for the theory of polar motion of the Earth. This angle is investigated for a broad spectrum of angular velocities of the inner core. To simplify the mathematical procedure, we model the dynamo field using an axisymmetric field generated by a system of electric currents within the outer core. The conductivity of the mantle is neglected. We find that the position of the dipole axis depends on the angular velocity of the inner core as well as on the distribution of the current system within the outer core. Coincidence of both axes can be reached if the angular velocity is high enough and if the current system is concentrated within a thin sheet near the outer core-inner core boundary.     

Modelling of Pc5 pulsation structure in the magnetosphere

Magnetohydrodynamic resonance theory is used to model the structure of the magnetospheric and ionospheric electric and magnetic fields associated with Pc5 geomagnetic pulsations. In this paper the variation of the fields across the invariant latitude of the resonance are computed. The results are combined with calculations of the variation along a field line to map the fields down to the ionosphere. In one case the results are compared with measurements obtained by the STARE auroral radar and show good agreement. The relationship between the width of the resonance region and ionospheric height-integrated Pedersen conductivity is computed and it is shown how auroral radar measurements of Pc5 oscillations could be used to determine ionospheric height-integrated Pedersen conductivity. It is pointed out that from these calculations it would be possible to identify the field line on which a satellite was located by comparing a Pc5 pulsation observed by the satellite, and the same pulsation observed by STARE.     

The role of the Hermanus Magnetic Observatory in geomagnetic field research

Geomagnetic field research carried out at the Hermanus Magnetic Observatory over the past decade is reviewed. An important aspect of this research has been the study of geomagnetic field variations, with particular emphasis on ULF geomagnetic pulsations. Features of geomagnetic pulsations which are unique to low latitude locations have been investigated, such as the cavity mode nature of low latitude Pi 2 pulsations and the role played by ionosphericO ⁺ ions in the field line resonances responsible for Pc 3 pulsations. A theoretical model has been developed which is able to account for the observed relationships between geomagnetic pulsations and oscillations in the frequency of HF radio waves traversing ionospheric paths. Other facets of the research have been geomagnetic field modelling, aimed at improving the accuracy and resolution of regional geomagnetic field models, and the development of improved geomagnetic activity indices.     

Mathematical theory of the ionospheric dynamo

A mathematical theory of the three-dimensional ionospheric dynamo is described, which is exact within the context of the assumption that the magnetic field lines are perfectly conducting.Explicit formulae are obtained for the electric fields, and hence for the current systems, which are excited by the dynamo e.m.f. associated with an arbitrary given wind field.They can also be applied to determine the fields arising from an external source of e.m.f., originating for example in the magnetosphere.The convergence properties of the solution are investigated for a simple symmetric wind field; it is found that inclusion of spherical harmonics up to degree 80 is sufficient to yield the features of the electrojet region, now fully coupled to the global current system.     

Aurora on Uranus: A Faraday disc dynamo mechanism

We propose a mechanism whereby the solar wind flowing past the magnetosphere of Uranus causes a Faraday disc dynamo topology to be established and power to be extracted from the kinetic energy of rotation of Uranus. An immediate consequence of this dynamo is the generation of Birkeland currents that flow in and out of the sunlit polar cap with the accompanying production of polar aurora. We calculate the power extracted from planetary rotation as a function of planetary dipole magnetic moment and the ionospheric conductivity of Uranus. For plausible values of ionospheric conductivity, the observed auroral power requires a magnetic moment corresponding to a surface equatorial field of the order of 4 Gauss, slightly larger than the value 1.8 Gauss given by the empirical “magnetic Bodes law”.     

Dissipative structures in the F-region of the equatorial ionosphere generated by Rayleigh-Taylor instability

The non-linear regime of electrostatic perturbations of the equatorial ionospheric F-region generated by Rayleigh-Taylor instability has been discussed, taking into account conductivity along magnetic field lines. A closed non-linear equation has been derived in the stationary limit for the polarization electric field potential. It coincides with the Karman equation of an ideal liquid. To solve the equation, the averaged variational Whitham method has been proposed. Some solutions localized along and across the geomagnetic field, B, as well as quasi-periodic solutions in the transverse direction, have been investigated. Non-linear longitudinal localization of perturbations has been shown to be due to electron-ion collisions.