Interplanetary Origin of Intense,Superintense and Extreme Geomagnetic Storms

We present a review on the interplanetary causes of intense geomagnetic storms (Dst≤−100 nT), that occurred during solar cycle 23 (1997–2005). It was reported that the most common interplanetary structures leading to the development of intense storms were: magnetic clouds, sheath fields, sheath fields followed by a magnetic cloud and corotating interaction regions at the leading fronts of high speed streams. However, the relative importance of each of those driving structures has been shown to vary with the solar cycle phase. Superintense storms (Dst≤−250 nT) have been also studied in more detail for solar cycle 23, confirming initial studies done about their main interplanetary causes. The storms are associated with magnetic clouds and sheath fields following interplanetary shocks, although they frequently involve consecutive and complex ICME structures. Concerning extreme storms (Dst≤−400 nT), due to the poor statistics of their occurrence during the space era, only some indications about their main interplanetary causes are known. For the most extreme events, we review the Carrington event and also discuss the distribution of historical and space era extreme events in the context of the sunspot and Gleissberg solar activity cycles, highlighting a discussion about the eventual occurrence of more Carrington-type storms.     

New Insights on Geomagnetic Storms from Model Simulations Using Multi-Spacecraft Data

The forecast of the terrestrial ring current as a major contributor to the stormtime Dst index and a predictor of geomagnetic storms is of central interest to ‘space weather’ programs. We thus discuss the dynamical coupling of the solar wind to the Earth's magnetosphere during several geomagnetic storms using our ring current-atmosphere interactions model and coordinated space-borne data sets. Our model calculates the temporal and spatial evolution of H⁺, O⁺, and He⁺ ion distribution functions considering time-dependent inflow from the magnetotail, adiabatic drifts, and outflow from the dayside magnetopause. Losses due to charge exchange, Coulomb collisions, and scattering by EMIC waves are included as well. As initial and boundary conditions we use complementary data sets from spacecraft located at key regions in the inner magnetosphere, Polar and the geosynchronous LANL satellites. We present recent model simulations of the stormtime ring current energization due to the enhanced large-scale convection electric field, which show the transition from an asymmetric to a symmetric ring current during the storm and challenge the standard theories of (a) substorm-driven, and (b) symmetric ring current. Near minimum Dst there is a factor of ∼ 10 variation in the intensity of the dominant ring current ion specie with magnetic local time, its energy density reaching maximum in the premidnight to postmidnight region. We find that the O⁺ content of the ring current increases after interplanetary shocks and reaches largest values near Dst minimum; ∼ 60% of the total ring current energy was carried by O⁺ during the main phase of the 15 July 2000 storm. The effects of magnetospheric convection and losses due to collisions and wave-particle interactions on the global ring current energy balance are calculated during different storm phases and intercompared. This revised version was published online in August 2006 with corrections to the Cover Date.     

Simulation of three-dimensional solar wind disturbances and resulting geomagnetic storms

A kinematic method of representing the three-dimensional solar wind flow is devised by taking into account qualitatively the stream-stream interaction which leads to the formation of a shock pair. Solar wind particles move radially away from the Sun, satisfying the frozen-magnetic field condition. The uniqueness of the present approach is that one can incorporate both theoretical and observational results by adjusting the parameters involved and that a self-consistent data set can be simulated. One can then infer the three-dimensional structure of the solar wind which is vital in understanding the interaction between the solar wind and the magnetosphere, and it is for this reason that the present kinematic method is devised. In the first part of this paper, the present kinematic method is described in detail by demonstrating that the following solar wind features can be simulated: (i) Variations of the solar wind quantities (such as the solar wind speed, the density and the IMF vector), associated with the solar rotation, at the Earth; (ii) the solar wind flow pattern in the meridian planes; (iii) the three-dimensional structure of the corotating interaction region (CIR); and (iv) the three-dimensional structure of the warped solar current sheet.In Section 2, the three-dimensional structure of solar wind disturbances are studied by introducing a flare-generated high speed stream into the two-stream model of the solar wind developed in Section 1. The treatment of the stream-stream interaction is generalized to deal with a flare-generated high speed stream, yielding a shock pair. The shock pair causes three-dimensional distortion of the solar current sheet as it propagates outward from the Sun. It is shown that a set of characteristic time variations of the solar wind speed, density, the interplanetary magnetic field magnitude B and angles (theta) and gf (phi) result at the time of the passage at the location of the Earth for a given set of flare conditions. These quantities allow us to compute the solar wind-magnetosphere energy coupling function . Time variations of the two geomagnetic indices AE and Dst are then estimated from . The simulated geomagnetic storms are compared with observed ones.In the third part, it is shown that recurrent geomagnetic storms can reasonably be reproduced, if fluctuating components of the interplanetary magnetic field (IMF) are superposed on the kinematic model of the solar wind developed in the first part. As an example, we simulate the fluctuating components by linearly polarized Alfvén waves and by random variations of the IMF angle (theta). Characteristics of the simulated and observed geomagnetic storms are discussed in terms of the simulated and observed AE and Dst indices. If the fluctuating components of the IMF can generally be identified as hydromagnetic waves, they may be an important cause for individual magnetospheric substorms, while the IMF magnitude B and the solar wind speed V modulate partially the intensity of magnetospheric substorms and storms.     

