SHIELDS project:Space Hazards Induced near Earth by Large, Dynamic Storms IC report
The main objective of SHIELDS project is to develop a new capability to understand, model, and predict the spacecraft Surface Charging Environment (SCE).
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The main objective of SHIELDS project is to develop a new capability to understand, model, and predict the spacecraft Surface Charging Environment (SCE).
We use a global hybrid code kinetic model to demonstrate how kinetic processes at the bow shock and within the foreshock can dramatically modify the solar wind just before its interaction with the magnetosphere. During periods of steady radial interplanetary magnetic field (IMF) orientation, the foreshock fills with a diffuse population of suprathermal ions. The ions generate cavities marked by enhanced temperatures, depressed densities, and diminished magnetic field strengths that convect antisunward into the bow shock with the solar wind flow. Tangential discontinuities marked by inward-pointing electric fields and normals transverse to the Sun-Earth line generate hot flow anomalies marked by hot tenuous plasmas bounded by outward propagating shocks. When the motional electric field in the magnetosheath points inward towards the Earth, a solitary bow shock appears. For typical IMF orientations, the solitary shocks should appear at poorly sampled high latitudes, but for strongly northward or southward IMF orientations the solitary shocks should appear on the flanks of the magnetosphere. Although quasi-perpendicular, solitary shocks should be marked by turbulent magnetosheath flows, often directed towards the Sun-Earth line, and abrupt spike-like enhancements in the density and magnetic field strength at the shock. Finally,we show how flux transfer events generated between parallel subsolar reconnection lines are destroyed upon encountering the magnetopause at latitudes above the cusp.
The injection region's formation, scale size, and propagation direction have been debated throughout the years, with new questions arising with increased plasma sheet observations by missions like Cluster and THEMIS. How do temporally and spatially smallscale injections relate to the larger injections historically observed at geosynchronous orbit? How to account for opposing propagation directions—earthward, tailward, and azimuthal—observed by different studies? To address these questions, we used a combination of multisatellite and ground–based observations to knit together a cohesive story explaining injection formation, propagation, and differing spatial scales and timescales. We used a case study to put statistics into context. First, fast earthward flows with embedded small–scale dipolarizing flux bundles transport both magnetic flux and energetic particles earthward, resulting in minutes–long injection signatures. Next, a large–scale injection propagates azimuthally and poleward/tailward, observed in situ as enhanced flux and on the ground in the riometer signal. The large–scale dipolarization propagates in a similar direction and speed as the large–scale electron injection. We suggest small–scale injections result from earthward–propagating, small–scale dipolarizing flux bundles, which rapidly contribute to the large–scale dipolarization. We suggest the large–scale dipolarization is the source of the large–scale electron injection region, such that as dipolarization expands, so does the injection. Here, the >90–keV ion flux increased and decreased with the plasma flow, which died at the satellites as global dipolarization engulfed them. We suggest the ion injection region at these energies in the plasma sheet is better organized by the plasma flow.
Three-dimensional particle-in-cell simulations of decaying whistler turbulence are carried out on a collisionless, homogeneous, magnetized, electron-ion plasma model. In addition, the simulations use an initial ensemble of relatively long wavelength whistler modes with a broad range of initial propagation directions with an initial electron beta β e = 0.05. The computations follow the temporal evolution of the fluctuations as they cascade into broadband turbulent spectra at shorter wavelengths. Three simulations correspond to successively larger simulation boxes and successively longer wavelengths of the initial fluctuations. The computations confirm previous results showing electron heating is preferentially parallel to the background magnetic field B o , and ion heating is preferentially perpendicular to B o . The new results here are that larger simulation boxes and longer initial whistler wavelengths yield weaker overall dissipation, consistent with linear dispersion theory predictions of decreased damping, stronger ion heating, consistent with a stronger ion Landau resonance, and weaker electron heating.
Here we present a fundamentally new type of space plasma spectrometer, the wide field of view plasma spectrometer, whose field of view is >1.25π ster using fewer resources than traditional methods. The enabling component is analogous to a pinhole camera with an electrostatic energy-angle filter at the image plane. Particle energy-per-charge is selected with a tunable bias voltage applied to the filter plate relative to the pinhole aperture plate. For a given bias voltage, charged particles from different directions are focused by different angles to different locations. Particles with appropriate locations and angles can transit the filter plate and are measured using a microchannel plate detector with a position-sensitive anode. Full energy and angle coverage are obtained using a single high-voltage power supply, resulting in considerable resource savings and allowing measurements at fast timescales. Lastly, we present laboratory prototype measurements and simulations demonstrating the instrument concept and discuss optimizations of the instrument design for application to space measurements.
On 7 January 2005 (Ap=40) prompt penetration electric field perturbations of opposite polarities were observed over Thumba and Jicamarca on a few occasions during 13:45–16:30 UT. However, the electric field was found to be eastward during 14:45–15:30 UT over both Thumba and Jicamarca contrary to the general expectation wherein opposite polarities are expected at nearly antipodal points. On closer scrutiny, three important observational features are noticed during 14:10–15:15 UT. First, during 14:10–14:45 UT, despite increasing southward interplanetary magnetic field (IMF) B z condition, the already westward electric field over Thumba weakened (less westward) while the eastward electric field over Jicamarca intensified (more eastward). Second, the electric field not only became anomalously eastward over Thumba but also got intensified further during 14:45–15:00 UT similar to Jicamarca. Third, during 15:00–15:15 UT, despite IMF B z remaining steadily southward, the eastward electric field continued to intensify over Thumba but weakened over Jicamarca. It is suggested that the changes in IMF B y component under southward IMF B z condition are responsible for skewing the ionospheric equipotential patterns over the dip equator in such a way that Thumba came into the same DP2 cell as that of Jicamarca leading to anomalous electric field variations. Magnetic field measurements along the Indian and Jicamarca longitude sectors and changes in high-latitude ionospheric convection patterns provide credence to this proposition. Therefore, in conclusion, the present investigation shows that the variations in IMF B y are fundamentally important to understand the prompt penetration effects over low latitudes.
Using the total radiation belt electron content calculated from Van Allen Probe phase space density, the time-dependent and global response of the outer radiation belt during storms is statistically studied. Using phase space density reduces the impacts of adiabatic changes in the main phase, allowing a separation of adiabatic and nonadiabatic effects and revealing a clear modality and repeatable sequence of events in storm time radiation belt electron dynamics. This sequence exhibits an important first adiabatic invariant (μ)-dependent behavior in the seed (150 MeV/G), relativistic (1,000 MeV/G), and ultrarelativistic (4,000 MeV/G) populations. The outer radiation belt statistically shows an initial phase dominated by loss followed by a second phase of rapid acceleration, while the seed population shows little loss and immediate enhancement. Lastly, the time sequence of the transition to the acceleration is also strongly μ dependent and occurs at low μ first, appearing to be repeatable from storm to storm.
The ion plasma sheet (~few hundred eV to ~few 10s keV) is usually dominated by H + ions. Here, changes in ion composition within the plasma sheet are explored both during individual events, and statistically during 54 calm-to-storm events and during 21 active-to-calm events. Ion composition data from the HOPE (Helium, Oxygen, Proton, Electron) instruments onboard Van Allen Probes satellites provide exceptional spatial and temporal resolution of the H + , O + , and He + ion fluxes in the plasma sheet. H+ shown to be the dominant ion in the plasma sheet in the calm-to-storm transition. However, the energy-flux of each ion changes in a quasi-linear manner during extended calm intervals. Heavy ions (O + and He + ) become increasingly important during such periods as charge-exchange reactions result in faster loss for H + than for O + or He + . Results confirm previous investigations showing that the ion composition of the plasma sheet can be largely understood (and predicted) during calm intervals from knowledge of: (a) the composition of previously injected plasma at the onset of calm conditions, and (b) use of simple drift-physics models combined with calculations of charge-exchange losses.
Understanding the physical processes that control the dynamics of energetic particles in the inner magnetosphere is important for both space-borne and ground-based assets essential to the modern society. The storm time distribution of ring current particles in the inner magnetosphere depends strongly on their transport in the evolving electric and magnetic fields along with particle acceleration and loss. Here, we investigated the ring current particle variations using observations and simulations. We compared the ion (H + , He + , and O + ) and electron flux and plasma pressure variations from Arase observations with the self-consistent inner magnetosphere model: Ring current Atmosphere interactions Model with Self Consistent magnetic field (RAM-SCB) during the 7–8 November 2017 geomagnetic storm. We investigated the contribution of the different species (ions and electrons) to the magnetic field deformation observed at ground magnetic stations (09°–45° MLat) using RAM-SCB simulations. The results show that the ions are the major contributor with ~88% and electrons contribute ~12% to the total ring current pressure. It is also found that the electron contribution is non-negligible (~18%) to the ring current in dawn-side during the main phase of the storm. Thus, the electron contribution to the storm time ring current is important and should not be neglected.
Not Available
Electron bulk energization in the diffusion region during anti-parallel symmetric reconnection entails two stages. First, the inflowing electrons are adiabatically trapped and energized by an ambipolar parallel electric field. Next, the electrons gain energy from the reconnection electric field as they undergo meandering motion. These collisionless mechanisms have been described previously, and they lead to highly structured electron velocity distributions. Furthermore, a simplified control-volume analysis gives estimates for how the net effective heating scales with the upstream plasma conditions in agreement with fully kinetic simulations and spacecraft observations.
The collisionless reconnection of two parallel flux ropes driven by both the coalescence and kink instabilities is examined using fully kinetic simulations in periodic and line-tied geometries. The three-dimensional reconnection rate is computed from the maximum of the quasi-potential, Ξ≡−∫E·dℓ, where the integral of the electric field is taken along the magnetic field lines across the system. In periodic simulations in which the kink mode is nearly suppressed, reconnection is driven by the coalescence instability, and the peak rate is within 3%–8% of comparable 2D simulations. When a strong kink growth is observed, the peak reconnection rate drops by 10%–25%, and there is a larger drop for lower guide field. With line-tied boundary conditions, the kink instability plays a key role in allowing the flux ropes to interact and partially reconnect. In this limit, the field lines with maximum quasi-potential are associated with a quasi-separatrix layer, and the electric field along these special field lines is supported predominantly by the divergence of the electron pressure tensor. Both of these features, along with the observed reconnection rate, are consistent with recent laboratory experiments on kinetic-scale flux ropes. In kinetic simulations, the non-gyrotropic pressure tensor terms contribute significantly more to the reconnecting electric field than do the gyrotropic terms, while contributions from the electron inertia are significant for field lines adjacent to the quasi-separatrix layer.
The rate of reconnection characterizes how quickly flux and mass can move into and out of the reconnection region. In the Terrestrial Reconnection EXperiment (TREX), the rate at which the antiparallel asymmetric reconnection occurs is modulated by the presence of a shock and a region of flux pileup in the high-density inflow. Simulations utilizing a generalized Harris-sheet geometry have tentatively shown agreement with TREX's measured reconnection rate scaling relative to system size, which is indicative of the transition from ion-coupled toward electron-only reconnection. Here, in this study, we present simulations tailored to reproduce the specific TREX geometry, which confirm both the reconnection rate scale and the shock jump conditions previously characterized experimentally in TREX. The simulations also establish an interplay between the reconnection layer and the Alfvénic expansions of the background plasma associated with the energization of the TREX drive coils; this interplay has not yet been experimentally observed.
We present the results of the first fully three-dimensional particle-in-cell simulations of decaying whistler turbulence in a magnetized, homogeneous, collisionless plasma in which both forward cascades to shorter wavelengths, and inverse cascades to longer wavelengths are allowed to proceed. For the electron beta β {sub e} = 0.10 initial value considered here, the early-time rate of inverse cascade is very much smaller than the rate of forward cascade, so that at late times the fluctuation energy in the regime of the inverse cascade is much weaker than that in the forward cascade regime. Similarly, the wavevector anisotropy in the inverse cascade regime is much weaker than that in the forward cascade regime.
Plasma turbulence consists of an ensemble of enhanced, broadband electromagnetic fluctuations, typically driven by multi-wave interactions which transfer energy in wavevector space via non- linear cascade processes. In addition, temperature anisotropy instabilities in collisionless plasmas are driven by quasi-linear wave–particle interactions which transfer particle kinetic energy to field fluctuation energy; the resulting enhanced fluctuations are typically narrowband in wavevector magnitude and direction. Whatever their sources, short-wavelength fluctuations are those at which charged particle kinetic, that is, velocity-space, properties are important; these are generally wavelengths of the order of or shorter than the ion inertial length or the thermal ion gyroradius. The purpose of this review is to summarize and interpret recent computational results concerning short-wavelength plasma turbulence, short-wavelength temperature anisotropy instabilities and relationships between the two phenomena.
We determine how the amplitude and phase of each wave field component of an electromagnetic (EM) wave is split between the fast and slow roots of the cold plasma dispersion relation as it propagates across a plasma boundary.
This report highlights progress made on plasma shear turbulence in order to showcase the Institutional Computing program's recent accomplishments.
The goals and highlights of this research are to: 1) Undertake two-fluid simulations of oblique tearing modes with $m_{i}/m_{e}$ = 100 and 1836, and $b_g$ ≤ 10 with further bridging the fluid and kinetic models and linear growth rates agree with kinetic theory for $b_g$ ≤ 1; 2) 3-D simulations of a doubly-periodic Harris Sheet with $m_{i}/m_{e}$ = 100 and $b_g$ = 2:5 that features current layer splits into multiple ~ de-scale diffusion regions, an amplification of current by 2x, and similar features to kinetic simulations; 3) Influence of boundary conditions: periodic vs. line-tied; and 4) Address how much of the kinetic nature of MR can be captured by a two-fluid model. (A fluid model requires dissipation.)