Magnetic Field and Plasma Asymmetries Between the Martian Quasi-Perpendicular and Quasi-Parallel Magnetosheaths

The Martian magnetosheath acts as a conduit for mass and energy transfer between the upstream solar wind and its induced magnetosphere. However, our understanding of its global properties remains limited. Using nine years of data from NASA’s Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we performed a quantitative statistical analysis to explore the spatial distribution of the magnetic fields, solar wind and planetary ions in the magnetosheath. We discovered significant asymmetries in the magnetic field, solar wind protons, and planetary ions between the quasi-perpendicular and quasi-parallel magnetosheaths. The asymmetries in the Martian magnetosheath exhibit both similarities and differences compared to those in the Earth’s and Venus’ magnetosheaths. These results indicate that the Martian magnetosheath is distinctly shaped by both shock geometry and planetary ions.

💡 Research Summary

This paper presents a comprehensive statistical study of magnetic field and plasma asymmetries in the Martian magnetosheath, focusing on the differences between quasi‑perpendicular (Q⊥) and quasi‑parallel (Q∥) shock geometries. Using nine years of observations from NASA’s MAVEN mission (December 2014–December 2023), the authors analyze magnetic field data from the MAG instrument, ion composition from STATIC, and solar‑wind proton measurements from SWIA. The analysis is performed in the Mars‑Solar‑Electric (MSE) coordinate system, but to treat both IMF Bx polarities uniformly the authors introduce a new Magnetosheath‑MSE (MMSE) system that flips the Y and Z axes when Bx is negative, ensuring that the Q∥ flank always lies in +Y_MMSE and the Q⊥ flank in –Y_MMSE.

Data are binned on a 0.1 Rₘ × 0.1 Rₘ grid in the MMSE X‑Y plane, limited to the magnetic equatorial band (|Z_MMSE| < 1 Rₘ) to reduce the influence of the heavy‑ion plume that dominates the ±E_sw hemispheres. Only bins containing more than 40 measurements and with |B_obs| ≤ 10 |B_model| (to exclude strong crustal‑field contamination) are retained. Each orbit’s measurements are normalized by the corresponding upstream solar‑wind conditions, and the magnetosheath region is defined using the bow‑shock (BS) and magnetic pile‑up boundary (MPB) models of Trotignon et al. (2006). The inner two‑thirds of the sheath are examined, and the sheath is further divided radially into three equal‑thickness layers.

The authors compute an angular coordinate θ between the position vector and +X_MMSE, assigning positive values to Q∥ and negative values to Q⊥. Data are rebinned in 15° overlapping angular sectors (50 % overlap) from –120° to +120°. For each sector they calculate median values of normalized magnetic field strength (|B|/|B_IMF|), proton density (n_H⁺/n_sw), proton speed (|v_H⁺|/|v_sw|), and heavy‑ion (O⁺, O²⁺) density. An asymmetry index A = 100 × (Q∥ – Q⊥)/(Q∥ + Q⊥) is then derived, with uncertainties estimated from the standard error of the mean.

Key findings:

1. Magnetic field strength – |B|/|B_IMF| is systematically larger in the Q⊥ sheath (average 3–5) than in the Q∥ sheath (average 2–4). The asymmetry index for magnetic field ranges from 3 % to 14 %, peaking near local noon where the shock compression is strongest. This behavior matches Rankine‑Hugoniot expectations: a quasi‑perpendicular shock amplifies the tangential magnetic component, while a quasi‑parallel shock leaves the normal component largely unchanged.

2. Proton density and velocity – Both n_H⁺/n_sw and |v_H⁺|/|v_sw| are higher on the Q∥ flank, especially between 60° and 90° from noon, yielding asymmetry values of 5 %–12 %. The Q∥ sheath corresponds to the foreshock region under typical Parker‑spiral IMF angles, where weaker magnetic turbulence allows solar‑wind protons to penetrate more deeply.

3. Heavy‑ion (O⁺, O²⁺) distribution – The heavy‑ion plume, aligned with the solar‑wind electric field (+E_sw), shows a pronounced enhancement on the Q∥ side. This reflects Mars’ extended exosphere and low gravity, which supply abundant planetary ions that are preferentially convected by the motional electric field.

4. Influence of shock geometry vs. planetary ions – The observed asymmetries arise from a combination of shock geometry (determining magnetic compression and foreshock extent) and the presence of planetary ions (modifying mass loading, pressure balance, and kinetic scales). The authors note that single‑fluid Rankine‑Hugoniot theory underestimates the magnetic amplification; a two‑fluid treatment (solar‑wind protons + planetary ions) is required, consistent with recent MAVEN shock‑crossing studies.

Comparisons with Earth and Venus reveal that while all three planets exhibit Q∥/Q⊥ asymmetries, Mars displays a distinct signature because planetary ions dominate the mass budget in its sheath. At Earth, electron temperature anisotropy and ion temperature dominate asymmetries; at Venus, induced magnetosphere effects lead to different spatial patterns.

The paper concludes that the Martian magnetosheath is uniquely shaped by both shock geometry and heavy‑ion loading, making it an excellent natural laboratory for studying kinetic plasma processes in mixed‑fluid environments. The authors recommend future work employing high‑resolution hybrid or multi‑fluid simulations to reproduce the observed asymmetry indices and to explore seasonal and solar‑cycle variations.