Search

Geospace plasma simulations have progressed toward more realistic descriptions of the solar wind-magnetosphere interaction from magnetohydrodynamic to hybrid ion-kinetic, such as the state-of-the-art Vlasiator model. Despite computational advances, electron scales have been out of reach in a global setting. eVlasiator, a novel Vlasiator submodule, shows for the first time how electromagnetic fields driven by global hybrid-ion kinetics influence electrons, resulting in kinetic signatures. We analyze simulated electron distributions associated with reconnection sites and compare them with Magnetospheric Multiscale (MMS) spacecraft observations. Comparison with MMS shows that key electron features, such as reconnection inflows, heated outflows, flat-top distributions, and bidirectional streaming, are in remarkable agreement. Thus, we show that many reconnection-related features can be reproduced despite strongly truncated electron physics and an ion-scale spatial resolution. Ion-scale dynamics and ion-driven magnetic fields are shown to be significantly responsible for the environment that produces electron dynamics observed by spacecraft in near-Earth plasmas.

Radial diffusion coefficients quantify non-adiabatic transport of energetic particles by electromagnetic field fluctuations in planetary radiation belts. Theoretically, radial diffusion occurs for an ensemble of particles that experience irreversible violation of their third adiabatic invariant, which is equivalent to a change in their Roederer L* parameter. Thus, the Roederer L* coordinate is the fundamental quantity from which radial diffusion coefficients can be computed. In this study, we present a methodology to calculate the Lagrangian derivative of L* from global magnetospheric simulations, and test it with an application to Vlasiator, a hybrid-Vlasov model of near-Earth space. We use a Hamiltonian formalism for particles confined to closed drift shells with conserved first and second adiabatic invariants to compute changes in the guiding center drift paths due to electric and magnetic field fluctuations. We investigate the feasibility of this methodology by computing the time derivative of L* for an equatorial ultrarelativistic electron population travelling along four guiding center drift paths in the outer radiation belt during a 5 minute portion of a Vlasiator simulation. Radial diffusion in this simulation is primarily driven by ultralow frequency waves in the Pc3 range (10-45 s period range) that are generated in the foreshock and transmitted through the magnetopause to the outer radiation belt environment. Our results show that an alternative methodology to compute detailed radial diffusion transport is now available and could form the basis for comparison studies between numerical and observational measurements of radial transport in the Earth's radiation belts.

Numerical simulations have grown to play a central role in modern sciences over the years. The ever-improving technology of supercomputers has made large and precise models available. However, this accuracy is often limited by the cost of computational resources. Lowering the simulation's spatial resolution in order to conserve resources can lead to key processes being unresolved. We have shown in a previous study how insufficient spatial resolution of the proton cyclotron instability leads to a misrepresentation of ion dynamics in hybrid-Vlasov simulations. This leads to larger than expected temperature anisotropy and loss-cone shaped velocity distribution functions. In this study, we present a sub-grid numerical model to introduce pitch-angle diffusion in a 3D Cartesian velocity space, at a spatial resolution where the relevant wave-particle interactions were previously not correctly resolved. We show that the method is successfully able to isotropize loss-cone shaped velocity distribution functions, and that this method could be applied to simulations in order to save computational resources and still correctly model wave-particle interactions.

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-ref

We study the interaction of upstream ultralow frequency (ULF) waves with collisionless shocks by analyzing the outputs of 11 2D local hybrid simulation runs. Our simulated shocks have Alfvenic Mach numbers between 4.29 and 7.42 and their theta BN angles are 15 degrees, 30 degrees, 45 degrees, and 50 degrees. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's MA and theta BN. The wavelength positively correlates with both parameters, with the dependence on theta BN being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on MA and theta BN. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely their wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties (theta BN, strength), so that different sections of the same ULF wavefront encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However, some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra, these features are weaker.

Collisionless shocks are some of the most efficient particle accelerators in heliospheric and astrophysical plasmas. Here we study and quantify ion acceleration at Earth's bow shock with observations from NASA's Magnetospheric Multiscale (MMS) satellites and in a global hybrid-Vlasov simulation. From the MMS observations, we find that quasiparallel shocks are more efficient at accelerating ions. There, up to 15% of the available energy goes into accelerating ions above 10 times their initial energy. Above a shock-normal angle of ~50°, essentially no energetic ions are observed downstream of the shock. We find that ion acceleration efficiency is significantly lower when the shock has a low Mach number (MA < 6) while there is little Mach number dependence for higher values. We also find that ion acceleration is lower on the flanks of the bow shock than at the subsolar point regardless of the Mach number. The observations show that a higher connection time of an upstream field line leads to somewhat higher acceleration efficiency. To complement the observations, we perform a global hybrid-Vlasov simulation with realistic solar-wind parameters with the shape and size of the bow shock. We find that the ion acceleration efficiency in the simulation shows good quantitative agreement with the MMS observations. With the combined approach of direct spacecraft observations, we quantify ion acceleration in a wide range of shock angles and Mach numbers.

