Articles | Volume 39, issue 3
Regular paper
24 Jun 2021
Regular paper |  | 24 Jun 2021

Magnetosheath plasma flow model around Mercury

Daniel Schmid, Yasuhito Narita, Ferdinand Plaschke, Martin Volwerk, Rumi Nakamura, and Wolfgang Baumjohann

The magnetosheath is defined as the plasma region between the bow shock, where the super-magnetosonic solar wind plasma is decelerated and heated, and the outer boundary of the intrinsic planetary magnetic field, the so-called magnetopause. Based on the Soucek–Escoubet magnetosheath flow model at the Earth, we present an analytical magnetosheath plasma flow model around Mercury. The model can be used to estimate the plasma flow magnitude and direction at any given point in the magnetosheath exclusively on the basis of the plasma parameters of the upstream solar wind. The model serves as a useful tool to trace the magnetosheath plasma along the streamline both in a forward sense (away from the shock) and a backward sense (toward the shock), offering the opportunity of studying the growth or damping rate of a particular wave mode or evolution of turbulence energy spectra along the streamline in view of upcoming arrival of BepiColombo at Mercury.

1 Introduction

The magnetosphere of a planet constitutes an obstacle to the super-magnetosonic solar wind. Upstream of the planet a bow shock emerges, because the interplanetary magnetic field (IMF), embedded in the solar wind, cannot simply penetrate the magnetosphere. At the bow shock the super-magnetosonic solar wind plasma is decelerated and heated. The region with the subsonic, heated plasma downstream of the bow shock is called magnetosheath. The magnetosheath plays an important role in the interaction between bow shock and magnetosphere as it conveys energy between the solar wind and the planetary magnetosphere.

One of the earliest magnetosheath plasma flow models is the hydrodynamic model introduced by Spreiter et al. (1966). Basically the model solves the gas-dynamic differential equations of an unmagnetized fluid around an obstacle, represented by the magnetosphere. It has successfully been tested against in situ spacecraft data (Song et al.1999; Stahara et al.1993; Spreiter and Alksne1968) and applied to model the magnetospheres of various planets in our solar system (see Stahara2002, for a review). A decisive drawback of this model, however, is the high complexity and computational demands to calculate numerically a set of differential equations.

To reduce the computational complexity, several analytical plasma flow models have alternatively been proposed (Russell et al.1983; Kallio and Koskinen2000; Romashets et al.2008). An analytical magnetosheath flow model, which has successfully been tested against spacecraft observations at Earth, has been implemented by Soucek and Escoubet (2012). This model is based on the magnetic field model developed by Kobel and Flückiger (1994) and later modified and extended by Génot et al. (2011) to obtain a magnetosheath plasma flow model. The essential advantage of this model is its compatibility with a wide range of bow shock and magnetopause models while retaining the simplicity and computational efficiency of the original magnetic field model. Furthermore, the model allows us to calculate the plasma flow velocity at any point in the magnetosheath using only the spacecraft position and solar wind parameter upstream of the bow shock.

In this work we follow the procedure proposed by Soucek and Escoubet (2012) and rescale their terrestrial magnetosheath flow model to the space environment at Mercury. First, we introduce the Hermean bow shock and magnetopause model, used to obtain the magnetosheath plasma flow model. Second, we revisit the magnetic field model of Kobel and Flückiger (1994) which Soucek and Escoubet (2012) used to determine the plasma velocity direction in the magnetosheath. Third, we extend the model by the Rankine–Hugoniot relations in a similar way as Génot et al. (2011) to determine the velocity magnitude downstream of the shock.

The aim of this paper is to provide a tool to estimate the plasma flow at a given point of spacecraft observation inside the Hermean magnetosheath on the basis of the solar wind conditions.

2 Bow shock and magnetopause model at Mercury

In the following we use an aberrated Mercury Solar Magnetospheric (MSM) coordinate system. This coordinate system is based on the Mercury Solar Orbital (MSO) coordinate system, but its origin is shifted northward by 479 km from the MSO origin to account for Mercury's dipole offset and rotated into the solar wind velocity direction. In the MSO coordinate system the XMSO axis points sunward, the YMSO axis points antiparallel to Mercury's orbital velocity and ZMSO=XMSO×YMSO completes the right-handed system. To compensate for the aberration of the solar wind direction due to the orbital motion of Mercury around the sun, the XMSO axis is rotated antiparallel to the solar wind flow velocity direction. In the MSM coordinate system the bow shock and magnetopause models are considered to be cylindrically symmetric around the XMSM axis, reducing the three dimensions {XMSM,YMSM,ZMSM} to two dimensions {XMSM,ρMSM} with ρ=YMSM2+ZMSM2.

