Hybrid-Vlasov simulation of auroral proton precipitation in the cusps : Comparison of northward and southward interplanetary magnetic field driving

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.Peer reviewe

Energy Flux Through the Magnetopause During Flux Transfer Events in Hybrid-Vlasov 2D Simulations

Solar wind-magnetosphere coupling drives magnetospheric dynamic phenomena by enabling energy exchange between magnetospheric and solar wind plasmas. In this study, we examine two-dimensional noon-midnight meridional plane simulation runs of the global hybrid-Vlasov code Vlasiator with southward interplanetary magnetic field driving. We compute the energy flux, which consists of the Poynting flux and hydrodynamic energy flux components, through the Earth's magnetopause during flux transfer events (FTEs). The results demonstrate the spatiotemporal variations of the energy flux along the magnetopause during an FTE, associating the FTE leading (trailing) edge with an energy injection into (escape from) the magnetosphere on the dayside. Furthermore, FTEs traveling along the magnetopause transport energy to the nightside magnetosphere. We identify the tail lobes as a primary entry region for solar wind energy into the magnetosphere, consistent with results from global magnetohydrodynamic simulations and observations.Peer reviewe

Asymmetries in the Earth's dayside magnetosheath : results from global hybrid-Vlasov simulations

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 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, 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.Peer reviewe

Non-locality of Earth's quasi-parallel bow shock : injection of thermal protons in a hybrid-Vlasov simulation

We study the interaction of solar wind protons with Earth's quasi-parallel bow shock using a hybrid-Vlasov simulation. We employ the global hybrid model Vlasiator to include effects due to bow shock curvature, tenuous upstream populations, and foreshock waves. We investigate the uncertainty of the position of the quasi-parallel bow shock as a function of several plasma properties and find that regions of non-locality or uncertainty of the shock position form and propagate away from the shock nose. Our results support the notion of upstream structures causing the patchwork reconstruction of the quasi-parallel shock front in a non-uniform manner. We propose a novel method for spacecraft data to be used to analyse this quasi-parallel reformation. We combine our hybrid-Vlasov results with test-particle studies and show that proton energization, which is required for injection, takes place throughout a larger shock transition zone. The energization of particles is found regardless of the instantaneous non-locality of the shock front, in agreement with it taking place over a larger region. Distortion of magnetic fields in front of and at the shock is shown to have a significant effect on proton injection. We additionally show that the density of suprathermal reflected particles upstream of the shock may not be a useful metric for the probability of injection at the shock, as foreshock dynamics and particle trapping appear to have a significant effect on energetic-particle accumulation at a given position in space. Our results have implications for statistical and spacecraft studies of the shock injection problem.Peer reviewe

Magnetosheath jet evolution as a function of lifetime : global hybrid-Vlasov simulations compared to MMS observations

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 Alfven Mach number M-A 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 M-A.Peer reviewe

Large-Scale Dune Aurora Event Investigation Combining Citizen Scientists' Photographs and Spacecraft Observations

Recently, citizen scientist photographs led to the discovery of a new auroral form called "the dune aurora" which exhibits parallel stripes of brighter emission in the green diffuse aurora at about 100 km altitude. This discovery raised several questions, such as (i) whether the dunes are associated with particle precipitation, (ii) whether their structure arises from spatial inhomogeneities in the precipitating fluxes or in the underlying neutral atmosphere, and (iii) whether they are the auroral manifestation of an atmospheric wave called a mesospheric bore. This study investigates a large-scale dune aurora event on 20 January 2016 above Northern Europe. The dunes were observed from Finland to Scotland, spanning over 1,500 km for at least 4 h. Spacecraft observations indicate that the dunes are associated with particle precipitation and reveal the presence of a temperature inversion layer below the mesopause during the event, creating suitable conditions for mesospheric bore formation. The analysis of a time lapse of pictures by a citizen scientist from Scotland leads to the estimate that, during this event, the dunes propagate toward the west-southwest direction at about 200 m s(-1), presumably indicating strong horizontal winds near the mesopause. These results show that citizen science and dune aurora studies can fill observational gaps and be powerful tools to investigate the least-known region of near-Earth space at altitudes near 100 km.Peer reviewe

Hybrid-Vlasov modelling of nightside auroral proton precipitation during southward interplanetary magnetic field conditions

Particle precipitation plays a key role in the coupling of the terrestrial magnetosphere and ionosphere by modifying the upper atmospheric conductivity and chemistry, driving field-aligned currents, and producing aurora. Yet quantitative observations of precipitating fluxes are limited, since ground-based instruments can only provide indirect measurements of precipitation, while particle telescopes aboard spacecraft merely enable point-like in situ observations with an inherently coarse time resolution above a given location. Further, orbit timescales generally prevent the analysis of whole events. On the other hand, global magnetospheric simulations can provide estimations of particle precipitation with a global view and higher time resolution. We present the first results of auroral (similar to 1-30 keV) proton precipitation estimation using the Vlasiator global hybrid-Vlasov model in a noon-midnight meridional plane simulation driven by steady solar wind with a southward interplanetary magnetic field. We first calculate the bounce loss-cone angle value at selected locations in the simulated nightside magnetosphere. Then, using the velocity distribution function representation of the proton population at those selected points, we study the population inside the loss cone. This enables the estimation of differential precipitating number fluxes as would be measured by a particle detector aboard a low-Earth-orbiting (LEO) spacecraft. The obtained differential flux values are in agreement with a well-established empirical model in the midnight sector, as are the integral energy flux and mean precipitating energy. We discuss the time evolution of the precipitation parameters derived in this manner in the global context of nightside magnetospheric activity in this simulation, and we find in particular that precipitation bursts ofPeer reviewe

