Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. We measure little heating across the shock (<10% of the ion kinetic energy) which is consistent with an absence of viscous dissipation.

A Fast Bow Shock Location Predictor-Estimator From 2D and 3D Analytical Models: Application to Mars and the MAVEN Mission

We present fast algorithms to automatically estimate the statistical position of the bow shock from spacecraft data, using existing analytical two-dimensional (2D) and three-dimensional (3D) models of the shock surface. We derive expressions of the standoff distances in 2D and 3D and of the normal to the bow shock at any given point on it. Two simple bow shock detection algorithms are constructed, one solely based on a geometrical predictor from existing models, the other using this predicted position to further refine it with the help of magnetometer data, an instrument flown on many planetary missions. Both empirical techniques are applicable to any planetary environment with a defined shock structure. Applied to the Martian environment and the NASA/MAVEN mission, the predicted shock position is on average within 0.15 planetary radius R p of the bow shock crossing. Using the predictor-corrector algorithm, this estimate is further refined to within a few minutes of the true crossing (≈0.05R p). Between 2014 and 2021, we detect 14,929 clear bow shock crossings, predominantly quasi-perpendicular. Thanks to 2D conic and 3D quadratic fits, we investigate the variability of the shock surface with respect to Mars Years (MY), solar longitude (Ls), and solar EUV flux levels. Although asymmetry in Y and Z Mars Solar Orbital coordinates is on average small, we show that for MY32 and MY35, Ls = [135°-225°] and high solar flux, it can become particularly noticeable, and is superimposed to the usual North-South asymmetry due in part to the presence of crustal magnetic fields.

Shocks in the Very Local Interstellar Medium

Large-scale disturbances generated by the Sun's dynamics first propagate through the heliosphere, influence the heliosphere's outer boundaries, and then traverse and modify the very local interstellar medium (VLISM). The existence of shocks in the VLISM was initially suggested by Voyager observations of the 2-3 kHz radio emissions in the heliosphere. A couple of decades later, both Voyagers crossed the definitive edge of our heliosphere and became the first ever spacecraft to sample interstellar space. Since Voyager 1's entrance into the VLISM, it sampled electron plasma oscillation events that indirectly measure the medium's density, increasing as it moves further away from the heliopause. Some of the observed electron oscillation events in the VLISM were associated with the local heliospheric shock waves. The observed VLISM shocks were very different than heliospheric shocks. They were very weak and broad, and the usual dissipation via wave-particle interactions could not explain their structure. Estimates of the dissipation associated with the collisionality show that collisions can determine the VLISM shock structure. According to theory and models, the existence of a bow shock or wave in front of our heliosphere is still an open question as there are no direct observations yet. This paper reviews the outstanding observations recently made by the Voyager 1 and 2 spacecraft, and our current understanding of the properties of shocks/waves in the VLISM. We present some of the most exciting open questions related to the VLISM and shock waves that should be addressed in the future.

Occurrence rate of ultra-low frequency waves in the foreshock of Mercury increases with heliocentric distance

Studies of Mercury's foreshock have analyzed in detail the properties of ultra-low frequency waves. However, an open question remains in regards to understanding favorable conditions for these planetary foreshocks waves. Here, we report that 0.05-0.41 Hz quasi-monochromatic waves are mostly present under quasi-radial and relatively low intensity Interplanetary Magnetic Field, based on 17 Mercury years of MESSENGER Magnetometer data. These conditions are consistent with larger foreshock size and reflection of solar wind protons, their most likely source. Consequently, we find that the wave occurrence rate increases with Mercury's heliocentric distance. Detection of these waves throughout Mercury's highly eccentric orbit suggests the conditions for backstreaming protons are potentially present for all of Mercury's heliocentric distances, despite the relatively low solar wind Alfvén Mach number regime. These results are relevant for planetary magnetospheres throughout the solar system, and the magnetospheres of exoplanets, and provide knowledge of particle acceleration mechanisms occurring inside foreshocks.

Electron Energy Partition across Interplanetary Shocks. I. Methodology and Data Product