Slavin et al. (2009) modeled the bow shock at Mercury by a conic section of the form

(1) ξ = ( x BS - x 0 ) 2 + ρ BS 2 = p ϵ 1 + ϵ cos ϕ ,

with x0 being the distance of the focus of the conic section from the dipole center along XMSM, ρBS=yBS2+zBS2 being the distance from the XMSM axis, p being the focal parameter and ϵ being the eccentricity. With the advent of the MESSENGER (MErcury Surface, Space ENvironment, GEophysics and Ranging;  Solomon et al.2007) spacecraft in an orbit around Mercury, it was possible to characterize the spatial location of the bow shock and magnetopause statistically. Winslow et al. (2013) determined that the best-fit parameters to the bow shock are given by x0=0.5 RM, p=2.75 RM and ϵ=1.04. With these parameters the extrapolated subsolar bow shock stand-off distance is RBS=1.9 RM (Mercury radii, 1 RM∼2440 km). For this work, it is advantageous to transform Eq. (1) into the origin of the MSM coordinate system with

(2) r BS = ( ξ cos ϕ + x 0 ) 2 + ( ξ sin ϕ ) 2 , θ = arccos ξ cos ϕ + x 0 r BS ,

where rBS is the distance from the dipole center to the bow shock, and θ is the angle between rBS and the XMSM axis. Figure 1 shows a schematic illustration of the parameters ξ, ϕ, rBS and θ, which are used in the formulation of the bow shock.

Korth et al. (2015) used the magnetopause model proposed from Shue et al. (1997) and found that the MESSENGER observations of magnetopause crossing are best fit by

(3) r MP = x MP 2 + ρ MP 2 = R MP 2 1 + cos θ α ,

with α=0.5 being the best-fit flaring parameter, and RMP=1.42 RM being the subsolar stand-off magnetopause distance.

Figure 1 shows the Slavin et al. (2009) bow shock model (S09-BS) and Korth et al. (2015) magnetopause model (K15-MP) evaluated from Eqs. (2) and (3), respectively.

Figure 1Schematic representation of the parameters used in the formulation of the bow shock and magnetopause in the MSM equatorial plane. The solid red line is the bow shock evaluated from Eq. (2) (S09-BS;  Slavin et al.2009). The solid black line represents the magnetopause determined by Eq. (3) (K15-MP;  Korth et al.2015). The dashed blue lines are the bow shock and magnetopause determined by Eq. (4) from the KF94 model (Kobel and Flückiger1994).


3 The KF94 magnetic field model

To obtain the magnetosheath plasma flow direction, we follow the procedure proposed by Soucek and Escoubet (2012) and use the magnetic field model developed by Kobel and Flückiger (1994). In the following we denote this model as KF94 model and mark all quantities pertaining to the KF94 model by a tilde, e.g., r̃. In the KF94 model the bow shock (BS) and magnetopause (MP) at Mercury are modeled by parabolic surfaces at a common focus with

(4) r ̃ { BS , MP } = - cos θ + cos 2 θ + 4 R { BS , MP } b { BS , MP } sin 2 θ 2 b { BS , MP } sin 2 θ ,

with bBS=1/(4RBS-2RMP) and bMP=1/(2RMP) defined by the stand-off distances R{BS,MP}.

Under the assumption that the IMF is parallel to the solar wind, the magnetic field lines of the KF94 model represent the flow lines of the solar wind and magnetosheath plasma. Using the magnetic field vector direction in the KF94 model, Soucek and Escoubet (2012) determined the flow velocity vector at a given position r=(x,ρ) by

(5) v ̃ x = v m ( C / 2 d - C / R MP ) , v ̃ ρ = v m ( C ρ / [ 2 d ( d + x - R MP / 2 ) ] ) ,

where vm corresponds to the flow velocity magnitude, d=|r-r0| is the difference between the given position in the magnetosheath and the parabolic surface focus r0=(RMP/2,0) and C=RMP(2RBS-RMP)/(2RBS-2RMP) is a constant defined by the bow shock and magnetopause stand-off distances.

