A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Electric circuits and their generators in the earth’s magnetosphere: The concept of electric circuits as applied to the first phase of the April 6, 2000, superstorm

{1} The first phase of the superstorm on April 6, 2000 was studied based on the analogy between systems of magnetospheric currents and wire electric currents. The conventional dataset supplemented with maps of ionospheric equivalent currents (ECs) and field-aligned currents (FACs) was also used. The application of this analogy made it possible to introduce spatial R.N inhomogeneities into FAC distributions in the two-dimensional ionosphere and three types of meridional current systems, MCS-0, MCS-1, and MCS-2, providing electric coupling of three Iijima and Potemra FAC Regions. This basis was used to describe the formation and observed dynamics of ionospheric auroral electrojets and three-dimensional current systems in a disturbed magnetosphere-ionosphere system. The results the modify known paradigms of the substorm current wedge (SCW). A new important fact was noted: simultaneously with the beginning of the disturbance expansion phase due to the stepwise growth in the dynamic pressure of the solar wind (SW), the stepwise growth in the area of polar cap and in the electromagnetic energy flux coming to the magnetosphere from the SW were observed.     

A new time-dependent ionosphere–magnetosphere coupling model: Comparison of field-aligned currents against ST5 observations

By using Tsyganenko's model for the magnetosphere's magnetic field, which links two hemispheres of the ionosphere, and adopting a practical boundary condition for the electric potential around the polar cap, we developed a new ionosphere–magnetosphere coupling model based on prairie view dynamo code (PVDC). The new model takes the variations in solar wind and interplanetary magnetic field, as well as the geomagnetic activity, into account. Rather than the previous version of PVDC that is useful only for quiet conditions, the new model enables to calculate the electric potential and currents in the ionosphere and the field-aligned current (FAC) off the ionosphere in quiet and disturbed times. Comparison of the calculated FAC with the measurements of Space Technology 5 (ST5) mission shows a good agreement.     

Observations of ionospheric flows and particle precipitation following a Sudden Commencement

On May 4, 1998, at 0227 UT an interplanetary shock crossed the WIND spacecraft, and half an hour later a Sudden Commencement occurred. Coinciding with the Sudden Commencement a rapid intensification of the flux of particle precipitation into the ionosphere was observed. Evidence is presented that the ionospheric electric fields were influenced by the associated dynamic variations of the ionospheric conductivities. Following the initial phase the ionospheric flow speeds increased rapidly over the next 20 min to more than 2000 m/s, in agreement with an increased effective coupling of the solar wind energy to the magnetosphere following the interplanetary shock that caused the Sudden Commencement. These strong flows were meandering in latitude, a type of plasma flow modulation that has been reported before to occur during Omega band events: a string of alternating field-aligned currents propagating eastward. The riometer absorption was found to be at a minimum in regions associated with outward directed field aligned currents. The riometer absorption regions (the regions of particle precipitation) were drifting with E × B drift speed of the ionospheric electrons.     

The effects of IMF and convection on thermal ion outflow in magnetosphere-ionosphere coupling

We study the influence of the interplanetary magnetic field (IMF) and convection electric field on the rate and destination of polar wind and other thermal (low-energy) ion outflows, and its resulting effects on magnetosphere–ionosphere coupling, using single-particle trajectory simulations in conjunction with ion velocity distribution measurements on Akebono and IMF and ionospheric convection data. We find that the ions preferentially feed the dusk sector of the plasma sheet when the IMF is duskward (B_y>0), and are more evenly distributed in the plasma sheet when the IMF is dawnward. The flow of oxygen ions originating from the noon or dusk sectors of the polar cap has a higher probability of reaching the magnetosphere and beyond compared with that from the dawn or midnight sectors, due to the increased centrifugal acceleration associated with the larger magnetic field curvature near noon and the increased convection electric field in the dusk sector. The flow is enhanced and confined to lower L-shells at times of strongly southward IMF, compared with that at times of northward IMF. The outflow rate to both the plasma sheet and the magnetotail correlates strongly with the ion temperature. As a result, the IMF and the convection electric fields affect both the overall magnitude and the detailed distribution of mass transfer from the ionosphere to the magnetosphere in magnetosphere–ionosphere coupling.     

