Search

Traditionally, data analysis and theory have been viewed as separate disciplines, each feeding into fundamentally different types of models. Modern deep learning technology is beginning to unify these two disciplines and will produce a new class of predictively powerful space weather models that combine the physical insights gained by data and theory. We call on NASA to invest in the research and infrastructure necessary for the heliophysics' community to take advantage of these advances., Comment: Heliophysics 2050 White Paper

We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms -- the ``combinatorial'' form related to the logarithm of the number of microstates per macrostate and the ``continuous'' form related to $f \ln f$, where $f$ is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of anti-parallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a position space entropy and a velocity space entropy, and we use this to investigate the nature of kinetic entropy transport during collisionless reconnection. We find the velocity space entropy of both electrons and ions increases in time due to plasma heating during magnetic reconnection, while the position space entropy decreases due to plasma compression. This project uses collisionless simulations, so it cannot address physical dissipation mechanisms; nonetheless, the infrastructure developed here should be useful for studies of collisional or weakly collisional heliospheric, planetary, and astrophysical systems. Beyond reconnection, the diagnostic is expected to be applicable to plasma turbulence and collisionless shocks., Comment: 20 pages, 7 figures

We compare a global high resolution resistive magnetohydrodynamics (MHD) simulation of Earth's magnetosphere with observations from the Magnetospheric Multiscale (MMS) constellation for a southward IMF magnetopause crossing during October 16, 2015 that was previously identified as an electron diffusion region (EDR) event. The simulation predicts a complex time-dependent magnetic topology consisting of multiple separators and flux ropes. Despite the topological complexity, the predicted distance between MMS and the primary separator is less than 0.5 Earth radii. These results suggest that global magnetic topology, rather than local magnetic geometry alone, determines the location of the electron diffusion region at the dayside magnetopause., Comment: 13 pages, 4 figures; Submitted to Journal of Geophysical Research

We present high resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall currents in the reconnection plane accelerate ions to the local Alfv\'en speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian "double magnetopause" structure, an Alfv\'enic ion jet extending across the upstream magnetopause and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves., Comment: 14 pages, 12 figures; presented at Geospace Environment Modeling (GEM) workshop (June, 2014) and Fall American Geophysical Union (AGU) meeting (December, 2014); submitted to Journal of Geophysical Research, December 2014

