Solar wind interaction with a planetary off-center mini-magnetosphere: A case study for Mars

This thesis asks a simple but nontrivial question: when magnetization is localized and asymmetric, does it help a planet retain its atmosphere, or can it open pathways that make loss easier?

---

Introduction

A long-standing idea in planetary habitability is that a global magnetic field acts as a protective shield, limiting atmospheric erosion by deflecting stellar-wind plasma [1,2]. At the same time, recent work has challenged this “magnetic umbrella” picture by arguing that magnetized planets can also open pathways for energy and momentum transfer into the upper atmosphere, potentially enhancing loss under some conditions [3,4,5]. This tension matters for exoplanets: if magnetization systematically helps or hurts atmospheric retention, it changes how we interpret which planets can plausibly maintain long-lived atmospheres and, by extension, meet one of the necessary ingredients for habitability.

To isolate the role of localized magnetization, we use Mars as a motivating case study: its localized crustal fields, weaker gravity, and nominal solar-wind conditions provide a unique natural laboratory to examine how magnetization alone reshapes the interaction and modulates the inferred planetary mass-flux proxy in our simulations. We compare three cases: a baseline with no dipole, and two “mini-magnetosphere” cases where the interplanetary magnetic field (IMF) is purely along the \(z\)-direction and oriented either southward or northward relative to the dipole axis, i.e. \((0,0,-1)\) or \((0,0,1)\).

---

Computational setup

We use the open-source PLUTO framework to evolve the three-dimensional compressible, resistive, magnetohydrodynamics (MHD) equations for a star-planet interaction problem. The solver advances conservation laws for mass, momentum, and total energy, together with the magnetic induction equation including magnetic diffusivity.

The plasma is treated with an adiabatic equation of state with polytropic index \(\gamma=5/3\), and we include a uniform magnetic diffusivity \(\eta=\times10^{9}\,\mathrm{m^2\,s^{-1}}\) motivated from [6] study.

The simulations are performed in a planet-fixed Cartesian coordinate system with the planet centered at \((0,0,0)\). The \(x\)-axis points from the planet against the pristine upstream stellar-wind flow, the \(z\)-axis points from the north magnetic pole toward the south magnetic pole (zero tilt by definition here), and the \(y\)-axis completes a right-handed system.

\[ x\in[-25R_M,\,100R_M],\quad y\in[-25R_M,\,25R_M],\quad z\in[-100R_M,\,100R_M], \qquad R_M=3389.5\,\mathrm{km} \]

The base mesh uses \([N_x,N_y,N_z]=[317,234,400]\) with non-uniform refinement. The finest region \([-3R_M,3R_M]\) is resolved at \(0.05R_M\) (about \(169.5\,\mathrm{km}\)), with coarser resolution farther out. Upstream conditions are injected at the left boundary (the \(y\)-\(z\) plane). The wind velocity is \([270,0,0]\,\mathrm{km\,s^{-1}}\), the temperature is \(4\times10^4\,\mathrm{K}\), and the density is set to four times the ambient density. The IMF has strength \(4\,\mathrm{nT}\) and is initialized purely along \(z\), oriented northward or southward depending on the run.

---

Mini-magnetosphere implementation

To mimic Mars’ highly inhomogeneous crustal magnetization, We introduce a localized, off-centered dipolar magnetic field in the southern hemisphere. This is motivated by observations and modeling showing that Mars’ strongest crustal sources are concentrated in the south, producing intense regional anomalies rather than a global dipole [7]. In the hybrid picture of the Mars-Sun interaction, strong crustal regions can generate localized closed-field structures and shift plasma boundaries, including magnetopause-like features above strong anomalies [8].

Concretely, the dipole is buried beneath the surface and offset to represent a southern-hemisphere anomaly. The goal is a minimal, tunable magnetic configuration that captures the key geometric consequence of crustal fields: a strong, localized magnetization, while keeping the system interpretable for controlled comparisons across IMF orientations. The above figure shows the initial magnetic field configuration.

---

Numerical simulations in action

To make the simulation less abstract, the animation above shows how the magnetic structure evolve and current sheets from the initial condition toward a quasi-steady interaction. Even in a single-fluid resistive-MHD model, you can see the bow-shock formation, field-line draping, and the way the localized southern anomaly reshapes the near-planet topology.

---

Results

Across the magnetized runs, a hybrid magnetic topology develops (above figure): solar-wind field lines drape around the planet and interact with the localized closed-field region associated with the mini-magnetosphere due to magnetic reconnection. In regions of strong shear and current density, the simulations show topological reconfiguration and weakening of magnetic structures in the wake.

Next, we investigated the effects of the dipolar field on the structure and location of plasma boundaries. Bow shock formed at a distance of roughly \(-1.52 R_M\), similar to [1,2], and displayed asymmetry (right panel in above figure) when mini-magnetosphere was introduced. However, the value was not significant. No magnetopuase boundary formed in our simulations, but the subsolar location of a similar region is estimated by identifying the respective signatures. The simulations exhibit the formation of a mini-magnetopause in the southern hemisphere. Besides, we observed the marsward displacement of this boundary during the periods of S-IMF.

Finally, we compared the atmospheric escape in no dople case against mini-magnetosphere when the IMF is directed in southern and northern directions. The above figure shows that the loss is maximum in the no-dipole case, favoring the role of mini-magnetosphere in the protection of the atmosphere. We also find that the higher atmospheric escape is happening in S-IMF direction relative to N-IMF, probably because the magnetic reconnection is favored in southward-directed IMF, which ushers the erosion of mini-magnetopause.

---

Discussion and conclusion

Overall, introducing a localized southern-hemisphere mini-magnetosphere changes the near-planet magnetic geometry and modestly alters the redistribution of flow and magnetic stresses in the interaction region. The magnetized runs show a more structured obstacle for the incident wind and a localized magnetopause-like shielding region above the anomaly, while the unmagnetized run more closely resembles a purely induced, draped-field interaction. These differences align with the broader hybrid magnetosphere picture motivated by strong crustal magnetic regions on Mars.

Atmospheric loss diagnostics must be interpreted carefully. The reported “escape” is best viewed as a localized planetary-plasma mass-flux proxy estimated within a single-fluid MHD model using tracers, not as a direct measurement of true ion outflow escape. The model does not include multi-species chemistry, kinetic escape physics, or the microphysical dissipation required to interpret field-line breaking as physical reconnection. Within this controlled framework, I find that a localized southern anomaly can act as a partial shield that reorganizes the interaction geometry and can modestly suppress the inferred local mass-flux proxy relative to the unmagnetized case, with sensitivity to IMF orientation.

The broader lesson is the one I care about: magnetization is not universally protective or universally harmful. The outcome depends on geometry, upstream forcing, and how dissipation in the modeled system regulates pathways for energy and momentum transfer into the upper atmosphere.

---

References

---

Project Details

Author: Sajal Gupta

Contributors: Dibyendu Nandy, Arnab Basak

Year: 2019-2020

Method: 3D resistive MHD (PLUTO), controlled comparisons across IMF orientations

Keywords: Mars, crustal magnetic fields, mini-magnetosphere, bow shock, draping, mass-flux proxy