Spontaneous substorms and ordered type of magnetospheric disturbances during the superstorm of November 20, 2003

The interval 0000-1400 UT of the superstorm of November 20, 2003, has been studies based on the ACE/WIND data and the MIT2 magnetogram inversion technique. The distributions of the electric potential and currents, field-aligned currents, and Joule heat in the ionosphere have been calculated. The variable magnetotail length and powers coming into the magnetosphere, ionosphere and ring current have been estimated. The selected superstorm intervals, when it became possible to identify the disturbance mode produced by the interaction between the variable solar wind dynamic pressure and IMF effects, have been described. Spontaneous substorms, two types of driven responses to changes in IMF or in the solar wind dynamic pressure (P _d ), zero events at simultaneous jumps of IMF and P _d , and a previously unknown mode of saturation of the ionospheric electric field at a redistribution of the energy coming into the magnetosphere between the ionosphere and ring current are among the selected modes.     

Interhemispheric comparison of TDIM E- and F-region ionospheres on 14 January 1988

The USU time-dependent ionospheric model (TDIM) simulated the northern (winter) and southern (summer) ionospheres as they responded to the changing solar wind and geomagnetic activity on 14 January 1988. This period began with moderately disturbed conditions, but as the IMF turned northward, the geomagnetic activity decreased. By 1400 UT, the IMF B_y component became strongly negative with B_z near zero; and eventually B_z turned southward. This began a period of intense activity as a magnetic storm developed. The magnetospheric electric field and auroral electron precipitation drivers for these simulations were obtained from the Naval Research Laboratories (NRL) Magnetohydrodynamic (MHD) magnetospheric simulation for this event.The F-region ionospheric simulations contrast the summer–winter hemispheres. Then, the difference in how the two hemispheres respond to the geomagnetic storm is related to the differences in magnetospheric energy deposition in the two hemispheres. This also emphasizes the role played by the E-region in the magnetosphere–ionosphere (M–I) coupling and subsequent lack of conjugacy in the two hemispheres. The F-region’s response to the changing geomagnetic conditions also demonstrates a striking lack of conjugacy. This manifests itself in a well-defined ionospheric morphology in the summer hemisphere and a highly irregular morphology in the winter hemisphere. These differences are found to be associated with the differences in the magnetospheric electric field input.     

Multimoment convecting flux tube model of the polar wind system with return current and microprocesses

Multimoment fluid simulation frameworks, which effectively account for anomalous transport due to microprocesses, combine best features of small-scale kinetic and global-scale MHD models. The most practical models of this type, 1D flux tube models, have been successfully used for realistic simulations of space plasmas including polar wind and magnetosphere–ionosphere coupling processes characterized by a wide range of temporal and spatial scales. Our earlier flux tube models with field-aligned current and microprocesses have been formulated for spatially stationary flux tubes. However, horizontal convection due to electric fields is an important aspect of the high-latitude ionosphere–polar wind system and typical time scales of the polar wind upflow are comparable to the transit time across the polar cap. To take into account this important feature we have added flux tube convection to our earlier model. Using typical convecting flux tube that starts outside auroral oval, then enters and leaves downward current region, it has been shown that anomalous transport effects due to current-driven microinstabilities significantly alter dynamics of several plasma moments and should be taken into account for an accurate interpretation and prediction of the observed data. Future applications of our new model have also been discussed.     

Influence of ionospheric flows on the effect of the planetary electric generator

Exact distributions of the electric field and current density in the vicinity of the rotating magnetized planet are determined within the planetary electric generator model. The model planetary plasma envelope is assumed to be inhomogeneous and consists of an atmosphere, an ionospheric layer rotating with respect to the atmosphere, and a magnetosphere located beyond the ionosphere. The model parameters, under which the influence of a large-scale flow of the ionospheric plasma on the electric field and current in the lower atmosphere is significant are determined. It is shown that the ionospheric superrotation reduces the electric field arising in the Earth’s lower atmosphere due to the planetary generator effect.     

Effects on the Ionosphere Due to Phenomena Occurring Below it

