# The Mie representation for Mercury's magnetic field

S. Toepfer, Y. Narita, K. -H. Glassmeier, D. Heyner, P. Kolhey, U.

Motschmann, Benoit Langlais

► To cite this version:

S. Toepfer, Y. Narita, K. -H. Glassmeier, D. Heyner, P. Kolhey, et al.. The Mie representation for Mercury's magnetic field. Earth Planets and Space, 2021, 73 (1), 10.1186/s40623-021-01386-4. hal-03284462

HAL Id: hal-03284462

https://hal.science/hal-03284462

Submitted on 27 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

FULL PAPER

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# The Mie representation for Mercury's magnetic field

S. Toepfer<fnref n="1" />*, Y. Narita²,³, K. -H. Glassmeier³, D. Heyner³, P. Kolhey³, U. Motschmann<fnref n="1" />,⁴ and B. Langlais⁵

**Abstract**

The parameterization of the magnetospheric field contribution, generated by currents flowing in the magnetosphere is of major importance for the analysis of Mercury's internal magnetic field. Using a combination of the Gauss and the Mie representation (toroidal–poloidal decomposition) for the parameterization of the magnetic field enables the analysis of magnetic field data measured in current carrying regions in the vicinity of Mercury. In view of the BepiColombo mission, the magnetic field resulting from the plasma interaction of Mercury with the solar wind is simulated with a hybrid simulation code and the internal Gauss coefficients for the dipole, quadrupole and octupole field are reconstructed from the data, evaluated along the prospective trajectories of the Mercury Planetary Orbiter (MPO) using Capon's method. Especially, it turns out that a high-precision determination of Mercury's octupole field is expectable from the future analysis of the magnetic field data measured by the magnetometer on board MPO. Furthermore, magnetic field data of the MESSENGER mission are analyzed and the reconstructed internal Gauss coefficients are in reasonable agreement with the results from more conventional methods such as the least-square fit.

**Keywords:** Mie representation, Poloidal and toroidal magnetic fields, Thin shell approximation, Gauss representation, Capon's method

**Introduction**

Characterization of Mercury's internal magnetic field is one of the primary goals of the BepiColombo mission (Benkhoff et al. 2010). The magnetic field in the vicinity of Mercury is composed of internal and external parts. The internal field originates in the dynamo-generated field, crustal remanent field and induction field, which are mainly dominated by the dipole, quadrupole and octupole field. The external field originates from currents flowing in the magnetosphere. For a clear separation of the internal magnetic field from the external field, each part of the magnetic field has to be modeled properly. Especially the parameterization of the external parts is of major importance since these parts

contribute a significant amount to the total field within Mercury's magnetosphere (Anderson et al. 2011).

Above the planetary surface, the internal part of the field is a current-free magnetic field which can be parameterized by the Gauss representation (Gauß 1839; Glassmeier and Tsurutani 2014). When only data in source-free regions without currents are analyzed, the Gauss representation also yields a parametrization for the external parts of the field (Olsen et al. 2010). But in the vicinity of Mercury significant currents are expected (Anderson et al. 2011, eg); and therefore, the Gauss representation is insufficient due to the current-generated magnetic field outside the planet. In this paper, we construct a parametric Mercury magnetic field model by extending the Gauss representation to the Gauss–Mie representation.

Previous studies parameterized the external parts by making use of a global magnetospheric model, such as the paraboloid model for Mercury's magnetosphere

<figure><img src="image_4.png" /></figure>

<footnote>*Correspondence: s.toepfer@tu-braunschweig.de</footnote>

<footnote n="1">¹ Institut für Theoretische Physik, Technische Universität Braunschweig,
Mendelssohnstraße 3, 38106 Braunschweig, Germany</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
(Alexeev et al. 2008), that has successfully been applied to the analysis of Mercury's internal and external magnetic field (Alexeev et al. 2010; Johnson et al. 2021), as well as the modified Tsyganenko model (Korth et al. 2004). As an alternative for these models the Mie representation is useful for decomposing the magnetic field in current carrying regions in the vicinity of Mercury (Backus 1986; Olsen 1997; Olsen et al. 2010).

The Mie representation, also known as toroidal-poloidal decomposition, is based on the solenoidality of the magnetic field and enables to decompose the field into its toroidal and poloidal parts (Backus 1986; Backus et al. 1996; Kazantsev and Kardakov 2019). This representation has successfully been applied to several problems in space plasma physics, especially for the analysis of the Earth's magnetosphere. For example, Engels and Olsen (1999) used the Mie representation for calculating the magnetic effect caused by a three dimensional current density model. Olsen (1997) applied the Mie representation to reconstruct ionospheric F-currents in the Earth's magnetosphere from MAGSAT data. Furthermore, Bayer et al. (2001) introduced the wavelet Mie representation of the magnetic field to calculate current densities from MAGSAT and CHAMP data. This approach has been expanded on the modeling of the Earth's magnetic field in terms of vector kernel functions (Mayer and Maier 2006). Kosik (1984) constructed a model for the Earth's magnetosphere based on the Mie representation.

In this work, the magnetic field in the vicinity of Mercury is parameterized by a combination of the Gauss and the Mie representation (hereafter, Gauss-Mie representation), based on the works of Backus (1986), Backus et al. (1996) and Olsen (1997), for analyzing Mercury's internal magnetic field. The desired internal Gauss coefficients for the dipole, quadrupole and octupole field are estimated with Capon's method (Capon 1969). The Capon method serves as a powerful inversion method for linear inverse problems and was used by, e.g., Motschmann et al. (1996), Glassmeier et al. (2001), Narita et al. (2006) and Narita (2012) to evaluate spatial spectra of space plasma waves. As shown by Toepfer et al. (2020a) the Capon method can also be used in a generalized way to compare actual measurements with theoretical models. Here, we expand on this method. First of all, the mathematical foundations of the Mie representation, as derived by Backus (1986) and Backus et al. (1996) are revisited. Afterwards, the resulting thin shell approximation is applied to simulated magnetic field data, that are evaluated along the future data points of the Mercury Planetary Orbiter (MPO). This enables the judgement of the expectable inversion results from the data of the magnetometer (Glassmeier et al. 2010; Heyner et al. 2020) on board MPO. Finally, the model is applied to MESSENGER data and the reconstructed Gauss

coefficients are compared with the results of former works by Anderson et al. (2012), Thébault et al. (2018) and Wardinski et al. (2019).

## The Mie representation

The solenoidality of the magnetic field $\underline{B}$ guarantees the existence of a vector potential $\underline{A}$, so that

$$ \underline{B} = \partial_x \times \underline{A} \qquad (1) $$

holds, where $\partial_x$ is the spatial derivative. Using spherical coordinates with radius $r \in [0, \infty)$, longitude $\lambda \in [0, 2\pi]$ and co-latitude angle $\theta \in [0, \pi]$ the vector potential can be separated into its component $\Psi_T \underline{r}$ parallel to $\underline{r} = r \mathbf{e}_r$, where $\mathbf{e}_r$ is the unit vector in radial direction, and its components $\underline{F} \times \underline{r}$ perpendicular to $\underline{r}$, yielding

$$ \underline{A} = \Psi_T \underline{r} + \underline{F} \times \underline{r} + \partial_x \varphi \qquad (2) $$

with the scalar functions $\Psi_T$ and $\varphi$ as well as a vector field $\underline{F}$ (Jacobs 1987; Krause and Rädler 1980). Because of $\partial_x \times \partial_x \varphi = 0$, the function $\varphi$ can be chosen properly so that $\partial_x \times \underline{F} = 0$ holds, and therefore

$$ \underline{F} = \partial_x \Psi_P \qquad (3) $$

with a scalar function $\Psi_P$ without changing the magnetic field (Jacobs 1987; Krause and Rädler 1980). The vector potential results in

$$ \underline{A} = \Psi_T \underline{r} + \partial_x \Psi_P \times \underline{r} \qquad (4) $$

$$ = \Psi_T \underline{r} + \partial_x \times (\Psi_P \underline{r}). \qquad (5) $$

Substituting $\underline{A}$ into Eq. (1) delivers

$$ \underline{B} = \partial_x \times (\Psi_T \underline{r}) + \partial_x \times [\partial_x \times (\Psi_P \underline{r})], \qquad (6) $$

which is called the Mie representation of the magnetic field (Backus 1986; Backus et al. 1996).

The first term on the right hand side in Eq. (6)

$$ \underline{B}_T = \partial_x \times (\Psi_T \underline{r}) \qquad (7) $$

is the toroidal part of the field and the second term

$$ \underline{B}_P = \partial_x \times [\partial_x \times (\Psi_P \underline{r})] \qquad (8) $$

is the poloidal part of $\underline{B}$.

From the definition of $\underline{B}_T$ and $\underline{B}_P$ it is clear that $\underline{B}_T$ is perpendicular to $\underline{B}_P$ and also perpendicular to $\underline{r}$. Therefore, the toroidal part of the field does not have a radial component. Furthermore, poloidal magnetic fields are generated by toroidal currents and vice versa:

$$ \partial_x \times \underline{B}_T = \partial_x \times [\partial_x \times (\underline{r} \Psi_T)] \qquad (9) $$

is a poloidal vector and
$$ \partial_{\underline{x}} \times \underline{B}_P = \partial_{\underline{x}} \times [\partial_{\underline{x}} \times [\partial_{\underline{x}} \times (r\Psi_P)]] = \partial_{\underline{x}} \times [(-\partial_{\underline{x}}^2 \Psi_P)r] \quad (10) $$

is a toroidal vector.

In curl-free regions where especially the poloidal current density $j_p$ vanishes, Ampère's law reads as follows:

$$ \mu_0 j_p = \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Psi_T r)] = -r \partial_{\underline{x}}^2 \Psi_T + \partial_{\underline{x}} \partial_r (r \Psi_T) = 0 \quad (11) $$

or equivalently

$$ \partial_{\underline{x}} \partial_r (r \Psi_T) = r \partial_{\underline{x}}^2 \Psi_T. \quad (12) $$

Therefore, the gradient of $\partial_r(r\Psi_T)$ has only a radial component, leaving us with

$$ \partial_{\lambda} \Psi_T = 0 = \partial_{\theta} \Psi_T, \quad (13) $$

so that the function $\Psi_T$ solely depends on the radial distance from the center ($\Psi_T = \Psi_T(r)$) and thus, the toroidal magnetic field vanishes

$$ \underline{B}_T = \partial_{\underline{x}} \times (\Psi_T r) = (\partial_{\underline{x}} \Psi_T) \times r = \partial_r \Psi_T e_r \times r = 0. \quad (14) $$

On the other hand, when the toroidal current density vanishes

$$ \mu_0 j_T = \partial_{\underline{x}} \times \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Psi_P r)] = \partial_{\underline{x}} \times [-r(\partial_{\underline{x}}^2 \Psi_P) e_r] \quad (15) $$

$$ = \frac{1}{r \sin(\theta)} \partial_{\varphi} [-r \partial_{\underline{x}}^2 \Psi_P] e_{\theta} - \frac{1}{r} \partial_{\theta} [-r \partial_{\underline{x}}^2 \Psi_P] e_{\varphi} \quad (16) $$

$$ = 0 \quad (17) $$

or equivalently

$$ \partial_{\varphi} \partial_{\underline{x}}^2 \Psi_P = 0 \quad (18) $$

and simultaneously