The storm-time ring current

In this paper I am reviewing recent advances and open disputes in the study of the terrestrial ring current, with emphasis on its storm-time dynamics. The ring current is carried by energetic charged particles flowing toroidally around the Earth, and creating a ring of westward electric current, centered at the equatorial plane and extending from geocentric distances of about 2 R _E to roughly 9 R _E. This current has a permanent existence due to the natural properties of charged particles in the geospace environment, yet its intensity is variable. It becomes more intense during global electromagnetic disturbances in the near-Earth space, which are known as space (or magnetic or geomagnetic) storms. Changes in this current are responsible for global decreases in the Earth's surface magnetic field, which is the defining feature of geomagnetic storms. The ring current is a critical element in understanding the onset and development of space weather disturbances in geospace. Ring current physics has long been driven by several paradigms, similarly to other disciplines of space physics: the solar origin paradigm, the substorm-driver paradigm, the large-scale symmetry paradigm, the charge-exchange decay paradigm. The paper addresses these paradigms through older and recent important investigations.     

Midlatitude Sporadic E. A Typical Paradigm of Atmosphere-Ionosphere Coupling

This paper provides a comprehensive update on sporadic E layers that is placed in the context of atmosphere-ionosphere coupling, exemplified here by the fundamental windshear theory processes that govern sporadic E layer formation and variability. Some basics of windshear theory are provided first, followed by a summary of key experimental results, their interpretation and physical understanding. The emphasis is placed on the wind shear control of the diurnal and sub-diurnal variability and altitude descent of sporadic E layers and the key role behind these properties of the diurnal and semidiurnal tides. Furthermore, the paper summarizes recent observations that establish a role also for the planetary waves in sporadic E layer occurrence and long-term variability. The possible mechanisms behind this interaction are examined and evidence is presented which shows that planetary waves affect sporadic E layers indirectly though the amplitude modulation of tides at lower altitudes in the MLT region. Only a brief mention is made about gravity wave effects on sporadic E, which apparently exist but cannot be as crucial in layer forming as thought in the past. There is now enough evidence to suggest that mid- and low-latitude sporadic E is not as “sporadic” as the name implies but a regularly occurring ionospheric phenomenon. This may suggest that the sporadic E layer physics can be incorporated in large-scale atmosphere-ionosphere coupling models.     

Why space physics needs to go beyond the MHD box