The terrestrial thermosphere and ionosphere form the most variable part of theEarth's atmosphere. Because our society depends on technological systems thatcan be affected by thermospheric and ionospheric phenomena, understanding,monitoring and ultimately forecasting the changes of the thermosphere–ionosphere system are of crucial importance to communications, navigation and the exploration of near-Earth space. The reason for the extreme variability of the thermosphere–ionosphere system isits rapid response to external forcing from various sources, i.e., thesolar ionizing flux, energetic charged particles and electric fields imposed via the interaction between the solar wind, magnetosphere and ionosphere, as well as coupling from below (meteorological influences) by the upward propagating, broad spectrum,internal atmospheric waves (planetary waves, tides, gravity waves) generated in thestratosphere and troposphere. Thunderstorms, typhoons, hurricanes, tornadoes andeven seismological events may also have observable consequences in the ionosphere.The release of trace gases due to human activity have the potential to cause changes inthe lower and the upper atmosphere.A brief overview is presented concerning the discoveries and experimentalresults that have confirmed that the ionosphere is subject to meteorologicalcontrol (especially for geomagnetic quiet conditions and for middle latitudes).D-region aeronomy, the winter anomaly of radiowave absorption, wave-liketravelling ionospheric disturbances, the non-zonality and regional peculiaritiesof lower thermospheric winds, sporadic-E occurrence and structure, spread-Fevents, the variability of ionospheric electron density profiles and Total ElectronContent, the variability of foF2, etc., should all be considered in connection withtropospheric and stratospheric processes. Ionospheric weather, as a part of spaceweather, (i.e., hour-to-hour and day-to-day variability of the ionospheric parameters)awaits explanation and prediction within the framework of the climatological, seasonal,and solar-cycle variations.     

A statistical study of magnetosphere–ionosphere coupling in the Lyon–Fedder–Mobarry global MHD model

The statistics of magnetosphere–ionosphere (MI) coupling derived from a two-month long run of the Lyon–Fedder–Mobarry (LFM) global simulation model are investigated. MI coupling characteristics such as polar cap potential and field-aligned current (FAC), downward Poynting flux and vorticity of ionospheric convection are compared with observed statistical averages and with results from the Weimer 05 empirical model. The comparisons for eight different IMF clock-angle orientations show that the LFM model produces reasonably accurate average distributions of the Region I and Region II currents. Both current systems have average amplitudes similar to those observed by the Iridium satellite constellation; however, the average LFM amplitudes are smaller by a factor of two compared with the values from the Weimer 05 model. The comparisons of polar cap potential show that the LFM model produces reasonable patterns of ionospheric convection, but the average cross polar cap potential (CPCP) is greater than the observed results by a factor of approximately 2 and greater than Weimer 05 by a factor of 1.5. The differences in convection in LFM results relative to the Weimer 05 model accounts for much of the difference in the Poynting flux patterns and integrated power produced by the two models. The comparisons of average ionospheric field-aligned vorticity show good agreement on the dayside; however, the LFM model gives higher nightside vorticity which may imply that the ionospheric conductance on the nightside is too small in the simulation.     

Earthquake-related Electric Field Changes Observed in the Ionosphere and Ground

The changes of the ionospheric electric field before and after four huge earthquakes, which include the Ms 8.7 earthquake of 2004 and the Ms 8.5 earthquake of 2005 in Sumatra of Indonesia, the Ms 8.0 Wenchuan earthquake of 2008 in China, the Ms 8.8 earthquake of 2010 in Chile, and their strong aftershocks are studied in this paper. The significant results revealed that the power spectral density of low-frequency electric field below 20 Hz in the ionosphere, a kind of electromagnetic radiation phenomena, increased abnormally before and after the earthquakes and partially corresponded to the increased power spectral density of the low-frequency geoelectric field in time. This research preliminarily indicates that the low-frequency electromagnetic radiation during the imminent stages before such earthquakes could be detected by the observation of the ionospheric electric field. However, the spatial, temporal, and intensive complexities of the electric field anomalies in the ionosphere before earthquakes have come in sight also.     

Penetration to the Earth’s surface of standing Alfvén waves excited by external currents in the ionosphere

The problem of boundary conditions for monochromatic Alfvén waves, excited in the magnetosphere by external currents in the ionospheric E-layer, is solved analytically. Waves with large azimuthal wave numbers m≫1 are considered. In our calculations, we used a model for the horizontally homogeneous ionosphere with an arbitrary inclination of geomagnetic field lines and a realistic height disribution of Alfvén velocity and conductivity tensor components. A relationship between such Alfvén waves on the upper ionospheric boundary with electromagnetic oscillations on the ground was detected, and the spatial structure of these oscillations determined.     