A global view of Pc3 wave activity in near-Earth space : Results from hybrid-Vlasov simulations

Ultra-low frequency (ULF) waves in the Pc3 range, with periods between 10-45 s, are routinely observed in Earth's dayside magnetosphere. They are thought to originate in the foreshock, which extends upstream of the quasi-parallel bow shock and is populated with shock-reflected particles. The foreshock is permeated with ULF waves generated by ion beam instabilities, most notably the "30-s " waves whose periods match those of the Pc3 waves and which are carried earthward by the solar wind flow. However, the global picture of Pc3 wave activity from the foreshock to the magnetosphere and its response to changing solar wind conditions is still poorly understood. In this study, we investigate the global distribution and properties of Pc3 waves across near-Earth space using global simulations performed with the hybrid-Vlasov model Vlasiator. The simulations enable us to study the waves in their global context, and compare their properties in the foreshock, magnetosheath and dayside magnetosphere, for different sets of upstream solar wind conditions. We find that in all three regions the Pc3 wave power peaks at higher frequencies when the interplanetary magnetic field (IMF) strength is larger, consistent with previous studies. The Pc3 wave power is significantly enhanced in all three regions for higher solar wind Alfven Mach number. As this parameter is known to affect the shock properties but has little impact inside the magnetosphere, this brings further support to the magnetospheric waves originating in the foreshock. Other parameters that are found to influence the foreshock wave power are the solar wind density and the IMF cone angle. Inside the magnetosphere, the wave power distribution depends strongly on the IMF orientation, which controls the foreshock position upstream of the bow shock. The wave power is largest when the angle between the IMF and the Sun-Earth line is smallest, suggesting that wave generation and transmission are most efficient in these conditions.Peer reviewe

Connection Between Foreshock Structures and the Generation of Magnetosheath Jets : Vlasiator Results

Earth’s magnetosheath consists of shocked solar wind plasma that has been compressed and slowed down at the Earth’s bow shock. Magnetosheath jets are pulses of enhanced dynamic pressure in the magnetosheath. Jets have been observed by numerous spacecraft missions, but their origin has remained unconfirmed, though several formation mechanisms have been suggested. In this study, we use a method for automatically identifying and tracking jets as well as foreshock compressive structures (FCSs) in four 2D runs of the global hybrid-Vlasov simulation Vlasiator. We find that up to 75% of magnetosheath jets are caused by FCSs impacting the bow shock. These jets propagate deeper into the magnetosheath than the remaining 25% of jets that are not caused by FCSs. We conduct a visual case study of one jet that was not caused by FCSs and find that the bow shock was not rippled before the formation of the jet.Earth's magnetosheath consists of shocked solar wind plasma that has been compressed and slowed down at the Earth's bow shock. Magnetosheath jets are pulses of enhanced dynamic pressure in the magnetosheath. Jets have been observed by numerous spacecraft missions, but their origin has remained unconfirmed, though several formation mechanisms have been suggested. In this study, we use a method for automatically identifying and tracking jets as well as foreshock compressive structures (FCSs) in four 2D runs of the global hybrid-Vlasov simulation Vlasiator. We find that up to 75% of magnetosheath jets are caused by FCSs impacting the bow shock. These jets propagate deeper into the magnetosheath than the remaining 25% of jets that are not caused by FCSs. We conduct a visual case study of one jet that was not caused by FCSs and find that the bow shock was not rippled before the formation of the jet. Plain Language Summary The space around Earth is filled with plasma, the fourth state of matter. Earth's magnetic field shields our planet from the stream of plasma coming from the Sun, the solar wind. The solar wind plasma is slowed down at the Earth's bow shock, before it flows against and around the Earth's magnetic field in the magnetosheath. Sometimes, pulses of high density or velocity can occur in the magnetosheath that have the potential to disturb the inner regions of near-Earth space where many spacecraft orbit. We call these pulses magnetosheath jets. Magnetosheath jets have been observed by many spacecraft over the past few decades, but how they form has remained unclear. In this study, we use the Vlasiator model to simulate plasma in near-Earth space and investigate the origins of magnetosheath jets. We find that the formation of up to 75% of these jets can be explained by compressive structures in the foreshock, a region populated by intense wave activity extending sunward of the quasi-parallel bow shock, where interplanetary magnetic field lines allow shock-reflected particles to travel back toward the Sun.Peer reviewe

Spatial filtering in a 6D hybrid-Vlasov scheme to alleviate adaptive mesh refinement artifacts : a case study with Vlasiator (versions 5.0, 5.1, and 5.2.1)

Numerical simulation models that are used to investigate the near-Earth space plasma environment require sophisticated methods and algorithms as well as high computational power. Vlasiator 5.0 is a hybrid-Vlasov plasma simulation code that is able to perform 6D (3D in ordinary space and 3D in velocity space) simulations using adaptive mesh refinement (AMR). In this work, we describe a side effect of using AMR in Vlasiator 5.0: the heterologous grid approach creates discontinuities due to the different grid resolution levels. These discontinuities cause spurious oscillations in the electromagnetic fields that alter the global results. We present and test a spatial filtering operator for alleviating this artifact without significantly increasing the computational overhead. We demonstrate the operator's use case in large 6D AMR simulations and evaluate its performance with different implementations.Peer reviewe