Magnetohydrodynamic (MHD) theory has been used in space physics for more than forty years, yet many important questions about space plasmas remain unanswered. We still do not understand how the solar wind is accelerated, how mass, momentum and energy are transported into the magnetosphere and what mechanisms initiate substorms. Questions have been raised from the beginning of the space era whether MHD theory can describe correctly space plasmas that are collisionless and rarely in thermal equilibrium. Ideal MHD fluids do not induce electromotive force, hence they lose the capability to interact electromagnetically. No currents and magnetic fields are generated, rendering ideal MHD theory not very useful for space plasmas. Observations from the plasma sheet are used as examples to show how collisionless plasmas behave. Interpreting these observations using MHD and ideal MHD concepts can lead to misleading conclusions. Notably, the bursty bulk flows (BBF) with large mean velocities left(〈 v 〉≥400 km s right) that have been interpreted previously as E×B flows are shown to involve much more complicated physics. The sources of these nonvanishing 〈 v 〉 events, while still not known, are intimately related to mechanisms that create large phase space gradients that include beams and acceleration of ions to MeV energies. The distributions of these nonvanishing 〈 v 〉 events are associated with large amplitude variations of the magnetic field at frequencies up to and exceeding the local Larmor frequency where MHD theory is not valid. Understanding collisionless plasma dynamics such as substorms in the plasma sheet requires the self-consistency that only kinetic theory can provide. Kinetic modeling is still undergoing continual development with many studies limited to one and two dimensions, but there is urgent need to improve these models as more and more data show kinetic physics is fundamentally important. Only then will we be able to make progress and obtain a correct picture of how collisionless plasmas work in space.     

Order out of Randomness: Self-Organization Processes in Astrophysics

Self-organization is a property of dissipative nonlinear processes that are governed by a global driving force and a local positive feedback mechanism, which creates regular geometric and/or temporal patterns, and decreases the entropy locally, in contrast to random processes. Here we investigate for the first time a comprehensive number of (17) self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous “order out of randomness”, during the evolution from an initially disordered system to an ordered quasi-stationary system, mostly by quasi-periodic limit-cycle dynamics, but also by harmonic (mechanical or gyromagnetic) resonances. The global driving force can be due to gravity, electromagnetic forces, mechanical forces (e.g., rotation or differential rotation), thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational (Balbus-Hawley) instability, the convective (Rayleigh-Bénard) instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or a loss-cone instability. Physical models of astrophysical self-organization processes require hydrodynamic, magneto-hydrodynamic (MHD), plasma, or N-body simulations. Analytical formulations of self-organizing systems generally involve coupled differential equations with limit-cycle solutions of the Lotka-Volterra or Hopf-bifurcation type.     

Satellite observations of power line harmonic radiation

In-situ spectral observations of power-line harmonic radiation (PLHR) are still quite rare and almost all the detailed characteristics have been derived from studies at Antarctic stations such as Siple and Halley, and their conjugates in North America. Because of the lack of more direct satellite evidence of PLHR and the difficulties in interpretation of morphological studies, such as those of Ariel 3 and 4, there is considerable controversy concerning the relative importance of PLHR and its contribution to wave-particle interactions (WPI) in the magnetosphere. The early Ariel 3 and 4 global surveys indicated that, in terms of true mean wave energy, there is no longitudinal localisation, the contribution of world-wide intense VLF emissions, associated with magnetic storms, being dominant. Also, the most intense wave emission, that of plasmaspheric hiss at ELF (< 1 kHz) exhibits little evidence of localisation. The PLHR phenomenon is most conspicuous by its persistence in quiet times (Kp ≤ 2+) at 45° < Λ < 55° over North America and its conjugate region, even though the less frequent strongest emissions, to which it gives rise in the summer, are located polewards at 3 < L < 5. In the northern winter, when spheric activity over both North America and its conjugate are low, there is a high occurrence of strong discrete emissions, which are more sharply localised than in the summer, on the NE industrial U.S.A. field line. The most recent Ariel 4 studies, particularly on the spheric wavefield over North America (using data from the Appleton Laboratory impulse counters) and on the character of the wavefield over the mainland and over the Atlantic immediately to the east (where the spheric contribution is similar) throw new light on the problem. It appears that the principal role of the PLHR may be to sustain duct structure and multihop propagation which is relatively much rarer over the Atlantic. Typical industrial PLHR consists of a series of narrow pulses at twice the mains frequency. It is suggested that these ‘artificial spherics’ may help to sustain the WPI and multihop duct structure. At L = 4, Yoshida et al. (1980) have shown that there is a strong, sharp maximum for WPIs originating in spherics, at f ? 3 kHz, in good agreement with Siple observations.     

Ring Current Dynamics