The magnetosphere–ionosphere system from the perspective of plasma circulation: A tutorial

This tutorial review examines the role of O⁺ in the dynamics of magnetosphere–ionosphere coupling. The life cycle of an O⁺ plasma element is considered as it circulates from the mid- to high-latitude ionosphere. Energization and diversion of the convecting plasma element into outflows involves Alfvénic turbulence at the low-altitude base of the cusp and plasmasheet boundary layer and in downward-current “pressure cookers.” Observational evidence indicating that O⁺ dominates the plasmasheet and ring current during extreme storm intervals is reviewed. The impacts of an O⁺-enriched plasma on solar wind–magnetosphere–ionosphere coupling are considered at both the micro and global scales. A synthesis of results from observation, theory and simulations suggests that the presence of O⁺ in the magnetosphere is both a disruptive and a moderating agent in maintaining the balance between dayside and nightside magnetic merging.     

Distribution of the field−aligned currents in the ionosphere: dawn–dusk asymmetry and its relation to the asymmetry between the two hemispheres

The dynamics of the system of field-aligned currents (FACs) and closing ionospheric Pedersen currents is considered with the use of original processing methods and the data from four substorms of different types. The total current system comprises of two parts. One is the well-known substorm current wedge (SCW) system, in which the zonal (westward ) current closes FACs in the R1 zone (region). The component 2 consists of two pairs of meridional currents flowing equatorward and poleward in the R1 region and creating regions R0 and R2 (according to the classification of Iijima and Potemra). It is shown that the total FAC of the disturbed magnetosphere–ionosphere system is dominated by the contribution of component 2, which contradicts the original version of the SCW model but is consistent with new data. The quantitative characteristics of the dawn–dusk asymmetry are determined for the FAC distribution in the ionosphere for each substorm. It is shown that the ratio of the average intensities of FACs in the regions R0 and R2 was I_R0/I_R2 ≥ 0.4, which contradicts the popular opinion that there are no FACs in the polar cap. In addition, a relatively rare event of the simultaneous start of the substorm explosive phase and the SSC caused by the dynamic impact of the solar wind when the polar cap expands (rather than compresses as usual) is considered.     

Earth’s Electromagnetic Environment

The natural spectrum of electromagnetic variations surrounding Earth extends across an enormous frequency range and is controlled by diverse physical processes. Electromagnetic (EM) induction studies make use of external field variations with frequencies ranging from the solar cycle which has been used for geomagnetic depth sounding through the 10\(^{-4}\)–10\(^4\) Hz frequency band widely used for magnetotelluric and audio-magnetotelluric studies. Above 10\(^4\) Hz, the EM spectrum is dominated by man-made signals. This review emphasizes electromagnetic sources at \(\sim\)1 Hz and higher, describing major differences in physical origin and structure of short- and long-period signals. The essential role of Earth’s internal magnetic field in defining the magnetosphere through its interactions with the solar wind and interplanetary magnetic field is briefly outlined. At its lower boundary, the magnetosphere is engaged in two-way interactions with the underlying ionosphere and neutral atmosphere. Extremely low-frequency (3 Hz–3 kHz) electromagnetic signals are generated in the form of sferics, lightning, and whistlers which can extend to frequencies as high as the VLF range (3–30 kHz).The roughly spherical dielectric cavity bounded by the ground and the ionosphere produces the Schumann resonance at around 8 Hz and its harmonics. A transverse resonance also occurs at 1.7–2.0 kHz arising from reflection off the variable height lower boundary of the ionosphere and exhibiting line splitting due to three-dimensional structure. Ground and satellite observations are discussed in the light of their contributions to understanding the global electric circuit and for EM induction studies.     

磁暴期间太阳风-磁层-电离层间的耦合过程

