Geoelectric Field Caused by Flux Transfer Events in an Ionosphere-Coupled Vlasiator Simulation

We report on the relationship between flux transfer events (FTEs) at Earth’s magnetopause and the geoelectric field that is induced near the FTEs’ magnetic footpoints. We study this system using the global hybrid-Vlasov code Vlasiator, which has recently been extended to model ionospheric physics. We also highlight the significance of 3D magnetic null points, which in our simulation can separate the FTEs into multiple flux ropes. Near the null points, the coiled FTE magnetic field lines are rerouted towards Earth, so that the magnetic footpoints are planted near the Region 1 ionospheric current system. The helicities of the flux ropes are organized by the y-component (GSE) of the magnetic field at the Earth’s magnetopause. This occurs in our simulation due to the absence of a y-component of the interplanetary magnetic field, which would normally supply the FTE guide field that determines the helicity. We observe Alfvenic, Earthward-flowing field-aligned currents generated near the magnetopause that correlate with the passage of FTEs nearby. These pulses of current coincide with the formation of rotational geoelectric field structures, that appear near the noon meridian and propagate around the auroral oval towards the nightside.

💡 Research Summary

This paper investigates how flux transfer events (FTEs) at Earth’s dayside magnetopause generate geoelectric fields that can induce geomagnetically induced currents (GICs) on the ground. The authors employ the global hybrid‑Vlasov code Vlasiator, which has recently been extended to include an electrostatic ionospheric model, thereby allowing a self‑consistent coupling between the magnetosphere and ionosphere. The simulation domain spans from –110 R_E to +50 R_E in the Sun‑Earth direction and ±58 R_E in the transverse directions, with adaptive mesh refinement achieving a finest cell size of 1000 km. A constant solar‑wind inflow (density 10 cm⁻³, velocity –750 km s⁻¹, temperature 5 × 10⁵ K) and a southward IMF (B = ‑5 nT) drive persistent dayside reconnection.

A key methodological advance is the identification of FTEs via three‑dimensional magnetic topology. The authors locate O‑type magnetic null lines (the cores of flux ropes) and X‑type null lines (reconnection sites) using the algorithm described in Alho et al. (2024). This enables the detection of multiple flux ropes within a single FTE, separated by three‑dimensional magnetic null points (nulls). Near these nulls, the coiled magnetic field lines are rerouted earthward, planting the footpoints close to the Region 1 ionospheric current system.

Because the simulation deliberately omits an IMF y‑component, the usual guide field that would set rope helicity is absent. Consequently, the helicities of the individual flux ropes become organized by the GSE y‑component of the magnetopause magnetic field, a novel result that confirms the theoretical expectation that IMF y controls rope twist.

The authors track field‑aligned currents (FACs) generated at the magnetopause. As an FTE passes, an Alfvénic pulse of FAC propagates earthward at the Alfvén speed, reaching the ionosphere after a modeled 4 s delay. The incoming FAC intensifies ionospheric Pedersen and Hall currents, producing a rotational geoelectric field structure. This structure first appears near the noon meridian and then drifts azimuthally around the auroral oval toward the nightside, closely resembling the “traveling convection vortices” reported in earlier MHD studies. The simulated geoelectric fields are derived from the time‑derivative of the magnetic perturbations, using a 1 s output cadence that is sufficient to resolve the high‑frequency components relevant for GIC generation.

The study highlights several implications for space‑weather forecasting. First, it demonstrates a direct causal chain from magnetopause reconnection (FTE) to ground‑level electric fields, mediated by Alfvénic FAC pulses and ionospheric current reconfiguration. Second, the presence of multiple flux ropes within an FTE and their helicity ordering suggest that the spatial pattern of ground‑induced electric fields can be more complex than the simple bipolar signatures traditionally assumed. Third, the rotational electric field patches could generate localized, high‑amplitude GICs even when the overall geomagnetic activity is modest.

Limitations are acknowledged. The ionospheric model uses a simplified conductance prescription based on precipitating ion moments, and the ground conductivity is represented by a uniform, isotropic layer. Realistic regional variations in crustal conductivity, as well as the contribution of the auroral electrojet system, would need to be incorporated for quantitative GIC predictions. Moreover, the absence of an IMF y‑component, while useful for isolating the helicity mechanism, does not capture the full range of solar‑wind conditions encountered in nature.

In conclusion, the ionosphere‑coupled Vlasiator simulation provides the first comprehensive, self‑consistent demonstration that FTEs can launch Alfvénic FAC pulses, reorganize ionospheric currents, and generate rotating geoelectric fields that propagate around the auroral oval. The work bridges a gap between magnetospheric reconnection physics and ground‑level space‑weather impacts, and it sets the stage for future studies that will integrate realistic IMF variability, high‑resolution ground conductivity models, and direct comparisons with simultaneous satellite‑ground observation campaigns.