Future Exploration of the Outer Heliosphere and Very Local Interstellar Medium by Interstellar Probe

A detailed overview of the knowledge gaps in our understanding of the heliospheric interaction with the largely unexplored Very Local Interstellar Medium (VLISM) are provided along with predictions of with the scientific discoveries that await. The new measurements required to make progress in this expanding frontier of space physics are discussed and include in-situ plasma and pick-up ion measurements throughout the heliosheath, direct sampling of the VLISM properties such as elemental and isotopic composition, densities, flows, and temperatures of neutral gas, dust and plasma, and remote energetic neutral atom (ENA) and Lyman-alpha (LYA) imaging from vantage points that can uniquely discern the heliospheric shape and bring new information on the interaction with interstellar hydrogen. The implementation of a pragmatic Interstellar Probe mission with a nominal design life to reach 375 Astronomical Units (au) with likely operation out to 550 au are reported as a result of a 4-year NASA funded mission study.

Magnetopause ripples going against the flow form azimuthally stationary surface waves

Surface waves process the turbulent disturbances which drive dynamics in many space, astrophysical and laboratory plasma systems, with the outer boundary of Earth’s magnetosphere, the magnetopause, providing an accessible environment to study them. Like waves on water, magnetopause surface waves are thought to travel in the direction of the driving solar wind, hence a paradigm in global magnetospheric dynamics of tailward propagation has been well-established. Here we show through multi-spacecraft observations, global simulations, and analytic theory that the lowest-frequency impulsively-excited magnetopause surface waves, with standing structure along the terrestrial magnetic field, propagate against the flow outside the boundary. Across a wide local time range (09–15h) the waves’ Poynting flux exactly balances the flow’s advective effect, leading to no net energy flux and thus stationary structure across the field also. Further down the equatorial flanks, however, advection dominates hence the waves travel downtail, seeding fluctuations at the resonant frequency which subsequently grow in amplitude via the Kelvin-Helmholtz instability and couple to magnetospheric body waves. This global response, contrary to the accepted paradigm, has implications on radiation belt, ionospheric, and auroral dynamics and potential applications to other dynamical systems.

The magnetopause surface waves (SW) that drive global plasma dynamics are thought, like waves on water, to travel with the driving solar wind. Here, the authors show that impulsively-excited SW, with standing structure along the geomagnetic field, are stationary by propagating against this flow.

Decay of Kelvin-Helmholtz Vortices at the Earth's Magnetopause Under Pure Southward IMF Conditions

At the Earth's low-latitude magnetopause, clear signatures of the Kelvin-Helmholtz (KH) waves have been frequently observed during periods of the northward interplanetary magnetic field (IMF), whereas these signatures have been much less frequently observed during the southward IMF. Here, we performed the first 3-D fully kinetic simulation of the magnetopause KH instability under the southward IMF condition. The simulation demonstrates that fast magnetic reconnection is induced at multiple locations along the vortex edge in an early nonlinear growth phase of the instability. The reconnection outflow jets significantly disrupt the flow of the nonlinear KH vortex, while the disrupted turbulent flow strongly bends and twists the reconnected field lines. The resulting coupling of the complex field and flow patterns within the magnetopause boundary layer leads to a quick decay of the vortex structure, which may explain the difference in the observation probability of KH waves between northward and southward IMF conditions.

Direct observations of a surface eigenmode of the dayside magnetopause

The abrupt boundary between a magnetosphere and the surrounding plasma, the magnetopause, has long been known to support surface waves. It was proposed that impulses acting on the boundary might lead to a trapping of these waves on the dayside by the ionosphere, resulting in a standing wave or eigenmode of the magnetopause surface. No direct observational evidence of this has been found to date and searches for indirect evidence have proved inconclusive, leading to speculation that this mechanism might not occur. By using fortuitous multipoint spacecraft observations during a rare isolated fast plasma jet impinging on the boundary, here we show that the resulting magnetopause motion and magnetospheric ultra-low frequency waves at well-defined frequencies are in agreement with and can only be explained by the magnetopause surface eigenmode. We therefore show through direct observations that this mechanism, which should impact upon the magnetospheric system globally, does in fact occur.

Surface waves on the boundary between a magnetosphere and the surrounding plasma might get trapped by the ionosphere forming an eigenmode. Here, Archer et al. show direct observations of this proposed mechanism at Earth’s magnetosphere by analyzing the response to an isolated fast plasma jet detected by the THEMIS satellites.