本文用行星际和地面磁场以及电离层资料,讨论了三次磁暴期间高、中纬电离层电场对太阳风和磁层内变化的响应。 分析表明,当IMF的B_x分量由北向转为南向时,太阳风驱动的磁层对流变化能直接反映出高纬电离层电位的变化。但持续南向的B_x再次增强时,太阳风输入的主要能量耗损于内磁层过程;电离层的响应表现为一个弛豫过程。当B_x由南转北时,环电流的消失对电离层的作用同样有弛豫的特点。此时,驱动电位已撤消,环电流是维持电离层电位的唯一外源。 本文用电路类比及简单模式法结论对上述几种实测情况进行了讨论。     

On the Leontovich boundary condition in geoelectromagnetism

The paper is dedicated to the 70th anniversary of the formulation of the Leontovich impedance boundary condition. This boundary condition is applied in geophysics for sounding of the lithosphere and diagnostics of the magnetosphere by electromagnetic waves. The use of the Leontovich condition is shown to open new possibilities for investigating the anharmonicity of MHD oscillations of the magnetosphere. The differential impedance equation, which is derived with the aid of the Leontovich boundary condition, can prove useful in the investigation of the electroconductivity of the lithosphere by the method of inductive sounding. It is shown within the framework of a simple model that the method of the Leontovich parabolic equation makes it possible to find a correction to the Leontovich boundary condition caused by the fact that, in general, the relation between the horizontal components of an electromagnetic field is not local. It is shown that the anharmonicity of MHD oscillations of the magnetosphere coupled with the nonlocality of the boundary condition on the Earth’s surface can lead to an apparent nonlinearity of the surface impedance calculated in accordance with the classical technique of magnetotelluric sounding.     

Method of characteristics in spherical geometry applied to a Harang-discontinuity situation

The method of characteristics for obtaining spatial distributions of ionospheric electrodynamic parameters from ground-based spatial observations of the ground magnetic disturbance and the ionospheric electric field is presented in spherical geometry. The method includes tools for separation of the external magnetic disturbance, its continuation to the ionosphere, and calculation of ionospheric equivalent currents. Based on these and the measured electric field distribution, the ionospheric Hall conductance is calculated as the primary output of the method. By estimating the Hall- to-Pedersen conductance ratio distribution, the remaining ionospheric electrodynamic parameters are inferred. The method does not assume = 0 to allow to study time-dependent situations. The application of this method to a Harang discontinuity (HD) situation on 27 October 1977, 17:39 UT, reveals the following: (1) The conductances at and north of the HD are clearly reduced as compared to the eastern electrojet region. (2) Plasma flow across the HD is observed, but almost all horizontal current is diverted into upward-flowing field-aligned currents (FACs) there. (3) The FACs connected to the Hall currents form a latitudinally aligned sheet with a magnitude peak between the electrically and magnetically defined HD, where break-up arcs are often observed. Their magnitude is larger than that of the more uniformly distributed FACs connected to the Pedersen currents. They also cause the southward shift of the magnetically defined HD with respect to the electrically defined one. (4) A tilt of the HD with respect to geomagnetic latitude as proposed by an earlier study on the same event, which used composite vector plot technique, and by statistical studies, is not observed in our single time-step analysis.Also at: Finnish Meteorological Institute, Geophysical Research, P.O. Box 503, FIN-00101 Helsinki, Finland.     

Characteristics of a stable are based on FAST and MIRACLE observations

A stable evening sector are is studied using observations from the FAST satellite at 1250 km altitude and the MIRACLE ground-based network, which contains all-sky cameras, coherent radars (STARE), and magnetometers. Both FAST and STARE observe a northward electric field region of about 200 km width and a field magnitude of about 50 mV/m southward of the arc, which is a typical signature for an evening-sector arc. The field-aligned current determined from FAST electron and magnetometer data are in rather good agreement within the arcs. Outside the arcs, the electron data misses the current carriers of the downward FAC probably because it is mainly carried by electrons of smaller energy than the instrument threshold. Studying the westward propagation speed of small undulations associated with the arc using the all-sky cameras gives a velocity of about 2 km//s. This speed is higher than the background ionospheric plasma speed (about 1 km//s), but it agrees rather well with the idea originally proposed by Davis that the undulations reflect an E × B motion in the acceleration region. The ground magnetograms indicate that the main current flows slightly south of the arc. Computing the ionospheric conductivity from FAST electron data and using the ground magnetograms to estimate the current yields an ionospheric electric field pattern, in rather good agreement with FAST results.