Explore the vast range of titles available.

Shocks are one of nature’s most powerful particle accelerators and have been connected to relativistic electron acceleration and cosmic rays. Upstream shock observations include wave generation, wave-particle interactions and magnetic compressive structures, while at the shock and downstream, particle acceleration, magnetic reconnection and plasma jets can be observed. Here, using Magnetospheric Multiscale (MMS) we show in-situ evidence of high-speed downstream flows (jets) generated at the Earth’s bow shock as a direct consequence of shock reformation. Jets are observed downstream due to a combined effect of upstream plasma wave evolution and an ongoing reformation cycle of the bow shock. This generation process can also be applicable to planetary and astrophysical plasmas where collisionless shocks are commonly found. Several mechanisms exist for formation of jets observed in Earth’s magnetosheath. Here, the authors show evidence of high-speed downstream flows generated at the Earth’s bow shock as a direct consequence of shock reformation, which is different than the proposed mechanisms.

Electron heating at collisionless shocks in space is a combination of adiabatic heating due to large‐scale electric and magnetic fields and non‐adiabatic scattering by high‐frequency fluctuations. The scales at which heating happens hints to what physical processes are taking place. In this letter, we study electron heating scales with data from the Magnetospheric Multiscale (MMS) spacecraft at Earth's quasi‐perpendicular bow shock. We utilize the tight tetrahedron formation and high‐resolution plasma measurements of MMS to directly measure the electron temperature gradient. From this, we reconstruct the electron temperature profile inside the shock ramp and find that the electron temperature increase takes place on ion or sub‐ion scales. Further, we use Liouville mapping to investigate the electron distributions through the ramp to estimate the deHoffmann‐Teller potential and electric field. We find that electron heating is highly non‐adiabatic at the high‐Mach number shocks studied here. Plain Language Summary Shock waves appear whenever a supersonic medium, such as a plasma, encounters an obstacle. The plasma, which consists of charged ions and free electrons, is heated by the shock wave through interactions with the electromagnetic fields. In this work, we investigate how electrons are heated at plasma shocks. A key parameter to electron heating is the thickness of the layer where the heating takes place. Here, we use observations from the four Magnetospheric Multiscale spacecraft that regularly cross the standing bow shock that forms when the supersonic plasma, known as the solar wind, encounters Earth's magnetic field. We find that the thickness of the shock is larger than previously reported and is on the scales where ion physics dominate. We also find that the electron heating deviates significantly from simple adiabatic heating. Key Points Using multipoint data from Magnetospheric Multiscale, we find that electron heating takes place on ion scales in the quasi‐perpendicular shock ramp We show that the time series of the temperature does not represent the spatial profile due to varying shock ramp speed Electron distributions in the ramp and downstream of the shock show that electrons are heated non‐adiabatically

Particle precipitation is a central aspect of space weather, as it strongly couples the magnetosphere and the ionosphere and can be responsible for radio signal disruption at high latitudes. We present the first hybrid-Vlasov simulations of proton precipitation in the polar cusps. We use two runs from the Vlasiator model to compare cusp proton precipitation fluxes during southward and northward interplanetary magnetic field (IMF) driving. The simulations reproduce well-known features of cusp precipitation, such as a reverse dispersion of precipitating proton energies, with proton energies increasing with increasing geomagnetic latitude under northward IMF driving, and a nonreversed dispersion under southward IMF driving. The cusp is also found more polewards in the northward IMF simulation than in the southward IMF simulation. In addition, we find that the bursty precipitation during southward IMF driving is associated with the transit of flux transfer events in the vicinity of the cusp. In the northward IMF simulation, dual lobe reconnection takes place. As a consequence, in addition to the high-latitude precipitation spot associated with the lobe reconnection from the same hemisphere, we observe lower-latitude precipitating protons which originate from the opposite hemisphere’s lobe reconnection site. The proton velocity distribution functions along the newly closed dayside magnetic field lines exhibit multiple proton beams travelling parallel and antiparallel to the magnetic field direction, which is consistent with previously reported observations with the Cluster spacecraft. In both runs, clear electromagnetic ion cyclotron waves are generated in the cusps and might further increase the calculated precipitating fluxes by scattering protons to the loss cone in the low-altitude cusp. Global kinetic simulations can improve the understanding of space weather by providing a detailed physical description of the entire near-Earth space and its internal couplings.

