Quasi-perpendicular Shock Research Articles

Dynamics of Plasma Turbulence at Earth’s Bow Shock and through the Magnetosheath

Abstract Earth’s magnetosheath can be treated as a natural laboratory to study turbulence development in confined space. The present study focuses on the characteristics of turbulent cascade downstream of the bow shock, where properties of turbulence are known to differ from those in the upstream solar wind. Characteristics of the turbulent spectrum are considered in two distinct points of the magnetosheath for two case studies. The analysis is based on high-resolution measurements of plasma parameters by the Spektr-R spacecraft and magnetic field data by the Themis/Arthemis mission. The measurements are performed for two distinct cases: in the dayside magnetosheath behind the quasi-perpendicular bow shock and in the nightside flank of the magnetosheath behind the quasi-parallel bow shock. The study focuses on the scales at which kinetic effects in plasma become significant and the turbulent spectrum is known to break. The analysis reveals that modification of the fluctuation spectrum at the bow shock is controlled by the distance of the measurement point from the bow shock’s nose. Also, performed statistical results suggest the influence of the large-scale parameters of the upstream solar wind and the type of the bow shock on the turbulent spectrum’s modification in the downstream region.

Solar Longitude Distribution of High-energy Proton Flares: Fluences and Spectra

Abstract The distribution of the longitudes of solar flares associated with the high-energy proton events called ground level events (GLEs) can be approximated by a Gaussian with a peak at ∼W60, with a full range from ∼E90 to ∼W150. The longitudes of flares associated with the top third (24 of 72) of GLEs in terms of their >430 MeV fluences (F 430) are primarily distributed over E20–W100 with a skew toward disk center. This 120° span in longitude is comparable to the latitudinal spans of powerful coronal mass ejections (CMEs) from limb flares. Only 5 of 24 strong GLEs are located within the W40–80 zone of good magnetic connection to Earth. GLEs with hard spectra, i.e., a spectral index SI30/200(= log(F 30/F 200)) < 1.5, also tend to avoid W40–80 source regions. Three-fourths of such events (16 of 21) arise in flares outside this range. The above tendencies favor a CME-driven shock source over a flare-resident acceleration process for high-energy solar protons. GLE spectra show a trend, with broad scatter, from hard spectra for events originating in eruptive flares beyond the west limb to soft spectra for GLEs with sources near central meridian. This behavior can be explained in terms of: (1) dominant near-Sun quasi-perpendicular shock acceleration of protons for far western (>W100) GLEs; (2) quasi-parallel shock acceleration for well-connected (W40–80) GLEs, and (3) proton acceleration/trapping at CME-driven bow shocks from central meridian (E20–W20) that strike the Earth.

Electron Acceleration in One-dimensional Nonrelativistic Quasi-perpendicular Collisionless Shocks

Abstract We study diffusive shock acceleration (DSA) of electrons in nonrelativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell simulations. By exploring the parameter space of sonic and Alfvénic Mach numbers we find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive nonresonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves, and undergo multiple cycles of shock-drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.

Shock waves in the magnetized cosmic web: the role of obliquity and cosmic ray acceleration

ABSTRACT Structure formation shocks are believed to be the largest accelerators of cosmic rays in the Universe. However, little is still known about their efficiency in accelerating relativistic electrons and protons as a function of their magnetization properties, i.e. of their magnetic field strength and topology. In this work, we analysed both uniform and adaptive mesh resolution simulations of large-scale structures with the magnetohydrodynamical grid code enzo, studying the dependence of shock obliquity with different realistic scenarios of cosmic magnetism. We found that shock obliquities are more often perpendicular than what would be expected from a random 3D distribution of vectors, and that this effect is particularly prominent in the proximity of filaments, due to the action of local shear motions. By coupling these results to recent works from particle-in-cell simulations, we estimated the flux of cosmic ray protons in galaxy clusters, and showed that in principle the riddle of the missed detection of hadronic γ-ray emission by the Fermi-LAT can be explained if only quasi-parallel shocks accelerate protons. On the other hand, for most of the cosmic web the acceleration of cosmic ray electrons is still allowed, due to the abundance of quasi-perpendicular shocks. We discuss quantitative differences between the analysed models of magnetization of cosmic structures, which become more significant at low cosmic overdensities.

On the Nature and Origin of Bipolar Electrostatic Structures in the Earth's Bow Shock