We demonstrate that the reconnection rate at the subsolar magnetopause is stronglycontrolled by the solar wind electric field and depends weakly on the local properties of the dissipationregion. Our approach is to match the solar wind and magnetospheric states in an internal boundarylayer described by the Cassak and Shay (2007, https://doi.org/10.1063/1.2795630) expression fortwo-dimensional asymmetric reconnection. Faraday's law along the Sun-Earth line determinesthe variation of the solar wind electric field from the bow shock to the magnetopause. While themagnetospheric plasma exerts some control over the reconnection rate, magnetic flux pileup in the sheathpartially compensates for any local reduction in the reconnection rate. For a fixed magnetospheric state,the reconnection rate is shown to be directly proportional to the solar wind electric field, thus explainingwhy the solar wind electric field correlates well with geomagnetic indices.

On 27 June 1996, the NASA Galileo spacecraft made humanity’s first flyby of Jupiter’s largest moon, Ganymede, discovering that it is the only moon known to possess an internally generated magnetic field. Resurrecting the original Galileo Plasma Subsystem (PLS) data analysis software, we processed the raw PLS data from G01 and for the first time present the properties of plasmas encountered. Entry into the magnetosphere of Ganymede occurred near the confluence of the magnetopause and plasma sheet. Reconnection-driven plasma flows were observed (consistent with an Earth-like Dungey cycle), which may be a result of reconnection in the plasma sheet, magnetopause, or might be Ganymede’s equivalent of a Low-Latitude Boundary Layer. Dropouts in plasma density combined with velocity perturbations afterward suggest that Galileo briefly crossed the cusps into closed magnetic field lines. Galileo then crossed the cusps, where field-aligned precipitating ions were observed flowing down into the surface, at a location consistent with observations by the Hubble Space Telescope. The density of plasma outflowing from Ganymede jumped an order of magnitude around closest approach over the north polar cap. The abrupt increase may be a result of crossing the cusp or may represent an altitude-dependent boundary such as an ionopause. More diffuse, warmer field-aligned outflows were observed in the lobes. Fluxes of particles near the moon on the nightside were significantly lower than on the dayside, possibly resulting from a diurnal cycle of the ionosphere and/or neutral atmosphere.

Secondary electrons are continuously generated via photoemission from sunlit spacecraft and instrument surfaces. These particles can subsequently contaminate low-energy channels of electron sensors. Spacecraft photoelectrons are measured at energies below that of a positive spacecraft potential and can be removed at the expense of energy resolution. However, fluxes of photoelectrons generated inside electron instruments are independent of spacecraft potential and must be fully characterized in order to correct electron data. Here we present observations of spacecraft and instrument photoelectron populations measured with the Dual Electron Spectrometers (DES) on NASA's Magnetospheric Multiscale (MMS) mission. We leverage observations from Earth's nightside plasma sheet taken during MMS commissioning and develop an empirical model of instrument photoelectrons. This model is used with DES velocity distribution functions to correct plasma moments and has been made publicly available on the MMS science data center for use by the scientific community.

Ion velocity distribution functions (VDFs) from a particle-in-cell simulation of asymmetric reconnection are investigated to reveal a two-scale structure of the ion diffusion region (IDR). Ions bouncing in the inner IDR are trapped mainly by the electric field normal to the current sheet (N direction), while those reaching the outer IDR are turned back mainly by the magnetic force. The resulting inner layer VDFs have counter-streaming populations along N with decreasing counter-streaming speeds away from the midplane while maintaining the out-of-plane speed, and the outer layer VDFs exhibit crescent shapes toward the out-of-plane direction. Observations of the above VDF features and the normal electric fields provide evidence for the two-scale meandering motion.

Although the effects of magnetic reconnection in magnetospheres can be observed at planetary scales, reconnection is initiated at electron scales in a plasma. Surrounding the electron diffusion region, there is an Ion-Decoupling Region (IDR) of the size of the ion length scales (inertial length and gyroradius). Reconnection at the Earths magnetopause often includes cold magnetospheric (few tens of eV), hot magnetospheric (10 keV), and magnetosheath (1 keV) ions, with different gyroradius length scales. We report observations of a subregion inside the IDR of the size of the cold ion population gyroradius (approx. 15 km) where the cold ions are demagnetized and accelerated parallel to the Hall electric field. Outside the subregion, cold ions follow the E x B motion together with electrons, while hot ions are demagnetized. We observe a sharp cold ion density gradient separating the two regions, which we identify as the cold and hot IDRs.

Measurements from the Magnetospheric Multiscale (MMS) mission are reported to show distinct features of electron energization and mixing in the diffusion region of the terrestrial magnetopause reconnection. At the ion jet and magnetic field reversals, distribution functions exhibiting signatures of accelerated meandering electrons are observed at an electron out-of-plane flow peak. The meandering signatures manifested as triangular and crescent structures are established features of the electron diffusion region (EDR). Effects of meandering electrons on the electric field normal to the reconnection layer are detected. Parallel acceleration and mixing of the inflowing electrons with exhaust electrons shape the exhaust flow pattern. In the EDR vicinity, the measured distribution functions indicate that locally, the electron energization and mixing physics is captured by two-dimensional reconnection, yet to account for the simultaneous four-point measurements, translational invariant in the third dimension must be violated on the ion-skin-depth scale.

The strength and time dependence of the electric field in a magnetopause diffusion region relate to the rate of magnetic reconnection between the solar wind and a planetary magnetic field. Here we use approximately 150 milliseconds measurements of energetic electrons from the Mercury Surface, Space Environment, GEochemistry, and Ranging (MESSENGER) spacecraft observed over Mercury's dayside polar cap boundary (PCB) to infer such small-scale changes in magnetic topology and reconnection rates. We provide the first direct measurement of open magnetic topology in flux transfer events at Mercury, structures thought to account for a significant portion of the open magnetic flux transport throughout the magnetosphere. In addition, variations in PCB latitude likely correspond to intermittent bursts of approximately 0.3 to 3 millivolts per meter reconnection electric fields separated by approximately 5 to10 seconds, resulting in average and peak normalized dayside reconnection rates of approximately 0.02 and approximately 0.2, respectively. These data demonstrate that structure in the magnetopause diffusion region at Mercury occurs at the smallest ion scales relevant to reconnection physics.

Each spacecraft of the recently launched magnetospheric multiscale MMS mission is equipped with Active Spacecraft Potential Control (ASPOC) Instruments, which control the spacecraft potential in order to reduce spacecraft charging effects. ASPOC typically reduces the spacecraft potential to a few volts. On several occasions during the commissioning phase of the mission, the ASPOC instruments were operating only on one spacecraft at a time. Taking advantage of such intervals, we derive photoelectron curves and also perform reconstructions of the uncontrolled spacecraft potential for the spacecraft with active control and estimate the electron plasma density during those periods. We also establish the criteria under which our methods can be applied.

We report ion velocity distribution functions (VDfs) observed by Magnetospheric Multiscale Mission (MMS) and present evidence for demagnetized ion Speiser motion during magnetopause reconnection. The demagnetization is observed in the vicinity of the X llne, as well as near the current sheet midlplane about tens of ion skin depths (d(sub 1)) away from the X line. Close to the X line before the outflow is built up, the VDFs are elongated, and the elongated part of VDFs rotates from the out-of-plane current direction toward the outflow directions downstream from the X line. Farther downstream, demagnetized ions exhibit a characteristic half-ring structure in the VDFs, as a result of the mixture of ions that have experienced different amounts of cyclotron turning around the magnetic field normal to the current sheet. Signatures of acceleration by electric fields are more pronounced in the VDFs near the X line than downstream.

Magnetic holes are ubiquitous in space plasmas, occurring in the solar wind, downstream of planetary bow shocks, and inside the magnetosphere. Recently, kinetic-scale magnetic holes have been observed near Earth's central plasma sheet. The Fast Plasma Investigation on NASA's Magnetospheric Multiscale (MMS) mission enables measurement of both ions and electrons with 2 orders of magnitude increased temporal resolution over previous magnetospheric instruments. Here we present data from MMS taken in Earth's nightside plasma sheet and use high-resolution particle and magnetometer data to characterize the structure of a subproton-scale magnetic hole. Electrons with gyroradii above the thermal gyroradius but below the current layer thickness carry a current sufficient to account for a 10-20 depression in magnetic field magnitude. These observations suggest that the size and magnetic depth of kinetic-scale magnetic holes is strongly dependent on the background plasma conditions.

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is insufficient with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.

The most common instrument for low energy plasmas consists of a top-hat electrostatic analyzer geometry coupled with a microchannel-plate (MCP)-based detection system. While the electrostatic optics for such sensors are readily simulated and parameterized during the laboratory calibration process, the detection system is often less well characterized. Furthermore, due to finite resources, for large sensor suites such as the Fast Plasma Investigation (FPI) on NASA's Magnetospheric Multiscale (MMS) mission, calibration data are increasingly sparse. Measurements must be interpolated and extrapolated to understand instrument behavior for untestable operating modes and yet sensor inter-calibration is critical to mission success. To characterize instruments from a minimal set of parameters we have developed the first comprehensive mathematical description of both sensor electrostatic optics and particle detection systems. We include effects of MCP efficiency, gain, scattering, capacitive crosstalk, and charge cloud spreading at the detector output. Our parameterization enables the interpolation and extrapolation of instrument response to all relevant particle energies, detector high voltage settings, and polar angles from a small set of calibration data. We apply this model to the 32 sensor heads in the Dual Electron Sensor (DES) and 32 sensor heads in the Dual Ion Sensor (DIS) instruments on the 4 MMS observatories and use least squares fitting of calibration data to extract all key instrument parameters. Parameters that will evolve in flight, namely MCP gain, will be determined daily through application of this model to specifically tailored in-flight calibration activities, providing a robust characterization of sensor suite performance throughout mission lifetime. Beyond FPI, our model provides a valuable framework for the simulation and evaluation of future detection system designs and can be used to maximize instrument understanding with minimal calibration resources.

Plasma measurements in space are becoming increasingly faster, higher resolution, and distributed over multiple instruments. As raw data generation rates can exceed available data transfer bandwidth, data compression is becoming a critical design component. Data compression has been a staple of imaging instruments for years, but only recently have plasma measurement designers become interested in high performance data compression. Missions will often use a simple lossless compression technique yielding compression ratios of approximately 2:1, however future missions may require compression ratios upwards of 10:1. This study aims to explore how a Discrete Wavelet Transform combined with a Bit Plane Encoder (DWT/BPE), implemented via a CCSDS standard, can be used effectively to compress count information common to plasma measurements to high compression ratios while maintaining little or no compression error. The compression ASIC used for the Fast Plasma Investigation (FPI) on board the Magnetospheric Multiscale mission (MMS) is used for this study. Plasma count data from multiple sources is examined: resampled data from previous missions, randomly generated data from distribution functions, and simulations of expected regimes. These are run through the compression routines with various parameters to yield the greatest possible compression ratio while maintaining little or no error, the latter indicates that fully lossless compression is obtained. Finally, recommendations are made for future missions as to what can be achieved when compressing plasma count data and how best to do so.

The MOST IDS team was tasked with focusing on two general areas: The first was to participate with the Fast Plasma Investigation (FPI) team in the development of virtual detectors that model the instrument responses of the MMS FPI sensors. The virtual instruments can be 'flown through' both simulation data (from magnetohydrodynamic, hybrid, and kinetic simulations) and Cluster and THEMIS spacecraft data. The goal is to determine signatures of magnetic reconnection expected during the MMS mission. Such signatures can serve as triggers for selection of burst mode downloads. The chapter contributed by the FPI team covers that effort in detail and, therefore, most of that work has not been included here. The second area of emphasis, and the one detailed in this chapter, was to build on past and present knowledge of magnetic reconnection and its physical signatures. Below we describe intensive analyses of Cluster and THEMIS data together with theoretical models and simulations that delineate the plasma signatures that surround sites of reconnection, including the effects of turbulence as well as the detailed kinetic signatures that indicate proximity to reconnection sites. In particular, we point out that particles are energized in several regions, not only at the actual site of reconnection.

We describe our experience using NVIDIA's CUDA (Compute Unified Device Architecture) C programming environment to implement a two-dimensional second-order MUSCL-Hancock ideal magnetohydrodynamics (MHD) solver on a GTX 480 Graphics Processing Unit (GPU). Taking a simple approach in which the MHD variables are stored exclusively in the global memory of the GTX 480 and accessed in a cache-friendly manner (without further optimizing memory access by, for example, staging data in the GPU's faster shared memory), we achieved a maximum speed-up of approx. = 126 for a sq 1024 grid relative to the sequential C code running on a single Intel Nehalem (2.8 GHz) core. This speedup is consistent with simple estimates based on the known floating point performance, memory throughput and parallel processing capacity of the GTX 480.

