The Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission is the fifth NASA Medium-class Explorer (MIDEX), launched on February 17, 2007 to determine the trigger and large-scale evolution of substorms. The mission employs five identical micro-satellites (hereafter termed “probes”) which line up along the Earth’s magnetotail to track the motion of particles, plasma and waves from one point to another and for the first time resolve space–time ambiguities in key regions of the magnetosphere on a global scale. The probes are equipped with comprehensive in-situ particles and fields instruments that measure the thermal and super-thermal ions and electrons, and electromagnetic fields from DC to beyond the electron cyclotron frequency in the regions of interest. The primary goal of THEMIS, which drove the mission design, is to elucidate which magnetotail process is responsible for substorm onset at the region where substorm auroras map (∼10 RE): (i) a local disruption of the plasma sheet current (current disruption) or (ii) the interaction of the current sheet with the rapid influx of plasma emanating from reconnection at ∼25 RE. However, the probes also traverse the radiation belts and the dayside magnetosphere, allowing THEMIS to address additional baseline objectives, namely: how the radiation belts are energized on time scales of 2–4 hours during the recovery phase of storms, and how the pristine solar wind’s interaction with upstream beams, waves and the bow shock affects Sun–Earth coupling. THEMIS’s open data policy, platform-independent dataset, open-source analysis software, automated plotting and dissemination of data within hours of receipt, dedicated ground-based observatory network and strong links to ancillary space-based and ground-based programs. promote a grass-roots integration of relevant NASA, NSF and international assets in the context of an international Heliophysics Observatory over the next decade. The mission has demonstrated spacecraft and mission design strategies ideal for Constellation-class missions and its science is complementary to Cluster and MMS. THEMIS, the first NASA micro-satellite constellation, is a technological pathfinder for future Sun-Earth Connections missions and a stepping stone towards understanding Space Weather.

The THEMIS Mission

In this review we will focus on a topic of fundamental importance for both plasma physics and astrophysics, namely the occurrence of large-amplitude low-frequency fluctuations of the fields that describe the plasma state. This subject will be treated within the context of the expanding solar wind and the most meaningful advances in this research field will be reported emphasizing the results obtained in the past decade or so. As a matter of fact, Ulysses’ high latitude observations and new numerical approaches to the problem, based on the dynamics of complex systems, brought new important insights which helped to better understand how turbulent fluctuations behave in the solar wind. In particular, numerical simulations within the realm of magnetohydrodynamic (MHD) turbulence theory unraveled what kind of physical mechanisms are at the basis of turbulence generation and energy transfer across the spectral domain of the fluctuations. In other words, the advances reached in these past years in the investigation of solar wind turbulence now offer a rather complete picture of the phenomenological aspect of the problem to be tentatively presented in a rather organic way.

/pdf/the-solar-wind-as-a-turbulence-laboratory-4mad4q4xdr.pdf

The Solar Wind as a Turbulence Laboratory

This paper presents a theoretical framework for understanding plasma turbulence in astrophysical plasmas. It is motivated by observations of electromagnetic and density fluctuations in the solar wind, interstellar medium and galaxy clusters, as well as by models of particle heating in accretion disks. All of these plasmas and many others have turbulent motions at weakly collisional and collisionless scales. The paper focuses on turbulence in a strong mean magnetic field. The key assumptions are that the turbulent fluctuations are small compared to the mean field, spatially anisotropic with respect to it and that their frequency is low compared to the ion cyclotron frequency. The turbulence is assumed to be forced at some system-specific outer scale. The energy injected at this scale has to be dissipated into heat, which ultimately cannot be accomplished without collisions. A kinetic cascade develops that brings the energy to collisional scales both in space and velocity. The nature of the kinetic cascade in various scale ranges depends on the physics of plasma fluctuations that exist there. There are four special scales that separate physically distinct regimes: the electron and ion gyroscales, the mean free path and the electron diffusion scale. In each of the scale ranges separated by these scales, the fully kinetic problem is systematically reduced to a more physically transparent and computationally tractable system of equations, which are derived in a rigorous way. In the inertial range above the ion gyroscale, the kinetic cascade separates into two parts: a cascade of Alfvenic fluctuations and a passive cascade of density and magnetic-field-strength fluctuations. The former are governed by the reduced magnetohydrodynamic (RMHD) equations at both the collisional and collisionless scales; the latter obey a linear kinetic equation along the (moving) field lines associated with the Alfvenic component (in the collisional limit, these compressive fluctuations become the slow and entropy modes of the conventional MHD). In the dissipation range below ion gyroscale, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like electron reduced magnetohydrodynamic (ERMHD) equations and a passive cascade of ion entropy fluctuations both in space and velocity. The latter cascade brings the energy of the inertial-range fluctuations that was Landau-damped at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for all of these cascades. The relationship between the theoretical models proposed in this paper and astrophysical applications and observations is discussed in detail.