We present a statistical analysis of large-amplitude bipolar electrostatic structures measured by the Magnetospheric Multiscale spacecraft in the Earth's bow shock. The analysis is based on 371 large-amplitude bipolar structures collected in nine supercritical quasi-perpendicular Earth's bow shock crossings. We have found that 361 of the bipolar structures have negative electrostatic potentials, and only 10 structures (> 1, a parameter range typical of collisionless shocks in the heliosphere and various astrophysical environments. This analysis indicates that the original mechanism of electron surfing acceleration involving electron phase space holes is not likely to be efficient in realistic collisionless shocks.

Constraining the coherence scale of the interstellar magnetic field using TeV gamma-ray observations of supernova remnants

ABSTRACT Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario, the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e. when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale λB around Vela Jr and SN1006 to $\lambda _B \simeq 13_{-4.3}^{+13}$ pc and $\lambda _B \gt 200_{-40}^{+50}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multifrequency spectra from the radio to TeV gamma-rays and match the observed gamma-ray morphology. Finally, to capture the propagation of an SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30 per cent for realistic parameters.

Nonstationary Quasiperpendicular Shock and Ion Reflection at Mars

AbstractCollisionless shocks in space plasma are regions of heating and acceleration of charged particles and dissipation of kinetic energy. These accelerated particles are the source of electromagnetic emissions from supernova remnants and other astrophysical structures. At high Mach numbers, shocks can be inherently nonstationary and exhibit modulated energy transfer and recurring plasma compression areas in the form of reformation. We use data from the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft to study reformation of the Martian bow shock, which has a relatively high curvature compared to that at Earth, and the upstream solar wind is often mass loaded with a population of pickup ions. We show evidence of ion reflection effects in reformation of a supercritical quasiperpendicular shock.

Suprathermal Electron Acceleration by a Quasi-perpendicular Shock: Simulations and Observations

Abstract The acceleration of suprathermal electrons in the solar wind is mainly associated with shocks driven by interplanetary coronal mass ejections (ICMEs). It is well known that the acceleration of electrons is much more efficient at quasi-perpendicular shocks than at quasi-parallel ones. Yang et al. studied the acceleration of suprathermal electrons with observations at a quasi-perpendicular ICME-driven shock event to claim the important role of shock-drift acceleration (SDA). Here, we perform test-particle simulations to study the acceleration of electrons in this event, by calculating the downstream electron intensity distribution for all energy channels assuming an initial distribution based on the average upstream intensities. Using simulations, we obtain the results similar to the observations from Yang et al. as follows. It is shown that the ratio of downstream to upstream intensities peaks at about 90° pitch angle. In addition, in each pitch angle direction the downstream electron energy spectral index is much larger than the theoretical index of diffusive shock acceleration. Furthermore, the estimated drift length is proportional to the electron energy but the drift time is almost energy independent. Finally, we use a theoretical model based on SDA to describe the drift length and drift time especially, to explain their energy dependence. These results indicate the importance of SDA in the acceleration of electrons by quasi-perpendicular shocks.

Collisionless Shocks as a Diagnostic Tool for Understanding Energetic Particle Transport in Space Plasmas

We study the transport of energetic particles accelerated at three different shock events observed in the solar wind by the ACE spacecraft. We consider particle propagation for a quasi-parallel, an oblique, and a quasi-perpendicular shock. The transport regime is deduced from the shape of the energetic particle profiles upstream of the shock, and for these events the profiles are well fitted by power-laws with slope $\beta$. This corresponds to a superdiffusive transport with the anomalous diffusion exponent $\alpha = 2 - \beta$ when $\beta < 1$, and to normal diffusion when $\beta\ge1$. We checked the resonant turbulence level upstream of the shocks, finding that this is statistically constant, so that the transport regime is not expected to change with the shock distance. For the three shocks under study, particle transport upstream of the shock is mostly superdiffusive, although the superdiffusive character appears to diminish with the increase of the shock normal angle $\theta_{\rm Bn}$. We discuss possible interpretations of these results.

Gamma-Ray and Neutrino Emissions due to Cosmic-Ray Protons Accelerated at Intracluster Shocks in Galaxy Clusters

Abstract We examine the cosmic-ray protons (CRp) accelerated at collisionless shocks in galaxy clusters using cosmological structure formation simulations. We find that in the intracluster medium (ICM) within the virial radius of simulated clusters, only ∼7% of shock kinetic energy flux is dissipated by the shocks that are expected to accelerate CRp—that is, supercritical, quasi-parallel (Q ∥) shocks with sonic Mach number M s ≥ 2.25. The rest is dissipated at subcritical shocks and quasi-perpendicular shocks, both of which may not accelerate CRp. Adopting the diffusive shock acceleration (DSA) model recently presented in Ryu et al., we quantify the DSA of CRp in simulated clusters. The average fraction of the shock kinetic energy transferred to CRp via DSA is assessed at ∼(1–2) × 10−4. We also examine the energization of CRp through reacceleration using a model based on the test-particle solution. Assuming that the ICM plasma passes through shocks three times on average through the history of the universe and that CRp are reaccelerated only at supercritical Q ∥-shocks, the CRp spectrum flattens by ∼0.05–0.1 in slope and the total amount of CRp energy increases by ∼40%–80% from reacceleration. We then estimate diffuse γ-ray and neutrino emissions, resulting from inelastic collisions between CRp and thermal protons. The predicted γ-ray emissions from simulated clusters lie mostly below the upper limits set by Fermi-LAT for observed clusters. The neutrino fluxes toward nearby clusters would be ≲10−4 of the IceCube flux at E ν = 1 PeV and ≲10−6 of the atmospheric neutrino flux in the energy range of E ν ≤ 1 TeV.

Spatial Distribution of Electromagnetic Waves near the Proton Cyclotron Frequency in ICME Sheath Regions Associated with Quasi-perpendicular Shocks: Wind Observations

Abstract Electromagnetic waves (EMWs) near the proton cyclotron frequency f cp are transverse left-handed (LH) or right-handed (RH) polarized waves, and are ubiquitous in the solar wind. However, the characteristics of these waves in the sheath regions of interplanetary coronal mass ejections (ICMEs) are poorly understood. Through a comprehensive survey of Wind magnetic field and plasma data using dynamic spectra and repeated filtering analyses, 700 EMW events (7.1% of the analysis time) are identified in the 62 ICME sheath regions associated with quasi-perpendicular shocks involved with a low shock Mach number M f and low upstream β 1. In the ICME sheath regions, outward (inward)-propagating LH (RH) EMWs have relatively higher counts and longer duration than inward (outward)-propagating LH (RH) EMWs in the plasma frame, consistent with previous STEREO observations. The spatial distributions of the magnetic field, plasma, and frequency parameters of EMWs are also presented in both spacecraft and plasma frames, especially the proton (alpha) temperature anisotropy , α abundance N α /N p, and normalized differential alpha-proton speed V d/V A. After removing the Doppler shift, 81.1% (59%) of all outward (inward)-propagating LH EMWs have a frequency below (above) 0.5f cp, while 68.3% (64%) of all outward (inward)-propagating RH EMWs have a frequency smaller (greater) than 0.5f cp. Further investigations of local plasma parameters reveal that different excitation mechanisms for EMWs are in different subregions of the ICME sheath regions. These results are helpful in understanding the important role of EMWs in the solar wind–ICME coupling process with different sheath regions.

Electron Firehose Instabilities in High-β Intracluster Shocks

Abstract The preacceleration of electrons through reflection and shock drift acceleration (SDA) is essential for the diffusive shock acceleration of nonthermal electrons in collisionless shocks. Previous studies suggested that, in weak quasi-perpendicular (Q ⊥) shocks in the high-β ( ) intracluster medium (ICM), the temperature anisotropy due to SDA-reflected electrons can drive the electron firehose instability, which excites oblique nonpropagating waves in the shock foot. In this paper, we investigate, through a linear analysis and particle-in-cell simulations, the firehose instabilities driven by an electron temperature anisotropy (ETAFI) and also by a drifting electron beam (EBFI) in β ∼ 100 ICM plasmas. The EBFI should be more relevant to describing the self-excitation of upstream waves in Q ⊥-shocks, since backstreaming electrons in the shock foot behave more like an electron beam rather than an anisotropic bi-Maxwellian population. We find that the basic properties of the two instabilities, such as the growth rate, γ, and the wavenumber of fast-growing oblique modes, are similar in the ICM environment, with one exception; while the waves excited by the ETAFI are nonpropagating (ω r = 0), those excited by the EBFI have a nonzero frequency ( ). However, the frequency is small with ω r &lt; γ. Thus, we conclude that the interpretation of previous studies for the nature of upstream waves based on the ETAFI remains valid in Q ⊥-shocks in the ICM.

A study on the dynamic spectral indices for SEP events on 2000 July 14 and 2005 January 20

We have studied the dynamic proton spectra for the two solar energetic particle (SEP) events on 2000 July 14 (hereafter GLE59) and 2005 January 20 (hereafter GLE69). The source locations of GLE59 and GLE69 are N22W07 and N12W58 respectively. Proton fluxes >30 MeV have been used to compute the dynamic spectral indices of the two SEP events. Our results show that spectral indices of the two SEP events increased more swiftly at early times, suggesting that the proton fluxes >30 MeV might be accelerated particularly by the concurrent flares at early times for the two SEP events. For the GLE69 with source location at N12W58, both flare site and shock nose are well connected with the Earth at the earliest time. However, only the particles accelerated by the shock driven by eastern flank of the CME can propagate along the interplanetary magnetic field line to the Earth after the flare. For the GLE59 with source location at N22W07, only the particles accelerated by the shock driven by western flank of the associated CME can reach the Earth after the flare. Our results also show that there was slightly more than one hour during which the proton spectra for GLE69 are softer than that for GLE59 after the flares, suggesting that the shock driven by eastern flank of the CME associated with GLE69 is weaker than the shock driven by the western flank of the CME associated with GLE59. The results support that quasi-perpendicular shock has stronger potential in accelerating particles than the quasi-parallel shock. These results also suggest that only a small part of the shock driven by western flank of the CME associated with the GLE59 is quasi-perpendicular.

Shock Drift Acceleration of Ions in an Interplanetary Shock Observed by MMS

Abstract An interplanetary (IP) shock wave was recorded crossing the Magnetospheric Multiscale constellation on 2018 January 8. Plasma measurements upstream of the shock indicate efficient proton acceleration in the IP shock ramp: 2–7 keV protons are observed upstream for about three minutes (∼8000 km) ahead of the IP shock ramp, outrunning the upstream waves. The differential energy flux of 2–7 keV protons decays slowly with distance from the ramp toward the upstream region (dropping by about half within 8 Earth radii from the ramp) and is lessened by a factor of about four in the downstream compared to the ramp (within a distance comparable to the gyroradius of ∼keV protons). Comparison with test-particle simulations has confirmed that the mechanism accelerating the solar wind protons and injecting them upstream is classical Shock Drift Acceleration (SDA). This example of observed proton acceleration by a low-Mach, quasi-perpendicular shock may be applicable to astrophysical contexts, such as supernova remnants or the acceleration of cosmic rays.

A computational study of the properties of the quasi-perpendicular fast forward shock event during solar maximum

MHD shock waves and discontinuities are often observed in the solar wind. We studied a strong shock wave event on 7 June 2014, during a solar maximum. The properties of forward shock are investigated and all physical parameters of interplanetary shock are analyzed. The upstream parameters, such as Alfven velocity and sound speed are calculated, and the angle ( $${\theta }$$ ) between the upstream magnetic field direction and the shock normal is estimated. To determine the propagation speed of the shock, a full numerical solution of shock adiabatic equation is carried out. The thermal pressure for electrons and ions are calculated and also the pressure ratio and the entropy change across the shock are estimated. The results show that the shock strength is about 3.88 with $${\theta }\approx \text{72 }^{\circ }$$ (a quasi-perpendicular case), and shock propagation velocity is about 678 km/s. The increase in the specific entropy is also evaluated.

Jets in the magnetosheath: IMF control of where they occur

Abstract. Magnetosheath jets are localized regions of plasma that move faster towards the Earth than the surrounding magnetosheath plasma. Due to their high velocities, they can cause indentations when colliding into the magnetopause and trigger processes such as magnetic reconnection and magnetopause surface waves. We statistically study the occurrence of these jets in the subsolar magnetosheath using measurements from the five Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft and OMNI solar wind data from 2008 to 2011. We present the observations in the BIMF–vSW plane and study the spatial distribution of jets during different interplanetary magnetic field (IMF) orientations. Jets occur downstream of the quasi-parallel bow shock approximately 9 times as often as downstream of the quasi-perpendicular shock, suggesting that foreshock processes are responsible for most jets. For an oblique IMF, with 30–60∘ cone angle, the occurrence increases monotonically from the quasi-perpendicular side to the quasi-parallel side. This study offers predictability for the numbers, locations, and magnetopause impact rates of jets observed during different IMF orientations, allowing us to better forecast the formation of these jets and their impact on the magnetosphere.

Collisionless Shocks Driven by Supersonic Plasma Flows with Self-Generated Magnetic Fields.

Collisionless shocks are ubiquitous in the Universe as a consequence of supersonic plasma flows sweeping through interstellar and intergalactic media. These shocks are the cause of many observed astrophysical phenomena, but details of shock structure and behavior remain controversial because of the lack of ways to study them experimentally. Laboratory experiments reported here, with astrophysically relevant plasma parameters, demonstrate for the first time the formation of a quasiperpendicular magnetized collisionless shock. In the upstream it is fringed by a filamented turbulent region, a rudiment for a secondary Weibel-driven shock. This turbulent structure is found responsible for electron acceleration to energies exceeding the average energy by two orders of magnitude.

Properties of Stream Interactions and Their Associated Shocks near 1.52 au: MAVEN Observations

Abstract The measurements from the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft are analyzed for the basic properties of solar wind stream interaction regions (SIRs) and their associated shocks at 1.52 au, as well as the evolution of SIRs from 1 to 1.52 au. A total of 149 SIRs are identified during the period from 2014 October to 2018 November, and 126 SIRs with high-quality data are selected for this study. The average occurrence rate of SIRs at 1.52 au is 36.3 yr−1, which is comparable to but slightly higher than that (32.4 yr−1) at 1 au, meaning that most SIRs are well formed at 1 au. The average duration of SIRs at 1.52 au is about 37.0 hr, comparable to that at 1 au (36.73 hr), indicating that SIRs have not yet expanded more rapidly as they are convected to 1.52 au. The maximum magnetic field strength and pressure of SIRs decrease significantly from 1 to 1.52 au. The shock association rates of SIRs increase from 20.3% to 33.3% or higher as SIRs evolve from 1 to 1.52 au. The forward shocks tend to occur twice more frequently than the reverse shocks. About 75% of shocks at 1.52 au are quasi-perpendicular shocks. The strength of the shocks becomes weaker and the average shock speed remains almost unchanged from 1 to 1.52 au. These results will help us understand the solar wind conditions at Mars and their potential impact on the Martian space environment.

Substorm Activity and Orientation of the Front of a Shock Wave of an Interplanetary Magnetic Cloud

The paper draws attention to the complex structure of the fast magnetic clouds of the solar wind, which, in addition to the actual cloud body, contain a turbulent transition layer (cloud sheath) with a large and irregular magnetic field following the shock wave. The orientation of the plane of the magnetic-cloud shock wave with respect to the interplanetary magnetic field modified by the shock wave propagating in the solar wind for 33 cases of the registration of fast magnetic clouds has been calculated. The dependence of the substorm activity in the auroral zone on the level of turbulent processes occurring in the sheath of magnetic clouds has been studied. It was taken into account that the turbulent phenomena in the sheath are largely determined by the orientation of the shock-wave plane with respect to the interplanetary magnetic field. It is shown that the level of magnetic activity in the auroral zone, which is characterized by the integral AL index, increases with a decrease in the angle between the direction normal to the shock-wave front and the vector of the interplanetary magnetic field. Thus, the most geoeffective are magnetic clouds with a quasi-parallel shock wave, and the least geoeffective are those with a quasi-perpendicular shock wave. It was concluded that the intensity of turbulent processes in the cloud sheath increases with a decrease in the magnetic field penetrating the sheath, which plays a stabilizing role for turbulent magnetohydrodynamic perturbations.

Electron Preacceleration in Weak Quasi-perpendicular Shocks in High-beta Intracluster Medium

Abstract Giant radio relics in the outskirts of galaxy clusters are known to be lit up by the relativistic electrons produced via diffusive shock acceleration (DSA) in shocks with low sonic Mach numbers, M s ≲ 3. The particle acceleration at these collisionless shocks critically depends on the kinetic plasma processes that govern the injection to DSA. Here, we study the preacceleration of suprathermal electrons in weak, quasi-perpendicular (Q ⊥) shocks in the hot, high-β (β = P gas/P B) intracluster medium (ICM) through two-dimensional particle-in-cell simulations. Guo et al. showed that, in high-β Q ⊥-shocks, some of the incoming electrons could be reflected upstream and gain energy via shock drift acceleration (SDA). The temperature anisotropy due to the SDA-energized electrons then induces the electron firehose instability (EFI), and oblique waves are generated, leading to a Fermi-like process and multiple cycles of SDA in the preshock region. We find that such electron preacceleration is effective only in shocks above a critical Mach number . This means that, in ICM plasmas, Q ⊥-shocks with M s ≲ 2.3 may not efficiently accelerate electrons. We also find that, even in Q ⊥-shocks with M s ≳ 2.3, electrons may not reach high enough energies to be injected to the full Fermi-I process of DSA, because long-wavelength waves are not developed via the EFI alone. Our results indicate that additional electron preaccelerations are required for DSA in ICM shocks, and the presence of fossil relativistic electrons in the shock upstream region may be necessary to explain observed radio relics.