This chapter reviews the current understanding of ring current dynamics. The terrestrial ring current is an electric current flowing toroidally around the Earth, centered at the equatorial plane and at altitudes of ∼10,000 to 60,000 km. Enhancements in this current are responsible for global decreases in the Earth’s surface magnetic field, which have been used to define geomagnetic storms. Intense geospace magnetic storms have severe effects on technological systems, such as disturbances or even permanent damage of telecommunication and navigation satellites, telecommunication cables, and power grids. The main carriers of the ring current are positive ions, with energies from ∼1 keV to a few hundred keV, which are trapped by the geomagnetic field and undergo an azimuthal drift. The ring current is formed by the injection of ions originating in the solar wind and the terrestrial ionosphere into the inner magnetosphere. The injection process involves electric fields, associated with enhanced magnetospheric convection and/or magnetospheric substorms. The quiescent ring current is carried mainly by protons of predominantly solar wind origin, while active processes in geospace tend to increase the abundance (both absolute and relative) of O⁺ ions, which are of ionospheric origin. During intense geospace magnetic storms, the O⁺ abundance increases dramatically. This increase has been observed to occur concurrently with the rapid intensification of the ring current in the storm main phase and to result in O⁺ dominance around storm maximum. This compositional change can affect several dynamic processes, such as species-and energy-dependent charge-exchange and wave-particle scattering loss.     

Numerical Magnetohydrodynamic (MHD) Modeling of Coronal Mass Ejections (CMEs)

The Coronal Mass Ejection (CME) is arguably the most important discovery of solar eruptive phenomena in the 20th century. It is now also recognized that CMEs have great impact on the Earth's environment by inducing geomagnetic storms. Thus, development of simulation models to understand the physical mechanisms of CME initiation and propagation has become a challenge in the solar MHD community. In this paper we shall summarize chronologically the development of the theoretical analyses, and the successes and failures of the numerical magnetohydrodynamic (MHD) simulations of coronal mass ejections (CMEs) during the past two decades. The chronological development of numerical simulation models and the evolution of the numerical methods to treat this class of problems are presented. The most appropriate way to model CMEs is to have (i) a realistic pre-event coronal atmosphere, and (ii) realistic driving mechanisms. Details of the progress and assessment of the theoretical and modeling efforts for the understanding of the physics of the CME initiation and propagation will be presented, and the numerical methods to construct these simulation models will be discussed.     

Interplanetary origin of geomagnetic storms

Around solar maximum, the dominant interplanetary phenomena causing intense magnetic storms (Dst<−100 nT) are the interplanetary manifestations of fast coronal mass ejections (CMEs). Two interplanetary structures are important for the development of storms, involving intense southward IMFs: the sheath region just behind the forward shock, and the CME ejecta itself. Whereas the initial phase of a storm is caused by the increase in plasma ram pressure associated with the increase in density and speed at and behind the shock (accompanied by a sudden impulse [SI] at Earth), the storm main phase is due to southward IMFs. If the fields are southward in both of the sheath and solar ejecta, two-step main phase storms can result and the storm intensity can be higher. The storm recovery phase begins when the IMF turns less southward, with delays of ≈1–2 hours, and has typically a decay time of 10 hours. For CMEs involving clouds the intensity of the core magnetic field and the amplitude of the speed of the cloud seems to be related, with a tendency that clouds which move at higher speeds also posses higher core magnetic field strengths, thus both contributing to the development of intense storms since those two parameters are important factors in genering the solar wind-magnetosphere coupling via the reconnection process. During solar minimum, high speed streams from coronal holes dominate the interplanetary medium activity. The high-density, low-speed streams associated with the heliospheric current sheet (HCS) plasma impinging upon the Earth's magnetosphere cause positive Dst values (storm initial phases if followed by main phases). In the absence of shocks, SIs are infrequent during this phase of the solar cycle. High-field regions called Corotating Interaction Regions (CIRs) are mainly created by the fast stream (emanating from a coronal hole) interaction with the HCS plasma sheet. However, because the B_z component is typically highly fluctuating within the CIRs, the main phases of the resultant magnetic storms typically have highly irregular profiles and are weaker. Storm recovery phases during this phase of the solar cycle are also quite different in that they can last from many days to weeks. The southward magnetic field (B_s) component of Alfvén waves in the high speed stream proper cause intermittent reconnection, intermittent substorm activity, and sporadic injections of plasma sheet energy into the outer portion of the ring current, prolonging its final decay to quiet day values. This continuous auroral activity is called High Intensity Long Duration Continuous AE Activity (HILDCAAs). Possible interplanetary mechanisms for the creation of very intense magnetic storms are discussed. We examine the effects of a combination of a long-duration southward sheath magnetic field, followed by a magnetic cloud B_s event. We also consider the effects of interplanetary shock events on the sheath plasma. Examination of profiles of very intense storms from 1957 to the present indicate that double, and sometimes triple, IMF B_s events are important causes of such events. We also discuss evidence that magnetic clouds with very intense core magnetic fields tend to have large velocities, thus implying large amplitude interplanetary electric fields that can drive very intense storms. Finally, we argue that a combination of complex interplanetary structures, involving in rare occasions the interplanetary manifestations of subsequent CMEs, can lead to extremely intense storms. This revised version was published online in June 2006 with corrections to the Cover Date.     

