Kinetic-scale magnetic turbulence and finite Larmor radius effects at Mercury

We use a nonstationary generalization of the higher-order structure function technique to investigate statistical properties of the magnetic field fluctuations recorded by MESSENGER spacecraft during its first flyby (01/14/2008) through the near Mercury’s space environment, with the emphasis on key boundary regions participating in the solar wind – magnetosphere interaction. Our analysis shows, for the first time, that kinetic-scale fluctuations play a significant role in the Mercury’s magnetosphere up to the largest resolvable time scale ~20 s imposed by the signal nonstationarity, suggesting that turbulence at this planet is largely controlled by finite Larmor radius effects. In particular, we report the presence of a highly turbulent and extended foreshock system filled with packets of ULF oscillations, broad-band intermittent fluctuations in the magnetosheath, ion-kinetic turbulence in the central plasma sheet of Mercury’s magnetotail, and kinetic-scale fluctuations in the inner current sheet encountered at the outbound (dawn-side) magnetopause. Overall, our measurements indicate that the Hermean magnetosphere, as well as the surrounding region, are strongly affected by non-MHD effects introduced by finite sizes of cyclotron orbits of the constituting ion species. Physical mechanisms of these effects and their potentially critical impact on the structure and dynamics of Mercury’s magnetic field remain to be understood.

💡 Research Summary

This paper presents a comprehensive analysis of magnetic field fluctuations measured by the MESSENGER spacecraft during its first flyby of Mercury on 14 January 2008, with a focus on kinetic‑scale turbulence and finite Larmor radius (FLR) effects. The authors employ a non‑stationary generalization of the higher‑order structure function (SF) technique, which is well suited for short, trend‑dominated data sets where traditional Fourier methods struggle. By calculating the q‑order structure functions Sₙ(τ)=⟨|δB(τ)|ⁿ⟩ over a range of time lags τ and extracting scaling exponents ζₙ from the power‑law relation Sₙ(τ)∝τ^{ζₙ}, they obtain a direct measure of the turbulence cascade. The second‑order exponent ζ₂ is linked to the spectral index β of the corresponding power spectrum via ζ₂=β−1 under a linear space‑time coupling assumption.

A key methodological innovation is the use of a continuous SF scalogram, which visualizes how scaling exponents evolve in time, allowing the authors to pinpoint transitions between distinct turbulence regimes even in the presence of strong non‑stationarity. The mapping from temporal scales τ to spatial wave numbers k is performed using the Taylor frozen‑flow hypothesis (k≈2π/(V₀τ)), where V₀ is the bulk plasma flow speed obtained from concurrent plasma measurements. This enables the identification of the ion‑gyro‑radius scale (kρᵢ≈1) and the associated crossover time τᵢ. From τᵢ, the ion gyro‑radius ρᵢ≈V₀τᵢ/2π and ion temperature Tᵢ≈mᵢ(V₀τᵢ/τ_{ci})² are estimated, providing physical parameters that are otherwise difficult to obtain in situ.

The flyby trajectory is divided into ten characteristic regions: unperturbed solar wind (SW1, SW2), four foreshock sub‑regions (FS1‑FS4), four magnetosheath sub‑regions (MS1‑MS4), the cross‑tail current sheet (CCS), the ion boundary layer (IBL), and a diamagnetic decrease (DD). In each region the authors compute SF exponents up to sixth order, revealing three generic scaling domains:

1. MHD regime (Region I) – Large scales exhibit β≈5/3 (ζ₂≈2/3), consistent with Kolmogorov‑type Alfvénic turbulence.

2. Ion‑kinetic regime (Region II) – At τ≈τᵢ (≈0.8–1.5 s depending on location) the spectra steepen to β≈2.3–2.5 (ζ₂≈1.3–1.5). This steepening reflects the breakdown of the fluid description as the ion gyro‑radius becomes comparable to the fluctuation scale, and FLR effects dominate the cascade.

3. Electron‑kinetic regime (Region III) – At still smaller scales (kρₑ≈1) the spectra would become even steeper, but the limited temporal resolution and non‑stationarity prevent a robust detection in the present data set.

Across all regions, higher‑order exponents ζₙ (n>2) display a positive dependence on n (∂ζₙ/∂n>0), indicating strong intermittency (“spikiness”) of the fluctuations. In some foreshock intervals flat ζₙ≈0 curves appear, signifying the presence of coherent structures such as shocks or wave packets embedded in the stochastic background.

Specific findings for each region include:

• Foreshock (FS1‑FS4): A highly turbulent environment filled with ultra‑low‑frequency (ULF) wave packets. The SF analysis reveals intermittent bursts superimposed on a background cascade, with ion gyro‑radii of 30–50 km and ion temperatures of 200–400 eV.

• Magnetosheath (MS1‑MS4): Consistent ion‑kinetic scaling (β≈2.0–2.3) throughout, with pronounced intermittency suggesting the operation of mirror, ion‑cyclotron, or other kinetic instabilities.

• Cross‑tail current sheet (CCS): The strongest ion‑kinetic signature is observed, τᵢ≈1.5 s, ρᵢ≈45 km, Tᵢ≈350 eV, indicating that kinetic turbulence dominates the plasma sheet dynamics.

• Ion boundary layer (IBL) and inner current sheet: ζ₂ approaches 1.4 and β exceeds 2.5, hinting at a possible transition toward electron‑scale turbulence.

The authors conclude that kinetic‑scale magnetic turbulence, governed by finite Larmor radius effects, is pervasive throughout Mercury’s near‑planet environment and extends up to the longest resolvable time scale (~20 s) limited only by data non‑stationarity. This challenges the conventional view that Mercury’s magnetosphere can be described primarily by MHD processes and underscores the necessity of incorporating kinetic physics into models of Mercury’s magnetic field structure, reconnection dynamics, and particle transport.

In summary, the study provides the first empirical evidence that ion‑kinetic turbulence is the dominant energy‑transfer mechanism in Mercury’s foreshock, magnetosheath, plasma sheet, and inner magnetopause. By introducing the continuous SF scalogram and demonstrating its robustness against non‑stationarity, the paper offers a powerful new tool for space‑plasma turbulence research, applicable to other planetary missions where data length and stationarity are limited.