Unveiling the Mixing and Transport Processes of Solar Wind and Planetary Ions in the Magnetopause Boundary Layer

TL;DR

This study probes KH-induced mixing and cross-boundary transport at the magnetopause using 3D global hybrid simulations with Earth as a representative magnetized planet. It introduces a kinetic framework and a mixing-rate metric MR to automatically delineate the KH-modulated magnetopause and to quantify solar wind entry and magnetospheric ion escape across a range of solar wind dynamic pressures and IMF clock angles. Key findings show that northward IMF concentrates transport near the equator with increasing flux as $P_d$ rises, while southward IMF, via reconnection coupled with KH, enhances solar wind injection and planetary escape near subsolar regions; By also modulates KH geometry and flux patterns. The results provide quantitative insight into cross-boundary mass and energy transfer, with implications for space weather, planetary atmospheres, and the evolution of magnetized plasmas in the heliosphere and beyond.

Abstract

Kelvin-Helmholtz (KH) vortices are widely observed in astrophysics and heliophysics, including at Jovian and terrestrial magnetopauses, the Martian sheath-ionosphere boundary, the heliopause, and within stellar accretion disks. These vortices play a critical role in transporting mass, momentum, and energy across boundary layers. Magnetized planets such as Earth exhibit a higher incidence of fully rolled-up, nonlinear KH vortices compared to non-magnetized planets like Mars. In contrast to previous magnetohydrodynamic (MHD) studies, this work adopts a kinetic point of view to quantify ion mixing rates using three-dimensional global hybrid simulations, with Earth as a representative case. This approach enables automated identification of the KH-modulated, corrugated magnetopause. For the first time, we provide a quantitative assessment of how solar wind conditions control solar wind entry and subsequent mixing with magnetospheric ions via KH waves. We find that under northward interplanetary magnetic field (IMF) conditions, the flux of particles crossing the dayside magnetopause increases with solar wind dynamic pressure and peaks in the KH region. Notably, the KH-modulated low-latitude boundary layer thins as the dynamic pressure increases. Under southward IMF conditions, coupled reconnection and KH structures further enhance solar wind injection and boost magnetospheric ion escape in the dayside, especially near the subsolar point where reconnection intensifies this exchange. These results also shed light on the evolution of space environments and mass transport at magnetized planets in the heliosphere and beyond.

Figures (8)

• Figure 1: Magnetopause Kelvin-Helmholtz waves under a northward IMF condition. (a) the magnetic field $\mathrm{B}$ ($\mathrm{nT}$), (b) the ion bulk velocity $V_\mathrm{x}$ ($\mathrm{km/s}$), (c) and (d) number densities of solar wind ions $N_\mathrm{sw}$ ($\mathrm{cm^{-3}}$) and magnetospheric ions $N_\mathrm{mp}$ ($\mathrm{cm^{-3}}$), respectively.
• Figure 2: The influence of solar wind dynamic pressure on ion mixing rates under the same northward IMF condition. The mixing rates are represented in the equatorial plane, where K-H fluctuations are believed to be most intense. From left to right, panels (a-e) demonstrate the trend in mixing rates at the magnetopause boundary layer under conditions corresponding to increasing dynamic pressures ($P$ from $<0.5\ \mathrm{nPa}$ to $>10\ \mathrm{nPa}$).
• Figure 3: The corresponding influence of solar wind dynamic pressure on the fluxes of entering solar wind ions (blue) and escaping magnetospheric ions (red) is shown on the magnetopause. The three-dimensional magnetopause surface is obtained from the peaks in mixing rates as depicted in Figure 2.
• Figure 4: Quantitative statistical results for the fluxes of entering solar wind ions (a) and escaping magnetospheric ions (b) traversing the magnetopause under different solar wind dynamic pressure conditions. The illustration features five concentric circles, each representing a distinct scenario, arranged from the innermost to the outermost layers, with increasing dynamic pressures ranging from $<0.5~\mathrm{nPa}$ to $>10~\mathrm{nPa}$, as indicated by red arrows. For each concentric circle, average flux measurements in six magnetopause sectors defined by north-south and dawn-dusk orientations are annotated within their respective sectors.
• Figure 5: The impact of IMF clock angle on ion mixing rates. The calculation of mixing rates and the statistics of particle flux are analogous to Figure 2. The perspective of these diagrams is from the Sun looking towards the Earth, with the horizontal axis representing the dawn-dusk direction and the vertical axis indicating geomagnetic north-south poles. Black clocks with red arrows indicate the corresponding IMF clock angle conditions for each case.