The Earth’s magnetosphere and its bow shock, which is formed by the interaction of the supersonic solar wind with the terrestrial magnetic field, constitute a rich natural laboratory enabling in situ investigations of universal plasma processes. Under suitable interplanetary magnetic field conditions, a foreshock with intense wave activity forms upstream of the bow shock. So-called 30 s waves, named after their typical period at Earth, are the dominant wave mode in the foreshock and play an important role in modulating the shape of the shock front and affect particle reflection at the shock. These waves are also observed inside the magnetosphere and down to the Earth’s surface, but how they are transmitted through the bow shock remains unknown. By combining state-of-the-art global numerical simulations and spacecraft observations, we demonstrate that the interaction of foreshock waves with the shock generates earthward-propagating, fast-mode waves, which reach the magnetosphere. These findings give crucial insight into the interaction of waves with collisionless shocks in general and their impact on the downstream medium.The Earth’s bow shock results from the interaction of the solar wind with the terrestrial magnetic field. With global numerical simulations and spacecraft observations, the transmission of fast magnetosonic waves through the bow shock is revealed.

Collisionless shocks can exhibit non‐stationary behavior even under steady upstream conditions, forming a complex transition region. Ion phase‐space holes, linked to shock self‐reformation and surface ripples, are a signature of this non‐stationarity. We statistically analyze their occurrence using 521 crossings of Earth's quasi‐perpendicular bow shock. Phase‐space holes appear in 65% of cases, though the actual rate may be higher as the holes may not be resolved during fast shock crossings. The occurrence rate peaks at 70% for shocks with Alfvén Mach numbers MA>7 ${M}_{A} > 7$. These findings suggest that Earth's quasi‐perpendicular bow shock is predominantly non‐stationary.

Magnetosheath jets are regions of high dynamic pressure, which can traverse from the bow shock towards the magnetopause. Recent modelling efforts, limited to a single jet and a single set of upstream conditions, have provided the first estimations about how the jet parameters behave as a function of position within the magnetosheath. Here we expand the earlier results by doing the first statistical investigation of the jet dimensions and parameters as a function of their lifetime within the magnetosheath. To verify the simulation behaviour, we first identify jets from Magnetosphere Multiscale (MMS) spacecraft data (6142 in total) and confirm the Vlasiator jet general behaviour using statistics of 924 simulated individual jets. We find that the jets in the simulation are in quantitative agreement with the observations, confirming earlier findings related to jets using Vlasiator. The jet density, dynamic pressure, and magnetic field intensity show a sharp jump at the bow shock, which decreases towards the magnetopause. The jets appear compressive and cooler than the magnetosheath at the bow shock, while during their propagation towards the magnetopause they thermalise. Further, the shape of the jets flatten as they progress through the magnetosheath. They are able to maintain their flow velocity and direction within the magnetosheath flow, and they end up preferentially to the side of the magnetosheath behind the quasi-parallel shock. Finally, we find that Vlasiator jets during low solar wind Alfvén Mach number MA are shorter in duration, smaller in their extent, and weaker in terms of dynamic pressure and magnetic field intensity as compared to the jets during high MA.

Kinetically driven plasma waves are fundamental for a description of the thermodynamical properties of the Earth's magnetosheath. The most commonly observed ion-scale instabilities are generated by temperature anisotropy of the ions, such as the mirror and proton cyclotron instabilities. We investigate here the spatial resolution dependence of the mirror and proton cyclotron instabilities in a global hybrid-Vlasov simulation using the Vlasiator model; we do this in order to find optimal resolutions and help future global hybrid-Vlasov simulations to save resources when investigating those instabilities in the magnetosheath. We compare the proton velocity distribution functions, power spectra and growth rates of the instabilities in a set of simulations with three different spatial resolutions but otherwise identical set-up. We find that the proton cyclotron instability is absent at the lowest resolution and that only the mirror instability remains, which leads to an increased temperature anisotropy in the simulation. We conclude that the proton cyclotron instability, its saturation and the reduction of the anisotropy to marginal levels are resolved at the highest spatial resolution. A further increase in resolution does not lead to a better description of the instability to an extent that would justify this increase at the cost of numerical resources in future simulations. We also find that spatial resolutions between 1.32 and 2.64 times the inertial length in the solar wind present acceptable limits for the resolution within which the velocity distribution functions resulting from the proton cyclotron instability are still bi-Maxwellian and reach marginal stability levels. Our results allow us to determine a range of spatial resolutions suitable for the modelling of the proton cyclotron and mirror instabilities and should be taken into consideration regarding the optimal grid spacing for the modelling of these two instabilities, within available computational resources.