$$ \partial_{\theta} \partial_{\underline{x}}^2 \Psi_P = 0 \quad (19) $$

the function $\partial_{\underline{x}}^2 \Psi_P$ solely depends on the radius. Therefore, the poloidal field

$$ \underline{B}_P = \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Psi_P r)] = -r \partial_{\underline{x}}^2 \Psi_P + \partial_{\underline{x}} \partial_r (r \Psi_P) \neq 0 \quad (20) $$

in general remains finite in current-free regions.

## Relation to the Gauss representation

When magnetic field data in curl-free regions (where $\partial_{\underline{x}} \times \underline{B} = 0$ holds) are analyzed, there exists a scalar potential $\Phi$ so that the magnetic field can be written as

$$ \underline{B} = -\partial_{\underline{x}} \Phi = -\partial_r \Phi e_r - \frac{1}{r} \partial_{\theta} \Phi e_{\theta} - \frac{1}{r \sin(\theta)} \partial_{\lambda} \Phi e_{\lambda}, \quad (21) $$

which is known as the Gauss representation of the magnetic field (Gauß 1839; Glassmeier and Tsurutani 2014). Simultanously, the Mie representation in due consideration that the toroidal magnetic field vanishes in curl-free regions is given by

$$ \underline{B} = \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Psi_P r)] \quad (22) $$

$$ = \frac{1}{r \sin(\theta)} \left[ -\partial_{\theta}(\sin(\theta)\partial_{\theta}\Psi_P) - \frac{1}{\sin(\theta)} \partial_{\lambda}^2 \Psi_P \right] e_r \quad (23) $$

$$ + \frac{1}{r} \partial_{\theta} \partial_{r} (r \Psi_{P}) e_{\theta} + \frac{1}{r \sin(\theta)} \partial_{\lambda} \partial_{r} (r \Psi_{P}) e_{\lambda} \quad (24) $$

Comparison of coefficients with Eq. (21) for $e_{\theta}$ shows that

$$ \frac{1}{r} \partial_r \partial_{\theta} (r \Psi_P) = - \frac{1}{r} \partial_{\theta} \Phi \quad (25) $$

and analogously for $e_{\lambda}$

$$ \frac{1}{r \sin(\theta)} \partial_r \partial_{\lambda} (r \Psi_P) = - \frac{1}{r \sin(\theta)} \partial_{\lambda} \Phi. \quad (26) $$

Consequently, when $\Psi_P$ is known, the scalar potential is given by

$$ \Phi = -\partial_r(r\Psi_P). \quad (27) $$

Comparison of the $e_r$-coefficients delivers

$$ \frac{1}{r \sin(\theta)} \left[ \partial_{\theta}(\sin(\theta)\partial_{\theta}\Psi_P) + \frac{1}{\sin(\theta)} \partial_{\lambda}^2 \Psi_P \right] = \partial_r \Phi \quad (28) $$

or equivalently

$$ \partial_S^2 \Psi_P = \frac{1}{r} \partial_r \Phi, \quad (29) $$

where

$$ \partial_S^2 = \frac{1}{r^2 \sin(\theta)} \partial_{\theta}(\sin(\theta)\partial_{\theta}) + \frac{1}{r^2 \sin^2(\theta)} \partial_{\lambda}^2 \quad (30) $$

is the angular part of the Laplacian. For a given $\Phi$, the function $\Psi_P$ can be found by solving Eqs. (27, 29) simultaneously.

As a consequence, the scalar function $\Psi_P$ and the scalar potential $\Phi$ are equivalent in curl-free regions and the Mie representation transists into the Gauss representation. Thus, the Gauss representation can be understood as a special case of the Mie representation.
## Parameterization of the magnetic field

Assuming that the conductivity of Mercury's mantle is negligibly small (like lunar regolith (Zharkova et al. 2020)), the planetary contribution to the field outside Mercury is purely poloidal. The currents flowing in the magnetosphere generate poloidal and toroidal magnetic fields that superpose with the curl-free planetary magnetic field. To be able to separate the planetary magnetic field out of the measured field and to parameterize it via the Gauss coefficients, a combined parametrization composed of the Mie and the Gauss representation (Gauss-Mie representation), which is based on the works of Backus (1986) and Olsen (1997) is necessary.

Suppose that the magnetic field in the vicinity of Mercury is measured within a spherical shell $S(a, c)$ with inner radius $a > R_M$, where $R_M$ indicates the radius of Mercury, outer radius $c$ and mean radius $b = \frac{1}{2}(a+c)$ as displayed in Fig. 1. The shell can be constructed independently of the orbit's geometry by conceptually covering the orbit of the spacecraft. Furthermore, the shell may include current carrying regions. Although the Mie representation enables us to analyze those currents, we focus on the analysis of Mercury's internal magnetic field.

Due to the underlying geometry, the space around Mercury can be decomposed into three disjoint radial zones:

* points in the region $r < a$ below the shell

* points in the region $r > c$ above the shell

* points in the region $a \le r \le c$ inside the shell layer

Making use of the superposition principle the total magnetic field $\underline{B}$ measured inside the shell layer ($a \le r \le c$) is a composition of the field $\underline{B}_{j \notin [a,c]}$ generated by currents flowing inside the shell and the field $\underline{B}_{j \in [a,c]}$ generated by currents flowing outside the shell. Again considering the underlying geometry, the second part can be divided into an internal part $\underline{B}^i$ resulting from currents flowing in the region $r < a$ and an external part $\underline{B}^e$ resulting from currents flowing in the region $r > c$, so that

$$ \underline{B}_{j \notin [a,c]} = \underline{B}^i + \underline{B}^e. \qquad (31) $$

As $\underline{B}^i$ and $\underline{B}^e$ have their sources beyond the shell they are purely poloidal and especially nonrotational within the shell. Thus, there exist scalar potentials $\Phi^i$ and $\Phi^e$ so that the field can be parameterized in the shell via the Gauss representation resulting in

$$ \underline{B}^i = - \partial_{\underline{x}} \Phi^i \qquad (32) $$

and

$$ \underline{B}^e = - \partial_{\underline{x}} \Phi^e, \qquad (33) $$

where the scalar potentials are given by (Gauß 1839; Glassmeier and Tsurutani 2014)

$$ \Phi^i = R_M \sum_{l=1}^{\infty} \sum_{m=0}^{l} \left(\frac{R_M}{r}\right)^{l+1} [g_l^m \cos(m\lambda) + h_l^m \sin(m\lambda)] P_l^m(\cos(\theta)) \qquad (34) $$

and

$$ \Phi^e = R_M \sum_{l=1}^{\infty} \sum_{m=0}^{l} \left(\frac{r}{R_M}\right)^l [q_l^m \cos(m\lambda) + s_l^m \sin(m\lambda)] P_l^m(\cos(\theta)). \qquad (35) $$

The expansion coefficients $g_l^m$ and $h_l^m$ are the internal Gauss coefficients, the coefficients $q_l^m$ and $s_l^m$ are the external Gauss coefficients and $P_l^m$ are the Schmidt normalized Legendre polynomials of degree $l$ and order $m$. Since Mercury's internal magnetic field is dominated by the internal dipole, quadrupole and octupole fields, the series expansions in Eqs. (34, 35) will be truncated at the degree $l=3$ for the practical application.

It should be noted that the internal field $\underline{B}^i$ is canonically described in a Mercury-Body-Fixed co-rotating coordinate system (MBF), whereas the external field $\underline{B}^e$ is canonically described in a Mercury-Solar-Orbital system (MSO) with the x-axis orientated towards the sun, the z-axis orientated parallel to the rotation axis, i.e. antiparallel to the internal dipole moment, and the y-axis completes the right-handed system (Heyner et al. 2020). Let $\underline{x}_{MSO} = (x_{MSO}, y_{MSO}, z_{MSO})^T$ define the MSO coordinate system and let $\underline{x}_{MBF} = (x_{MBF}, y_{MBF}, z_{MBF})^T$ be the coordinates of the co-rotating MBF system. Then, the internal parts of the field are given by

$$ \underline{B}^i(\underline{x}_{MBF}) = -\partial_{\underline{x}_{MBF}} \Phi^i(\underline{x}_{MBF}) = \underline{H}^i(\underline{x}_{MBF}) \underline{g}^i, \qquad (36) $$

whereas the external parts are described in the MSO system, i.e.

$$ \underline{B}^e(\underline{x}_{MSO}) = -\partial_{\underline{x}_{MSO}} \Phi^e(\underline{x}_{MSO}) = \underline{H}^e(\underline{x}_{MSO}) \underline{g}^e, \qquad (37) $$

where the terms of the series expansion are arranged into the matrices $\underline{H}^i$ and $\underline{H}^e$ and the corresponding Gauss coefficients are summarized into the vectors $\underline{g}^i$ and $\underline{g}^e$.

The co-rotating MBF system can be transformed into the MSO system via

$$ \underline{x}_{MSO} = \underline{H} \underline{x}_{MBF}, \qquad (38) $$
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Cross section of the spherical shell S(a, c) with inner radius a, outer radius c and mean radius b = (a + c)/2. The inner sphere with radius R_M symbolizes Mercury</figcaption></figure>

where A describes a rotation matrix around the z-axis depending on the angular velocity of Mercury's self-rotation measured within the MSO system.

For the practical application it is convenient to describe both parts of the field in one coordinate system, for example the MSO system. The transformed data are given by

$$ \underline{B}(\underline{x}_{\text{MSO}}) = \underline{A}\underline{B}^i(\underline{x}_{\text{MBF}}) + \underline{B}^e(\underline{x}_{\text{MSO}}) \quad (39) $$

$$ = \underline{A}\underline{H}^i(\underline{x}_{\text{MBF}})\underline{g}^i + \underline{H}^e(\underline{x}_{\text{MSO}})\underline{g}^e \quad (40) $$

in the MSO system.

Since for the first validation the model will be applied to simulated magentic field data, it is useful to match the coordinate system of the parametrization with the coordinate system of the simulation. Therefore, in the following all parts of the magnetic field are described in a Mercury-Body-Fixed anti-solar orientated coordinate system (MASO) with coordinates $\underline{x} = (x, y, z)^T$, where the x-axis is orientated towards the nightside of Mercury (away from the sun), the z-axis is orientated antiparallel to the internal dipole moment and the y-axis completes the right-handed system, so that

