Simulations of solar wind - magnetosheath - magnetopause interactions

This thesis investigates interactions between solar wind and the magnetosphere of the Earth using two global magnetosphericsimulation models, GUMICS-4 and Vlasiator, which are both developed in Finland. The main topic of the thesis is magnetic reconnection at the dayside magnetopause, its drivers and global effects. Magnetosheath mirror mode waves and their evolution, identification and impacts on the local reconnection rates at the magnetopause are also discussed. This thesis consists of four peer-reviewed papers and an introductory part. GUMICS-4 is a magnetohydrodynamic model solving plasma as a single magnetized fluid. Vlasiator is the world’s first global magnetospheric hybrid-Vlasov simulation model, which solves the motion of ions by describing them as velocity distribution functions, whereas electrons are described as a charge neutralizing fluid. Vlasiator is able to solve ion scale physics in a global scale simulation. However, it is computationally heavy and the global simulations are currently describing Earth’s magnetosphere only in two spatial dimensions, whereas the velocity space is three dimensional. This thesis shows that magnetic reconnection at the dayside magnetopause is controlled by several factors. The impact of dipole tilt angle and sunward component of the interplanetary magnetic field on magnetopause reconnection is investigated with a set of GUMICS-4 simulations. Using Vlasiator simulations, this thesis shows that local reconnection rate is highly variable even during steady solar wind and correlates well with an analytical model for 2D asymmetric reconnection. It is also shown that the local reconnection rate is affected by local variations in the magnetosheath plasma. Fluctuations in the magnetosheath parameters near X-lines are partly generated by mirror mode waves that are observed to grow in the quasi-perpendicular magnetosheath. These results show that that the local reconnection rate at the X-lines is affected not only by the fluctuations in the inflow parameters but also by reconnection at nearby X-lines. Outflow from stronger X-lines pushes against the weaker ones and might ultimately suppress reconnection in the weaker X-lines. Magnetic islands, 2D representations of FTEs, form between X-lines in the Vlasiator simulations. FTEs propagate along the dayside magnetopause driving bow waves in the magnetosheath. The bow waves propagate upstream all the way to the bow shock causing bulges in the shock, from which solar wind particles can reflect back to the solar wind causing local foreshocks. The overall conclusion of this thesis is that the ion scale kinetic physics is important to accurately model the solar wind – magnetosheath – magnetopause interactions. Vlasiator results show a strong scale-coupling between ion and global scales: global scale phenomena have an impact on the local physics and the local phenomena may have unexpected impacts on the global dynamics of the magnetosphere. Neglecting the global scales in local ion scale simulations and vice versa may therefore lead to incomplete description of the solar wind – magnetosphere interactions

On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks

In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.Peer reviewe

Using mutual information to investigate non-linear correlation between AE index, ULF Pc5 wave activity and electron precipitation

In this study, we use mutual information from information theory to investigate non-linear correlation between geomagnetic activity indicated by auroral electrojet (AE) index with both the global ultra low frequency (ULF) Pc5 wave power and medium energy (>= 30 keV) electron precipitation at the central outer radiation belt. To investigate the energy and magnetic local time (MLT) dependence of the non-linearity, we calculate the mutual information and Pearson correlation coefficient separately for three different energy ranges (30-100 keV, 100-300 keV and >= 300 keV) and four different MLT sectors (0-6, 6-12, 12-18, 18-24). We compare results from 2 years 2004 and 2007 representing geomagnetically more active and less active years, respectively. The correlation analysis between the AE index and electron precipitation shows a clear MLT and energy dependence in both active and quiet conditions. In the two lowest energy ranges of the medium energy electrons (30-100 keV and 100-300 keV) both non-linear correlation and Pearson correlation indicate strong dependence with the AE index in the dawn sector. The linear dependence indicated by the Pearson correlation coefficient decreases from dawn to dusk while the change in the non-linear correlation is smaller indicating an increase in the non-linearity from dawn to dusk. The non-linearity between the AE index and electron precipitation is larger at all MLT sectors except MLTs 6-12 during geomagnetically more active year when larger amount of the activity is driven by interplanetary coronal mass ejections (ICMEs) compared to lower activity year with high speed stream (HSS) and stream interaction region (SIR) driven activity. These results indicate that the processes leading to electron precipitation become more non-linear in the dusk and during geomagnetically more active times when the activity is driven by ICMEs. The non-linearity between the AE index and global ULF Pc5 activity is relatively low and seems not to be affected by the difference in the geomagnetic activity during the 2 years studied.Peer reviewe

Reconnection rates and X line motion at the magnetopause : Global 2D-3V hybrid-Vlasov simulation results

We present results from a first study of the local reconnection rate and reconnection site motion in a 2D-3V global magnetospheric self-consistent hybrid-Vlasov simulation with due southward interplanetary magnetic field. We observe magnetic reconnection at multiple locations at the dayside magnetopause and the existence of magnetic islands, which are the 2-D representations of flux transfer events. The reconnection locations (the X lines) propagate over significant distances along the magnetopause, and reconnection does not reach a steady state. We calculate the reconnection rate at the location of the X lines and find a good correlation with an analytical model of local 2-D asymmetric reconnection. We find that despite the solar wind conditions being constant, the reconnection rate and location of the X lines are highly variable. These variations are caused by magnetosheath fluctuations, the effects of neighboring X lines, and the motion of passing magnetic islands.Peer reviewe

On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks

In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role

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

Observatories of the Solar Corona and Active Regions (OSCAR)

Coronal Mass Ejections (CMEs) and Corotating Interaction Regions (CIRs) are major sources of magnetic storms on Earth and are therefore considered to be the most dangerous space weather events. The Observatories of Solar Corona and Active Regions (OSCAR) mission is designed to identify the 3D structure of coronal loops and to study the trigger mechanisms of CMEs in solar Active Regions (ARs) as well as their evolution and propagation processes in the inner heliosphere. It also aims to provide monitoring and forecasting of geo- effective CMEs and CIRs. OSCAR would contribute to significant advancements in the field of solar physics, improvements of the current CME prediction models, and provide data for reliable space weather forecasting. These objectives are achieved by utilising two spacecraft with identical instrumentation, located at a heliocentric orbital distance of 1 AU from the Sun. The spacecraft will be separated by an angle of 68° to provide optimum stereoscopic view of the solar corona. We study the feasibility of such a mission and propose a preliminary design for OSCAR

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