/pdf/astrophysical-gyrokinetics-kinetic-and-fluid-turbulent-48o03f94s0.pdf

Astrophysical Gyrokinetics: Kinetic and Fluid Turbulent Cascades in Magnetized Weakly Collisional Plasmas

Kinetic plasma physics of the solar corona and solar wind are reviewed with emphasis on the theoretical understanding of the in situ measurements of solar wind particles and waves, as well as on the remote-sensing observations of the solar corona made by means of ultraviolet spectroscopy and imaging. In order to explain coronal and interplanetary heating, the microphysics of the dissipation of various forms of mechanical, electric and magnetic energy at small scales (e.g., contained in plasma waves, turbulences or non-uniform flows) must be addressed. We therefore scrutinise the basic assumptions underlying the classical transport theory and the related collisional heating rates, and also describe alternatives associated with wave-particle interactions. We elucidate the kinetic aspects of heating the solar corona and interplanetary plasma through Landau- and cyclotron-resonant damping of plasma waves, and analyse in detail wave absorption and micro instabilities. Important aspects (virtues and limitations) of fluid models, either single- and multi-species or magnetohydrodynamic and multi-moment models, for coronal heating and solar wind acceleration are critically discussed. Also, kinetic model results which were recently obtained by numerically solving the Vlasov‐Boltzmann equation in a coronal funnel and hole are presented. Promising areas and perspectives for future research are outlined finally.

/pdf/kinetic-physics-of-the-solar-corona-and-solar-wind-bduwadfu6p.pdf

Kinetic Physics of the Solar Corona and Solar Wind

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 Bz 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 (Bs) component of Alfven 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 Bs 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 Bs 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.

Interplanetary origin of geomagnetic storms

Twenty one differences between CME-driven geomagnetic storms and CIR-driven geomagnetic storms are tabulated. (CME-driven includes driving by CME sheaths, by magnetic clouds, and by ejecta; CIR-driven includes driving by the associated recurring high-speed streams.) These differences involve the bow shock, the magnetosheath, the radiation belts, the ring current, the aurora, the Earth's plasma sheet, magnetospheric convection, ULF pulsations, spacecraft charging in the magnetosphere, and the saturation of the polar cap potential. CME-driven storms are brief, have denser plasma sheets, have strong ring currents and Dst, have solar energetic particle events, and can produce great auroras and dangerous geomagnetically induced currents; CIR-driven storms are of longer duration, have hotter plasmas and stronger spacecraft charging, and produce high fluxes of relativistic electrons. Further, the magnetosphere is more likely to be preconditioned with dense plasmas prior to CIR-driven storms than it is prior to CME-driven storms. CME-driven storms pose more of a problem for Earth-based electrical systems; CIR-driven storms pose more of a problem for space-based assets.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2005JA011447

Differences between CME‐driven storms and CIR‐driven storms

The coupling of the plasma sheet to the solar wind is studied statistically using measurements from various satellite pairs: one satellite in the solar wind and one in either the magnetotail central plasma sheet or the near-Earth plasma sheet. It is found that the properties of the plasma sheet are highly correlated with the properties of the solar wind: specifically that (1) the density of the plasma sheet is strongly correlated with the density of the solar wind, (2) the temperature of the plasma sheet is strongly correlated with the velocity of the solar wind, (3) the particle pressure and total pressure of the plasma sheet are strongly correlated with the ram pressure of the solar wind, (4) B y in the plasma sheet is strongly correlated with the B y in the solar wind, (5) B z in the plasma sheet is weakly correlated with the B z in the solar wind, (6) E y in the plasma sheet is weakly correlated with the E y in the solar wind, and (7) plasma sheet earthward-tailward flow velocity is weakly anticorrelated with the solar wind velocity. After removing these solar wind dependencies, the dependencies of the properties of the plasma sheet on geomagnetic activity are reduced and changed. The time lags between the solar wind density and the plasma sheet density are investigated statistically and on a case-by-case basis; it is found that solar wind material reaches the midtail plasma sheet in ∼2 hours, reaches the near-Earth nightside plasma sheet in ∼4 hours, and reaches the dayside plasma sheet in ∼15 hours. The pathway for mass transport into the plasma sheet is explored; it is suggested that solar wind material is transported rapidly up the tail to the near-Earth region via the plasma sheet boundary layer, followed by poloidal transport via E x B convection in the near-Earth (low β) plasma sheet or followed by eddy diffusion in the not-too-distant (high β) plasma sheet. Analogies are drawn between the high Reynolds number plasma sheet and a high Reynolds number fluid wake.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/97JA02986