4 The magnetosheath plasma flow model around Mercury

To obtain the magnetosheath plasma flow at a specific point (r=(x,ρ)) in the magnetosheath at Mercury, we evaluate the plasma flow direction first and then determine the magnitude of the velocity vector from the Rankine–Hugoniot relation across the bow shock.

A magnetic field model is used to describe the plasma flow here. The reason for this is explained as follows. Assumptions are made such that there are no proton sinks or sources in the magnetosheath. Strictly speaking, this assumption is only weakly justified because of the neutral particles such as the hydrogen corona and the sodium exosphere around Mercury. A stationary case is taken in the fluid picture. The continuity equation reads then as

(6) ( n U ) = 0 ,

where n and U are the density and velocity of protons, respectively. Now we make an analogy such that

(7) ( n U ) B ,

holds, and Eq. (6) is equivalent to the divergence-free condition of magnetic field:

(8) B = 0 .

For a more detailed discussion, see e.g., Génot et al. (2011).

4.1 Plasma flow direction

Following the procedure proposed by Soucek and Escoubet (2012), we rescale the plasma flow direction from the KF94 model to Mercury's space environment as follows:

  1. As a first step, we calculate the angle θ=arccos(x/x2+ρ2) between r and the XMSM axis.

  2. Then we estimate the fractional distance, , of r between the bow shock and magnetopause from Eqs. (2) and (3) with

    (9) F = r ( θ ) - r BS ( θ ) r BS ( θ ) - r MP ( θ ) .
  3. Now we change into the KF94 model and calculate r̃BS(θ) and r̃MP(θ) from Eq. (4) with the angle θ. Note that the stand-off distances (RBS and RMP) and the focus (r0=(RMP/2,0)) in Eq. (4) are the best-fit values from Eqs. (2) and (3).

  4. In a next step we determine the fractional position within the magnetosheath in the KF94 model with

    (10) r ̃ ( θ ) = F [ r ̃ BS ( θ ) - r ̃ MP ( θ ) ] + r ̃ BS ( θ ) ,

    according to Eq. (9).

  5. Using Eq. (5) we evaluate the KF94 flow velocity vector, ṽ for ṽ=(ṽx,ṽρ) at the position r̃(θ). Note that the velocity magnitude vm is determined in a later step.

  6. With the obtained flow velocity vector, v, we are able to estimate the new position of an adjacent point along the same flow line r̃=r̃+ṽΔt, by choosing an infinitesimally small time increment Δt.

  7. Next we determine the angle between the new position and the XMSM axis, θ, and the fractional distance inside the KF94 magnetosheath, , using Eq. (9).

  8. Finally we transform the new position r̃(θ) back from the KF94 model to the original MSM reference frame where the magnetosheath is confined by Eqs. (2) and (3). The new position, r, inside this magnetosheath is then given by

    (11) r = F [ r BS ( θ ) - r MP ( θ ) ] + r BS ( θ ) ,

    and thus the plasma flow direction can be determined by v=(r-r)/Δt.

Applying recursively this procedure (steps 1–8) yields the plasma flow line within Mercury's magnetosheath. In Fig. 2 five examples of flow lines are shown.

Figure 2Schematic representation of the flow lines (red) in the MSM equatorial plane. The solid black lines are the bow shock and magnetopause evaluated from Eqs. (2) and (3), respectively.


4.2 Plasma flow magnitude

