Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heliosphere”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Exploring the Outer Heliosphere Through ENA Observations from IMAP

The Interstellar Mapping and Acceleration Probe (IMAP) mission is poised to revolutionize our understanding of two of the most central issues in heliophysics today: the nature of the interaction of the solar wind with the local interstellar medium (LISM), and the acceleration of energetic particles throughout the heliosphere. Here we focus on the anticipated contributions to this endeavor by measurement of hydrogen energetic neutral atoms (ENAs) that emanate from the solar wind–LISM interface region. ENA imaging allows us to remotely sense the structure and plasma characteristics of the heliosheath and the LISM just beyond the heliopause, and how these regions evolve in response to solar activity. These instruments, IMAP-Lo, IMAP-Hi, and IMAP-Ultra, cover energies from 10 eV to 300 keV with unprecedented angular resolution, energy resolution, and detection sensitivity. In this paper, we present the IMAP mission science specific to the study of the heliosphere/LISM interaction. We review what we have learned about the heliosphere interaction based on prior observations and theoretical studies. We then present the outstanding science questions that have consequently arisen and how IMAP plans to resolve them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spectral Properties of Globally Distributed ENA Fluxes across Diverse Regions of the Heliosphere

This study analyzes energetic neutral atom (ENA) spectral properties across distinct regions of globally distributed flux (GDF) sky maps, using Interstellar Boundary Explorer data from a full solar cycle, corrected for time dispersion. By time-shifting the data to the heliosheath using GDF source distances from D. B. Reisenfeld et al., we achieve a more accurate representation of heliosheath GDF energy spectra. We quantify ENA spectral characteristics, heliosheath line-of-sight-integrated proton pressure, and heliosheath proton temperature, comparing these to solar wind properties at 1 au and interplanetary scintillation-derived solar wind data. Our findings show that the spectral index is generally anticorrelated with heliosheath proton temperature and pressure, except in the central tail, where a partial positive correlation is observed. The lowest spectral index values occur when high-latitude heliosheath regions are dominated by fast solar wind from polar coronal holes. The south pole exhibits the flattest energy spectra due to plasma heating from both fast solar wind and a late-2014 pressure pulse. The central tail shows shorter variability (5–6 yr) for spectral index and heliosheath proton temperature, while proton pressure follows the 11 yr solar cycle. Most spectral shapes exhibit a “knee” distribution, peaking during solar maximum, with an “ankle” shape observed only at the south pole during solar cycle transitions. Asymmetry in proton pressure in the lobes is driven by the draping effect of the local interstellar magnetic field. This study provides insights into the energetic properties of GDF across the heliosphere, enhancing our understanding of the heliospheric environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Occurrence of magnetic reconnection in the heliospheric current sheet

We have analyzed the solar wind properties associated with a comparable number of the heliospheric current sheet (HCS) crossing events by Parker Solar Probe (PSP) ranging from 0.07 to 0.3 au and Advanced Composition Explorer (ACE) at 0.99 au. Nearly all PSP events (7 out of 8) show signatures of magnetic reconnection, which are more frequent than the ACE events (5 out of 8) that show reconnection. Because the HCS reconnection events have occurred in a variety of wind speeds and plasma conditions, for each event, we propose to define an approximate aspect ratio (width/length) of HCS as the ratio between the absolute HCS width (derived from observations) and the distance traveled by Alfvén waves over the propagation time of the solar wind. We find that the aspect ratio defined in such a way tends to be smaller than 0.01 for most reconnecting events, and becomes much larger than 0.01 for non-reconnecting events. This analysis also explains the different occurrence rates of reconnection observed by PSP and ACE. Potential consequences of magnetic reconnection at the HCS are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Global heliospheric termination shock strength in the solar–interstellar interaction

A heliospheric termination shock (HTS) surrounds our Solar System at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before encountering the interstellar medium. HTS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the global properties of the HTS beyond in situ measurements from Voyager in only two directions of the sky. Here we fill this gap by extracting the HTS strength using particle-in-cell, test particle and magnetohydrodynamic simulations, constrained by Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS-accelerated particles. Our results reveal there is a higher compression near the poles during solar minimum compared with solar maximum due to the higher Mach number flow. North–south asymmetries arise from the disparate evolution of the polar coronal holes, while minimum compression near the flanks is probably due to SW slowing from mass loading over a greater distance to the HTS. The results imply a strong connection between the HTS strength and the SW and interstellar medium dynamics.

