Using a 3D multispecies magnetohydrodynamic model, we investigated the effect of the solar wind dynamic pressure (Pd) with different densities and velocities on the subsolar standoff distance (r0) of the Martian magnetic pileup boundary (MPB). We fixed the solar maximum condition, the strongest crustal field located in the dayside region, and the Parker spiral interplanetary magnetic field at Mars. We simulated 35 cases with a Pd range of 0.1494 to 7.323 nPa (solar wind number density n ∈ [1, 9] cm−3, and solar wind velocity V ∈ [−258, −1344] km s−1). The main results are as follows. (1) r0 decreases with increasing Pd according to the power-law relations. For the same Pd, a higher solar wind velocity (lower density) results in a larger r0 of the Martian MPB. (2) A higher solar wind density leads to a lower ratio of the compressed magnetic field strength to the crustal field strength and a larger plasma β under the same Pd. This indicates that the thermal pressure at the Martian MPB plays a significant role for the compressed magnetic field. Because the magnetic pileup process is stronger for a higher solar wind velocity, the magnetic pressure at the Martian MPB is increased. As a result, the thermal pressure decreases and r0 of the Martian MPB increases. (3) We present a new formula of r0 with the parameters of the solar wind dynamic pressure, number density, and velocity.

Using global magnetohydrodynamic simulations, we construct a 3D parametric model of the Martian magnetic pileup boundary (MPB). This model employs a modified parabola function defined by four parameters. The effects of the solar wind dynamic pressure, the solar wind densities and velocities, and the intensity and orientation of the interplanetary magnetic field (IMF) are examined using 267 simulation cases. The results from our parametric model show that (1) the MPB moves closer to Mars when the upstream solar wind dynamic pressure (Pd) increases, the subsolar standoff distance decreases and the flaring degree of the Martian MPB increases with an increasing Pd according to the power-law relations. For the same Pd, a higher solar wind velocity (a lower density) leads to a farther location of the MPB from Mars, along with a larger flaring degree, which is explained by the higher solar wind convection electric fields and a stronger magnetic pileup process under these conditions. (2) Larger Y or Z components of the IMF, BY or BZ, result in a thicker pileup region and a farther MPB location from Mars, as well as a decrease in the flaring degree. The radial IMF component, BX, has little effect on the geometry of the MPB. (3) In most of the simulations used to derive the current parametric model, the strongest Martian crustal magnetic field is located on the dayside. However, for a larger value of the southward IMF, the Martian MPB is located farther away in the northern hemisphere instead of the southern hemisphere. The north-south asymmetry of the Martian MPB with the southern hemisphere being farther away is observed for other IMF directions. We suggest that the magnetic reconnection of the southward IMF with the crustal field that occurs at middle latitudes of the southern hemisphere results in different magnetic field topologies and the closer location of the MPB under these conditions. Our model results show a relatively good agreement with the previous empirical and theoretical models.

The Martian bow shock (BS) is generated with the mass-loading and magnetic pileup processes when the solar wind interacts with the Martian ionosphere. In this vein, the interplanetary magnetic field (IMF) frozen in the solar wind can affect the location of the Martian BS, which is less reported. Based on the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we manually identify 10,283 BS crossings during a period of the gradually declining solar cycle phase (2014 October–2020 December) and investigate the effects of the intensity and orientation of the IMF on the Martian BS. In the Mars Solar Orbital coordinate system, our results show the following: (1) The Martian BS, including the subsolar and flank regions, linearly moves away from Mars when the IMF intensity increases, which confirms the theoretical and the MHD simulation results. (2) Under the radial IMF condition, we first demonstrate that the subsolar and flank regions of the Martian BS are situated closer to Mars compared to other IMF situations. This might be caused by the weaker magnetic pileup process and the “low-pressure magnetosheath” model under the radial IMF condition. (3) Moreover, the cross section of the Martian BS is elongated in the north–south direction when the Y component of the IMF is dominant, which is on account of the fast magnetosonic speed effect and verifies the elongation phenomenon of the terrestrial BS. The IMF intensity and orientation effects cannot be ignored and should be considered in future models of the Martian BS.