Hybrid-Vlasov modeling of three-dimensional dayside magnetopause reconnection

Dayside magnetic reconnection at the magnetopause, which is a major driver of space weather, is studied for the first time in a three-dimensional (3D) realistic setup using a hybrid-Vlasov kinetic model. A noon-midnight meridional plane simulation is extended in the dawn-dusk direction to cover 7 Earth radii. The southward interplanetary magnetic field causes magnetic reconnection to occur at the subsolar magnetopause. Perturbations arising from kinetic instabilities in the magnetosheath appear to modulate the reconnection. Its characteristics are consistent with multiple, bursty, and patchy magnetopause reconnection. It is shown that the kinetic behavior of the plasma, as simulated by the model, has consequences on the applicability of methods such as the four-field junction to identify and analyze magnetic reconnection in 3D kinetic simulations.Peer reviewe

A possible source mechanism for magnetotail current sheet flapping

The origin of the flapping motions of the current sheet in the Earth's magnetotail is one of the most interesting questions of magnetospheric dynamics yet to be solved. We have used a polar plane simulation from the global hybrid-Vlasov model Vlasiator to study the characteristics and source of current sheet flapping in the center of the magnetotail. The characteristics of the simulated signatures agree with observations reported in the literature. The flapping is initiated by a hemispherically asymmetric magnetopause perturbation, created by subsolar magnetopause reconnection, that is capable of displacing the tail current sheet from its nominal position. The current sheet displacement propagates downtail at the same pace as the driving magnetopause perturbation. The initial current sheet displacement launches a standing magnetosonic wave within the tail resonance cavity. The travel time of the wave within the local cavity determines the period of the subsequent flapping signatures. Compression of the tail lobes due to added flux affects the cross-sectional width of the resonance cavity as well as the magnetosonic speed within the cavity. These in turn modify the wave travel time and flapping period. The compression of the resonance cavity may also provide additional energy to the standing wave, which may lead to strengthening of the flapping signature. It may be possible that the suggested mechanism could act as a source of kink-like waves that have been observed to be emitted from the center of the tail and to propagate toward the dawn and dusk flanks.Peer reviewe

Ion Acceleration by Flux Transfer Events in the Terrestrial Magnetosheath

We report ion acceleration by flux transfer events in the terrestrial magnetosheath in a global two-dimensional hybrid-Vlasov polar plane simulation of Earth's solar wind interaction. In the model we find that propagating flux transfer events created in magnetic reconnection at the dayside magnetopause drive fast-mode bow waves in the magnetosheath, which accelerate ions in the shocked solar wind flow. The acceleration at the bow waves is caused by a shock drift-like acceleration process under stationary solar wind and interplanetary magnetic field upstream conditions. Thus, the energization is not externally driven but results from plasma dynamics within the magnetosheath. Energetic proton populations reach the energy of 30 keV, and their velocity distributions resemble time-energy dispersive ion injections observed by the Cluster spacecraft in the magnetosheath.Peer reviewe

Tail reconnection in the global magnetospheric context : Vlasiator first results

The key dynamics of the magnetotail have been researched for decades and have been associated with either three-dimensional (3-D) plasma instabilities and/or magnetic reconnection. We apply a global hybrid-Vlasov code, Vlasiator, to simulate reconnection self-consistently in the ion kinetic scales in the noon-midnight meridional plane, including both dayside and nightside reconnection regions within the same simulation box. Our simulation represents a numerical experiment, which turns off the 3-D instabilities but models ion-scale reconnection physically accurately in 2-D. We demonstrate that many known tail dynamics are present in the simulation without a full description of 3-D instabilities or without the detailed description of the electrons. While multiple reconnection sites can coexist in the plasma sheet, one reconnection point can start a global reconfiguration process, in which magnetic field lines become detached and a plasmoid is released. As the simulation run features temporally steady solar wind input, this global reconfiguration is not associated with sudden changes in the solar wind. Further, we show that lobe density variations originating from dayside reconnection may play an important role in stabilising tail reconnection.Peer reviewe

Evidence for transient, local ion foreshocks caused by dayside magnetopause reconnection

We present a scenario resulting in time-dependent behaviour of the bow shock and transient, local ion reflection under unchanging solar wind conditions. Dayside magnetopause reconnection produces flux transfer events driving fast-mode wave fronts in the magnetosheath. These fronts push out the bow shock surface due to their increased downstream pressure. The resulting bow shock deformations lead to a configuration favourable to localized ion reflection and thus the formation of transient, travelling foreshock-like field-aligned ion beams. This is identified in two-dimensional global magnetospheric hybrid-Vlasov simulations of the Earth's magnetosphere performed using the Vlasiator model (http://vlasiator.fmi.fi). We also present observational data showing the occurrence of dayside reconnection and flux transfer events at the same time as Geotail observations of transient foreshock-like field-aligned ion beams. The spacecraft is located well upstream of the fore-shock edge and the bow shock, during a steady southward interplanetary magnetic field and in the absence of any solar wind or interplanetary magnetic field perturbations. This indicates the formation of such localized ion foreshocks

Hybrid-Vlasov simulation of soft X-ray emissions at the Earth’s dayside magnetospheric boundaries

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

ULF Wave Transmission Across Collisionless Shocks : 2.5D Local Hybrid Simulations

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

Ion Acceleration by Flux Transfer Events in the Terrestrial Magnetosheath

We report ion acceleration by flux transfer events in the terrestrial magnetosheath in a global two‐dimensional hybrid‐Vlasov polar plane simulation of Earth's solar wind interaction. In the model we find that propagating flux transfer events created in magnetic reconnection at the dayside magnetopause drive fast‐mode bow waves in the magnetosheath, which accelerate ions in the shocked solar wind flow. The acceleration at the bow waves is caused by a shock drift‐like acceleration process under stationary solar wind and interplanetary magnetic field upstream conditions. Thus, the energization is not externally driven but results from plasma dynamics within the magnetosheath. Energetic proton populations reach the energy of 30 keV, and their velocity distributions resemble time‐energy dispersive ion injections observed by the Cluster spacecraft in the magnetosheath. </p

Ion Acceleration Efficiency at the Earth's Bow Shock : Observations and Simulation Results

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 similar to 50 degrees, 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 (M ( A ) < 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.Peer reviewe

Spacecraft measurements constraining the spatial extent of a magnetopause reconnection X line

Multispacecraft measurements from the Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission are used to probe the spatial extent of an X line at the dayside magnetopause. A case study from 21 April 2014 is presented where two THEMIS spacecraft have a near-simultaneous encounter with the equatorial dayside magnetopause separated by 3.9 Earth radii. Both spacecraft observe similar steady inflow conditions with southward interplanetary magnetic field and a high magnetic shear angle at the magnetopause (133 degrees) boundary. One spacecraft observes clear fluid and kinetic signatures of active magnetic reconnection, while the other spacecraft does not observe reconnection. The predicted location of reconnection across the magnetopause is found using several theoretical models and a Block Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) MHD simulation. Each model predicts a continuous X line passing close to the two spacecraft, suggesting both would observe reconnection, if active. Using the constraints of the multipoint measurements, the extent or length L of the reconnection is estimated to be 2.4L <5.2h in local time or 6L<14 R-E.Peer reviewe