Astrophysical plasmas↗
DOE OSTI2026

Energy Evolution from the Chromosphere to the Heliosphere in the 2021 October 28 Solar Eruption

We perform a detailed study of the energetics for a well-observed solar eruption and flare that occurred on 2021 October 28. This event included a GOES class X1.0 flare, a global extreme-UV (EUV) wave, and a coronal mass ejection (CME) that reached speeds of >2000 km s −1 . The event was observed from a variety of spacecraft in NASA’s Heliophysics System Observatory, including multiple missions near Earth, STEREO-A off the Sun-Earth line, and Solar Orbiter, near the Sun-Earth line at about 0.8 au. Using remote sensing, in situ observations, and in some cases scaling laws based on previous observations, we characterize the following quantities: free magnetic energy, energy in nonthermal electrons, energy in nonthermal ions, bolometric energy, energy deposited in the chromosphere, thermal energy radiated in the flare loops, energy dissipated by the EUV wave, CME kinetic and gravitational potential energy, CME energy flux in the heliosphere, and the energy partition in the CME shock. We find that the total energy released during the event is consistent with estimates of the pre-event stored magnetic energy, and the CME kinetic + potential energy dominates the energy partition.

solar coronal mass ejection shocks↗
DOE OSTI2026

The Role of Magnetic Reconnection in Energizing Protons and Heavier Ions at the Heliospheric Current Sheet

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions, indicating possible energization by magnetic reconnection up to 10 s–100 s keV nucleon −1 . Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff E max consistent with in situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as E max ∝ (Q/M) α with α ∼ 0.8–1.1, approximately consistent with PSP measurements in the broader range α ∼ 0.6–1.7. In the limit where ions are injected at the same energy per nucleon, α can be as low as ∼0.3. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