Analyses of 15,314 electron velocity distribution functions (VDFs) within ±2 hr of 52 interplanetary (IP) shocks observed by the Wind spacecraft near 1 au are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar VDF, while both the halo and beam/strahl components were best fit to bi-kappa VDF. This is the first statistical study to show that the core electron distribution is better fit to a self-similar VDF than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The ranges of values defined by the lower and upper quartiles for the kappa exponents are κ _ec ~ 5.40-10.2 for the core, κ _eh ~ 3.58-5.34 for the halo, and κ _eb ~ 3.40-5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents is s _ec ~ 2.00-2.04, and those of asymmetric bi-self-similar core exponents are p _ec ~ 2.20-4.00 for the parallel exponent and q _ec ~ 2.00-2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.

Direct evidence of nonstationary collisionless shocks in space plasmas

Collisionless shocks are ubiquitous throughout the universe: around stars, supernova remnants, active galactic nuclei, binary systems, comets, and planets. Key information is carried by electromagnetic emissions from particles accelerated by high Mach number collisionless shocks. These shocks are intrinsically nonstationary, and the characteristic physical scales responsible for particle acceleration remain unknown. Quantifying these scales is crucial, as it affects the fundamental process of redistributing upstream plasma kinetic energy into other degrees of freedom-particularly electron thermalization. Direct in situ measurements of nonstationary shock dynamics have not been reported. Thus, the model that best describes this process has remained unknown. Here, we present direct evidence demonstrating that the transition to nonstationarity is associated with electron-scale field structures inside the shock ramp.

Revisiting the structure of low-Mach number, low-beta, quasi-perpendicular shocks

A study of the structure of 145 low-Mach number (M ≤ 3), low-beta (β ≤ 1), quasi-perpendicular interplanetary collisionless shock waves observed by the Wind spacecraft has provided strong evidence that these shocks have large-amplitude whistler precursors. The common occurrence and large amplitudes of the precursors raise doubts about the standard assumption that such shocks can be classified as laminar structures. This directly contradicts standard models. In 113 of the 145 shocks (~78%), we observe clear evidence of magnetosonic-whistler precursor fluctuations with frequencies ~0.1-7 Hz. We find no dependence on the upstream plasma beta, or any other shock parameter, for the presence or absence of precursors. The majority (~66%) of the precursors propagate at ≤45° with respect to the upstream average magnetic field and most (~87%) propagate ≥30° from the shock normal vector. Further, most (~79%) of the waves propagate at least 20° from the coplanarity plane. The peak-to-peak wave amplitudes (δB _pk-pk) are large with a range of maximum values for the 113 precursors of ~0.4-13 nT with an average of ~2 nT. When we normalize the wave amplitudes to the upstream averaged magnetic field and the shock ramp amplitude, we find average values of ~40% and ~220%, respectively.

Proton-helium spectral anomaly as a signature of cosmic ray accelerator

The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants must hinge on the full consistency of acceleration theory with the observations; direct proof is impossible because of CR-orbit scrambling. Recent observations indicate deviations between helium and proton CR rigidity spectra inconsistent with the theory. By considering an initial (injection) phase of the diffusive shock acceleration, where elemental similarity does not apply, we demonstrate that the spectral difference is, in fact, a unique signature of the acceleration mechanism. Collisionless shocks inject more He(2+) when they are stronger and so produce harder He(2+) spectra. The injection bias is due to Alfvén waves driven by the more abundant protons, so the He(2+) ions are harder to trap by these waves. By fitting the p/He ratio to the PAMELA data, we bolster the diffusive shock acceleration case for resolving the century-old mystery of CR origin.

Dispersive nature of high mach number collisionless plasma shocks: Poynting flux of oblique whistler waves

Whistler wave trains are observed in the foot region of high Mach number quasiperpendicular shocks. The waves are oblique with respect to the ambient magnetic field as well as the shock normal. The Poynting flux of the waves is directed upstream in the shock normal frame starting from the ramp of the shock. This suggests that the waves are an integral part of the shock structure with the dispersive shock as the source of the waves. These observations lead to the conclusion that the shock ramp structure of supercritical high Mach number shocks is formed as a balance of dispersion and nonlinearity.