To evaluate the magnetosheath plasma velocity magnitude, vm, we apply the Rankine–Hugoniot (RH) equations, which relate the upstream (u) with the downstream (d) plasma conditions. The downstream plasma flow velocity directly behind the bow shock, vd, is derived by the following procedure:

  1. From the given spacecraft position in the magnetosheath, r=(x,ρ), we trace the flow line back to the bow shock. Thereto we iteratively apply steps (1)–(8) from above, with reversed increments r̃=r̃-ṽΔt in step (6), until the bow shock is reached (F=0). Then we calculate the angle θBS between the XMSM axis and the bow shock intersection at (xBS,ρBS) with θBS=arccos(xBS/xBS2+ρBS2).

  2. In a next step we determine the bow shock tangent t^ and normal n^ unit vectors where the back-traced flow line intersects the bow shock. For any point along the bow shock, the normal, n, and tangent, t, vector can easily be computed by

    (12) t = d r BS d θ cos θ - r BS sin θ e ^ x + d r BS d θ sin θ + r BS cos θ e ^ ρ , n = d r BS d θ sin θ + r BS cos θ e ^ x - d r BS d θ cos θ + r BS sin θ e ^ ρ ,

    where drBSdθ is numerically calculated from two consecutive points along the bow shock given by Eq. (2). The shock frame of reference is then defined by the normalized normal and tangent vector at θBS.

  3. In the shock reference frame the RH equations can be combined to determine the downstream velocity vector component parallel vnd and perpendicular vtd to the shock normal (see e.g., Génot2008):

    (13) v n d = v n u 1 R , v t d = v n u tan θ Vn + 1 - 1 / R M A u 2 / ( R cos 2 θ Bn ) - 1 tan θ Bn ,

    where vnu is the upstream velocity vector component parallel to the shock normal, R=ρd/ρu is the compression ratio between the upstream and downstream mass density, θVn=arctan(vtu/vnu) is the angle between the upstream velocity vector and the shock normal, θBn=arctan(Btu/Bnu) is the angle between upstream magnetic field vector and shock normal, and MAu=vnuμ0ρuBnu is the upstream Alfvén Mach number.

  4. In Eq. (13) all parameters pertain to the upstream side, except the compression ratio . However, can also be expressed by exclusively upstream parameters with (see e.g., Anderson1963)

    (14) ( M A u 2 - R ) 2 ( γ β u R + M A u 2 cos 2 θ Bn [ ( γ - 1 ) R - ( γ + 1 ) ] ) + R M A u 2 sin 2 θ Bn ( [ γ + ( 2 - γ ) R ] M A u 2 + R [ ( γ - 1 ) R - ( γ + 1 ) ] ) = 0

    where βu=(2μ0pu)/Bu2 is the ratio of the upstream thermal to magnetic pressures, and γ is the polytropic index which is typically assumed to be γ=5/3. Equation (14) is equivalent to Eq. (2.43) in Anderson (1963) with different notations. The solution is given in the form of compression ratio as a function of the shock angle θBn. Solutions exist for a compression ratio in the range 1R4 when using Eq. (14). Another class of solutions also exists for the expansion (ℛ<1) with a decrease of entropy from the upstream onto the downstream side. The latter case is not physically relevant and is not considered here. The upper limit of compression ratio (ℛ=4) corresponds to the limit of high Alfvén Mach number (MA→∞) under a polytropic index of γ=5/3.

  5. By solving Eq. (14) for the downstream velocity magnitude vd=vnd2+vtd2 is therefore entirely determined by only the upstream plasma parameters.

Since the detailed density profile along the flow line is unknown, we assume in a first approximation a constant plasma density and thus a constant velocity magnitude along the flow line. Therefore, the velocity magnitude at a given point r directly corresponds to the velocity magnitude downstream of the shock, and vm=vd. Although this assumption has the tendency to underestimate the velocity close to the magnetopause, it yields satisfactory results in a first approach (Génot et al.2011).

The entire procedure from above is implemented in an IDL computer program which can be retrieved from OSF (Schmid2020). The program is designed to evaluate the plasma flow velocity vector at a given observation point of a spacecraft inside the Hermean magnetosheath exclusively on the basis of the upstream solar wind conditions. As the solar wind input parameters, we use the solar wind propagation model of Tao et al. (2005), which is modified by the orbital motion of Mercury. The model is a one-dimensional magnetohydrodynamic model and takes the OMNI dataset as input to compute propagation at all solar system bodies including Mercury. The correction for the orbital motion is achieved by adding the solar wind velocity vector, VSW, (which is radially away from the Sun) and the orbital motion velocity vector of Mercury, Vmercury, with V=VSW+Vmercury. To obtain Vmercury, we use the dataset provided from the Navigation and Ancillary Information Facility (NAIF; Acton1996), which provide the distance between Mercury and Sun, D, and absolute velocity of Mercury, Vmercury. To determine aberration velocity vector, we first calculate the aberration angle ϕ on the basis of Mercury's elliptical orbit with

(15) ϕ = arctan b / a 2 - b 2 sin arccos ( ( 1 - p / D ) / e ) ,

where a is the semimajor axis, b the semiminor axis, p the semi-latus rectum and e the eccentricity of the ellipse. With the aberration angle ϕ and the consideration whether Mercury moves towards or away from the sun, we subsequently obtain the aberration velocity vector Vmercury,x=±Vmercurysin(ϕ) and Vmercury,y=±Vmercurycos(ϕ). The transformation due to this abberation effect is made by applying a two-dimensional rotation matrix to the spatial coordinates (spanning the x and ρ coordinates):