Murtas, Giulia [West Virginia Univ., Morgantown, W↗
DOE OSTI2026

Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

Shadows cast in the cosmic ray (CR) muon sky by the Sun were located using muon data from the MINOS far detector in Northern Minnesota. The shadows were observed independently across three time periods; near solar minimum, near solar maximum, and over the entire 13 year span of the data. A distribution of muon positions for each shadow was then sampled to simulate CR motions near the Sun using the Parker spiral model of the Heliospheric Magnetic Field (HMF) and a dipole model of the Geomagnetic Field (GMF). The resulting particle distributions were then compared to their position with respect to the Sun. Results show that the Parker spiral model is most consistent with the solar minimum shadow and least consistent with the solar maximum shadow, as expected. The simple Parker spiral is more consistent with the data for a harder CR spectrum than is actually present, indicating the need for a more detailed HMF model. Plausible modifications to the Parker spiral model which would affect the overall shift of the Sun's CR shadow are discussed.

Adamson, P. [Fermilab]↗
DOE OSTI2025

Simultaneous Observation of Ion-scale Wave Packets with Opposite Polarizations and Their Implications on the Generation Region in the Inner Heliosphere

This paper reports a dispersion analysis of two wave packets simultaneously observed near the local proton gyrofrequency by the Parker Solar Probe. The observed wave event exhibits clear two-banded wave packets both propagating along the magnetic field, characterized by left-handed (L-mode) and right-handed (R-mode) polarizations simultaneously. By incorporating the Doppler shift effect into a linear dispersion analysis, we find two possible scenarios that explain these simultaneous opposite polarizations: (1) Two inherently L-mode waves in the plasma frame, propagate parallel and antiparallel to the solar wind velocity, with similar wave frequencies and wave numbers. The polarization of the antiparallel propagating wave reverses as it moves sunward in the plasma frame while still comoving with the solar wind in the stationary frame. This reversal manifests the polarization of the wave as an R-mode in the spacecraft frame. (2) Simultaneous L-mode and R-mode waves propagate parallel to the solar wind velocity, with different wave frequencies and wave numbers. Concurrent proton observations during the wave event reveal a dominant anisotropic ($T$⟂/$T$ ∥ > 1) core distribution with a drifting beam population. Estimation of the linear growth rate for both L-mode and R-mode waves suggests that both scenarios are plausible, indicating that the observation is near the wave-generation region. We explore the potential impact of these simultaneous waves on solar wind heating and scattering effects, hypothesizing that such waves might enhance efficiency compared to waves with a single wave packet, contingent upon the statistical significance of such waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Interstellar Mapping And Acceleration Probe High Energy (IMAP-Hi) Neutral Atom Imager

The IMAP-Hi Energetic Neutral Atom (ENA) Imager on NASA’s Interstellar Mapping and Acceleration Probe (IMAP) mission (McComas et al. 2018a, 2025) is designed to measure ENAs from the global interaction between the heliosphere and the local interstellar medium (LISM). These ENAs are initially plasma ions of solar wind origin that are neutralized by charge exchange with the cold neutral atoms of LISM that freely flow through the heliosphere-LISM interaction region. IMAP-Hi consists of two identical single-pixel sensors, each covering the ENA spectral range from 0.44 keV to 15.6 keV over nine contiguous energy passbands and having an approximately conical field-of-view (FOV) of 4.1o full width at half maximum (FWHM). The Hi-45 sensor points 45o relative to the spacecraft spin axis from the antisunward direction; each spacecraft spin, it measures ENA intensity over a circular swath with half-cone angle 45o centered on the ecliptic plane. The Hi-90 sensor points 90o relative to the spin axis; each spacecraft spin, it measures ENA intensity over a great circle in the sky, sampling both the north and south ecliptic poles. As the IMAP spin vector is re-pointed daily toward the Sun, the ecliptic longitude of the swaths moves daily by ∼1o such that a full sky map is acquired by Hi-90 every six months and a complete low latitude (−45o to +45o) map is acquired by Hi-45 annually. The IMAP-Hi sensor design has direct heritage from the IBEX-Hi imager on the Interstellar Boundary Explorer (IBEX) mission, with substantial improvements in energy range, energy resolution, angular resolution, signal-to-noise ratio, and, for ecliptic latitudes within ±45o, temporal resolution and exposure time. The global ENA maps acquired by IMAP-Hi partially overlap in energy and viewing with the ENA maps acquired by the IMAP-Lo and IMAP-Ultra ENA imagers, which we combine to answer fundamental questions about the structure and dynamics of the interaction of the heliosphere and the LISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling Particle Acceleration and Release from Solar Eruptions

Determining the relative contribution of solar flares versus coronal mass ejections in large solar energetic particle (SEP) events is a long-standing problem. Flare-accelerated particles may travel through complex magnetic fields in the eruption region and escape into interplanetary space, thereby contributing to large SEP events. The process by which flare accelerated particles are released into the heliosphere is poorly understood and yet is critical to advancing our understanding of SEPs. In this work, we address the release problem by solving the focused transport equation in the context of a 2.5D ARMS magnetohydrodynamic simulation of a breakout coronal mass ejection (CME)/flare event. We find that particles accelerated by flare reconnection can be released into interplanetary space through interchange reconnection between closed and open field lines. These particles can contribute directly to SEP events and may become an important seed population for further acceleration by CME-driven shocks. Additionally, we find that the energetic particle fluxes in the inner heliosphere remain elevated for an extended period, allowing them to contribute to SEP acceleration by subsequent CMEs. This study represents the first direct particle modeling of how flare-accelerated particles can contribute to major SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Puzzling Variation of Gamma Rays from the Sun over the Solar Cycle Revealed with Fermi-LAT

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic disk component, and (2) a leptonic extended component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11 yr solar cycle, being highest during solar minimum and lowest during solar maximum, as it varies with the cosmic-ray flux. No study has yet analyzed the flux variation of each component over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 yr of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near Earth. We find that the flux variation of the disk anticorrelates with the sunspot number and correlates with cosmic-ray protons, as expected, confirming its emission mechanism. In contrast, the extended component exhibits a more complex variation: despite an initial anticorrelation with the sunspot number, we find neither anticorrelation with the sunspot number nor correlation with cosmic-ray electrons over the full 15 yr period. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and may differ for electrons and protons or, alternatively, that there is an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.

Acharyya, A. [Univ. of Southern Denmark, Odense (D↗
DOE OSTI2025

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗
DOE OSTI2025

On Alfvénic turbulence of solar wind streams observed by Solar Orbiter during March 2022 perihelion and their source regions

It has been recently accepted that the standard classification of the solar wind solely according to flow speed is outdated, and particular interest has been devoted to the study of the origin and evolution of so-called Alfvénic slow solar wind streams and to what extent such streams resemble or differ from fast wind. In March 2022, Solar Orbiter completed its first nominal phase perihelion passage. During this interval, it observed several Alfvénic streams, allowing for characterization of fluctuations in three slow wind intervals (AS1-AS3) and comparison with a fast wind stream (F) at almost the same heliocentric distance. This work makes use of Solar Orbiter plasma parameters from the Solar Wind Analyzer (SWA) and magnetic field measurements from the magnetometer (MAG). The magnetic connectivity to the solar sources of selected solar wind intervals was reconstructed using a ballistic extrapolation based on measured solar wind speed down to the (spherical) source surface at 2.5 R s below which a potential field extrapolation was used to map back to the Sun. The source regions were identified using SDO/AIA observations. A spectral analysis of in situ measured magnetic field and velocity fluctuations was performed to characterize correlations, Alfvénicity, normalized cross-helicity, and residual energy in the frequency domain as well as intermittency of the fluctuations and spectral energy transfer rate estimated via mixed third-order moments. A machine learning technique was used to separate proton core, proton beam, and alpha particles and to study v − b correlations for the different ion populations in order to evaluate the role played by each population in determining the Alfvénic content of solar wind fluctuations. The comparison between fast wind and Alfvénic slow wind intervals highlights the differences between the two solar wind regimes: The fast wind is characterized by larger amplitude fluctuations, and magnetic and velocity fluctuations are closer to equipartition of energy. In fact the Alfvénic slow wind streams appear to be on a spectrum of wind types, with AS1, originating from open field lines neighboring active regions and displaying similarities with the fast wind in terms of fluctuation amplitude and turbulence characteristics, but not with respect to the alpha particles and proton beams. The other two slow streams differed both in their sources as well as plasma characteristics, with AS2 coming from the expansion of a narrow coronal hole corridor and AS3 from a region straddling a pseudostreamer. The latter displayed the coldest and highest density but the slowest stream with the smallest fluctuation amplitude and greatest magnetic energy excess. It also showed the largest scatter in proton beam speeds and the greatest difference in speed between proton beam and alpha particles. This study shows how the old fast–slow solar wind dichotomy, already called into question by the observations of slower Alfvénic solar wind streams, should further be refined, as the Alfvénic slow wind, originating in different solar wind regions, show significant differences in density, temperature, and proton and alpha-particle properties in the inner heliosphere. The observations presented here provide the starting point for a better understanding of the origin and evolution of different solar wind streams as well as the evolving turbulence contained within.

magnetohydrodynamics (MHD)↗
DOE OSTI2026

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Scaler Rates from the Pierre Auger Observatory: A New Proxy of Solar Activity

The modulation of low-energy galactic cosmic rays reflects interplanetary magnetic field variations and can provide useful information on solar activity. An array of ground-surface detectors can reveal the secondary particles, which originate from the interaction of cosmic rays with the atmosphere. In this work, we present an investigation of the low-threshold rate (scaler) time series recorded in 16 yr of operation by the Pierre Auger Observatory surface detectors in Malargüe, Argentina. Through an advanced spectral analysis, we detected highly statistically significant variations in the time series with periods ranging from the decadal to the daily scale. We investigate their origin, revealing a direct connection with solar variability. Thanks to their intrinsic very low noise level, the Auger scalers allow a thorough and detailed investigation of the galactic cosmic-ray flux variations in the heliosphere at different timescales and can, therefore, be considered a new proxy of solar variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Electron-scale Magnetic Holes Generation Driven by Whistler-to-Bernstein Mode Conversion in Fully Kinetic Plasma Turbulence

Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