Earth's magnetosheath forms the interface between our planetary magnetosphere and the supermagnetosonic solar wind flow. The magnetosheath is bounded by the bow shock separating it from the upstream solar wind, and by the magnetopause, which is the outer boundary of the magnetosphere. An important factor determining the nature of the bow shock and the magnetosheath is the orientation of the interplanetary magnetic field (IMF). Where the bow shock normal direction is quasi-parallel to the IMF, solar wind particles reflected at the bow shock can stream back toward the Sun, forming the turbulent foreshock upstream of the shock (e.g., Eastwood et al., 2005;Hoppe et al., 1981). Consequently, the magnetosheath behind the quasi-parallel bow shock is more turbulent (e.g., Dimmock et al., 2014;Gutynska et al., 2015) than the quasi-perpendicular magnetosheath. In recent years, one of the most important avenues in magnetosheath research has been the study of magnetosheath jets (e.g., Plaschke et al., 2018), which are regions of high dynamic pressure. Magnetosheath jets are widely associated with the quasi-parallel magnetosheath

Abstract. 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.

Abstract. Magnetosheath jets are regions of high dynamic pressure, which 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 making 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 Multi-Scale (MMS) spacecraft data (6142 in total) and confirm the Vlasiator jet general behaviour using a statistics of 924 simulated individual jets. We find that the jets in the simulation are in excellent 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. 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.

Abstract. 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.

We present methods and algorithms that allow the Vlasiator code to run global, three-dimensional hybrid-Vlasov simulations of Earth's entire magnetosphere. The key ingredients that make Vlasov simulations at magnetospheric scales possible are the sparse velocity space implementation and spatial adaptive mesh refinement. We outline the algorithmic improvement of the semi-Lagrangian solver for six-dimensional phase space quantities, discuss the coupling of Vlasov and Maxwell equations' solvers in a refined mesh, and provide performance figures from simulation test runs that demonstrate the scalability of this simulation system to full magnetospheric runs.

Abstract. 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.

Abstract. 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 accurately the polarity of the asymmetries, 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. We investigate how the asymmetries change when the angle between the IMF and the Sun-Earth line is reduced and when the Alfven 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 Alfven Mach number leads to a reduced magnetic field asymmetry and a decrease in the variability of the magnetosheath density and velocity, the latter likely due to weaker foreshock processes. Our results highlight the strong impact of the foreshock on global magnetosheath properties, in particular on the magnetosheath density, which is extremely sensitive to transient foreshock processes.

The precipitation of charged particles from the magnetosphere into the ionosphere is one of the crucial coupling mechanisms between these two regions of geospace and is associated with multiple space weather effects, such as global navigation satellite system signal disruption and geomagnetically induced currents at ground level. While precipitating particle fluxes have been measured by numerous spacecraft missions over the past decades, it often remains difficult to obtain global precipitation patterns with a good time resolution during a substorm. Numerical simulations can help to bridge this gap and improve the understanding of mechanisms leading to particle precipitation at high latitudes through the global view they offer on the near-Earth space system. We present the first results on auroral (0.5–50 keV) proton precipitation within a 3-dimensional simulation of the Vlasiator hybrid-Vlasov model. The run is driven by southward interplanetary magnetic field conditions with constant solar wind parameters. We find that, on the dayside, cusp proton precipitation exhibits the expected energy–latitude dispersion and takes place in the form of successive bursts associated with the transit of flux transfer events formed through dayside magnetopause reconnection. On the nightside, the precipitation takes place within the expected range of geomagnetic latitudes, and it appears clearly that the precipitating particle injection is taking place within a narrow magnetic local time span, associated with fast Earthward plasma flows in the near-Earth magnetotail. Finally, the simulated precipitating fluxes are compared to observations from Defense Meteorological Satellite Program spacecraft during driving conditions similar to those in the simulation and are found to be in good agreement with the measurements.