(16) x y = cos θ a - sin θ a sin θ a cos θ a x ρ ,

where the aberration angle θa is given by the radial solar wind velocity, VSW, and the apparent solar wind velocity, V with θa=arccosVSWV. In Fig. 3 the results of the model are shown for the average solar wind plasma parameters during the entire MESSENGER operation service between 2011 and 2015. After modifying the solar wind velocity vector of the Tao et al. (2005) model by the orbital motion of Mercury (Acton1996), the average input solar wind plasma parameters for our model are density of nu40cm-3, temperature of T≈18 eV, flow speed of |Vu|-400km/s, magnetic field magnitude of |Bu|20nT with the radial component Br≈18 nT (ignoring the sign) and the tangential component Bt≈16 nT (ignoring the sign). The mean values of density, temperature, flow speed and magnetic field are valid for nearly 1500 d of observations of MESSENGER (confirmed by one of the reviewers). It is worthwhile to note that one finds an angle of 23 from the Tao dataset, which is consistent to the spiral angle at Mercury at an average position of 0.4 astronomical units (AU) from the Sun for a solar wind speed of 400km/s. The Bx component is computed from the By component using the Parker spiral field in the Tao model. The Alfvén Mach number in the upstream region is MA=vu/VA5.8 (with an Alfvén speed of VA69km/s), and the plasma parameter beta (upstream) is β=2μ0nkBT/B20.72 (where μ0 is the permeability of free space, and kB is the Boltzmann constant) in our setup. Color coded is the obtained velocity magnitude vm. Additionally, the back-traced flow line from a virtual spacecraft located at xMSM=-3RM and ρMSM=3 RM is plotted (green line). At the bow shock intersection the calculated shock normal n^ and tangent t^ are illustrated in thin black lines. The thick black line downstream of the bow shock shows the velocity vector direction evaluated from the RH relations, which is in good agreement with the streamline direction determined by the KF94 model.

Figure 3Color-coded flow speed calculated from Eq. (13) with the averaged upstream parameters from the Tao et al. (2005) solar wind propagation model at Mercury between 2011 and 2015. Also plotted is a schematic representation of back-traced flow line (green) from a virtual spacecraft (S/C) to the bow shock (BS) with the respective shock normal n^ and tangent t^ obtained from Eq. (12). The thick black line downstream of the BS is the downstream velocity vector direction determined by Eq. (13). Alfvén Mach number is 5.8, and plasma beta is 0.79.


At the virtual spacecraft position the model predicts a magnetosheath plasma flow velocity of vx-200km/s and vρ17km/s.

A naive picture of the compression by a factor of 4 (in the limit of high Mach number) is not realistic, because the interplanetary magnetic field around Mercury reaches a magnitude of 20 to 50 nT, and the Alfvén Mach number is correspondingly smaller than that around the Earth by 20 to 50, respectively. A picture of the constant density and a reduced flow speed to 1/4 of the solar wind speed at the magnetopause (along the streamline tangential to the magnetopause) is not valid, either, since the adiabatic expansion breaks down, and the model is not applicable to the flow in the subsolar region and at the magnetopause.

5 Discussion and conclusions

Here we present the first analytical magnetosheath plasma flow model for the space environment around Mercury. The model is based on the magnetosheath model by Soucek and Escoubet (2012), which has successfully been tested against spacecraft observations at Earth. The proposed model is relatively simple to implement and provides the possibility to trace the flow lines inside the Hermean magnetosheath.

The model presented in this paper is generally robust and easy to implement for its analytic expression using upstream parameters. It can help to determine the local plasma conditions of a spacecraft in the magnetosheath exclusively on the basis of the upstream solar wind parameters. Two applications are in mind in view of the BepiColombo mission, where the Mercury Magnetospheric Orbiter (MMO also referred to as Mio) will probe Mercury's magnetosheath and solar wind with unprecedented fast measurements of the particle distribution functions: (1) the Tátrallyay method to observationally determine the growth rate or damping rate of specific mode such as the mirror mode along the streamline (Tátrallyay and Erdös2002; Tátrallyay et al.2008) and (2) the Guicking method to observationally track the spectral evolution of turbulent fluctuations along the streamline (Guicking et al.2010, 2012).