The Space Physics Environment Data Analysis System (SPEDAS)

With the advent of the Heliophysics/Geospace System Observatory (H/GSO), a complement of multi-spacecraft missions and ground-based observatories to study the space environment, data retrieval, analysis, and visualization of space physics data can be daunting. The Space Physics Environment Data Analysis System (SPEDAS), a grass-roots software development platform (www.spedas.org), is now officially supported by NASA Heliophysics as part of its data environment infrastructure. It serves more than a dozen space missions and ground observatories and can integrate the full complement of past and upcoming space physics missions with minimal resources, following clear, simple, and well-proven guidelines. Free, modular and configurable to the needs of individual missions, it works in both command-line (ideal for experienced users) and Graphical User Interface (GUI) mode (reducing the learning curve for first-time users). Both options have "crib-sheets," user-command sequences in ASCII format that can facilitate record-and-repeat actions, especially for complex operations and plotting. Crib-sheets enhance scientific interactions, as users can move rapidly and accurately from exchanges of technical information on data processing to efficient discussions regarding data interpretation and science. SPEDAS can readily query and ingest all International Solar Terrestrial Physics (ISTP)-compatible products from the Space Physics Data Facility (SPDF), enabling access to a vast collection of historic and current mission data. The planned incorporation of Heliophysics Application Programmer's Interface (HAPI) standards will facilitate data ingestion from distributed datasets that adhere to these standards. Although SPEDAS is currently Interactive Data Language (IDL)-based (and interfaces to Java-based tools such as Autoplot), efforts are under-way to expand it further to work with python (first as an interface tool and potentially even receiving an under-the-hood replacement). We review the SPEDAS development history, goals, and current implementation. We explain its "modes of use" with examples geared for users and outline its technical implementation and requirements with software developers in mind. We also describe SPEDAS personnel and software management, interfaces with other organizations, resources and support structure available to the community, and future development plans.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (10.1007/s11214-018-0576-4) contains supplementary material, which is available to authorized users.

How Accurately Can We Measure the Reconnection Rate M for the MMS Diffusion Region Event of 11 July 2017?

We investigate the accuracy with which the reconnection electric field E M can be determined from in situ plasma data. We study the magnetotail electron diffusion region observed by National Aeronautics and Space Administration's Magnetospheric Multiscale (MMS) on 11 July 2017 at 22:34 UT and focus on the very large errors in E M that result from errors in an L M N boundary normal coordinate system. We determine several L M N coordinates for this MMS event using several different methods. We use these M axes to estimate E M. We find some consensus that the reconnection rate was roughly E M = 3.2 ± 0.6 mV/m, which corresponds to a normalized reconnection rate of 0.18 ± 0.035. Minimum variance analysis of the electron velocity (MVA-v e), MVA of E, minimization of Faraday residue, and an adjusted version of the maximum directional derivative of the magnetic field (MDD-B) technique all produce reasonably similar coordinate axes. We use virtual MMS data from a particle-in-cell simulation of this event to estimate the errors in the coordinate axes and reconnection rate associated with MVA-v e and MDD-B. The L and M directions are most reliably determined by MVA-v e when the spacecraft observes a clear electron jet reversal. When the magnetic field data have errors as small as 0.5% of the background field strength, the M direction obtained by MDD-B technique may be off by as much as 35°. The normal direction is most accurately obtained by MDD-B. Overall, we find that these techniques were able to identify E M from the virtual data within error bars ≥20%.

Simultaneous Remote Observations of Intense Reconnection Effects by DMSP and MMS Spacecraft During a Storm Time Substorm

During a magnetic storm on 23 June 2015, several very intense substorms took place, with signatures observed by multiple spacecraft including DMSP and Magnetospheric Multiscale (MMS). At the time of interest, DMSP F18 crossed inbound through a poleward expanding auroral bulge boundary at 23.5 h magnetic local time (MLT), while MMS was located duskward of 22 h MLT during an inward crossing of the expanding plasma sheet boundary. The two spacecraft observed a consistent set of signatures as they simultaneously crossed the reconnection separatrix layer during this very intense reconnection event. These include (1) energy dispersion of the energetic ions and electrons traveling earthward, accompanied with high electron energies in the vicinity of the separatrix; (2) energy dispersion of polar rain electrons, with a high-energy cutoff; and (3) intense inward convection of the magnetic field lines at the MMS location. The high temporal resolution measurements by MMS provide unprecedented observations of the outermost electron boundary layer. We discuss the relevance of the energy dispersion of the electrons, and their pitch angle distribution, to the spatial and temporal evolution of the boundary layer. The results indicate that the underlying magnetotail magnetic reconnection process was an intrinsically impulsive and the active X-line was located relatively close to the Earth, approximately at 16-18 R_E.

