In this paper, we present a detailed study of the effects of the interplanetary magnetic field (IMF) strength on the foreshock properties at small and large scales. Two simulation runs performed with ...the hybrid‐Vlasov code Vlasiator with identical setup but with different IMF strengths, namely, 5 and 10 nT, are compared. We find that the bow shock position and shape are roughly identical in both runs, due to the quasi‐radial IMF orientation, in agreement with previous magnetohydrodynamic simulations and theory. Foreshock waves develop in a broader region in the higher IMF strength run, which we attribute to the larger growth rate of the waves. The velocity of field‐aligned beams remains essentially the same, but their density is generally lower when the IMF strength increases, due to the lower Mach number. Also, we identify in the regular IMF strength run ridges of suprathermal ions which disappear at higher IMF strength. These structures may be a new signature of the foreshock compressional boundary. The foreshock wave field is structured over smaller scales in higher IMF conditions, due to both the period of the foreshock waves and the transverse extent of the wave fronts being smaller. While the foreshock is mostly permeated by monochromatic waves at typical IMF strength, we find that magnetosonic waves at different frequencies coexist in the other run. They are generated by multiple beams of suprathermal ions, while only a single beam is observed at typical IMF strength. The consequences of these differences for solar wind‐magnetosphere coupling are discussed.
Plain Language Summary
Our solar system is filled with a stream of particles escaping from the Sun, called the solar wind. The Earth is shielded from these particles by its magnetic field, which creates a magnetic bubble around our planet, the magnetosphere. Because the solar wind flow is supersonic, a bow shock forms in front of the magnetosphere to slow it down. The outermost region of the near‐Earth space is called the foreshock. It is a very turbulent region, filled with particles reflected off the Earth's bow shock, and with a variety of magnetic waves. These waves can be transmitted inside the magnetosphere and create disturbances in the magnetic field on the Earth's surface. In this work, we use supercomputer simulations to study how the foreshock changes when the solar magnetic field, carried by the solar wind, intensifies. This happens in particular during solar storms, which create stormy space weather at Earth and can have adverse consequences on, for example, spacecraft electronics and power grids. We find that the foreshock properties are very different during these events compared to normal conditions and that these changes may have consequences in the regions closer to Earth.
Key Points
We study the effects of the interplanetary magnetic field strength on the foreshock properties using hybrid‐Vlasov simulations
The foreshock is structured over smaller scales for higher magnetic field strength; waves are less monochromatic and grow in a larger region
Suprathermal ion density decreases at higher field strength, and we identify a possible new signature of the foreshock compressional boundary
Context. Solar energetic particles observed in association with coronal mass ejections (CMEs) are produced by the CME-driven shock waves. The acceleration of particles is considered to be due to ...diffusive shock acceleration (DSA). Aims. We aim at a better understanding of DSA in the case of quasi-parallel shocks, in which self-generated turbulence in the shock vicinity plays a key role. Methods. We have developed and applied a new Monte Carlo simulation code for acceleration of protons in parallel coronal shocks. The code performs a self-consistent calculation of resonant interactions of particles with Alfvén waves based on the quasi-linear theory. In contrast to the existing Monte Carlo codes of DSA, the new code features the full quasi-linear resonance condition of particle pitch-angle scattering. This allows us to take anisotropy of particle pitch-angle scattering into account, while the older codes implement an approximate resonance condition leading to isotropic scattering. We performed simulations with the new code and with an old code, applying the same initial and boundary conditions, and have compared the results provided by both codes with each other, and with the predictions of the steady-state theory. Results. We have found that anisotropic pitch-angle scattering leads to less efficient acceleration of particles than isotropic. However, extrapolations to particle injection rates higher than those we were able to use suggest the capability of DSA to produce relativistic particles. The particle and wave distributions in the foreshock as well as their time evolution, provided by our new simulation code, are significantly different from the previous results and from the steady-state theory. Specifically, the mean free path in the simulations with the new code is increasing with energy, in contrast to the theoretical result.
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.
Flux transfer events (FTEs) are transient magnetic flux ropes at Earth's dayside magnetopause formed due to magnetic reconnection. As they move across the magnetopause surface, they can generate ...disturbances in the ultralow frequency (ULF) range, which then propagate into the magnetosphere. This study provides evidence of ULF waves in the Pc2 wave frequency range (>0.1 Hz) caused by FTEs during dayside reconnection using a global 3D hybrid‐Vlasov simulation (Vlasiator). These waves resulted from FTE formation and propagation at the magnetopause are particularly associated with large, rapidly moving FTEs. The wave power is stronger in the morning than afternoon, showing local time asymmetry. In the pre and postnoon equatorial regions, significant poloidal and toroidal components are present alongside the compressional component. The noon sector, with fewer FTEs, has lower wave power and limited magnetospheric propagation.
Plain Language Summary
The Earth's magnetosphere is a dynamic region shaped by the interplay between the solar wind and Earth's magnetic field. This interaction occurs at the boundary of the magnetosphere (magnetopause) through a process known as magnetic reconnection, giving rise to Flux Transfer Events (FTEs), which are magnetic structures that carry flux and energy into the magnetosphere. These FTEs form either in sudden bursts, patchy patterns or in a continuous, and relatively stable way making the magnetopause surface dynamic. As the FTEs move along the boundary of the magnetosphere, they create compressed regions and lead to wave generation that can extend into the magnetosphere. The study uses an advanced 3D hybrid‐Vlasov simulation model to analyze waves originated from FTE formation and propagation at the magnetopause. We find that rapidly moving and large FTEs have a significant impact on the magnetopause, leading to the generation of ULF waves with frequency above 0.1 Hz. This shows first direct evidence supporting previous theoretical speculations regarding the ability of FTEs to generate waves near the magnetopause.
Key Points
Dayside Pc2 waves (>0.1 Hz) have been detected in a 3D hybrid‐Vlasov simulation
These waves exhibit lower intensity within the magnetosphere at noon, compared to the prenoon and postnoon sectors
Pc2 waves observed in the simulation are associated with largest and fast moving flux transfer events initiated by subsolar reconnection
Context. The source of high-energy protons (above ~500 MeV) responsible for ground level enhancements (GLEs) remains an open question in solar physics. One of the candidates is a shock wave driven by ...a coronal mass ejection, which is thought to accelerate particles via diffusive-shock acceleration. Aims. We perform physics-based simulations of proton acceleration using information on the shock and ambient plasma parameters derived from the observation of a real GLE event. We analyse the simulation results to find out which of the parameters are significant in controlling the acceleration efficiency and to get a better understanding of the conditions under which the shock can produce relativistic protons. Methods. We use the results of the recently developed technique to determine the shock and ambient plasma parameters, applied to the 17 May 2012 GLE event, and carry out proton acceleration simulations with the Coronal Shock Acceleration (CSA) model. Results. We performed proton acceleration simulations for nine individual magnetic field lines characterised by various plasma conditions. Analysis of the simulation results shows that the acceleration efficiency of the shock, i.e. its ability to accelerate particles to high energies, tends to be higher for those shock portions that are characterised by higher values of the scattering-centre compression ratio rc and/or the fast-mode Mach number MFM. At the same time, the acceleration efficiency can be strengthened by enhanced plasma density in the corresponding flux tube. The simulations show that protons can be accelerated to GLE energies in the shock portions characterised by the highest values of rc. Analysis of the delays between the flare onset and the production times of protons of 1 GV rigidity for different field lines in our simulations, and a subsequent comparison of those with the observed values indicate a possibility that quasi-perpendicular portions of the shock play the main role in producing relativistic protons.
ABSTRACT The event-averaged charge state of heavy ion solar energetic particles (SEPs), measured at 1 au from the Sun, typically increases with the ions' kinetic energy. The origin of this behavior ...has been ascribed to processes taking place within the acceleration region. In this paper we study the propagation through interplanetary space of SEP Fe ions, injected near the Sun with a variety of charge states that are uniformly distributed in energy, by means of a 3D test particle model. In our simulations, due to gradient and curvature drifts associated with the Parker spiral magnetic field, ions of different charge propagate with very different efficiencies to an observer that is not magnetically well connected to the source region. As a result we find that, for many observer locations, the 1 au event-averaged charge state , as obtained from our model, displays an increase with particle energy E, in qualitative agreement with spacecraft observations. We conclude that drift-associated propagation is a possible explanation for the observed distribution of versus E in SEP events, and that the distribution measured in interplanetary space cannot be taken to represent that at injection.
The foreshock, extending upstream of Earth's bow shock, is a region of intense electromagnetic wave activity and nonlinear phenomena, which can have global effects on geospace. It is also the first ...geophysical region encountered by solar wind disturbances journeying toward Earth. Here, we present the first observations of considerable modifications of the foreshock wave field during extreme events of solar origin called magnetic clouds. Cluster's multispacecraft data reveal that the typical quasi‐monochromatic foreshock waves can be completely replaced by a superposition of waves each with shorter correlation lengths. Global numerical simulations further confirm that the foreshock wave field is more intricate and organized at smaller scales. Ion measurements suggest that changes in shock‐reflected particle properties may cause these modifications of the wave field. This state of the foreshock is encountered only during extreme events at Earth, but intense magnetic fields are typical close to the Sun or other stars.
Plain Language Summary
Solar storms are giant clouds of particles ejected from the Sun into space during solar eruptions. When solar storms are directed toward Earth, they can cause large disturbances in near‐Earth space, for example, disrupting communications or damaging spacecraft electronics. Understanding in detail what happens when solar storms reach Earth is crucial to mitigate their effects. Using measurements from the Cluster spacecraft, we investigate how solar storms modify the properties of the very first region of near‐Earth space they encounter when journeying toward Earth. This region, called the foreshock, extends ahead of the protective bubble formed by the Earth's magnetic field, the magnetosphere. The foreshock is home to intense electromagnetic waves, and disturbances in this region can perturb the Earth's magnetosphere. Our study reveals that solar storms modify profoundly the foreshock, resulting in a more complex wave activity. Global numerical simulations performed with the Vlasiator code confirm our findings. These changes could affect the regions of space closer to Earth, for example, in modifying the wave properties or the amount of solar particles entering the Earth's magnetosphere. This needs to be taken into account to better anticipate the effects of solar storms at Earth.
Key Points
When reaching geospace, magnetic clouds modify significantly the properties of the first geophysical region they encounter, the foreshock
Typical quasi‐monochromatic foreshock waves are replaced by a superposition of waves at different periods with a shorter transverse extent
Multiple field‐aligned beams observed during one event suggest a link between the multiple wave periods and the suprathermal ion properties
Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause, the magnetosheath and the polar cusps by ...remote sensing techniques. The Solar wind–Magnetosphere–Ionosphere Link Explorer (SMILE) and Lunar Environment heliospheric X-ray Imager (LEXI) missions aim to obtain soft X-ray images of near-Earth space thanks to their Soft X-ray Imager (SXI) instruments. While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission, the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma. To investigate the possible signatures of ion-kinetic-scale processes in soft X-ray images, we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model. The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field. We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations, providing meridional and equatorial views. We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirror mode wave structures in the magnetosheath and flux transfer events (FTEs) at the magnetopause. Our results suggest that, although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers, mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images. For instance, a local increase by 30% in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12% enhancement in the line-of-sight- and time integrated soft X-ray emissivity originating from this region. Likewise, a proton density increase by 14% in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%. These are likely conservative estimates, given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind. These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale, transient phenomena occurring on the dayside.
Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional ...approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution.
Aims. We model the transport of GLE-energy solar protons using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. We derive 1 AU observables and compare the simulation results with data from PAMELA.
Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data.
Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the position of the observer and is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, resulting in a rollover in the spectrum.