At the moment a comparison with plasma data is technically not feasible for our model. Due to the limited particle measurements on board MESSENGER, it is not possible to obtain the plasma parameters properly in the solar wind and magnetosheath, that is, by covering the full velocity distributions and to compare with the model velocities. Above all, the plasma instrument is located behind the heat shield and has just a limited field of view. Due to this fact, the majority (thermal core part) of solar wind particles cannot be detected. Comparison with the numerical simulations would be another possibility to test for the model, but a quantitative comparison remains a challenge for the reason that there are large discrepancies in the density and flow velocity among various simulation models on the dayside from the subsolar point to the northern terminator (Aizawa et al.2021).

An approximation that the magnetic field is more aligned with the solar wind flow direction is more justified at Mercury than at the Earth because of the Parker spiral nature. One of the possible consequences of our assumption is that the magnetic field magnitude would change or evolve in the same sense as that of the plasma (or particle number flux) in the magnetosheath. The correlation between the magnetic field and that of the particle flux in the magnetosheath would ideally be tested against the plasma and magnetic field data on the arrival of BepiColombo at Mercury. The change in the number density can be interpreted as the change in the cross-sectional area of a flux tube across which the plasma streams. The change in the flow velocity can then be compared with that from the adiabatic expansion and that from the measurement.

The following remarks are drawn as scientific message. First, Fig. 1 visually demonstrates that different models predict different shapes of the tail and magnetosheath, which is an overlooked issue in the Mercury magnetosphere community. Second, Fig. 2 shows that the flow lines near the subsolar region (Sun-to-planet line if neglecting the planetary orbital motion) expand abruptly so that the adiabatic expansion may break down. In particular, the adiabatic expansion plays an important role in predicting the flow in the magnetosheath. Logical continuation of our model construction would be to evaluate also the density and velocity profile along the flow line and to test under which conditions the adiabatic expansion breaks down.

Although the proposed model has a good performance overall for a wide range of upstream conditions, the accuracy strongly depends on the used bow shock and magnetopause model. Here we utilize the bow shock and magnetopause model from Slavin et al. (2009) and Korth et al. (2015), which were adopted from MESSENGER boundary crossing observations.

The presented model is cylindrically symmetric around the XMSM axis. In reality, however, non-radial IMF conditions will lead to a spatially asymmetric magnetosheath (Nishino et al.2008; Dimmock and Nykyri2013; Dimmock et al.2016). On the quasi-perpendicular side, where the shock-normal angles θBn are greater than 45, the magnetosheath is known to be thicker with larger plasma flow velocities than on the quasi-parallel side, where θBn<45. Such asymmetries cannot be reproduced by the simple model presented here but should be addressed in future work.

Furthermore, the method used to determine the flow velocity magnitude can possibly be improved. Here we assumed a constant plasma density and velocity along the flow line which has the tendency to underestimate the plasma velocity in regions with lower densities, e.g., close to the magnetopause. Génot et al. (2011) proposed a simple ad hoc model of a plasma density profile which has been implemented by Soucek and Escoubet (2012). While this ad hoc density model showed good correspondence with in situ spacecraft plasma observation at Earth, the solar wind and magnetospheric conditions at other planets can be very different (like at Mercury) and thus might give a worse prediction (Soucek and Escoubet2012). Our model inherits the properties from the Soucek–Escoubet model by scaling the Kobel–Flückinger model of the near-Earth environment: (1) time stationary flow and (2) axially symmetric around the axis of (apparent) solar wind penetrating the planet. Assumption of time stationary flow may break down when the change in the solar wind state is not negligible. Assumption of the axisymmetric magnetosheath may also break down when the magnetopause location is not symmetric between the northern and the southern hemisphere (in particular, in the tail region).

At this stage we decided not to include an ad hoc density profile, also because it can hardly be tested due to the limited plasma observations around Mercury. The assumption of constant density implies a constant velocity along a given flow line. The velocity profile may vary considerably from that estimated in the earlier models such as the Spreiter model (Spreiter et al.1966), the Genot model (Génot et al.2011), and the Soucek–Escoubet model (Soucek and Escoubet2012). As mentioned above, the constant density will likely underestimate the propagation timing in the magnetosheath. In future work such a density profile should be evaluated and included.

Code availability