This paper describes properties and behavior of magnetic reconnection and flux transfer events (FTEs) on the dayside magnetopause using the global hybrid-Vlasov code Vlasiator. We investigate two simulation runs with and without a sunward (positive)B-x component of the interplanetary magnetic field (IMF) when the IMF is southward. The runs are two-dimensional in real space in the noon-midnight meridional (polar) plane and three-dimensional in velocity space. Solar wind input parameters are identical in the two simulations with the exception that the IMF is purely southward in one but tilted 45 degrees toward the Sun in the other. In the purely southward case (i.e., without B-x) the magnitude of the magnetos heath magnetic field component tangential to the magnetopause is larger than in the run with a sunward tilt. This is because the shock normal is perpendicular to the IMF at the equatorial plane, whereas in the other run the shock configuration is oblique and a smaller fraction of the total IMF strength is compressed at the shock crossing. Hence, the measured average and maximum reconnection rate are larger in the purely southward run. The run with tilted IMF also exhibits a north-south asymmetry in the tangential magnetic field caused by the different angle between the IMF and the bow shock normal north and south of the equator. Greater north-south asymmetries are seen in the FTE occurrence rate, size, and velocity as well; FTEs moving toward the Southern Hemisphere are larger in size and observed less frequently than FTEs in the Northern Hemisphere.

Today, the near‐Earth space is facing a paradigm change as the number of new spacecraft is literally skyrocketing. Increasing numbers of small satellites threaten the sustainable use of space, as without removal, space debris will eventually make certain critical orbits unusable. A central factor affecting small spacecraft health and leading to debris is the radiation environment, which is unpredictable due to an incomplete understanding of the near‐Earth radiation environment itself and its variability driven by the solar wind and outer magnetosphere. This paper presents the FORESAIL‐1 nanosatellite mission, having two scientific and one technological objectives. The first scientific objective is to measure the energy and flux of energetic particle loss to the atmosphere with a representative energy and pitch angle resolution over a wide range of magnetic local times. To pave the way to novel model‐in situ data comparisons, we also show preliminary results on precipitating electron fluxes obtained with the new global hybrid‐Vlasov simulation Vlasiator. The second scientific objective of the FORESAIL‐1 mission is to measure energetic neutral atoms of solar origin. The solar energetic neutral atom flux has the potential to contribute importantly to the knowledge of solar eruption energy budget estimations. The technological objective is to demonstrate a satellite deorbiting technology, and for the first time, make an orbit maneuver with a propellantless nanosatellite. FORESAIL‐1 will demonstrate the potential for nanosatellites to make important scientific contributions as well as promote the sustainable utilization of space by using a cost‐efficient deorbiting technology.

We use a global hybrid-Vlasov simulation for the magnetosphere, Vlasiator, to investigate magnetosheath high-speed jets. Unlike many other hybrid-kinetic simulations, Vlasiator includes an unscaled geomagnetic dipole, indicating that the simulation spatial and temporal dimensions can be given in SI units without scaling. Thus, for the first time, this allows investigating the magnetosheath jet properties and comparing them directly with the observed jets within the Earth's magnetosheath. In the run shown in this paper, the interplanetary magnetic field (IMF) cone angle is 30°, and a foreshock develops upstream of the quasi-parallel magnetosheath. We visually detect a structure with high dynamic pressure propagating from the bow shock through the magnetosheath. The structure is confirmed as a jet using three different criteria, which have been adopted in previous observational studies. We compare these criteria against the simulation results. We find that the magnetosheath jet is an elongated structure extending earthward from the bow shock by ∼ 2.6 RE, while its size perpendicular to the direction of propagation is ∼ 0.5R/E. We also investigate the jet evolution and find that the jet originates due to the interaction of the bow shock with a high-dynamic-pressure structure that reproduces observational features associated with a short, large-amplitude magnetic structure (SLAMS). The simulation shows that magnetosheath jets can develop also under steady IMF, as inferred by observational studies. To our knowledge, this paper therefore shows the first global kinetic simulation of a magnetosheath jet, which is in accordance with three observational jet criteria and is caused by a SLAMS advecting towards the bow shock., Export Date: 22 October 2018; Article; Correspondence Address: Palmroth, M.; Department of Physics, University of HelsinkiFinland; email: minna.palmroth@helsinki.fi; Funding details: IA104416; Funding details: 704681; Funding details: CSC, China Scholarship Council; Funding details: NNX17AI45G; Funding details: 309937; Funding details: 267144; Funding details: 312351; Funding details: 682068-PRESTISSIMO, ERC, European Research Council; Funding details: 200141-QuESpace, ERC, European Research Council; Funding text: Acknowledgements. We acknowledge the European Research Council for Starting grant 200141-QuESpace, with which Vlasiator (http://helsinki.fi/vlasiator; last access: 4 September 2018) was developed, and Consolidator grant 682068-PRESTISSIMO awarded to further develop Vlasiator and use it for scientific investigations. We gratefully also acknowledge the Finnish Centre of Excellence in Research of Sustainable Space (Academy of Finland grant numbers 312351, 267144, and 309937). Primož Kajdicˇ’s work was supported by DGAPA/PAPIIT grant IA104416. The CSC – IT Center for Science in Finland is acknowledged for the Grand Challenge award leading to the results shown in here. We acknowledge valuable discussions within the International Space Science Institute (ISSI) team 350, called “Jets downstream of collisionless shocks”, led by Ferdinand Plaschke and Heli Hietala. Lucile Turc acknowledges Marie Sklodowska-Curie grant 704681. Heli Hietala was supported by the Turku Collegium for Science and Medicine and NASA NNX17AI45G. We thank Jonas Suni for producing data for the figures in the revised version. QC 20181112