Structure, force balance, and topology of Earth's magnetopause

The magnetopause deflects the solar wind plasma and confines Earth's magnetic field. We combine measurements made by the four spacecraft of the Magnetospheric Multiscale mission to demonstrate how the plasma and magnetic forces at the boundary affect the interaction between the shocked solar wind and Earth's magnetosphere. We compare these forces with the plasma pressure and examine the electron distribution function. We find that the magnetopause has sublayers with thickness comparable to the ion scale. Small pockets of low magnetic field strength, small radius of curvature, and high electric current mark the electron diffusion region. The flow of electrons, parallel and antiparallel to the magnetic field, reveals a complex topology with the creation of magnetic ropes at the boundary.

Global observations of magnetospheric high-m poloidal waves during the 22 June 2015 magnetic storm

We report global observations of high-m poloidal waves during the recovery phase of the 22 June 2015 magnetic storm from a constellation of widely spaced satellites of five missions including Magnetospheric Multiscale (MMS), Van Allen Probes, Time History of Events and Macroscale Interactions during Substorm (THEMIS), Cluster, and Geostationary Operational Environmental Satellites (GOES). The combined observations demonstrate the global spatial extent of storm time poloidal waves. MMS observations confirm high azimuthal wave numbers (m ~ 100). Mode identification indicates the waves are associated with the second harmonic of field line resonances. The wave frequencies exhibit a decreasing trend as L increases, distinguishing them from the single-frequency global poloidal modes normally observed during quiet times. Detailed examination of the instantaneous frequency reveals discrete spatial structures with step-like frequency changes along L. Each discrete L shell has a steady wave frequency and spans about 1 R_E , suggesting that there exist a discrete number of drift-bounce resonance regions across L shells during storm times.

Rippled Quasiperpendicular Shock Observed by the Magnetospheric Multiscale Spacecraft

Collisionless shock nonstationarity arising from microscale physics influences shock structure and particle acceleration mechanisms. Nonstationarity has been difficult to quantify due to the small spatial and temporal scales. We use the closely spaced (subgyroscale), high-time-resolution measurements from one rapid crossing of Earth's quasiperpendicular bow shock by the Magnetospheric Multiscale (MMS) spacecraft to compare competing nonstationarity processes. Using MMS's high-cadence kinetic plasma measurements, we show that the shock exhibits nonstationarity in the form of ripples.

A comparative study of dipolarization fronts at MMS and Cluster

We present a statistical study of dipolarization fronts (DFs), using magnetic field data from MMS and Cluster, at radial distances below 12 R_E and 20 R_E , respectively. Assuming that the DFs have a semicircular cross section and are propelled by the magnetic tension force, we used multispacecraft observations to determine the DF velocities. About three quarters of the DFs propagate earthward and about one quarter tailward. Generally, MMS is in a more dipolar magnetic field region and observes larger-amplitude DFs than Cluster. The major findings obtained in this study are as follows: (1) At MMS ∼57 % of the DFs move faster than 150 km/s, while at Cluster only ∼35 %, indicating a variable flux transport rate inside the flow-braking region. (2) Larger DF velocities correspond to higher B_z  values directly ahead of the DFs. We interpret this as a snow plow-like phenomenon, resulting from a higher magnetic flux pileup ahead of DFs with higher velocities.

Transient, small-scale field-aligned currents in the plasma sheet boundary layer during storm time substorms

We report on field-aligned current observations by the four Magnetospheric Multiscale (MMS) spacecraft near the plasma sheet boundary layer (PSBL) during two major substorms on 23 June 2015. Small-scale field-aligned currents were found embedded in fluctuating PSBL flux tubes near the separatrix region. We resolve, for the first time, short-lived earthward (downward) intense field-aligned current sheets with thicknesses of a few tens of kilometers, which are well below the ion scale, on flux tubes moving equatorward/earthward during outward plasma sheet expansion. They coincide with upward field-aligned electron beams with energies of a few hundred eV. These electrons are most likely due to acceleration associated with a reconnection jet or high-energy ion beam-produced disturbances. The observations highlight coupling of multiscale processes in PSBL as a consequence of magnetotail reconnection.