The foreshock located upstream of Earth's bow shock hosts a wide variety of phenomena related to the reflection of solar wind particles from the bow shock and the subsequent formation of ultra-low-frequency (ULF) waves. In this work, we investigate foreshock cavitons, which are transient structures resulting from the non-linear evolution of ULF waves, and spontaneous hot flow anomalies (SHFAs), which are thought to evolve from cavitons as they accumulate suprathermal ions while being carried to the bow shock by the solar wind. Using the global hybrid-Vlasov simulation model Vlasiator, we have conducted a statistical study in which we track the motion of individual cavitons and SHFAs in order to examine their properties and evolution. In our simulation run where the interplanetary magnetic field (IMF) is directed at a sunward–southward angle of 45∘, continuous formation of cavitons is found up to ∼11 Earth radii (RE) from the bow shock (along the IMF direction), and caviton-to-SHFA evolution takes place within ∼2 RE from the shock. A third of the cavitons in our run evolve into SHFAs, and we find a comparable amount of SHFAs forming independently near the bow shock. We compare the properties of cavitons and SHFAs to prior spacecraft observations and simulations, finding good agreement. We also investigate the variation of the properties as a function of position in the foreshock, showing that transients close to the bow shock are associated with larger depletions in the plasma density and magnetic field magnitude, along with larger increases in the plasma temperature and the level of bulk flow deflection. Our measurements of the propagation velocities of cavitons and SHFAs agree with earlier studies, showing that the transients propagate sunward in the solar wind rest frame. We show that SHFAs have a greater solar wind rest frame propagation speed than cavitons, which is related to an increase in the magnetosonic speed near the bow shock.

Bounded by the bow shock and the magnetopause, the magnetosheath forms the interface between solar wind and magnetospheric plasmas and regulates solar wind–magnetosphere coupling. Previous works have revealed pronounced dawn–dusk asymmetries in the magnetosheath properties. The dependence of these asymmetries on the upstream parameters remains however largely unknown. One of the main sources of these asymmetries is the bow shock configuration, which is typically quasi-parallel on the dawn side and quasi-perpendicular on the dusk side of the terrestrial magnetosheath because of the Parker spiral orientation of the interplanetary magnetic field (IMF) at Earth. Most of these previous studies rely on collections of spacecraft measurements associated with a wide range of upstream conditions which are processed in order to obtain average values of the magnetosheath parameters. In this work, we use a different approach and quantify the magnetosheath asymmetries in global hybrid-Vlasov simulations performed with the Vlasiator model. We concentrate on three parameters: the magnetic field strength, the plasma density, and the flow velocity. We find that the Vlasiator model reproduces the polarity of the asymmetries accurately but that their level tends to be higher than in spacecraft measurements, probably because the magnetosheath parameters are obtained from a single set of upstream conditions in the simulation, making the asymmetries more prominent. A set of three runs with different upstream conditions allows us to investigate for the first time how the asymmetries change when the angle between the IMF and the Sun–Earth line is reduced and when the Alfvén Mach number decreases. We find that a more radial IMF results in a stronger magnetic field asymmetry and a larger variability of the magnetosheath density. In contrast, a lower Alfvén Mach number leads to a reduced magnetic field asymmetry and a decrease in the variability of the magnetosheath density, the latter likely due to weaker foreshock processes. Our results highlight the strong impact of the quasi-parallel shock and its associated foreshock on global magnetosheath properties, in particular on the magnetosheath density, which is extremely sensitive to transient quasi-parallel shock processes, even with the perfectly steady upstream conditions in our simulations. This could explain the large variability of the density asymmetry levels obtained from spacecraft measurements in previous studies.