$$ \underline{x} = \begin{bmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end{bmatrix} \underline{x}_{\text{MSO}}. \quad (41) $$

As already mentioned in the introduction ("Introduction" section) there is no current-free shell-like region around Mercury (Olsen et al. 2010). The currents flowing in the shell generate toroidal $\underline{B}_T^{sh}$ and poloidal $\underline{B}_P^{sh}$ magnetic fields which superpose with $\underline{B}^i$ and $\underline{B}^e$. Thus, the total measured field within the shell is composed of four parts given by

$$ \underline{B} = \underline{B}_{\underline{l} \notin [a,c]} + \underline{B}_{\underline{l} \in [a,c]} = \underline{B}^i + \underline{B}^e + \underline{B}_T^{sh} + \underline{B}_P^{sh} \quad (42) $$

$$ = -\partial_{\underline{x}} \Phi^i - \partial_{\underline{x}} \Phi^e + \partial_{\underline{x}} \times \left[ \partial_{\underline{x}} \times \left( \underline{r} \Psi_P^{sh} \right) \right] + \partial_{\underline{x}} \times \left( \underline{r} \Psi_T^{sh} \right), \quad (43) $$

where each part of the field is described either by a scalar potential $\Phi^i$, $\Phi^e$ or a scalar function $\Psi_P^{sh}$, $\Psi_T^{sh}$.

The scalar potentials $\Phi^i$ and $\Phi^e$ are already parameterized by the Gauss coefficients. In the following, a proper parameterization for the scalar functions $\Psi_T^{sh}$ and $\Psi_P^{sh}$ is required. Because of the underlying spherical geometry it is straightforward to expand the functions into spherical harmonics

$$ \Psi_T^{sh} = \sum_{l=1}^{\infty} \sum_{m=0}^{l} [a_l^m(r) \cos(m\lambda) + b_l^m(r) \sin(m\lambda)] P_l^m(\cos(\theta)) \quad (44) $$

and

$$ \Psi_P^{sh} = R_M \sum_{l=1}^{\infty} \sum_{m=0}^{l} [c_l^m(r) \cos(m\lambda) + d_l^m(r) \sin(m\lambda)] P_l^m(\cos(\theta)), \quad (45) $$

where $a_l^m(r)$, $b_l^m(r)$, $c_l^m(r)$ and $d_l^m(r)$ are the expansion coefficients which in general depend on the radius r and again $P_l^m$ are the Schmidt normalized Legendre polynomials. Since the toroidal and poloidal fields can be locally generated by currents flowing in the shell the radial dependences of the fields and the expansion coefficients, respectively, are unknown.

### Series expansion of the coefficients

Since the exact radial dependence of the expansion coefficients $a_l^m(r)$ and $b_l^m(r)$ is unknown, it is useful to expand these functions into a Taylor series in the vicinity of the mean radius b of the shell. Within this series expansion it is advisable not to incorporate the effect of all components of the poloidal current density to the toroidal magnetic field at once. Here, we first concentrate on the radial component of the current density and consider the horizontal components in higher orders of the Taylor series.

The toroidal magnetic field $\underline{B}_T^{sh}$ is generated by poloidal currents $\underline{j}_p$ (cf. Eq. (11)). Ampère's law yields
$$ \mu_0 j_P = \partial_x \times \underline{B}_T^{sh} = \partial_x \times \left[ \partial_x \times \left( r \Psi_T^{sh} \right) \right] \quad (46) $$

$$ = \frac{1}{r \sin(\theta)} \left[ -\partial_{\theta} \left( \sin(\theta) \partial_{\theta} \Psi_T^{sh} \right) - \frac{1}{\sin(\theta)} \partial_{\lambda}^2 \Psi_T^{sh} \right] e_r \quad (47) $$

$$ + \frac{1}{r} \partial_{\theta} \partial_{r} (r \Psi_{T}^{sh}) e_{\theta} + \frac{1}{r \sin(\theta)} \partial_{\lambda} \partial_{r} (r \Psi_{T}^{sh}) e_{\lambda}. \quad (48) $$

The components of the horizontal $e_{\theta}$- and $e_{\lambda}$-direction are proportional to $\partial_r (r \Psi_T^{sh})$. Therefore, the ansatz

$$ a_l^m(r) = \frac{R_M}{r} \left[ a_l^m + a_l^{m'} \rho + \mathcal{O}(\rho^2) \right] \quad (49) $$

and

$$ b_l^m(r) = \frac{R_M}{r} \left[ b_l^m + b_l^{m'} \rho + \mathcal{O}(\rho^2) \right], \quad (50) $$

where $\rho = \frac{r-b}{R_M}$ and $a_l^m, a_l^{m'}, b_l^m, b_l^{m'}$ are constants for each pair of $l$ and $m$, is utilized. In the first order of the Taylor series expansion in the vicinity of the mean radius $b$, where $\Psi_T^{sh} \sim \frac{1}{r}$, the horizontal components of $j_P$ vanish and only the contributions of the radial currents driving the toroidal magnetic field are considered (Olsen 1997). Using higher orders of the Taylor series, also the contributions of the horizontal components of $j_P$ to the toroidal magnetic field in the vicinity of the mean radius $b$ can be incorporated.

The scalar function of the toroidal magnetic field results in

$$ \Psi_T^{sh} = \frac{R_M}{r} \sum_{l=1}^{\infty} \sum_{m=0}^{l} \left[ \alpha_l^m + \alpha_l^{m'} \rho + \mathcal{O}(\rho^2) \right] P_l^m(\cos(\theta)), \quad (51) $$

where $\alpha_l^m = a_l^m \cos(m\lambda) + b_l^m \sin(m\lambda)$ and $\alpha_l^{m'} = a_l^{m'} \cos(m\lambda) + b_l^{m'} \sin(m\lambda)$ (Olsen 1997). Thereby, each order of the Taylor series is linked with an additional set of expansion coefficients $a_l^m, b_l^m, a_l^{m'}, b_l^{m'}$ and so on which can be reconstructed from the data in analogy to the Gauss coefficients.

From a mathematical point of view the scalar function $\Psi_P^{sh}$ of the poloidal magnetic field $\underline{B}_P^{sh}$ can be parametrized analogously to the toroidal counterpart. But within the reconstruction procedure the poloidal fields that are generated by toroidal currents flowing inside the shell cannot be distinguished from the internally and externally driven poloidal fields, since these fields follow the same topological structure. But when the half thickness of the shell, defined by $h = (c-a)/2$ is smaller than the length scale on which the toroidal currents change in radial direction, the shell is called a thin shell and the scalar function $\Psi_P^{sh}$ of the poloidal field $\underline{B}_P^{sh}$ vanishes within

this thin shell approximation (Backus 1986; Backus et al. 1996) as illustrated in the following section.

### The thin shell approximation

The thin shell approximation (Backus 1986; Backus et al. 1996) finally allows the separation of the poloidal field into its internal and external contributions.

Conferring to Eq. (6), the Mie representation for the magnetic field in the whole space $\mathbb{R}^3$ is given by

$$ \underline{B} = \partial_{\underline{x}} \times (\Psi_T \underline{r}) + \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Psi_P \underline{r})]. \quad (52) $$

Following Ampère's law the current density $j$ is also sole- noidal and can as well be parameterized via the Mie rep- resentation resulting in

$$ \mu_0 j = \partial_{\underline{x}} \times (\Gamma_T \underline{r}) + \partial_{\underline{x}} \times [\partial_{\underline{x}} \times (\Gamma_P \underline{r})] \quad (53) $$

with related scalar functions $\Gamma_T$ and $\Gamma_P$.

Since the poloidal part of the current density corre- sponds with the curl of the toroidal magnetic field, the comparision of Eq. (53) with Eq. (9) shows that the scalar function $\Psi_T$ and $\Gamma_P$ are the same

$$ \Psi_T = \Gamma_P. \quad (54) $$

Analogously, the toroidal part of the current density cor- responds with the curl of the poloidal magnetic field and the comparision of Eq. (53) with Eq. (10) shows that the functions $\Psi_P$ and $\Gamma_T$ are related via

$$ \partial_{\underline{x}}^2 \Psi_P = -\Gamma_T, \quad (55) $$

so that the function $\Psi_P$ is given by the Green's function method

$$ \Psi_P(\underline{r}) = \frac{1}{4\pi} \int_{\mathbb{R}^3} \frac{\Gamma_T(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3\underline{r}'. \quad (56) $$

Due to the underlying geometry the toroidal current density $j_T$ flowing in the whole space can be written as the sum of the toroidal currents $j_T^i$ flowing in the region $r < a$, the toroidal currents $j_T^e$ flowing in the region $r > c$ and the toroidal currents $j_T^{sh}$ flowing inside the spherical shell in the region $a \le r \le c$, so that

$$ j_T^i = j_T^i + j_T^e + j_T^{sh}. \quad (57) $$

Thereby, the Mie representation of each part is given by

$$ j_T^i = \partial_{\underline{x}} \times (\Gamma_T^i \underline{r}) \quad (58) $$

$$ j_T^e = \partial_{\underline{x}} \times (\Gamma_T^e \underline{r}) \quad (59) $$
$$j_{-T}^{sh} = \partial_{\underline{x}} \times (\Gamma_{T}^{sh} \underline{r}) \qquad (60)$$

with $\Gamma_T^i = \Gamma_T \chi_{[r<a]}$, $\Gamma_T^{sh} = \Gamma_T \chi_{[a<r<c]}$, where

$$\chi_I(x) = \begin{cases} 1, & x \in I \\ 0, & x \notin I \end{cases} \qquad (61)$$

is the indicator function of the interval $I$. Using this segmentation the scalar function of the poloidal field can be rewritten as

$$\Psi_P(\underline{r}) = \frac{1}{4\pi} \int_{V(r'<a)} \frac{\Gamma_T^i(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3 r' \\ + \frac{1}{4\pi} \int_{V(r'>a)} \frac{\Gamma_T^e(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3 r' \\ + \frac{1}{4\pi} \int_{S(a,c)} \frac{\Gamma_T^{sh}(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3 r' \qquad (62)$$

Thus, the part of the scalar function that corresponds to the poloidal magnetic field which is generated inside the shell is given by

$$\Psi_P^{sh}(\underline{r}) = \frac{1}{4\pi} \int_{S(a,c)} \frac{\Gamma_T^{sh}(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3 r'. \qquad (63)$$

Since $2b = a + c$ and $2h = c - a$ the bounds of integration can be rewritten as $a = b - h$ and $c = b + h$, so that

$$\Psi_P^{sh}(\underline{r}) = \frac{1}{4\pi} \int_{S(b-h,b+h)} \frac{\Gamma_T^{sh}(\underline{r}')}{|\underline{r} - \underline{r}'|} d^3 r'. \qquad (64)$$

Analogously to the scalar function $\Psi_T$, the function $\Gamma_T$ can be expanded into a Taylor series in the vicinity of the mean radius $b$, resulting in

$$\Gamma_T^{sh}(\underline{r}') = \sum_{n=0}^{\infty} \frac{1}{n!} \partial_{r'}^n \Gamma_T^{sh}|_{r'=b} (r' - b)^n. \qquad (65)$$

$$|\Psi_P^{sh}(\underline{r})| \le \frac{1}{4\pi} \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{h}{L}\right)^n \max_{\theta', \lambda'} \{|\Gamma_T^{sh}(b)|\} \int_{S(b-h,b+h)} \frac{1}{|\underline{r} - \underline{r}'|} d^3 r'. \qquad (73)$$

Substituting the Taylor series into the function $\Psi_P^{sh}$ delivers

$$\Psi_P^{sh}(\underline{r}) = \frac{1}{4\pi} \sum_{n=0}^{\infty} \frac{1}{n!} \int_{S(b-h,b+h)} \partial_{r'}^n \Gamma_T^{sh}|_{r'=b} \frac{(r' - b)^n}{|\underline{r} - \underline{r}'|} d^3 r'. \qquad (66)$$

For the further evaluation of the integral it is assumed that the derivatives of the toroidal currents with respect to $r$ are bounded, i.e., there exists a constant $L > 0$, so that

$$L^n |\partial_{r'}^n j_{-T}| \le |j_{-T}| \qquad (67)$$

for $n \in \mathbb{N}$. Thus, $L$ represents the length scale on which the toroidal currents change in radial direction. From

$$\begin{aligned} j_{-T} &= \partial_{\underline{x}} \times (\Gamma_T \underline{r}) \\ &= \frac{1}{\sin(\theta)} \partial_{\lambda} \Gamma_T e_{\theta} - \partial_{\theta} \Gamma_T e_{\lambda} \\ &= \left( \frac{1}{\sin(\theta)} e_{\theta} \partial_{\lambda} - e_{\lambda} \partial_{\theta} \right) \Gamma_T \end{aligned} \qquad (68)$$

and therefore

$$\partial_{r'}^n j_{-T} = \left( \frac{1}{\sin(\theta)} e_{\theta} \partial_{\lambda} - e_{\lambda} \partial_{\theta} \right) \partial_{r'}^n \Gamma_T \qquad (69)$$

it follows that

$$L^n |\partial_{r'}^n \Gamma_T^{sh}| \le |\Gamma_T^{sh}| \qquad (70)$$

or equivalently

$$|\partial_{r'}^n \Gamma_T^{sh}| \le \frac{|\Gamma_T^{sh}|}{L^n}. \qquad (71)$$

Since $r' \in [b-h, b+h]$, each summand within the Taylor series can be estimated upwards via

$$|\partial_{r'}^n \Gamma_T^{sh}(b)| \cdot |r' - b|^n \le \frac{|\Gamma_T^{sh}(b)|}{L^n} \cdot h^n \qquad (72)$$

delivering for Eq. (66)
The integral in Eq. (73) may be evaluated for any $\underline{r} \in \mathbb{R}^3$ but we restrict to $\underline{r}$ inside the shell as only there the magnetic field is measured. Then, the remaining integral results in

density. In contrast to the toroidal magnetic field, the poloidal magnetic field is also measurable in current-free region. Thus, the poloidal field is a superposition of fields generated by currents flowing inside and outside

$$ \int_{S(b-h,b+h)} \frac{1}{|\underline{r} - \underline{r}'|} d^3\underline{r}' = \frac{4\pi}{3} \frac{1}{r} [r^3 - (b-h)^3] + 2\pi [(b+h)^2 - r^2], \quad (74) $$

utilizing that the coordinate system can be chosen properly so that $\theta'$ defines the angle between $\underline{r}$ and $\underline{r}'$. Therefore

the shell as well as currents flowing within the shell. This superposition is verified in Eq. (62). Therefore, the amount of the poloidal field generated by currents

$$ |\Psi_P^{sh}(\underline{r})| \le \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{\underline{h}}{L}\right)^n \max_{\theta', \lambda'} \left\{ |\Gamma_T^{sh}(b)| \right\} \left\{ \frac{1}{3} \left[ r^2 - \frac{(b-h)^3}{r} \right] + \frac{1}{2} [(b+h)^2 - r^2] \right\}. \quad (75) $$

For all $r \in [b-h, b+h]$, the function

$$ f(r) = \frac{1}{3} \left[ r^2 - \frac{(b-h)^3}{r} \right] + \frac{1}{2} [(b+h)^2 - r^2] \quad (76) $$

is non-negative and reaches its maximum value at $r = b-h$ with the related function value $f(b-h) = 2bh$. Therefore, an upper bound for $\Psi_P^{sh}$ can finally be estimated as

$$ \begin{aligned} |\Psi_P^{sh}(\underline{r})| & \le 2bh \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{\underline{h}}{L}\right)^n \max_{\theta', \lambda'} \{|\Gamma_T^{sh}(b)|\} \\ & = 2bL \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{\underline{h}}{L}\right)^{n+1} \max_{\theta', \lambda'} \{|\Gamma_T^{sh}(b)|\}. \end{aligned} \quad (77) $$

When $h \ll L$, the spherical shell is called a thin shell (Backus 1986; Backus et al. 1996) and the scalar function $\Psi_P^{sh}$ of the poloidal magnetic field $\underline{B}_P^{sh}$ vanishes. The scalar function $\Psi_T^{sh}$, however, remains finite for all $h$ as shown in Appendix A. Thus, if the shell may be regarded as thin, then the contribution of the toroidal currents in this shell to the poloidal magnetic field may be neglected. The poloidal magnetic field is mainly driven by currents beyond the shell. The contribution of the poloidal currents in this shell to the toroidal magnetic field may not be neglected.

From a first point of view the thin shell approximation is not an intuitive approximation. Considering the above presented nature of poloidal and toroidal magnetic fields it can be understood as follows:

The toroidal magnetic field only exists within current carrying regions (cf. Eq. 14) and thus, it is measurable only within these regions. Therefore, the spatial extent of the regions where the poloidal currents flow does not influence the strength of the toroidal magnetic field. It solely depends on the strength of the poloidal current

flowing within the shell has to be compared with the amount of the internal/external contributions. Furthermore, the poloidal field does not solely depend on the strength of the toroidal current density, since for the evaluation of the integrals also the volume where the current density flows is vital. Thus, a small toroidal current density that flows within a large volume outside/ inside the shell can have a larger contribution to the field measured within the shell than a stronger current flowing within the thin shell.

## Application of the thin shell approximation

The thin shell approximation is applied to parameterize the magnetic field in the vicinity of Mercury by the Gauss-Mie representation to reconstruct Mercury's internal magnetic field. The internal and external scalar potentials $\Phi^i$ and $\Phi^e$ are expanded into spherical harmonics up to third degree and order, representing the internal and external dipole, quadrupole and octupole fields. The scalar function $\Psi_T^{sh}$ of the toroidal magnetic field $\underline{B}_T^{sh}$ is expanded into spherical harmonics up to second degree and order and additionally into a first order Taylor series for the radius. The scalar function $\Psi_P^{sh}$ of the poloidal field $\underline{B}_P^{sh}$ is neglected within the thin shell approximation. Therefore, the total field is parameterized by 46 expansion coefficients, i.e. Gauss internal dipole ($g_1^0, g_1^1, h_1^1$), Gauss internal quadrupole ($g_2^0, g_2^1, h_2^1, g_2^2, h_2^2$), Gauss internal octupole ($g_3^0, g_3^1, h_3^1, g_3^2, h_3^2, g_3^3, h_3^3$), Gauss external dipole ($q_1^0, q_1^1, s_1^1$), Gauss external quadrupole ($q_2^0, q_2^1, s_2^1, q_2^2, s_2^2$), Gauss external octupole ($q_3^0, q_3^1, h_3^1, q_3^2, s_3^2, h_3^3, s_3^3$), toroidal coefficients ($a_1^0, a_1^1, b_1^1, a_2^0, a_2^1, b_2^1, a_2^2, b_2^2, a_1'^0, a_1'^1, b_1'^1, a_2'^0, a_2'^1, b_2'^1, a_2'^2, b_2'^2$). These 46 coefficients are estimated with Capon's method (Capon 1969; Toepfer et al. 2020a, b). The Capon method and the underlying model are tested against simulated data and MESSEN-GER in situ data around Mercury.
## Hybrid simulation of Mercury's magnetosphere

For the first application of the thin shell approximation simulated magnetic field data are analyzed. The magnetic field resulting from the plasma interaction of Mercury with the solar wind is simulated with the hybrid code AIKEF (Müller et al. 2011), that has successfully been applied to several problems in Mercury's plasma interaction, (Exner et al. 2018, 2020, e.g.). The internal Gauss coefficients $g_1^0 = -190$ nT (dipole field), $g_2^0 = -78$ nT (quadrupole field) and $g_3^0 = -20$ nT (octupole field) (Anderson et al. 2012; Thébault et al. 2018; Wardinski et al. 2019) are implemented in the simulation code. The interplanetary magnetic field with a magnitude of $B_{IMF} = 20$ nT is orientated along the vector $(x, y, z)^T = (0.0, 0.43, 0.9)^T$ in the MASO frame. The solar wind velocity of $v_{sw} = 400$ km/s points along the x-axis (away from the Sun) and the solar wind proton density number was chosen to $n_{sw} = 30$ cm$^{-3}$. The resulting magnitude of the magnetic field and the corresponding current density in the x-z-plane are displayed in Figs. 2, 3.

The internal dipole field dominates the geometry of Mercury's magnetosphere. Yet the quadrupole field in terms of the apparently shifted dipole field is visible. The influence of the octupole field is not clearly noticeable on the field lines in the figure. Furthermore, the distribution of the simulated current density shows that there exist no completely current-free region around Mercury.

## Reconstruction of the Gauss coefficients from simulated data

For the reconstruction of the internal Gauss coefficients implemented in the simulation code, first of all, magnetic field data at a distance of $0.2 R_M$ from Mercury's surface are evaluated. The data are retrieved along meridional circular orbits around Mercury. The orbital plane is rotated about the rotation axis (z-axis) from $-50°$ (afternoon/post-midnight sector) over $0°$ (noon/midnight, x-z-plane) to $50°$ (morning/pre-midnight sector). For this synthetically generated ideal case in terms of the thin shell approximation, the spherical shell that covers the circular orbits has a vanishing thickness $h = 0$ so that the application of the thin shell approximation is surely valid. The reconstructed internal Gauss coefficients are listed in Table 1. The optimal diagonal loading parameter which determines how the data are weighted within the application of Capon's method results in $\sigma_{opt.} \approx 1000$ nT (Toepfer et al. 2020b).

The deviation $|g_{int} - g_{C}^{int}|$ between the reconstructed internal coefficients $g_C^{int}$ and the implemented internal coefficients $g_{int}^{\text{re}}$ results in 4.0 nT, i.e. $|g_{int}^{\text{re}} - g_{C}^{int}|/|g_{int}^{\text{re}}| \approx 1.9\%$ and thus, the implemented

coefficients are reconstructed from the data with good precision.

Concerning the BepiColombo mission the generated circular orbits are idealized cases which are not realizable in practice. To investigate the applicability of the thin shell approximation for elliptical orbits we analyze the magnetic field data along the prospective orbits of MPO. The orbits are generated in analogy to the circular orbits, i.e., with the same longitudinal extent. Along this trajectories the distances of the data points from Mercury's surface vary from $0.12 R_M$ up to $0.6 R_M$ resulting in a shell thickness of $2h \approx 0.48 R_M$. Although the shell is now much thicker the thin shell approximation works successfully. The reconstructed Gauss coefficients are displayed in Table 2. The optimal diagonal loading parameter for the application of Capon's method results in $\sigma_{opt.} \approx 1000$ nT.

The deviation between the reconstructed and the implemented internal coefficients results in 4.1 nT, i.e. 1.9% and thus, these values agree with the coefficients reconstructed from the data evaluated along the circular orbits. Since the data are evaluated along the MPO orbits, it is expectable that also Mercury's internal octupole field will be reconstructed with good precision from the data of the magnetometer on board MPO.

It should be noted that the extension of the underlying model by the parameterization of the external parts of the field using the Gauss-Mie representation improves the results significantly. When only the internal parts

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Simulated magnitude of the current density *j* in the x-z-plane. The grey circle of radius 1 *R*<sub>M</sub> symbolizes Mercury</figcaption></figure>

of the field *B*<sup>*i*</sup> are considered in the model, the devia-
tion between the implemented and the reconstructed
coefficients results in 29.7 nT or 14.4%, respectively (cf.
Table 2). Additional parameterization of the external
poloidal fields *B*<sup>*e*</sup> using the external scalar potential *Φ*<sup>*e*</sup>
yields a deviation of 10.0 nT or 4.9%, respectively. Thus,
for the analysis of Mercury's internal magnetic field, the
Gauss-Mie representation is a suitable alternative to the
application of global magnetospheric models.

**Validity of the thin shell approximation**

The coefficients reconstructed from the data evalu-
ated along the planned *MPO* trajectories are basically
in agreement with that from the data along the *circular*
*orbits* although the shell covering the MPO orbits has a
finite thickness.

To investigate the limits of the thin shell approximation
within the hybrid simulation, the thickness of the shell is
increased incrementally. Referring to Eq. (67), the length
scale *L* is estimated via radial variation of the current
distribution

$$
L = \frac{|\dot{j}_H|}{|\partial_r \dot{j}_H|} \qquad (78)
$$

for the horizontal currents $\dot{j}_H = \dot{j} - (\dot{j} \cdot \mathbf{e}_r) \mathbf{e}_r$ which are
used as a proxy for the toroidal currents at each point
along the orbit. Since $L$ represents a local quantity that
varies for each point along the orbit, whereas the half

<table><caption>**Table 1** Implemented and reconstructed Gauss coefficients for the dipole, quadrupole and octupole field resulting from the simulated data along the circular orbits at a distance of 0.2 *R*<sub>M</sub> from Mercury's surface</caption>
  <thead>
    <tr>
      <th>Gauss coefficient</th>
      <th>Input in nT</th>
      <th>Output Capon in nT</th>
    </tr>
  </thead>
  <tbody>
    <tr>
      <td>g<sup>0</sup><sub>1</sub></td>
      <td>-190.0</td>
      <td>-191.2</td>
    </tr>
    <tr>
      <td>g<sup>1</sup><sub>1</sub></td>
      <td>0.0</td>
      <td>-1.9</td>
    </tr>
    <tr>
      <td>h<sup>1</sup><sub>1</sub></td>
      <td>0.0</td>
      <td>-1.7</td>
    </tr>
    <tr>
      <td>g<sup>0</sup><sub>2</sub></td>
      <td>-78.0</td>
      <td>-76.0</td>
    </tr>
    <tr>
      <td>g<sup>1</sup><sub>2</sub></td>
      <td>0.0</td>
      <td>1.3</td>
    </tr>
    <tr>
      <td>h<sup>1</sup><sub>2</sub></td>
      <td>0.0</td>
      <td>-0.5</td>
    </tr>
    <tr>
      <td>g<sup>2</sup><sub>2</sub></td>
      <td>0.0</td>
      <td>-0.5</td>
    </tr>
    <tr>
      <td>h<sup>2</sup><sub>2</sub></td>
      <td>0.0</td>
      <td>-0.2</td>
    </tr>
    <tr>
      <td>g<sup>0</sup><sub>3</sub></td>
      <td>-20.0</td>
      <td>-20.3</td>
    </tr>
    <tr>
      <td>g<sup>1</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>0.5</td>
    </tr>
    <tr>
      <td>h<sup>1</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>0.4</td>
    </tr>
    <tr>
      <td>g<sup>2</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>-0.1</td>
    </tr>
    <tr>
      <td>h<sup>2</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>0.9</td>
    </tr>
    <tr>
      <td>g<sup>3</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>1.0</td>
    </tr>
    <tr>
      <td>h<sup>3</sup><sub>3</sub></td>
      <td>0.0</td>
      <td>-0.2</td>
    </tr>
  </tbody>
</table>

thickness *h* is a global quantity, the set of resulting length
scales is averaged over the number of points along the
orbit resulting in the mean length scale ⟨*L*⟩. It should be
noted that also the poloidal currents have horizontal
components. Thus, the estimation of the length scale *L*
with the horizontal currents *j*<sub>**H**</sub> as a proxy for the toroidal
currents is not exact, but it is sufficient for a qualitative
discussion.

The deviations between the ideal coefficients imple-
mented in the simulation and the reconstructed coeffi-
cients for varying values of *h*/⟨*L*⟩ are displayed in Table 3.

When *h* approaches ⟨*L*⟩ the deviations are greater
than 12 nT. This deviation is of the same order when the
parameterization is restricted to the Gauss representa-
tion. But for the data points along the MPO orbits, where
*h*/⟨*L*⟩ ≈ 0.47, the application of thin shell approximation
is justified.

Since a shell of thickness 2h ≈ 0.48 R<sub>M</sub> is called thin,
the name *thin shell* is misleading, although this naming
has been adopted within the literature. Conferring to the
global length scale of 1 R<sub>M</sub> the shell of thickness 0.48 R<sub>M</sub>
is not thin, but compared with the current length scale
*L* it is. Therefore, the term *thin* has to be understood in a
mathematical sense.
<table><caption>**Table 2** Implemented and reconstructed Gauss coefficients for the dipole, quadrupole and octupole field resulting from the simulated data along the future MPO orbits. For comparison in the last columns the reconstructed Gauss coefficients resulting from the sole parameterization of the internal parts $B^1$ (Gauss internal, 15 coefficients considered within the model) as well as the coefficients resulting from the parameterization of the internal $B^1$ and external poloidal fields $B^2$ (Gauss internal external, 30 coefficients considered within the model) are presented</caption><thead><tr><th rowspan="2">Gauss<br>coefficient</th><th rowspan="2">Input in nT</th><th colspan="3">Output capon in nT</th></tr><tr><th>Gauss-Mie</th><th>Gauss internal</th><th>Gauss<br>internal<br>external</th></tr></thead><tbody><tr><td>$g_1^0$</td><td>-190.0</td><td>-190.9</td><td>-188.3</td><td>-190.2</td></tr><tr><td>$g_1^1$</td><td>0.0</td><td>-2.1</td><td>0.8</td><td>-0.8</td></tr><tr><td>$h_1^1$</td><td>0.0</td><td>1.7</td><td>-18.6</td><td>-5.7</td></tr><tr><td>$g_2^0$</td><td>-78.0</td><td>-76.1</td><td>-78.5</td><td>-76.0</td></tr><tr><td>$g_2^1$</td><td>0.0</td><td>0.7</td><td>5.6</td><td>3.3</td></tr><tr><td>$h_2^1$</td><td>0.0</td><td>0.0</td><td>0.7</td><td>-0.1</td></tr><tr><td>$g_2^2$</td><td>0.0</td><td>0.3</td><td>0.5</td><td>0.4</td></tr><tr><td>$h_2^2$</td><td>0.0</td><td>0.4</td><td>4.1</td><td>4.5</td></tr><tr><td>$g_3^0$</td><td>-20.0</td><td>-19.5</td><td>-14.3</td><td>-20.8</td></tr><tr><td>$g_3^1$</td><td>0.0</td><td>0.6</td><td>-0.5</td><td>0.4</td></tr><tr><td>$h_3^1$</td><td>0.0</td><td>-0.1</td><td>0.9</td><td>-0.2</td></tr><tr><td>$g_3^2$</td><td>0.0</td><td>-1.9</td><td>15.2</td><td>-2.1</td></tr><tr><td>$h_3^2$</td><td>0.0</td><td>0.2</td><td>0.4</td><td>-0.2</td></tr><tr><td>$g_3^3$</td><td>0.0</td><td>0.8</td><td>-3.0</td><td>-0.4</td></tr><tr><td>$h_3^3$</td><td>0.0</td><td>-0.4</td><td>14.4</td><td>5.1</td></tr><tr><td>$q_1^0$</td><td>-</td><td>-34.0</td><td>-</td><td>-35.9</td></tr><tr><td>$q_1^1$</td><td>-</td><td>-8.3</td><td>-</td><td>-7.2</td></tr><tr><td>$s_1^1$</td><td>-</td><td>-6.4</td><td>-</td><td>-9.1</td></tr></tbody></table>

## Reconstruction of the Gauss coefficients from MESSENGER data

The Gauss-Mie representation has successfully been validated for the simulated data. For the reconstruction of the Gauss coefficients from the MESSENGER data only data points in the northern hemisphere within a distance of $r \le 1.5 R_M$ and $x > -0.4 R_M$ from Mercury's surface can be considered because the orbits do not cover the southern hemisphere properly.

For the reconstruction of Mercury's internal magnetic field, the data from nine pairs of MESSENGER orbits with different orientations of the periapsis between 10. August 2012 and 14. July 2014 are analyzed. The combination of the orbits improves the model resolution compared to the analysis of single orbits. A discussion of the model resolution, as provided by Connerney (1981), can be found in the Appendix B. As a proof of concept, only

<table><caption>**Table 3** Deviation between the implemented and the reconstructed Gauss coefficients for varying ratios of the half thickness $h$ to the mean current length scale $\langle L \rangle$</caption><thead><tr><th>$h$ in $R_M$</th><th>$\langle L \rangle$ in $R_M$</th><th>$h/\langle L \rangle$</th><th>Deviation in nT</th><th>Deviation in %</th></tr></thead><tbody><tr><td>0.24</td><td>0.51</td><td>0.47</td><td>4.1</td><td>1.9</td></tr><tr><td>0.32</td><td>0.49</td><td>0.65</td><td>4.8</td><td>2.3</td></tr><tr><td>0.43</td><td>0.46</td><td>0.94</td><td>12.8</td><td>6.2</td></tr><tr><td>0.62</td><td>0.42</td><td>1.5</td><td>13.5</td><td>7.0</td></tr><tr><td>1.09</td><td>0.37</td><td>2.9</td><td>16.8</td><td>8.1</td></tr></tbody></table>

a small subset of the whole MESSENGER data set is analyzed. In the case of a small and noisy data set the performance of the estimator can be improved by separating the data set into several subsets (Meir 1994). Therefore, the Gauss coefficients are reconstructed for each pair of orbits and the results are averaged over the nine pairs. Outliers were not included within the averaging. The resulting mean values are listed in Table 4. The reconstructed coefficients for each pair of orbits and the standard deviations of the mean values are listed in Table 5 of Appendix C.

The mean values of the reconstructed internal Gauss coefficients and the external dipole coefficients are in feasible agreement with the values provided by Anderson et al. (2012), Thébault et al. (2018) and Wardinski et al. (2019). Furthermore, the related standard deviations for each coefficient are within the range of the variations resulting from the time varying magnetic field (Wardinski et al. 2019). To classify the coefficients reconstructed from the simulated magnetic field data (Table 2) in terms of that resulting from the MESSENGER data, the simulated data are furthermore evaluated along the MESSENGER trajectories. The reconstructed Gauss coefficients are listed in Table 6 of Appendix C. These coefficients are in agreement with the results presented in Table 4. Since the reconstructed coefficients resulting from the data evaluated along the MPO trajectories are in better agreement with the implemented coefficients than the coefficients resulting from the data evaluated along the MESSENGER trajectories, the restriction of the data points to the northern hemisphere and the related degradation of the model resolution influences the results significantly (Heyner et al. 2020, cf). Furthermore, it should be noted that the analysis of the length scales of the current densities (cf. "Validity of the thin shell approximation" section) cannot be performed for the MESSENGER orbits, because the current densities in the vicinity of the orbits are unknown.
<table><caption>**Table 4** Gauss coefficients for the internal dipole, quadrupole and octupole field and for the external dipole field reconstructed from MESSENGER data. In the last column the ranges of Gauss coefficients reconstructed from MESSENGER data by Anderson et al. (2012), Thébault et al. (2018) and Wardinski et al. (2019) are shown</caption><thead><tr><th>Gauss coefficient</th><th>Output capon in nT</th><th>Ranges of former works in nT</th></tr></thead><tbody><tr><td>g<sup>0</sup><sub>1</sub></td><td>-175.9</td><td>-215.8 to -190.0</td></tr><tr><td>g<sup>1</sup><sub>1</sub></td><td>-1.7</td><td>-2.9 to -0.9</td></tr><tr><td>h<sup>1</sup><sub>1</sub></td><td>4.5</td><td>0.8 to -2.7</td></tr><tr><td>g<sup>0</sup><sub>2</sub></td><td>-82.0</td><td>-83.2 to -57.7</td></tr><tr><td>g<sup>1</sup><sub>2</sub></td><td>0.8</td><td>-1.5 to -3.4</td></tr><tr><td>h<sup>1</sup><sub>2</sub></td><td>1.0</td><td>-1.4 to -0.0</td></tr><tr><td>g<sup>2</sup><sub>2</sub></td><td>2.9</td><td>-7.0 to -1.4</td></tr><tr><td>h<sup>2</sup><sub>2</sub></td><td>1.7</td><td>-3.3 to -0.4</td></tr><tr><td>g<sup>0</sup><sub>3</sub></td><td>-19.1</td><td>-36.7 to -15.7</td></tr><tr><td>g<sup>1</sup><sub>3</sub></td><td>1.1</td><td>1.8 to -4.1</td></tr><tr><td>h<sup>1</sup><sub>3</sub></td><td>0.4</td><td>0.3 to -0.8</td></tr><tr><td>g<sup>2</sup><sub>3</sub></td><td>-0.5</td><td>-1.5 to -9.2</td></tr><tr><td>h<sup>2</sup><sub>3</sub></td><td>-5.9</td><td>0.9 to -2.6</td></tr><tr><td>g<sup>3</sup><sub>3</sub></td><td>-5.8</td><td>-2.5 to -1.4</td></tr><tr><td>h<sup>3</sup><sub>3</sub></td><td>-3.3</td><td>0.1 to -0.3</td></tr><tr><td>q<sup>0</sup><sub>1</sub></td><td>-36.6</td><td>-39.7 to -23.2</td></tr><tr><td>q<sup>1</sup><sub>1</sub></td><td>-7.4</td><td>-0.2 to -0.7</td></tr><tr><td>s<sup>1</sup><sub>1</sub></td><td>0.8</td><td>-0.1 to -1.5</td></tr></tbody></table>

## Summary and outlook

In the vicinity of Mercury no completely current-free region is present. Therefore, the Gauss representation does not yield a proper parametrization of Mercury's magnetospheric field. Extension of the Gauss representation to the Gauss-Mie representation allows a more complete characterization of Mercury's internal and magnetospheric field.

For the parameterization of the magnetic field the orbit where the magnetic field is measured, is conceptually covered by a spherical shell. Due to the underlying geometry the total measured magnetic field is a superposition of internal and external poloidal fields generated by toroidal currents flowing outside the spherical shell as well as toroidal and poloidal fields generated by currents flowing within the shell. Thereby, each component of the field is either described by a scalar potential ($\Phi^i, \Phi^e$) or a scalar function ($\Psi_p^{sh}, \Psi_T^{sh}$). These potentials and functions can be expanded into spherical harmonics. When the thickness of the spherical shell is smaller than the length scale on which the toroidal current density changes in radial direction the shell is called a thin shell. Then the poloidal field generated by currents flowing inside the shell is negligible compared to the poloidal field generated by currents flowing beyond the shell,

whereas the toroidal field remains finite. In the case of the planned MPO orbits, the thin shell approximation is a reasonable choice.

For the application of the thin shell approximation the internal Gauss coefficients for the dipole, quadrupole and octupole field are implemented in the simulation code AIKEF and the magnetic field data resulting from the plasma interaction of Mercury with the solar wind are simulated in the vicinity of Mercury. The data are evaluated along the planned MPO orbits and the 46 expansion coefficients, describing the internal, external and the toroidal field are reconstructed with Capon's method. Since the reconstructed internal Gauss coefficients are in good agreement with that implemented in the simulation code, the parameterization of the magnetic field using the Gauss-Mie representation is a suitable alternative to the application of global magnetospheric models. Even the implemented Gauss coefficient for the octupole field of Mercury can be reconstructed accurately and therefore, it is expectable that Mercury's internal octupole field will be reconstructed with high precision from the magnetometer data on board MPO. Thus, it is worthwhile to investigate the analysis of higher multipoles, such as hexadecapole, from the data along the MPO trajectories.

Furthermore, the thin shell approximation is applied to reconstruct Mercury's internal magnetic field from the data of the MESSENGER mission. The results are in reasonable agreement with former works. Since only the data points in the northern hemisphere are vital for the analysis of the MESSENGER data, the symmetrically distributed MPO orbits will deliver a better model resolution than the MESSENGER orbits.

Concerning the BepiColombo mission the combination of the Gauss representation with the Mie representation is a useful model for the analysis of Mercury's internal magnetic field. As the BepiColombo mission consists of two elements, the planetary orbiter and the magnetospheric orbiter, measurements of any gradients of Mercury's magnetic field are possible which may lead to further improvements of the methods presented here. Besides the analysis of the internal magnetic field, the reconstructed coefficients for the toroidal magnetic field can be used for calculating poloidal current systems, e.g. field aligned currents, in the vicinity of the orbit where the data are evaluated.
## Appendix A: The thin shell approximation for the toroidal magnetic field

As derived in "The thin shell approximation" section, the scalar function of the poloidal magnetic field can be estimated upwards as

$$|\Psi_P^{sh}(r)| \le 2bL \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{h}{L}\right)^{n+1} \max_{\theta', \lambda'} \{|\Gamma_T^{sh}(b)|\}. \quad (79)$$

When $h \ll L$, the function $\Psi_P^{sh}$ vanishes. Now the influence of the thin shell approximation for the toroidal magnetic field $\mathbf{B}_T^{sh}$ is discussed.

The scalar function $\Psi_T^{sh}$ of the toroidal magnetic field $\mathbf{B}_T^{sh}$ can be expanded analogously into a Taylor series in the vicinity of the mean radius $b$

$$\Psi_T^{sh}(\mathbf{r}) = \Gamma_P^{sh}(\mathbf{r}) = \sum_{n=0}^{\infty} \frac{1}{n!} \partial_r^n \Gamma_P^{sh}(\mathbf{r})|_{\mathbf{r}=b}(r-b)^n, \quad (80)$$

using that the scalar functions $\Psi_T^{sh}$ and $\Gamma_P^{sh}$ are the same (cf. Eq. 54).

If we furthermore assume that also the radial derivatives of the poloidal current density are bounded so that

$$|\partial_r^n \Gamma_P^{sh}| \le \frac{|\Gamma_P^{sh}|}{L^n}, \quad (81)$$

the function $\Psi_T^{sh}(r)$ can be estimated upwards as

$$\begin{aligned} |\Psi_T^{sh}(\mathbf{r})| &\le \sum_{n=0}^{\infty} \frac{1}{n!} |\partial_r^n \Gamma_P^{sh}(b)| \cdot |r-b|^n \\ &\le \sum_{n=0}^{\infty} \frac{1}{n!} \frac{|\Gamma_P^{sh}(b)|}{L^n} |r-b|^n. \end{aligned} \quad (82)$$

Using that $\Psi_T^{sh}(\mathbf{r})$ is solely evaluated within the spherical shell, where $r \in [b-h, b+h]$ and therefore $|r-b| \le h$ results in

$$|\Psi_T^{sh}(\mathbf{r})| \le \sum_{n=0}^{\infty} \frac{1}{n!} \left(\frac{h}{L}\right)^n |\Gamma_P^{sh}(b)|. \quad (83)$$

When $h \ll L$, the function $\Psi_T^{sh}$ remains finite.

## Appendix B: Model resolution

For the analysis of Mercury's internal magnetic field from the MESSENGER data only data points in the northern hemisphere can be considered. This limitation impairs the condition number $\kappa$ of the shape matrix $\underline{H}$ from $\kappa \approx 226$ for the data points along the MPO orbits to $\kappa \approx 10^8$ for the data points along a single MESSENGER orbit, where the shape matrix describes the spatial distribution of the data

points. The condition number is defined as the ratio of the largest and the smallest singular value of the shape matrix. Therefore, some of the 46 singular values, corresponding to the 46 expansion coefficients, of the shape matrix have to be dropped within the numerical calculation to improve the condition number. Within the low rank approximation (Eckart and Young 1936) only $k$ singular values, where $k \le 46$, can be considered and therefore, the shape matrix $\underline{H}$ is approximated by a shape matrix $\underline{H}_k$ which has a lower rank and a lower condition number.

Capon's filter matrix

$$\underline{w}^{\dagger} = [\underline{H}^{\dagger} \underline{M}^{-1} \underline{H}]^{-1} \underline{H}^{\dagger} \underline{M}^{-1} \quad (84)$$

which is the key parameter for calculating Capon's estimator

$$\underline{g}_C = \underline{w}^{\dagger} \underline{B} \quad (85)$$

fulfills the distortionless constraint

$$\underline{w}^{\dagger} \underline{H} = \underline{I}, \quad (86)$$

where $\underline{I}$ is the identity matrix and $\underline{M} = (\underline{B} \circ \underline{B})$ describes the data covariance matrix of the magnetic field measurements $\underline{B}$ (Toepfer et al. 2020a, b). As a consequence of the low rank approximation Capon's filter matrix is modified to

$$\underline{w}^{\dagger}_{k} = [\underline{H}^{\dagger}_{k} \underline{M}^{-1} \underline{H}_{k}]^{-1} \underline{H}^{\dagger}_{k} \underline{M}^{-1} \quad (87)$$

so that

$$\underline{R} = \underline{w}^{\dagger}_{k} \underline{H} \neq \underline{I}. \quad (88)$$

The matrix $\underline{R} = \underline{w}^{\dagger}_{k} \underline{H}$ is called the model resolution matrix. Because of

$$\underline{g}_C^k = \underline{R} \underline{g}, \quad (89)$$

where $g_C^k$ denotes the estimator resulting from the consideration of $k$ singular values and $g$ is the ideal coefficient vector implemented in the simulation code, the diagonal elements of $\underline{R}$ identify the resolution of each coefficient (Connerney 1981). When $\underline{R} = \underline{I}$ each coefficient is resolved for 100%. If the resolution is smaller than 100%, there exist model parameter covariances.

The more singular values are considered within the estimation, the better the model resolution becomes, whereas the condition number of the shape matrix increases and thus, the influence of measurement errors increases.

To achieve a compromise between the resolution and the condition number, the coefficients are estimated
for different numbers of singular values and the change of the coefficients resulting from *k* singular values to the coefficients resulting from *k* − 1 singular values is regarded. For the final computation the maximum number of singular values is chosen from which the reconstructed coefficients are almost constant. In Fig. 4

<figure><img src="image_8.png" /></figure>

the procedure is exemplarily illustrated for the coefficient $g_1^0$, which corresponds to the diagonal element $R_{11}$ by analyzing the data of one pair of MESSENGER orbits.

Changing the number of considered singular values incrementally from 46 to 37 the resulting coefficient changes significantly. For $k \le 37$ the values are almost constant and therefore, 37 singular values, corresponding to a model resolution of 88% for the coefficient $g_1^0$ are considered within the estimation.

**Appendix C: Tables: Reconstructed Gauss coefficients from MESSENGER orbits**
<table><caption>**Table 5** Gauss coefficients for the dipole, quadrupole and octupole field reconstructed from selected pairs of MESSENGER orbits</caption><thead><tr><th>Orbits</th><th>g<sup>0</sup><sub>1</sub></th><th>g<sup>1</sup><sub>1</sub></th><th>h<sup>1</sup><sub>1</sub></th><th>g<sup>0</sup><sub>2</sub></th><th>g<sup>1</sup><sub>2</sub></th><th>h<sup>1</sup><sub>2</sub></th><th>g<sup>2</sup><sub>2</sub></th><th>h<sup>2</sup><sub>2</sub></th><th>g<sup>0</sup><sub>3</sub></th><th>g<sup>1</sup><sub>3</sub></th><th>h<sup>1</sup><sub>3</sub></th><th>g<sup>2</sup><sub>3</sub></th><th>h<sup>2</sup><sub>3</sub></th><th>g<sup>3</sup><sub>3</sub></th><th>h<sup>3</sup><sub>3</sub></th><th>q<sup>0</sup><sub>1</sub></th><th>q<sup>1</sup><sub>1</sub></th><th>s<sup>1</sup><sub>1</sub></th></tr></thead><tbody><tr><td>2013-T245, 2014-T21</td><td>-197.4</td><td>2.0</td><td>-2.0</td><td>-71.6</td><td>-10.1</td><td>21.1</td><td>-40.1</td><td>-8.5</td><td>-6.7</td><td>9.0</td><td>1.7</td><td>9.5</td><td>-7.3</td><td>23.5</td><td>80.2</td><td>-65.2</td><td>-47.8</td><td>-13.7</td></tr><tr><td>2014-T143, 2014-T195</td><td>-190.0</td><td>18.0</td><td>15.3</td><td>-66.7</td><td>16.9</td><td>-7.6</td><td>5.2</td><td>2.7</td><td>-11.5</td><td>-10.9</td><td>-24.5</td><td>4.0</td><td>-11.4</td><td>-43.4</td><td>-21.9</td><td>-24.2</td><td>-12.1</td><td>-21.9</td></tr><tr><td>2013-T330, 2014-T107</td><td>-158.5</td><td>-15.1</td><td>-1.6</td><td>-96.6</td><td>17.0</td><td>-2.2</td><td>-6.4</td><td>2.1</td><td>-10.0</td><td>-1.8</td><td>12.5</td><td>-6.4</td><td>-8.2</td><td>3.8</td><td>7.1</td><td>5.6</td><td>2.7</td><td>3.3</td></tr><tr><td>2013-T330, 2014-T195</td><td>-164.1</td><td>-6.5</td><td>6.6</td><td>-81.2</td><td>12.7</td><td>6.3</td><td>12.0</td><td>-52.2</td><td>-13.6</td><td>-5.9</td><td>-14.1</td><td>-6.3</td><td>3.6</td><td>-32.5</td><td>21.5</td><td>-5.8</td><td>-16.2</td><td>14.3</td></tr><tr><td>2014-T107, 2014-T143</td><td>-155.6</td><td>-2.6</td><td>5.0</td><td>-81.7</td><td>22.8</td><td>-8.0</td><td>13.5</td><td>30.7</td><td>-20.8</td><td>-8.9</td><td>-3.6</td><td>-20.3</td><td>-18.8</td><td>-32.1</td><td>-19.9</td><td>-24.5</td><td>-0.8</td><td>16.8</td></tr><tr><td>2012-T223, 2012-T311</td><td>-190.4</td><td>0.2</td><td>0.5</td><td>-76.0</td><td>73.0</td><td>11.9</td><td>-30.0</td><td>-82.1</td><td>-17.4</td><td>-65.6</td><td>-68.5</td><td>-41.4</td><td>88.9</td><td>-178.1</td><td>53.2</td><td>9.6</td><td>41.1</td><td>-41.5</td></tr><tr><td>2013-T330, 2014-T21</td><td>-163.3</td><td>-9.0</td><td>7.9</td><td>-80.6</td><td>-8.5</td><td>5.6</td><td>-9.8</td><td>-27.2</td><td>-1.8</td><td>3.5</td><td>-11.5</td><td>14.4</td><td>-39.1</td><td>-22.6</td><td>35.5</td><td>11.4</td><td>-27.4</td><td>-5.5</td></tr><tr><td>2014-T107, 2013-T245</td><td>-162.3</td><td>-4.7</td><td>-27.2</td><td>-96.6</td><td>-9.5</td><td>16.0</td><td>-29.6</td><td>32.0</td><td>-27.5</td><td>10.7</td><td>31.2</td><td>-21.7</td><td>-13.2</td><td>25.1</td><td>35.5</td><td>-61.5</td><td>-1.7</td><td>-38.4</td></tr><tr><td>2014-T195</td><td>-163.7</td><td>8.7</td><td>-17.3</td><td>-87.0</td><td>-12.8</td><td>-21.2</td><td>-22.1</td><td>10.5</td><td>-32.9</td><td>12.8</td><td>17.3</td><td>-18.4</td><td>13.7</td><td>53.2</td><td>-41.4</td><td>-38.2</td><td>-16.2</td><td>12.3</td></tr><tr><td>Mean value</td><td>-175.9</td><td>-1.7</td><td>4.5</td><td>-82.0</td><td>0.8</td><td>1.0</td><td>2.9</td><td>1.7</td><td>-19.1</td><td>1.1</td><td>0.4</td><td>-0.5</td><td>-5.9</td><td>-5.8</td><td>-3.3</td><td>-36.6</td><td>-7.4</td><td>0.8</td></tr><tr><td>Standard deviation</td><td>15.8</td><td>5.9</td><td>6.2</td><td>10.2</td><td>13.9</td><td>8.2</td><td>10.6</td><td>7.8</td><td>8.5</td><td>9.2</td><td>12.7</td><td>12.1</td><td>11.0</td><td>26.8</td><td>21.2</td><td>23.2</td><td>8.4</td><td>15.0</td></tr></tbody></table>
<table><caption>**Table 6** Gauss coefficients for the dipole, quadrupole and octupole field reconstructed from the simulated data evaluated at selected pairs of MESSENGER orbits</caption><thead><tr><th>Orbits</th><th>g<sup>0</sup><sub>1</sub></th><th>g<sup>1</sup><sub>1</sub></th><th>h<sup>1</sup><sub>1</sub></th><th>g<sup>0</sup><sub>2</sub></th><th>g<sup>1</sup><sub>2</sub></th><th>h<sup>1</sup><sub>2</sub></th><th>g<sup>2</sup><sub>2</sub></th><th>h<sup>2</sup><sub>2</sub></th><th>g<sup>0</sup><sub>3</sub></th><th>g<sup>1</sup><sub>3</sub></th><th>h<sup>1</sup><sub>3</sub></th><th>g<sup>2</sup><sub>3</sub></th><th>h<sup>2</sup><sub>3</sub></th><th>g<sup>3</sup><sub>3</sub></th><th>h<sup>3</sup><sub>3</sub></th><th>g<sup>0</sup><sub>1</sub></th><th>q<sup>1</sup><sub>1</sub></th><th>s<sup>1</sup><sub>1</sub></th></tr></thead><tbody><tr><td>2013-T245, 2014-T21</td><td>-136.4</td><td>9.9</td><td>3.7</td><td>-103.8</td><td>-16.0</td><td>5.2</td><td>-16.3</td><td>-4.0</td><td>-24.4</td><td>7.6</td><td>-6.8</td><td>-5.7</td><td>-5.1</td><td>10.6</td><td>10.6</td><td>-39.2</td><td>-13.0</td><td>-11.1</td></tr><tr><td>2014-T143, 2014-T195</td><td>-179.6</td><td>-10.5</td><td>-4.4</td><td>-93.7</td><td>21.6</td><td>7.4</td><td>-3.6</td><td>-15.0</td><td>-12.2</td><td>-10.7</td><td>-10.9</td><td>3.7</td><td>5.7</td><td>-13.9</td><td>2.1</td><td>-31.9</td><td>-27.5</td><td>4.3</td></tr><tr><td>2013-T330, 2014-T107</td><td>-148.9</td><td>-15.1</td><td>5.6</td><td>-105.7</td><td>15.4</td><td>7.6</td><td>3.7</td><td>-34.3</td><td>-15.2</td><td>-5.8</td><td>-8.4</td><td>-6.1</td><td>6.9</td><td>-17.1</td><td>11.7</td><td>-20.2</td><td>-22.8</td><td>-5.9</td></tr><tr><td>2013-T330, 2014-T195</td><td>-149.3</td><td>-12.6</td><td>3.7</td><td>-99.4</td><td>10.4</td><td>4.1</td><td>7.0</td><td>-31.0</td><td>-18.2</td><td>-3.5</td><td>-4.1</td><td>-6.8</td><td>2.3</td><td>-16.9</td><td>17.0</td><td>-12.8</td><td>-18.7</td><td>-6.8</td></tr><tr><td>2014-T107, 2014-T143</td><td>-190.4</td><td>-21.2</td><td>-8.6</td><td>-97.0</td><td>25.1</td><td>6.4</td><td>-15.8</td><td>-17.3</td><td>-8.1</td><td>-11.3</td><td>-10.0</td><td>12.3</td><td>7.4</td><td>-10.0</td><td>13.1</td><td>-43.2</td><td>-36.7</td><td>-2.5</td></tr><tr><td>2012-T223, 2012-T311</td><td>-150.1</td><td>16.3</td><td>-1.8</td><td>-113.9</td><td>27.9</td><td>2.7</td><td>24.8</td><td>-39.4</td><td>-28.8</td><td>-44.8</td><td>-4.4</td><td>-57.4</td><td>77.5</td><td>-79.4</td><td>-16.4</td><td>-41.5</td><td>-9.8</td><td>-35.4</td></tr><tr><td>2013-T330, 2014-T21</td><td>-105.0</td><td>-14.8</td><td>7.5</td><td>-131.5</td><td>6.4</td><td>2.1</td><td>3.1</td><td>-17.1</td><td>-13.9</td><td>-1.6</td><td>0.0</td><td>-17.1</td><td>-3.2</td><td>-8.5</td><td>9.0</td><td>-6.6</td><td>-7.8</td><td>-22.6</td></tr><tr><td>2014-T107, 2013-T245</td><td>-144.5</td><td>-28.4</td><td>1.6</td><td>-117.6</td><td>36.5</td><td>21.9</td><td>1.7</td><td>-10.5</td><td>-8.7</td><td>-13.7</td><td>-18.8</td><td>-5.7</td><td>2.3</td><td>-25.1</td><td>6.6</td><td>-4.6</td><td>-37.9</td><td>6.8</td></tr><tr><td>2014-T195</td><td>-131.5</td><td>6.8</td><td>-6.3</td><td>-114.8</td><td>-51.6</td><td>-11.7</td><td>-1.9</td><td>0.4</td><td>-19.0</td><td>46.0</td><td>-9.7</td><td>-0.1</td><td>-41.7</td><td>82.1</td><td>-7.4</td><td>-53.3</td><td>-3.4</td><td>51.7</td></tr><tr><td>Mean value</td><td>-163.7</td><td>1.1</td><td>0.1</td><td>-104.0</td><td>7.6</td><td>5.1</td><td>1.7</td><td>-10.6</td><td>-18.8</td><td>-5.6</td><td>-6.8</td><td>-3.5</td><td>2.3</td><td>-9.3</td><td>6.5</td><td>-31.5</td><td>-14.7</td><td>-2.5</td></tr><tr><td>Standard deviation</td><td>19.9</td><td>14.0</td><td>5.6</td><td>8.1</td><td>14.3</td><td>2.2</td><td>3.9</td><td>7.4</td><td>5.9</td><td>7.3</td><td>3.7</td><td>4.3</td><td>4.9</td><td>10.4</td><td>7.1</td><td>12.4</td><td>8.6</td><td>6.9</td></tr></tbody></table>
## Acknowledgements

The authors are grateful for stimulating discussions and helpful suggestions by Alexander Schwenke and the technical support by Willi Exner.

## Authors' contributions

All authors contributed conception and design of the study; ST and UM wrote the first draft of the manuscript; all authors contributed to manuscript revision. All authors read and approved the final manuscript.

## Funding

Open Access funding enabled and organized by Projekt DEAL. We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the Technische Universität Braunschweig. The work by Y. Narita is supported by the Austrian Space Applications Programme at the Austrian Research Promotion Agency under contract 865967. D. Heyner and K.-H. Glassmeier were supported by the German Ministerium für Wirtschaft und Energie and the German Zentrum für Luft- und Raumfahrt under contract 50 QW1501.

## Availability of data and materials

Simulation data can be provided upon request.

## Declarations

### Competing interests

The authors declare that they have no conflict of interest.

## Author details

¹ Institut für Theoretische Physik, Technische Universität Braunschweig, Mен delssohnstraße 3, 38106 Braunschweig, Germany. ² Space Research Institute, Austrian Academy of Sciences, Schmiedlstraße 6, 8042 Graz, Austria. ³ Institut für Geophysik und extraterrestrische Physik, Technische Universität Braunschweig, Mен delssohnstraße 3, 38106 Braunschweig, Germany. ⁴ Deutsches Zentrum für Luft- und Raumfahrt, Institut für Planetenforschung, Rutherfordstraße 2, 12489 Berlin, Germany. ⁵ Laboratoire de Planétologie et Géodynamique, UMR 6112, CNRS, Université de Nantes, Université d'Angers, 2 Rue de la Houssinière, 44000 Nantes, France.

Received: 18 December 2020 Accepted: 20 February 2021

Published online: 08 March 2021

## References

Alexeev II, Belenkaya ES, Bobrovnikov SYu, Slavin JA, Sarantos M (2008) Parabola model of Mercury's magnetosphere. J Geophys Res. 113:A12210. https://doi.org/10.1029/2008JA013368

Alexeev II, Belenkaya ES, Slavin JA, Korth H, Anderson BJ, Baker DN, Boardsen SA, Johnson CL, Purucker ME, Sarantos M, Solomon SC (2010) Mercury's magnetospheric magnetic field after the first two MESSENGER flybys. Icarus 209:23-39. https://doi.org/10.1016/j.icarus.2010.01.024

Anderson BJ, Johnson CL, Korth H, Purucker ME, Winslow RM, Slavin JA, Solomon SC, McNutt RL Jr, Raines JM, Zurbuchen TH (2011) The global magnetic field of Mercury from MESSENGER orbital observations. Science 333:1859-1862. https://doi.org/10.1126/science.1211001

Anderson BJ, Johnson CL, Korth H, Winslow RM, Borovsky JE, Purucker ME et al (2012) Low-degree structure in Mercury's planetary magnetic field. J Geophys Res. 117:E00L12. https://doi.org/10.1029/2012JE004159

Backus G (1986) Poloidal and toroidal fields in geomagnetic field modeling. Rev Geophys. 24:75-109. https://doi.org/10.1029/RG024i001p00075

Backus G, Parker R, Constable C (1996) Foundations of Geomagnetism. Cambridge University Press, Cambridge

Bayer M, Freeden W, Maier T (2001) A vector wavelet approach to iono- and magnetospheric geomagnetic satellite data. J Atmos Solar-Terrestr Phys. 63:581-597. https://doi.org/10.1016/S1364-6826(00)00234-0

Benkhoff J, van Casteren J, Hayakawa H, Fujimoto M, Laakso H, Novara M, Ferri P, Middleton HR, Ziethe R (2010) BepiColombo-comprehensive exploration of mercury: mission overview and science goals. Planet Space Sci. 85(1-2):2-20. https://doi.org/10.1016/j.pss.2009.09.020

Capon J (1969) High resolution frequency-wavenumber spectrum analysis. Proc IEEE 57:1408-1418. https://doi.org/10.1109/PROC.1969.7278

Connerney JEP (1981) The magnetic field of Jupiter: a generalized inverse approach. J Geophys Res. 86:7679-7693. https://doi.org/10.1029/JA086iA09p07679

Eckart C, Young G (1936) The approximation of one matrix by another of lower rank. Psychometrika 1:211-218. https://doi.org/10.1007/BF02288367

Engels U, Olsen N (1999) Computation of magnetic fields within source regions of ionospheric and magnetospheric currents. J Atmos Solar-Terrestr Phys. 60:1585-1592. https://doi.org/10.1016/S1364-6826(98)00094-7

Exner W, Heyner D, Liuzzo L, Motschmann U, Shiota D, Kusano K, Shibayama T (2018) Coronal mass ejection hits mercury: A.I.K.E.F. hybrid-code results compared to MESSENGER data. Planet Space Sci 153:89-99. https://doi.org/10.1016/j.pss.2017.12.016

Exner W, Simon S, Heyner D, Motschmann U (2020) Influence of Mercury's exosphere on the structure of the magnetosphere. J Geophys Res. 125:e27691. https://doi.org/10.1029/2019JA027691

Gauß CF (1839) Allgemeine Theorie des Erdmagnetismus: Resultate aus den Beobachtungen des magnetischen Vereins im Jahre 1838, edited by: Gauss, C. F. and Weber, W., 1-57, Weidmannsche Buchhandlung, Leipzig, 1839

Glassmeier K-H, Motschmann U, Dunlop M, Balogh A, Acuña MH, Carr C, Musmann G, Fornacçon K-H, Schweda K, Vogt J, Georgescu E, Buchert S (2001) Cluster as a wave telescope-first results from the fluxgate magnetometer. Ann Geophys. 19:1439-1447. https://doi.org/10.5194/angeo-19-1439-2001

Glassmeier K-H, Auster H-U, Heyner D, Okraska F, Carr C, Berghofer G, Anderson BJ et al (2010) The fluxgate magnetometer of the BepiColombo Mercury Planetary Orbiter. Planet Space Sci. 58:287-299. https://doi.org/10.1016/j.pss.2008.06.018

Glassmeier K-H, Tsurutani BT (2014) Carlfried gauss-general theory of terrestrial magnetism-a revised translation of the German text. Hist Geo Space Sci. 5:11-62. https://doi.org/10.5194/hgss-5-11-2014

Heyner D, Auster H-U, Fornacon K-H, Carr C, Richter I, Mith JZD, Kolhey P et al (2020) The BepiColombo planetary magnetometer MPO-MAG: What can we learn from the Hermean magnetic field? Space Sci Rev.

Jacobs JA (1987) Geomagnetism 2. Press Acad, London

Johnson CL, Purucker ME, Korth H, Anderson BJ, Winslow RM, Al Asad MMH, Slavin JA, Alexeev II, Phillips RJ, Zuber MT, Solomon SC (2012) Messenger observations of Mercury's magnetic field structure. J Geophys Res. 117:0014. https://doi.org/10.1029/2012JE004217

Kazantsev SG, Kardakov VB (2019) Poloidal-toroidal decomposition of solenoidal vector fields in the ball. J Applied Industr Math. 13:480-499. https://doi.org/10.1134/S1990478919030098

Korth H, Anderson BJ, Acuña Mario H, Slavin JA, Tsyganenko NA, Solomon SC, McNutt RL (2004) Determination of the properties of Mercury's magnetic field by the messenger mission. Planet Space Sci. 52:733-746. https://doi.org/10.1016/j.pss.2003.12.008

Kosic JK (1984) Quantitative magnetospheric magnetic field modelling with toroidal and poloidal vector fields. Planet Space Sci. 32:965-974. https://doi.org/10.1016/0032-0633(84)90053-9

Krause F, Rädler K-H (1980) Mean-field magnetohydrodynamics and dynamo theory. Pergamon Press, Oxford

Mayer C, Maier T (2006) Separating inner and outer Earth's magnetic field from CHAMP satellite measurements by means of vector scaling functions and wavelets. Geophys J Int 167(3):1188-1203. https://doi.org/10.1111/j.1365-246X.2006.03199.x

Meir R (1994) Bias, variance and the combination of least squares estimators, MIT Press, Cambridge, MA, USA, Proceedings of the 7th International Conference on Neural Information Processing Systems, 295-302

Motschmann U, Woodward TI, Glassmeier K-H, Southwood DJ, Pinçon J-L (1996) Wavelength and direction filtering by magnetic measurements at satellite arrays: generalized minimum variance analysis. J Geophys Res. 101:4961-4966. https://doi.org/10.1029/95JA03471

Müller J, Simon S, Motschmann U, Schüle J, Glassmeier K-H, Pringle GJ (2011) A.I.K.E.F.: Adaptive hybrid model for space plasma simulations. Comp Phys Comm 182:946-966. https://doi.org/10.1016/j.cpc.2010.12.033

Narita Y, Glassmeier K-H, Treumann RA (2006) Wave-number spectra and intermittency in the terrestrial Foreshock region. Phys Rev Lett. 97:191101. https://doi.org/10.1103/PhysRevLett.97.191101

Narita Y (2012) Plasma Turbulence in the Solar System. Springer-Verlag, Berlin Heidelberg. https://doi.org/10.1007/978-3-642-25667-7
Olsen N (1997) Ionospheric F currents at middle and low latitudes estimated from Magsat data. J Geophys Res. 102:4569–4576. https://doi.org/10.1029/96JA02949

Olsen N, Glassmeier K-H, Jia X (2010) Separation of the magnetic field into external and internal parts. Space Sci Rev 152:135–157. https://doi.org/10.1007/s11214-009-9563-0

Thébault E, Langlais B, Oliveira JS, Amit H, Leclercq L (2018) A time-averaged regional model of the Hermean magnetic field. Phys Earth Planet Interior 276:93–105. https://doi.org/10.1016/j.pepi.2017.07.001

Toepfer S, Narita Y, Heyner D, Motschmann U (2020a) The Capon method for Mercury's magnetic field analysis. Front Phys. 8:249. https://doi.org/10.3389/fphy.2020.00249

Toepfer S, Narita Y, Heyner D, Kolhey P, Motschmann U (2020b) Mathematical foundation of Capon's method for planetary magnetic field analysis.

Geosci Instrum Method Data Syst. 9:471–481. https://doi.org/10.5194/gi-9-471-2020

Wardinski I, Langlais B, Thébault E (2019) Correlated time-varying magnetic field and the core size of Mercury. J Geophys Res. 124:2178–2197. https://doi.org/10.1029/2018JE005835

Zharkova AY, Kreslavsky MA, Head JW, Kokhanov AA (2020) Regolith textures on Mercury: comparison with the Moon. Icarus 351:113945. https://doi.org/10.1016/j.icarus.2020.113945

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com