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 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 degrees, 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 the 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.
Foreshock cavitons and spontaneous hot flow anomalies: A
statistical study with a global hybrid-Vlasov simulation
