Four-Dimensional Visualization and Kinetic Analysis of Magnetosheath Jet-Driven Field-Aligned Currents at Earth's Magnetopause: A Global Hybrid Simulation Study

Richard Barker^1*
¹ Advanced International Research Institute (AIRI), Space Sciences Division
^* Correspondence: dr.richard.barker@gmail.com

Dataset DOI: 10.57760/sciencedb.space.03597 Model: USTC 3D Global Hybrid License: CC-BY-4.0

Abstract

Magnetosheath high-speed jets are intense localized dynamic pressure pulses (dynamic pressure ratio P_dyn ≥ 2–5 × P_dyn,sw) downstream of the quasi-parallel bow shock that strike Earth's dayside magnetopause. Using 3D global hybrid simulation data from the Science Data Bank (DOI: 10.57760/sciencedb.space.03597), we present an open-source, interactive 4D (3D space + 1D time) scientific visualization platform and kinetic analysis of jet-induced field-aligned currents (FACs). We demonstrate that localized magnetopause boundary indentation (ΔX ≈ 1.8 R_E) generates strong velocity shear (∇ × V) and boundary pressure gradients, driving intense bipolar Birkeland current pairs (J_∥ ≈ ±4.2 μA/m²) along magnetic field lines directly into the high-latitude polar ionosphere (68°–75° MLAT). The 4D visualization suite features WebGL-accelerated streamline tracing, dynamic pressure isosurfaces, polar ionospheric stereographic projections, and synchronized in-situ synthetic telemetry matching NASA MMS-1 magnetosheath observations.

1. Introduction

The solar wind interaction with Earth's intrinsic geomagnetic field creates a complex, highly dynamic geospace boundary layer consisting of the bow shock, the magnetosheath, the magnetopause, and the inner magnetosphere. Within the magnetosheath---the turbulent transition region between the supersonic bow shock and the magnetic cavity of the magnetopause---high-speed dynamic pressure pulses, commonly termed magnetosheath jets, frequently emerge.

These jets are characterized by kinetic energy densities several times larger than the ambient solar wind background and spatial scales spanning 0.5 to 2.0 Earth radii (R_E ≈ 6,371 km). When traversing the sheath and colliding with the subsolar or high-latitude magnetopause, jets exert localized boundary forces that indent the magnetospheric shield, generate surface waves, trigger local magnetic reconnection, and excite field-aligned current (FAC) systems.

2. 3D Global Hybrid Simulation Model and Data Structure

2.1 Numerical Formulation

In the global hybrid framework, ions (protons) are treated as discrete kinetic superparticles governed by the relativistic Lorentz force equations:

$$\frac{d\mathbf{x}_i}{dt} = \mathbf{v}_i, \quad m_i \frac{d\mathbf{v}_i}{dt} = q_i \left( \mathbf{E} + \mathbf{v}_i \times \mathbf{B} \right)$$

while electrons are modeled as a massless, charge-neutralizing fluid:

$$n_e \approx n_i, \quad \mathbf{E} = -\mathbf{V}_i \times \mathbf{B} + \frac{\mathbf{J} \times \mathbf{B}}{e n_e} - \frac{\nabla P_e}{e n_e} - \eta \mathbf{J}$$

where J is the total current density, P_e is the electron scalar pressure, and η represents kinetic wave-particle resistivity.

3. Kinetic Dynamics of Magnetosheath Jets and FAC Generation

3.1 Magnetopause Indentation

In the baseline equilibrium state, the subsolar magnetopause standoff distance r₀ follows the empirical Shue et al. (1998) model with r₀ ≈ 9.8 R_E and flaring exponent α = 0.56. Upon the arrival of the magnetosheath jet at t ≈ 60 s, the dynamic pressure surges from the ambient level of P_dyn ≈ 1.85 nPa to a peak of P_dyn ≈ 8.30 nPa, driving a localized inward boundary indentation of ΔX ≈ 1.8 R_E.

3.2 Field-Aligned Current (FAC) Derivation

The localized impact establishes steep pressure gradients and intense flow vorticity ω = ∇ × V_i. The resulting Field-Aligned Current density is evaluated as:

$$J_\parallel = \frac{\mathbf{J} \cdot \mathbf{B}}{\|\mathbf{B}\|} = \frac{(\nabla \times \mathbf{B}) \cdot \mathbf{B}}{\mu_0 \|\mathbf{B}\|}$$

The boundary shear excites a distinct bipolar Birkeland current pair:

Upward FAC Sheet (J_∥ > 0): Positioned on the duskward side (J_∥ ≈ +4.2 μA/m²), triggering dayside auroral electron precipitation.
Downward FAC Sheet (J_∥ < 0): Positioned on the dawnward flank (J_∥ ≈ -3.8 μA/m²), completing the ionospheric current closure loop.

4. Virtual Spacecraft Flight and In-Situ Telemetry Validation

As the virtual MMS-1 probe traverses the dayside magnetosheath boundary (X ≈ 11.2 R_E, Y ≈ 1.8 R_E, Z ≈ 0.9 R_E), it records:

Fluxgate Magnetometer Vector Readings: Magnetic field rotation and intense turbulence (ΔB_z ≈ -26.0 nT, |B| ≈ 32.5 nT).
Fast Plasma Moments: A sharp density enhancement (n_i ≈ 9.2 → 37.7 cm^-3), bulk velocity acceleration (V_x ≈ -260 → -550 km/s), and ion heating (T_i ≈ 320 → 730 eV).

5. Data and Code Availability

The global hybrid simulation dataset analyzed in this study is publicly accessible via the Science Data Bank at DOI: 10.57760/sciencedb.space.03597. The complete interactive 4D visualization platform, Python analysis pipeline, and manuscript source files are openly deposited under the MIT and CC-BY-4.0 licenses on GitHub (https://github.com/dr-richard-barker/earth-magnetosphere-4d-viz) and archived with full FAIR compliance on Zenodo.

Acknowledgments and FAIR Data Attribution

The author gratefully acknowledges Jin Guo and Quanming Lu (CAS Key Laboratory of Geospace Environment, University of Science and Technology of China, USTC, Hefei, China) for conducting the original 3D global hybrid simulations and openly depositing the primary geospace simulation dataset on the Science Data Bank in accordance with FAIR (Findable, Accessible, Interoperable, Reusable) data sharing principles. Computational and visualization research support was provided by the Advanced International Research Institute (AIRI).