Absolute intensity measurements in the extreme ultraviolet spectrum of solar radiation

Experimental results and problems of absolute intensity measurements of solar electromagnetic radiation in the extreme ultraviolet (EUV) and soft X-ray region of the spectrum (designated cumulatively as XUV for brevity) are reviewed. The numerous practical problems are divided in two major areas, (a) general problems of heterochromatic absolute XUV spectrophotometry in the laboratory and (b) specific problems characteristic of requirements of solar physics, the physics of planetary atmospheres, and existing restrictions of space technology. Within the first area (a) emphasis is placed on recent progress toward justified reliance on ionization detectors without necessary connection to source standards. For the second area (b), emphasis is placed on the immediate need to have existing exploratory observations followed by a new phase of more systematic experiments of increased accuracy.     

Electric Fields and Magnetic Fields in the Plasmasphere: A Perspective From CLUSTER and IMAGE

The electric field and magnetic field are basic quantities in the plasmasphere measured since the 1960s. In this review, we first recall conventional wisdom and remaining problems from ground-based whistler measurements. Then we show scientific results from Cluster and Image, which are specifically made possible by newly introduced features on these spacecraft, as follows. 1. In situ electric field measurements using artificial electron beams are successfully used to identify electric fields originating from various sources. 2. Global electric fields are derived from sequences of plasmaspheric images, revealing how the inner magnetospheric electric field responds to the southward interplanetary magnetic fields and storms/substorms. 3. Understanding of sub-auroral polarization stream (SAPS) or sub-auroral ion drifts (SAID) are advanced through analysis of a combination of magnetospheric and ionospheric measurements from Cluster, Image, and DMSP. 4. Data from multiple spacecraft have been used to estimate magnetic gradients for the first time.     

Nonlinear problems of physics of the geomagnetic pulsations

The current state of research involving manifestations of nonlinearity in geomagnetic pulsations is reviewed. Traditionally, the attention of researchers was focused on the effects of resonant interaction of geomagnetic pulsations with small groups of energetic particles, which actually means the study of the quasi-linear relaxation of radiation belt ions, the modulation of auroral electron fluxes, etc. The present review concentrates on the problem of the nonlinear effect influence of pulsations on the backgroud (cold) plasma and on the geomagnetic field. This kind of interaction results in a significant modification of the plasma distribution in the magnetosphere. Self-consistent wave structures—solitons and vortices may occur as well. Such nonlinear effects contribute to physics of geomagnetic pulsations and are also of fundamental importance for general physics. Another set of more narrow problems considered in the review, is related to phenomenological modeling of fluctuational and critical phenomena in the magnetosphere. The essence of our approach is to present the magnetosphere as a black box, whose properties should be determined by the statistical characteristics of its output signals. This approach to phenomenology can be a useful supplement to the methods of microscopic modeling aimed at detecting nonlinear manifestations of geomagnetic pulsations.     

Gravity Probe B Data Analysis