Geoeffective jets impacting the magnetopause are very common

The subsolar magnetosheath is penetrated by transient enhancements in dynamic pressure. These enhancements, also called high-speed jets, can propagate to the magnetopause, causing large-amplitude yet localized boundary indentations on impact. Possible downstream consequences of these impacts are, e.g., local magnetopause reconnection, impulsive penetration of magnetosheath plasma into the magnetosphere, inner magnetospheric and boundary surface waves, drop outs and other variations in radiation belt electron populations, ionospheric flow enhancements, and magnetic field variations observed on the ground. Consequently, jets can be geoeffective. The extend of their geoeffectiveness is influenced by the amount of mass, momentum, and energy they transport, i.e., by how large they are. Their overall importance in the framework of solar wind-magnetosphere coupling is determined by how often jets of geoeffective size hit the dayside magnetopause. In this paper, we calculate such jet impact rates for the first time. From a large data set of Time History of Events and Macroscale Interactions during Substorms (THEMIS) multispacecraft jet observations, we find distributions of scale sizes perpendicular and parallel to the direction of jet propagation. They are well modeled by an exponential function with characteristic scales of 1.34R_E (perpendicular) and 0.71R_E (parallel direction), respectively. Using the distribution of perpendicular scale sizes, we derive an impact rate of jets with cross-sectional diameters larger than 2R_E on a reference area of about 100RE2 of the subsolar magnetopause. That rate is about 3 per hour in general, and about 9 per hour under low interplanetary magnetic field cone angle conditions (<30°), which are favorable for jet occurrence in the subsolar magnetosheath.

Two states of magnetotail dipolarization fronts: A statistical study

We study the ion density and temperature in the predipolarization and postdipolarization plasma sheets in the Earth's magnetotail using 9 years (2001-2009) of Cluster data. For our study we selected cases when Cluster observed dipolarization fronts (DFs) with an earthward plasma flow greater than 150km/s. We perform a statistical study of the temperature and density variations during the DF crossings. Earlier studies concluded that on average, the temperature increases while the density decreases across the DF. Our statistical results show a more diverse picture: While ∼54% of the DFs follow this pattern (category A), for ∼28% the temperature decreases while the density increases across the DF (category B). We found an overall decrease in thermal pressure for category A DFs with a more pronounced decrease at the DF duskside, while DFs of category B showed no clear pattern in the pressure change. Both categories are associated with earthward plasma flows but with some difference: (1) category A flows are faster than category B flows, (2) the observations indicate that category B flows are directed perpendicular to the current in the near-Earth current sheet while category A flows are tilted slightly duskward from this direction, and (3) the background Bz of category B is higher than that of category A. Based on these results, we hypothesize that after reconnection takes place, a bursty bulk flow emerges with category A characteristics, and as it travels earthward, it further evolves into category B characteristics, which is in a more dipolarized region with slower plasma flow (closer to the flow-braking region).

On the generation of magnetosheath high-speed jets by bow shock ripples

[1]The terrestrial magnetosheath is embedded with coherent high-speed jets of about 1R_E in scale, predominantly during quasi-radial interplanetary magnetic field (IMF). When these high dynamic pressure (P_dyn) jets hit the magnetopause, they cause large indentations and further magnetospheric effects. The source of these jets has remained controversial. One of the proposed mechanisms is based on ripples of the quasi-parallel bow shock. In this paper, we combine for the first time, 4 years of subsolar magnetosheath observations from the Time History of Events and Macroscale Interactions during Substorms mission and corresponding NASA/OMNI solar wind conditions with model calculations of a rippled bow shock. Concentrating on the magnetosheath close to the shock during intervals when the angle between the IMF and the Sun-Earth line was small, we find that (1) 97% of the observed jets can be produced by local ripples of the shock under the observed upstream conditions; (2) the coherent jets form a significant fraction of the high P_dyn tail of the magnetosheath flow distribution; (3) the magnetosheath P_dyn distribution matches the flow from a bow shock with ripples that have a dominant amplitude to wavelength ratio of about 9% (∼0.1R_E/1R_E) and are present ∼12% of the time at any given location.