An IDL program to evaluate plasma flow velocity vector in Mercury's magnetosheath from solar wind parameters of the Tao solar wind propagation model can be retrieved from OSF: (last access: 25 May 2021) (Schmid2020).

Data availability

The plasma data of the heliospheric Tao model are open-access data and can be retrieved on the AMDA website ( (last access: 25 November 2020), Centre de Données de la Physique des Plasmas (CDPP), 2018) via the WorkSpace Explorer: DataBase/Solar Wind Propagation Models/Tao Model/SW Input OMNI (Tao et al.2005). The orbital motion data of Mercury are provided by the Navigation and Ancillary Information Facility (NAIF) and can be retrieved on the NAIF website under (last access: 6 May 2020) (Acton1996).

Author contributions

DS initiated this study, collected the data and implemented the method. FP, YN, MV and WB helped with evaluating the article.

Competing interests

The authors declare that they have no conflict of interest.

Financial support

This research has been supported by the Österreichische Forschungsförderungsgesellschaft (grant no. 865967).

Review statement

This paper was edited by Anna Milillo and reviewed by two anonymous referees.


Acton, C. H.: Ancillary data services of NASA's Navigation and Ancillary Information Facility, Planet. Space Sci., 44, 65–70,, 1996. a, b, c

Anderson, J. E.: Magnetohydrodynamic shock waves, Cambridge: MIT Press, 1963. a, b

Aizawa, S., Griton, L. S., Fatemi, S., Exner, W., Deca, J., Pantellini, F., Yagi, M., Heyner, D., Génot, V., André, N., Amaya, J., Murakami, G., Beigbeder, L., Gangloff, M., Bouchemit, M., Budnik, E., and Usui, H.: Cross-comparison of global simulation models applied to Mercury's dayside magnetosphere, Planet. Space Sci., 198, 105176,, 2021. a

Dimmock, A. P. and Nykyri,K.: The Statistical Mapping of Magnetosheath Plasma Properties Based on THEMIS Measurements in the Magnetosheath Interplanetary Medium Reference Frame, J. Geophys. Res.-Space, 118, 4963–4976,, 2013. a

Dimmock, A. P., Nykyri, K., Osmane, A., and Pulkkinen, T. I.: Statistical mapping of ULF Pc3 velocity fluctuations in the Earth's dayside magnetosheath as a function of solar wind conditions, Adv. Spac. Res., 58, 196–207,, 2016. a

Génot, V.: Mirror and Firehose Instabilities in the Heliosheath, Astrophys. J., 687, 119–122,, 2008. a

Génot, V., Broussillou, L., Budnik, E., Hellinger, P., Trávníček, P. M., Lucek, E., and Dandouras, I.: Timing mirror structures observed by Cluster with a magnetosheath flow model, Ann. Geophys., 29, 1849–1860,, 2011. a, b, c, d, e, f

Guicking, L., Glassmeier, K.-H., Auster, H.-U., Delva, M., Motschmann, U., Narita, Y., and Zhang, T. L.: Low-frequency magnetic field fluctuations in Venus' solar wind interaction region: Venus Express observations, Ann. Geophys., 28, 951–967,, 2010. a

Guicking, L., Glassmeier, K.-H., Auster, H.-U., Narita, Y., and Kleindienst, G.: Low-frequency magnetic field fluctuations in Earth's plasma environment observed by THEMIS, Ann. Geophys., 30, 1271–1283,, 2012. a

Kallio, E. J. and Koskinen, H. E. J.: A semiempirical magnetosheath model to analyze the solar wind-magnetosphere interaction, J. Geophys. Res.-Space, 105, 469–479,, 2000. a

Kobel, E. and Flückiger, E. O.: A model of the steady state magnetic field in the magnetosheath, J. Geophys. Res.-Space, 99, 617–622,, 1994. a, b, c, d

Korth, H., Tsyganenko, N. A., Johnson, C. L., Philpott, L. C., Anderson, B. J., Al Asad, M. M., Solomon, S. C., and McNutt Jr., R. L.: Modular model for mercury's magnetospheric magnetic field confined within the average observed magnetopause, J. Geophys. Res.-Space, 120, 4503–4518,, 2015. a, b, c, d

Nishino, M. N., Fujimoto, M., Phan, T.-D., Mukai, T., Saito, Y., Kuznetsova, M. M., and Rastätter, L.: Anomalous Flow Deflection at Earth's Low-Alfvén-Mach-Number Bow Shock, Phys. Res. Lett., 101, 065003,, 2008. a