The driving of the plasma sheet by the solar wind

Twelve electron-acceleration mechanisms and 10 generator mechanisms for auroral arcs are examined, and a characteristic auroral-arc thickness is worked out for each mechanism except one. The arc thicknesses are then mapped down to the ionosphere along the terrestrial magnetic-field lines; near the Earth, a dipole magnetic-field model is used, and farther from the Earth, the mapping includes the effects of magnetic-field-line draping. The 21 theoretical models all predict auroral-arc thicknesses that are at least an order of magnitude wider than the optically observed arcs. As an alternative explanation of the observed narrow auroral arcs, the acceleration of ionospheric electrons to produce airglow in electrical discharge mechanisms appears improbable. Also unsuccessfully explored is the possibility that the observed narrow auroral arcs are caused by interference effects when Alfven waves reflect off the ionosphere. Suggestions are made for future ground-based auroral-arc measurements.

Auroral arc thicknesses as predicted by various theories

[1] It is argued here that the inner heliosphere is filled with a network of entangled magnetic flux tubes and that the flux tubes are fossil structures that originate at the solar surface. 65,860 flux tubes are collected from seven years of measurements with the ACE spacecraft at 1 AU by spotting the flux tube walls with large changes in the magnetic field direction and the vector flow velocity. The tube walls are associated with large changes in the ion entropy density and the alpha-to-proton ratio. The median size of the flux tubes at 1 AU is 4.4 x 10 5 km. The tubes are larger in slow wind than in fast wind. The tubes are on-average aligned with the Parker spiral, with a large spread in orientations. This large spread may be caused by slight misalignments of tubes in the corona. The flux tubes map to granule and supergranule sizes on the Sun. The amounts of magnetic flux in the tubes at 1 AU correspond to the amounts of magnetic flux in field concentrations in the magnetic carpet. It is argued that the flux tubes do not reconnect during the ∼ 100-h advection to 1 AU owing to the expansion of the solar wind. The flux tube texture impacts the flow properties of the solar wind, turbulence in the solar wind, energetic-particle propagation in the inner heliosphere, and the driving of the Earth's magnetosphere. A method for using measurements of the flux tube walls for the remote sensing of magnetic field dynamics in the magnetic carpet is suggested.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2007JA012684

Flux tube texture of the solar wind: Strands of the magnetic carpet at 1 AU?

The bulk flows and magnetic-field fluctuations of the plasma sheet are investigated using single-point measurements from the ISEE-2 Fast Plasma Experiment and fluxgate magnetometer. Ten several-hour-long intervals of continuous data (with 3 s and 12 s time resolution) are analysed. The plasma-sheet flow appears to be strongly ‘turbulent’ (i.e. the flow is dominated by fluctuations that are unpredictable, with rms velocities[Gt ]mean velocities and with field fluctuations≈mean fields). The flow velocities are typically sub-Alfvenic. The flow-velocity probability distribution P(v) is constructed, and is found to be well fitted by exponential functions. Autocorrelation functions [Ascr ](τ) are constructed, and the autocorrelation times τcorr for the flow velocities are found to be about 2 min. From the flow measurements, an estimate of the mixing length in the plasma sheet is produced, yielding Lmix≈2 Earth radii; correspondingly, the plasma-sheet material appears to be well mixed in density and temperature. An eddy viscosity for the plasma sheet is also estimated. Power spectra, which are constructed from the v(t) and B(t) time series, have portions that are power laws with spectral indices that are near the range of those expected for turbulence theories. The plasma sheet may provide a laboratory for the study of turbulence in parameter regimes different from that of solar-wind turbulence: the plasma sheet is a β[Gt ]1, hot-ion plasma, and the turbulence may be strongly driven rather than well developed. The turbulent nature of the flow and the disordered nature of the magnetic field have implications for the transport of plasma-sheet material, for the penetration of the solar-wind electric field into the plasma sheet, and for the calculation of particle orbits in the magnetotail.

Joseph E. Borovsky

Papers

Differences between CME‐driven storms and CIR‐driven storms

The driving of the plasma sheet by the solar wind

Auroral arc thicknesses as predicted by various theories

Flux tube texture of the solar wind: Strands of the magnetic carpet at 1 AU?

The Earth's plasma sheet as a laboratory for flow turbulence in high-β MHD