We report our findings comparing the geometric factor (GF) as determined from simulations and laboratory measurements of the new Dual Electron Spectrometer (DES) being developed at NASA Goddard Space Flight Center as part of the Fast Plasma Investigation on NASA's Magnetospheric Multiscale mission. Particle simulations are increasingly playing an essential role in the design and calibration of electrostatic analyzers, facilitating the identification and mitigation of the many sources of systematic error present in laboratory calibration. While equations for laboratory measurement of the Geometric Factpr (GF) have been described in the literature, these are not directly applicable to simulation since the two are carried out under substantially different assumptions and conditions, making direct comparison very challenging. Starting from first principles, we derive generalized expressions for the determination of the GF in simulation and laboratory, and discuss how we have estimated errors in both cases. Finally, we apply these equations to the new DES instrument and show that the results agree within errors. Thus we show that the techniques presented here will produce consistent results between laboratory and simulation, and present the first description of the performance of the new DES instrument in the literature.

We present our experience porting a Hall MHD code to a Graphics Processing Unit (GPU). The code is a 2nd order accurate MUSCL-Hancock scheme which makes use of an HLL Riemann solver to compute numerical fluxes and second-order finite differences to compute the Hall contribution to the electric field. The divergence of the magnetic field is controlled with Dedner?s hyperbolic divergence cleaning method. Preliminary benchmark tests indicate a speedup (relative to a single Nehalem core) of 58x for a double precision calculation. We discuss scaling issues which arise when distributing work across multiple GPUs in a CPU-GPU cluster.

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Recent global MHD simulations of Flux Transfer Events (FTE's)[Dorelli and Bhattachar ee, JGR, 114, 2009] demonstrate that they are topologically complex flux ropes which extend large distances away from the subsolar magnetopause. Thus, FTE's represent a significant perturbation to the magnetopause magnetic field topology. Specifically, prior to the formation of the first FTE of the simulation, the dayside magnetopause has a relatively simple topology consisting of a single magnetic separator draping northward (even for southward IMF) over the dayside magnetopause (in a manner similar to the "overdraping" which occurs under northward IMF conditions). After the first FTE forms, however, the magnetopause becomes topologically more complex, consisting of multiple separators and braided flux domains. One consequence of this topological complexity is that open flux tubes have direct access to the magnetosphere near the subsolar region (such direct access is not possible during steady separator reconnection). Since the magnetopause topology affects the mapping of the solar wind electric field to the magnetosphere, an interesting question arises: What is the effect of FTE generation on the global magnetospheric convection pattern? In this talk, we address this question by exploring the differences between magnetospheric and ionospheric convection before and after the formation of the first FTE of the simulation.