This is the first of five connected papers detailing progress on the Gravity Probe B (GP-B) Relativity Mission. GP-B, launched 20 April 2004, is a landmark physics experiment in space to test two fundamental predictions of Einstein’s general relativity theory, the geodetic and frame-dragging effects, by means of cryogenic gyroscopes in Earth orbit. Data collection began 28 August 2004 and science operations were completed 29 September 2005. The data analysis has proven deeper than expected as a result of two mutually reinforcing complications in gyroscope performance: (1) a changing polhode path affecting the calibration of the gyroscope scale factor C _g against the aberration of starlight and (2) two larger than expected manifestations of a Newtonian gyro torque due to patch potentials on the rotor and housing. In earlier papers, we reported two methods, ‘geometric’ and ‘algebraic’, for identifying and removing the first Newtonian effect (‘misalignment torque’), and also a preliminary method of treating the second (‘roll-polhode resonance torque’). Central to the progress in both torque modeling and C _g determination has been an extended effort on “Trapped Flux Mapping” commenced in November 2006. A turning point came in August 2008 when it became possible to include a detailed history of the resonance torques into the computation. The East-West (frame-dragging) effect is now plainly visible in the processed data. The current statistical uncertainty from an analysis of 155 days of data is 5.4 marc-s/yr (～14% of the predicted effect), though it must be emphasized that this is a preliminary result requiring rigorous investigation of systematics by methods discussed in the accompanying paper by Muhlfelder et al. A covariance analysis incorporating models of the patch effect torques indicates that a 3–5% determination of frame-dragging is possible with more complete, computationally intensive data analysis.     

Physics of Magnetospheric Variability

Many widely used methods for describing and understanding the magnetosphere are based on balance conditions for quasi-static equilibrium (this is particularly true of the classical theory of magnetosphere/ionosphere coupling, which in addition presupposes the equilibrium to be stable); they may therefore be of limited applicability for dealing with time-variable phenomena as well as for determining cause-effect relations. The large-scale variability of the magnetosphere can be produced both by changing external (solar-wind) conditions and by non-equilibrium internal dynamics. Its developments are governed by the basic equations of physics, especially Maxwell’s equations combined with the unique constraints of large-scale plasma; the requirement of charge quasi-neutrality constrains the electric field to be determined by plasma dynamics (generalized Ohm’s law) and the electric current to match the existing curl of the magnetic field. The structure and dynamics of the ionosphere/magnetosphere/solar-wind system can then be described in terms of three interrelated processes: (1) stress equilibrium and disequilibrium, (2) magnetic flux transport, (3) energy conversion and dissipation. This provides a framework for a unified formulation of settled as well as of controversial issues concerning, e.g., magnetospheric substorms and magnetic storms.     

Modelling of the electromagnetic field in the interplanetary space and in the Earth's magnetosphere

A magnetohydrodynamic model of the solar wind flow is constructed using a kinematic approach. It is shown that a phenomenological conductivity of the solar wind plasma plays a key role in the forming of the interplanetary magnetic field (IMF) component normal to the ecliptic plane. This component is mostly important for the magnetospheric dynamics which is controlled by the solar wind electric field. A simple analytical solution for the problem of the solar wind flow past the magnetosphere is presented. In this approach the magnetopause and the Earth's bow shock are approximated by the paraboloids of revolution. Superposition of the effects of the bulk solar wind plasma motion and the magnetic field diffusion results in an incomplete screening of the IMF by the magnetopause. It is shown that the normal to the magnetopause component of the solar wind magnetic field and the tangential component of the electric field penetrated into the magnetosphere are determined by the quarter square of the magnetic Reynolds number. In final, a dynamic model of the magnetospheric magnetic field is constructed. This model can describe the magnetosphere in the course of the severe magnetic storm. The conditions under which the magnetospheric magnetic flux structure is unstable and can drive the magnetospheric substorm are discussed. The model calculations are compared with the observational data for September 24–26, 1998 magnetic storm (Dst _min=−205 nT) and substorm occurred at 02:30 UT on January 10, 1997. This revised version was published online in August 2006 with corrections to the Cover Date.