Romashets, E. P., Poedts, S., and Vandas, M.: Modeling of the magnetic field in the magnetosheath region, J. Geophys. Res.-Space, 113, A2,, 2008. a

Russell, C. T., Luhmann, J. G., Odera, T. J., and Stuart, W. F.: The rate of occurrence of dayside Pc 3,4 pulsations: The L-value dependence of the IMF cone angle effect, Geophys. Res. Lett., 10, 663–666,, 1983. a

Schmid, D.: Mercury Streamline Model: IDL code, OSF, available at: (last access: 25 May 2021), 2020. a, b

Shue, J.-H., Chao, J. K., Fu, H. C., Russell, C. T., Song, P., Khurana, K. K., and Singer,H. J.: A new functional form to study the solar wind control of the magnetopause size and shape, J. Geophys. Res.-Space, 102, 9497–9511,, 1997. a

Slavin, J. A., Anderson, B. J., Zurbuchen, T. H., Baker, D. N., Krimigis, S. M., Acuña, M. H., Benna, M., Boardsen, S. A., Gloeckler, G., Gold, R. E., Ho, G. C., Korth, H., McNutt Jr., R. L., Raines, J. M., Sarantos, M., Schriver, D., Solomon, S. C., and Trávnícek, P.: Messenger observations of Mercury's magnetosphere during northward IMF, Geophys. Res. Lett., 36, 2,, 2009. a, b, c, d

Solomon, S. C., McNutt Jr., R. L., Gold, R. E., and Domingue, D. L.: MESSENGER Mission Overview, Space Sci. Rev., 131, 3–39,, 2007. a

Song, P., Russell, C. T., Zhang, X. X., Stahara, S. S., Spreiter, J. R., and Gombosi, T. I.: On the processes in the terrestrial magnetosheath: 2. Case study, J. Geophys. Res.-Space, 104, 357–373,, 1999. a

Soucek, J. and Escoubet, C. P.: Predictive model of magnetosheath plasma flow and its validation against Cluster and THEMIS data, Ann. Geophys., 30, 973–982,, 2012. a, b, c, d, e, f, g, h, i, j

Spreiter, J. and Alksne, A.: Comparison of theoretical predictions of the flow and magnetic fields exterior to the magnetosphere with the observations of pioneer 6, Planet Space Sci., 16, 971–979,, 1968. a

Spreiter, J., Summers, A., and Alksne, A.: Hydromagnetic flow around the magnetosphere, Planet Space Sci., 14, 223–253,, 1966. a, b

Stahara, S. S.: Adventures in the magnetosheath: two decades of modeling and planetary applications of the Spreiter magnetosheath model, Planet Space Sci., 50, 421–442,, 2002.  a

Stahara, S. S., Rachiele, R. R., Molvik, G. A., and Spreiter, J. R.: Development of a preliminary solar wind transport magnetosheath forecast model, NASA STI/Recon Technical Report N., 1993. a

Tátrallyay, M. and Erdös, G.: The evolution of mirror mode fluctuations in the terrestrial magnetosheath, Planet Space Sci., 50, 593–599,, 2002. a

Tátrallyay, M., Erdös, G., Balogh, A., and Dandouras, I.: The evolution of mirror type magnetic fluctuations in the magnetosheath based on multipoint observations, Adv. Space Res., 41, 1537–1544,, 2008. a

Tao, C., Kataoka, R., Fukunishi, H., Takahashi, Y., and Yokoyama T.: Magnetic field variations in the jovian magnetotail induced by solar wind dynamic pressure enhancements, J. Geophys. Res.-Space, 110, A11,, 2005. a, b, c, d

Winslow, R. M., Anderson, B., Johnson, C., Slavin, J., Korth, H., Purucker, M., Baker, D. N., and Solomon, S.: Mercury's magnetopause and bow shock from MESSENGER Magnetometer observations, J. Geophys. Res.-Space, 118, 2213–2227,, 2013. a

Short summary
In this work we present the first analytical magnetosheath plasma flow model for the space environment around Mercury. The proposed model is relatively simple to implement and provides the possibility to trace the flow lines inside the Hermean magnetosheath. It can help to determine the the local plasma conditions of a spacecraft in the magnetosheath exclusively on the basis of the upstream solar wind parameters.