# MoMo: a new empirical model of the Mars ionospheric total electron content based on Mars Express MARSIS data

Nicolas Bergeot, Olivier Witasse, Sébastien Le Maistre, Pierre-Louis Blelly, Wlodek Kofman, Kerstin Peter, Véronique Dehant, Jean-Marie Chevalier

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

## ► To cite this version:

Nicolas Bergeot, Olivier Witasse, Sébastien Le Maistre, Pierre-Louis Blelly, Wlodek Kofman, et al. MoMo: a new empirical model of the Mars ionospheric total electron content based on Mars Express MARSIS data. Journal of Space Weather and Space Climate, 2019, 9, pp.A36. 10.1051/swsc/2019035. insu-02938684

HAL Id: insu-02938684

https://insu.hal.science/insu-02938684

Submitted on 15 Sep 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# MoMo: a new empirical model of the Mars ionospheric total electron content based on Mars Express MARSIS data

Nicolas Bergeot<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*, Olivier Witasse<fnref n="4" />, Sébastien Le Maistre<fnref n="2" />, Pierre-Louis Blelly<fnref n="5" />, Wlodek Kofman<fnref n="6" />,<fnref n="7" />, Kerstin Peter<fnref n="8" />, Véronique Dehant<fnref n="2" />,<fnref n="3" />, and Jean-Marie Chevalier<fnref n="1" />,<fnref n="2" />

Received 27 February 2019 / Accepted 9 September 2019

**Abstract**—Aims: Several scientific landers and rovers have reached the Martian surface since the 1970s. Communication between the asset (i.e., lander or rover) and Mars orbiters or Earth antennas uses radio signals in UHF to X-band frequencies passing through the Mars' ionosphere. It is consequently necessary to take into account electron density variation in the Mars' ionosphere to correct the refraction of the signal transmitted.

**Methods:** We developed a new empirical model of the Mars' ionosphere called MoMo. It is based on the large database of Total Electron Content (TEC) derived from the subsurface mode of the Mars Express MARSIS radar. The model provides vertical TEC as a function of solar zenith angle, solar activity, solar longitude and location. For validation, the model is compared with Mars Express radio occultation data as well as with the numerical model IPIM (IRAP Plasmasphere-Ionosphere Model).

**Results:** We discussed the output of the model in terms of climatology behaviour of the Mars' ionosphere. The output of MoMo is then used to quantify the impact of the Martian ionosphere for radio-science experiments. From our results, the effect is of the order of $10^{-3}$ mm s$^{-1}$ in Doppler observables especially around sunrise and sunset. Consequently, this new model could be used to support the data analysis of any radio-science experiment and especially for present InSight RISE and futur ExoMars LARA instruments aiming at better understand the deep-interior of Mars.

## 1 Introduction

Radio-science experiments are used to obtain motions of spacecraft or asset (e.g., lander and rover) around or on planets or moons in space. One usually performs Doppler measurements alone on the radio signal between the Earth and the asset in order to reconstruct its position and velocity variations. Ranging is also used for obtaining ephemerides of planets or moons. For a spacecraft or asset on its surface at Mars, dedicated radio frequency bands for space communications are determined by the International Telecommunication Union (ITU). The frequencies used are in UHF (0.3–3.0 GHz), S-band (2.0–4.0 GHz), X-band (8.0–12.0 GHz), and in the future, Ka-band (26.5–40 GHz). At present most of the spacecraft

use UHF for inter-satellite links and X-band for Earth-Mars links. Previously, landers such as Viking and or Pathfinder used S-band.

Concerning recent and future missions, the NASA lander InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport), which landed successfully in November 2018, conducts a radio-science experiment called RISE (Rotation and Interior Structure Experiment, Folkner et al., 2018). The lander signal in X-band is routinely tracked from Earth, with the goal of measuring the rotational movements of the red planet, in order to study its deep interior. A very similar and complementary radio-science experiment called LaRa (Lander Radio science, Dehant et al., 2009, 2011, 2019), was also designed to refine the rotation model of Mars and to further constrain its interior. LaRa is part of the scientific payload of the Russian surface platform of the

<footnote n="1">¹ Solar-Terrestrial Centre of Excellence (STCE), Space Pole, 1180 Brussels, Belgium</footnote>

<footnote n="2">² Royal Observatory of Belgium Brussels, 1180 Uccle, Belgium</footnote>

<footnote n="3">³ Earth and Climate, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium</footnote>

<footnote n="4">⁴ European Space Agency, ESTEC-Scientific Support Office, ESTEC, 2201 AZ Noordwijk, The Netherlands</footnote>

<footnote n="5">⁵ Institut de Recherche en Astrophysique et Planétologie (IRAP), 31400 Toulouse, France</footnote>

<footnote n="6">⁶ Université Grenoble Alpes, CNRS, CNES, IPAG, 38400 Grenoble, France</footnote>

<footnote n="7">⁷ Space Research Centre, PAS, 00 716 Warsaw, Poland</footnote>

<footnote n="8">⁸ Rheinisches Institut für Umweltforschung an der Universität zu Köln, Abt. Planetenforschung, 50931 Köln, Germany</footnote>

<footnote>*Corresponding author: nicolas.bergeot@oma.be</footnote>
ESA-Roscosmos ExoMars 2020 mission. The X-band measure-
ments provided by these two missions have unprecedented
accuracy (and low power level in the case of LaRa). In partic-
ular, at least centimeter accuracy in the position determined
from the phase delay is required. Doppler measurements from
RISE and LaRa instruments are at a precision as high as
0.02 mm s⁻¹ at 60 s integration time corresponding to the
instrument precision requirement. However, the impact of a
liquid core on Martian Orientation Parameters is lower than
0.004 mm s⁻¹ in X-band (Yseboodt et al., 2017) and could
be detected from data time-series. To reach this accuracy, the
ionospheres of Mars and Earth as well as the interplanetary
plasma contributions need to be corrected for along the line
of sight. Because the ionosphere and plasma effects depend
on frequency, an adequate combination of radio signals on
two frequencies allows scientists to obtain corrected radio links
(see Mars Express radio science experiment in Pätzold et al.,
2004). However, landers have usually only one frequency
disabling the correction for these effects. This is also the case
for the radio-science experiments RISE and Lara which use only
one S-band frequency. It is thus important to take into account
the ionospheric effect on radio signals using external models to
correct for ionospheric and plasma effects.

Ionosphere effects for the Earth can be modelled from, e.g., GNSS observations (e.g., Klobuchar et al., 1996) and are nowadays well corrected for most radio-wave based applications. Solar plasma is another source of degradation of the radio-science data quality. Pretty low when Mars elongation is above ~15°, the solar plasma noise strongly increases for smaller angles. Data acquired close to solar conjunctions are usually excluded from radio-science analysis or they are properly weighted (e.g., Zuber et al., 2007) when combined with data acquired at larger angle. Concerning the Mars' ionosphere, several models exist. Let us cite the IPIM model (IRAP Plasmasphere-Ionosphere Model, Marchaudon & Bleilly, 2015), a state of the art numerical model, which is heavy to run for routine processes. Another existing model is the recent empirical model called NeMars (Sánchez-Cano et al., 2013) constructed from MARSIS data. This model is based on data from a special mode of the Mars Express radar (active ionospheric sounding), which does not probe the entire ionospheric layer (Sánchez-Cano et al., 2012; Morgan et al., 2013). Peter et al. (2014) published a 1-dimensional model of the ionosphere, used to interpret the Mars Express radio-occultation data. Recently, the Mars International Reference Ionosphere (MIRI-2018, Mendillo et al., 2018) provides a new module for Total Electron Content (TEC) using more than 120,000 values from the Mars Reconnaissance Orbiter/SHAllow RADar (MRO/SHARAD). The TEC output available on the web (http://sirius.bu.edu/miri/miri.php) is driven by Solar Zenith Angle (SZA) and solar index only. The output TEC is the integrated value for the electron density for altitudes comprise between 80 and 270 km while the ionopause (i.e., upper limit of the ionospheric plasma) is generally considered to be more than 400 km (e.g., Morel et al., 2004).

The aim of the present paper is to provide a new model of
the Mars' ionosphere with few coefficients, simple to use
but taking into account for areophysical considerations
(i.e., absorption of the solar EUV in a terrestrial atmosphere,
seasons, hemispheric dichotomies) and adapted to correct data
for any experiment using radio signals. In order to build such

tool, we used existing radar data derived from Mars Express
TEC measurements between the surface up to 900 km. Our
model is called model of Mars' ionosphere (MoMo) and can
be used to predict the Martian ionospheric TEC in a simple
way. The goal is to provide corrections to be applied for any
investigations using radio signal probing the Mars' ionosphere
at a given location, at a given time and for a given solar activity
level.

The first section of the paper describes the data used to build
our model, the model itself and the validation of our approach
with independent radio occultation data and the IPIM model.
The second section describes the climatology behaviour of the
Mars' ionosphere using our model, and provides some applica-
tions to radio-science experiments.

## 2 Data and model

This section is divided in three sub-sections. First we
describe the MARSIS data used. Secondly, we describe the
strategies and equations used to build the model and finally
we validate the model for different SZA and solar activity levels.

## 2.1 Data and method

We analysed the ionospheric data from 6,870,739 echoes (6850 tracks) from the Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) radar sounder in subsurface mode. The data set consists of vertical Total Electron Content (vTEC, i.e., the vertical integral of the electron density from the satellite to the ground, expressed in TEC units with 1 TECu = $10^{16} e^{-} m^{-2}$) for the period 19-06-2005 to 01-02-2014 (Safaeinili et al., 2007; Mouginot et al., 2008). Only data with acceptable signal to noise ratio threshold (i.e., SNR is greater than 15 dB, Grima & Kofman, 2008) where used. The data from year 2005 were also excluded after quality visual check (i.e., too much noisy and unrealistic vTEC values). Finally, the data set used consists of 5,308,594 echoes (more than 75% of the entire data set). The main limitation is that the MARSIS radar sounder in subsurface mode operates primarily at SZA higher than 60°. Indeed, although the SZA of the whole data set varies between 4° and 171°, a large part of the measurements (~85%) are acquired between 60° and 120° around the solar terminator (i.e., when SZA ≈ 108°). Sánchez-Cano et al. (2015) carried out a critical analysis of different vTEC data sets and processing techniques and recommend to consider an error bar between 0.3 and 1 TECu for SZA lower than 75°. Thus, we paid attention to validate our model at different ranges of SZA.

The first driver used to model the variation of the vTEC is
the EUV emission from the Sun accounting for the photo-
ionization for a given solar activity level. The index widely used
to model the vTEC on Earth is the daily integrated emission
from the solar disc at 2800 MHz (10.7 cm wavelength called
F10.7) recorded routinely since 1947 (Covington, 1969) for
1 AU distance. The F10.7 is expressed in solar flux unit (sfu,
with 1 sfu = $10^{-22}$ Wm$^{-2}$ Hz$^{-1}$). Nowadays, an additional
index called F10.7P (i.e., the mean of daily F10.7 of the day
of interest and its average over the previous 81 days), is consid-
ered as the new proxy which better represents the solar EUV
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Daily solar indices during the 23rd and 24th solar cycles in Solar Flux Unit (sfu). Daily F10.7 flux from Penticton radio telescope observations (black) and derived F10.7P (blue). The red dots are the F10.7P used for Mars. The black box represents the MARSIS data period used in this study.</figcaption></figure>

fluxes causing the photo-ionisation of the Earth neutral atmosphere when comparing with satellite SEM/SOHO EUV data (Liu et al., 2006). While this new index is nowadays used to predict the Earth ionosphere activity (Liu & Chen, 2009; Bilitza et al., 2011; Bergeot et al., 2013), the Mars' ionosphere is generally studied using the F10.7 (Lillis et al., 2010; Mendillo et al., 2013) corresponding to F10.7 corrected from the Sun to the Mars distance ($d_{S-M}$) and defined as $\overline{\text{F10.7}} = \frac{F_{10.7}}{d_{S-M}^{1/2}}$.

The $\overline{\text{F10.7P}}$ (Fig. 1) is used in this study as index to describe the solar EUV emission at Mars. Starting from the F10.7 from the Penticton radio telescope, we derive the $\overline{\text{F10.7P}}$ using the relation:

$$ \overline{\text{F10.7P}} = \frac{\overline{\text{F10.7}} + \overline{\text{F10.7A}}}{2}, \quad (1) $$

where $\overline{\text{F10.7A}}$ is the 81-previous days mean of daily $\overline{\text{F10.7}}$.

We made some tests taking into account for the time taken for the Sun to rotate from facing Earth to facing Mars. This delay, which varies between -14 and +14 days depending on the configuration, has a negligible effect on the resulting $\overline{\text{F10.7}}$ (lower than 3 sfu) implying an error lower than 4% on the output of our model. Thus, for more simplicity, we neglected this time delay/advance to estimate the $\overline{\text{F10.7}}$ as it is more simple for potential users of our model.

For the period of MARSIS data used in this study (see Fig. 1), the F10.7P (resp. $\overline{\text{F10.7P}}$) ranges from 65 sfu (resp. 24 sfu) to 205 sfu (resp. 77 sfu). Therefore, we considered in the present study a F10.7P (i.e., at 1 AU) of 70 sfu for low, 120 sfu for moderate and 170 sfu for high solar activity levels.

The other parameter taken into account for our model is the Solar Longitude (Ls), which characterises the position of the planet around the Sun, and thus the seasons. Regarding the repartition time and space domain with respect to Ls, SZA and $\overline{\text{F10.7P}}$ of the MARSIS data set, two different sets of Ls are considered: from mid-spring to mid-autumn comprising the summer (Ls in [45°; 225°]) and from mid-autumn to

mid-summer comprising the winter season (Ls in [225°; 45°]) in the Northern hemisphere. This allows taking into account the variations in the thermospheric reservoir (for the ionosphere) between summer and winter.

The last parameter considered is the inter-hemispheric dichotomy of Mars. To highlight the general distribution of the vTEC on Mars, the mean of the vTEC values for SZA lower than 75° for the entire period (2006–2014) on a 10° × 10° grid is shown in Figure 2.

There is a clear dichotomy between the Northern and Southern vTEC general level, with higher vTEC values in the South (~0.6 TECu; Min: 0.4 TECu; Max: 0.8 TECu) compared to the North (~0.5 TECu; Min: 0.4 TECu; Max: 0.6 TECu). Additionally, the standard deviation of the distribution (top left Fig. 2) is higher at low latitudes (i.e., close to the equator) compared to mid- and high latitudes. This indicates a relative stability of the mean vTEC at high latitudes and consequently the reliability of higher values in the Southern hemisphere compared to the North. This is due to the effects on the ionosphere of the crustal magnetic field (e.g., Flynn et al., 2017) which are stronger in the southern hemisphere (Acuña et al., 2001). To support our conclusion on the hemispherical dichotomy, we also considered the mean for the SZA and the $\overline{\text{F10.7P}}$ (Figure 2 middle and right). From the global mean repartition of these quantities, there is no clear dependence of the mean vTEC with respect to the SZA and $\overline{\text{F10.7P}}$ neither latitude nor longitude. In our model, we take into consideration this dichotomy by analysing separately the Northern and the Southern hemispheres.

## 2.2 Empirical model of Mars' ionosphere

The main driver of the evolution of the vTEC at a given location is the solar radiation whose penetration in the atmosphere depends on the SZA. The SZA as the main parameter to constrain the model. In order to fill lack of data at low SZA, we considered the formulation for the Chapman grazing
<figure><img src="image_4.png" /><figcaption>Fig. 2. vTEC, SZA and F10.7P distribution on Mars for SZA lower than 75°. Bottom: mean values every 10° in longitude and latitude for the entire period (2006–2014) for vTEC (left), SZA (middle) and F10.7P (right). Top: standard deviation of the mean at 1σ.</figcaption></figure>

incidence integral noted ch(Xp, SZA). The use of this function allows the determination of the absorption of the solar EUV in a terrestrial atmosphere, and consequently the ionisation rate, over a wide range of solar illuminations (Chapman, 1931). The simplified function is given by (Eq. 9 in Smith & Smith, 1972):

$$ch(X_p, SZA) = \sqrt{\frac{\pi}{2} X_p} \times f(SZA), \quad (2)$$

where Xp is the distance of the estimated ionospheric thin shell layer with respect to the centre of Mars (3392 + 140 km) divided by the atmospheric scale height (~15 km), SZA the Solar Zenith Angle and f(SZA) the functions defined in Smith & Smith, 1972. Using this formulation, the function is divided in two ranges of SZA: SZA ≤ 90° (see Eq. (12) in Smith & Smith, 1972) and SZA > 90° (see Eq. 15 in the same reference). The use of this formulation allows us to extrapolate the model for SZA lower than 60° where there is a lack of data, taking into account for the physical processes defined in Chapman (1931) to describe the absorption of Solar radiation in a planetary atmosphere.

Among all the tests we made, our best theoretical vTEC (in terms of differences between the predicted values and the data in input) is modelled according to:

$$v\text{TEC}(\text{SZA}, \text{Hem}, \text{Ls}, \overline{\text{F10.7P}}) = \alpha_1 + \frac{\alpha_2}{\sqrt{ch(X_p, \text{SZA})}}, \quad (3)$$

where ch(Xp, SZA) is the function given in equation (2), α1 and α2 coefficients are parameters estimated to fit the observations by minimizing the squared differences between the modelled and observed vTEC from MARSIS data in a least-square sense.

At this stage, the vTEC is predicted for a given range of F10.7P index taking into account for the heliocentric distance of Mars (see Eq. (1), at a given SZA, in a given hemisphere (Hem) and for a given Ls range (Ls = [45°; 225°] or Ls = [225°; 45°]). For that, different tools and data from SPICE/NAIF (Acton, 1996; Acton et al., 2018) were used to retrieve information such as solar longitude, heliocentric distance or SZA at a given location and time.

Figure 3 shows the estimated α1 and α2 coefficients for (1) each hemisphere, (2) for the two Ls ranges defined above and (3) for different ranges of F10.7P (from 20 to 80 sfu with a step of 5 sfu).

As first order approximation, the α1 coefficients can be considered as constant with respect to the F10.7P (see Fig. 3 left panels), while the α2 coefficients are linearly dependent on F10.7P (see Fig. 3 right panels). Therefore, to make MoMo applicable for any solar activity level and heliocentric distance, equation (3) is refined as follows:

$$\text{TEC}(\text{SZA}, \text{Hem}, \text{Ls}, \overline{\text{F10.7P}}) = \text{mean}(\alpha_1) + \frac{\beta_1 + \beta_2 \times \overline{\text{F10.7P}}}{\sqrt{ch(X_p, \text{SZA})}}, \quad (4)$$

where mean(α1) is the mean of the α1 and α2 is replaced by a linear equation parameterised by β1 and β2 and defined as α2 = β1 + β2 × F10.7P.

Finally, MoMo consists of 12 coefficients: three coefficients for each hemisphere and solar longitude range. The mean of the differences between the model and the observations is 0.01 ± 0.06 TECu (Figure 4) with only ~9% of the absolute value of differences larger than 0.1 TECu between MoMo outputs and MARSIS observations. The coefficients to be used in
<figure><img src="image_5.png" /></figure>

**Fig. 3.** Coefficients $\alpha_1$ and $\alpha_2$ estimated from the least-square adjustment of equation (3). Top: for the Northern hemisphere; Bottom: for the Southern hemisphere. The colours are related to the Martian seasons with $45^\circ \le L_s < 225^\circ$ for the red curve and $225^\circ \le L_s < 45^\circ$ for the green curve. The lines are the estimated $\alpha_1$ and $\alpha_2$ coefficients with respect to the $F10.7P$. The crosses are the coefficients obtained for the different $F10.7P$ index while the lines are the values used in the model.

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

**Fig. 4.** Histograms of the differences between MoMo output and MARSIS data set used to constrain the model.

equation (4) and the RMS of the model are available in Table 1. The model (subroutine *sub_momo.f*) and the $F10.7P$ index are available at https://lara.oma.be/marsatmo/iono/momo.html.

## 2.3 MoMo model validation

To validate the model, the output from MoMo is compared with a completely independent data set, consisting of 250 electron density profiles of the dayside ionosphere (SZA [50°–85°]) extracted from Mars Express radio occultation experiments MaRs (Peter et al., 2014). The mean differences between MoMo and MaRs data is $0.01 \pm 0.05$ TECu for SZA greater than $50^\circ$ (Fig. 5), which is as the same level as the differences with MARSIS data. The mean of the differences in the North is $0.00 \pm 0.04$ TECu, while the differences in the South are $0.03 \pm 0.06$ TECu. This could be due to more intense ionospheric variability in the South not well represented in MoMo which aims at representing the climatology of the vTEC. Furthermore, these differences are lower than 10% of the vTEC values. Consequently, our model reflects well the dayside ionosphere vTEC especially for mid- to low-solar elevations (i.e., larger than $50^\circ$ of SZA).

MoMo is also compared with the IPIM numerical model (Marchaudon & Blelly, 2015) which includes physical processes in the ionosphere (e.g., electrodynamics, energetics, transport, photo-chemistry, etc). One of the outputs of this model are the density profiles which we integrated with respect to the altitude to compare with MoMo runs (Fig. 6). We consider here three runs for IPIM: case 1 corresponding to location in the Southern hemisphere during Northern Winter and low solar activity (22-06-2009, $F10.7 = 70$ sfu, i.e., $F10.7P = 34.8$ sfu, latitude $= 45^\circ$ S, $L_s = 289.6^\circ$); case 2 for a location in the Northern hemisphere during Northern Spring and intermediate solar activity (18-08-2013, $F10.7 = 129$ sfu; $F10.7P = 47.9$ sfu, latitude $= 45^\circ$ N, $L_s = 8.6^\circ$); case 3 for a location in the Northern hemisphere during Northern Autumn and intermediate solar activity (24-10-2012, $F10.7 = 129$ sfu; $F10.7P = 61.8$ sfu, latitude $= 20^\circ$ N, $L_s = 194.0^\circ$).

The shapes of the vTEC values of the two models are very consistent with differences between $-0.12$ and $0.06$ TECu (mean of $0.01 \pm 0.03$ TECu) showing the efficiency of the MoMo estimations even at low SZA. Despite MoMo's simplicity, the behaviour of the Mars' ionosphere is well predicted. This comforts us in the reliability of MoMo outputs even at low SZA. Presently, no data are available for SZA lower than $30^\circ$. This is the reason why we considered physical model (Smith & Smith, 1972) to extrapolate our vTEC at low SZA. In the future, satellites data with SZA $= 0^\circ$ will be necessary to validate MoMo and all the models available.

## 3 Model applications

In this second section, MoMo outputs are used for different applications. The first sub-section is dedicated to the climatology of the Mars' ionosphere while the second sub-section deals with the impact of the Mars TEC on radio-science experiments.

### 3.1 vTEC climatology on Mars

To access the climatological behaviour of the Martian ionospheric, MoMo is used to predict the Martian vTEC for different seasons, hemispheres and solar activity levels. Figure 7 presents the daily variations (for $0^\circ \le SZA \le 110^\circ$) of the
<table><caption>**Table 1.** Coefficients for the MoMo model to predict the vTEC in the Mars' ionosphere. These coefficients are given for different ranges of solar longitude (Ls) and the two hemispheres. The number of data per range of SZA and the associated RMS of the differences between the vTEC predicted by the model and the observed vTEC are reported in the last columns.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th rowspan="2">Mean(α<sub>1</sub>)</th><th rowspan="2">β<sub>1</sub></th><th rowspan="2">β<sub>2</sub></th><th colspan="2">#Data</th><th colspan="2">RMS (TECu)</th></tr><tr><th>10<sup>-2</sup></th><th>10<sup>-1</sup></th><th>SZA &lt; 75°</th><th>SZA ≥ 75°</th><th>SZA &lt; 75°</th><th>SZA ≥ 75°</th></tr></thead><tbody><tr><td rowspan="2">North</td><td>45° ≤ Ls &lt; 225°</td><td>3.284</td><td>2.624</td><td>1.564</td><td>159,997</td><td>983,540</td><td>0.10</td><td>0.05</td></tr><tr><td>225° ≤ Ls &lt; 45°</td><td>3.004</td><td>2.950</td><td>1.439</td><td>258,079</td><td>1,070,191</td><td>0.12</td><td>0.05</td></tr><tr><td rowspan="2">South</td><td>45° ≤ Ls &lt; 225°</td><td>2.473</td><td>5.521</td><td>0.964</td><td>217,032</td><td>1,280,671</td><td>0.09</td><td>0.05</td></tr><tr><td>225° ≤ Ls &lt; 45°</td><td>3.577</td><td>-0.222</td><td>2.287</td><td>220,611</td><td>1,118,474</td><td>0.11</td><td>0.05</td></tr></tbody></table>

MoMo, model of Mars ionosphere; RMS, root mean square; SZA, solar Zenith angle; TECu, total electron content units; vTEC, vertical total electron content.

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** MoMo empirical model results for the Northern (left) and Southern (middle) hemispheres. The colours correspond to the different seasons (i.e., Ls); the solid, dashed and dotted lines represent high, moderate and low solar activity level respectively. Also represented are the radio occultation data (Peter et al., 2014) for different solar activity levels: high (triangle), moderate (square) and low (circle). On the right is the scatter plot showing the MoMo prediction model versus the vTEC from the Mars Express radio occultation data for the two hemispheres.</figcaption></figure>

vTEC at different locations (top Northern hemisphere, bottom Southern Hemisphere corresponding to 20° N and 20° S respectively), different solar activity level (F10.7P at Earth distance of 70, 120 and 170 sfu) and different seasons. The seasonal differences are given by the two ranges of Ls used to construct the model but also in the F10.7P index (Eq. (1)) which is dependent on the heliocentric distance, and consequently on the Ls variations during Mars sidereal year.

At first glance, it appears in Figure 7 that: (1) The maximum of the vTEC (0.9 TECu during low solar activity and 2 TECu during high solar activity) is observed in the Southern hemisphere in Summer. This can be explained by the fact that during this season, the planet is at perihelion, with its rotation axis oriented to receive more solar flux in the South than the North; (2) The night side of the Mars' ionosphere (SZA ≈ 110°) is constant and independent from the solar level activity with a mean value of 0.03 ± 0.01 TECu. This value is consistent with the conclusion from Mendillo et al. (2013); (3) The Northern and Southern hemispheres (NH and SH) are behaving more or less in the same way with an increase phase in vTEC values from Autumn to Winter and then a decreasing phase from Spring to Summer resulting from the fact that the Mars' ionosphere is mainly driven by the solar flux angle (i.e., SZA), the solar activity and the heliocentric distance. The peaks in Spring/NH and Autumn/SH compared to the next season confirm the recent results of Sánchez-Cano et al. (2018),

<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Comparison between MoMo empirical model and IPIM numerical model. The three cases correspond to different locations, solar activity levels and seasons (see description in the text). In red are the MoMo outputs while in blue are the IPIM outputs.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** Seasonal climatological behaviour of the vTEC on Mars for Northern (top, 20° N) and Southern (bottom, 20° S) hemispheres. For each season and hemisphere, each curve corresponds to the daily variation of the vTEC (i.e., SZA from 110° to 0°). The colours correspond to different solar cycle activity levels (F10.7P = 70, 120 and 170 sfu).</figcaption></figure>

who identified a coupling between the lower and upper atmospheres for this range of Ls. The authors concluded that these peaks in the vTEC at these periods are due to an increase in the thermospheric density and changes in ion composition.

Consequently, MoMo shows a Mars' ionosphere variability mainly driven by the solar illumination, activity, and the seasons as expected, with amplitude variations of the vTEC over an entire day lower than ~2 TECu.

## 3.2 Application to radio-science experiments

MoMo is now used to estimate the impact of the ionosphere on radio-science experiments between Mars surface assets (e.g., lander, rover, ...) and Mars orbiters or Earth antennas. For that purpose, we produce vTEC maps for Mars with MoMo. Four of those maps, one per season, are shown in Figure 8 for a specific case of noon local time at the meridian origin with a maximum solar activity level of 170 sfu at 1 AU.

These maps are then employed to estimate the slant TEC (sTEC, see Fig. 9) in the line of sight of an external object, and to quantify the impact on different wave propagation parameters: phase delay, and Doppler shift. Different radio frequencies from UHF to X-band are used nowadays for Mars for telecommunications and radio-science experiments. Communications between surface assets and Mars orbiters are commonly performed in UHF, while direct-to-Earth radio tracking is now performed in X-band (S-band was used for Viking landers). In both cases, the electromagnetic waves interact with the sTEC present in the Mars' ionosphere.

We consider two hypothetical surface assets positioned at S10° and N10° in latitudes and 0° in longitude (see crosses in Fig. 8). The transmitting/receiving antennas on Earth (or on-board a Mars orbiting spacecraft) are assumed to be at

different elevation angles seen from the asset: 20° and 55° which correspond to the typical elevation range for present and future radio-science operation of RISE (Folkner et al., 2018) and LaRa (Dehant et al., 2019). We also consider an elevation angle of 90° (i.e., at the zenith of the asset). The azimuth of the external object is chosen to be 90° E, i.e., the asset Eastern direction. Only the winter solstice (Ls = 270°) is considered when the vTEC values are maximum. Additionally, the solar activity level is fixed to be in moderate phase (F10.7P = 120 sfu). The frequencies considered are within UHF, S- and X-bands (0.4, 2.0 and 8.0 GHz, respectively). The vTEC from MoMo is projected into slant TEC (sTEC) along the line-of-sight (LOS) at the Ionospheric Pierce Point (IPP) considering a thin shell layer at 140 km altitude corresponding to the expected altitude of the maximum electron density (Morel et al., 2004, Fig. 8).

Finally the phase delay is estimated (i.e., time delay of the phase, see Tables 2 and 3 in Ho et al., 2002) in meters by using the relation $I = \frac{40.3 \times s_{TEC} \times 10^{16}}{f^2}$ with f the frequency considered. We also estimate the instant Doppler shift due to the Mars' ionosphere TEC change defines as $\Delta f_{\text{iono}} = \frac{40.3 \times d(\text{sTEC} \times 10^{16})}{c \times f} \frac{d(\text{sTEC} \times 10^{16})}{dt}$ expressed in Hz (mHz is used here) with c the speed of light. Finally, the velocity change due to the Mars' Ionosphere is obtained using $\Delta v_{\text{iono}} = \frac{c \times \Delta f_{\text{iono}}}{2 \times f}$ expressed in m s⁻¹ (mm s⁻¹ in this study) and assuming a two-way radio-link. These different quantities are given with respect to the Solar Local Time (SLT) at the asset location.

As shown in Figure 10, the maximum phase delay is obtained at 12:20 SLT with a stronger effect at lower elevation and for an asset situated in the Southern hemisphere. In that
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Global maps of Mars vTEC. The different maps correspond to the different outputs of the MoMo model for different seasons (in the North hemisphere) at local noon at the meridian origin. The F10.7P is fixed for a high solar activity level (170 sfu) corresponding to F10.7P of 70 sfu for spring, 61 sfu for summer, 79 sfu for autumn and 88 sfu for winter. The black crosses are the location of two assets on the Mars surface (N10° and S10° in latitude; 0° in longitude) used in Figures 9 and 10.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Projection used to retrieve the sTEC in a line-of-sight from the MoMo vTEC predictions. IPP stands for Ionospheric Pierce Point. *e* is the elevation angle of the external object with respect to the asset.</figcaption></figure>

case, the values in phase delay reach more than 8 m in UHF, 32 cm in S-band and 2 cm in X-band. At 90° of elevation (i.e., the Earth or the satellite is at the zenith of the asset) the effect decreases to 3 m, 11 cm and <0.7 cm in the Northern hemisphere. Mars ephemeris is mainly constrained by the range measurements of in-situ spacecraft, including those acquired in the late seventies by the Viking landers 1 and 2. Together, these two landers provided about 1000 range data points in S-band (e.g., Kuchynka et al., 2014), which must be affected by the Martian ionosphere at the tens of centimetre level. This is only 3–10 times below the current accuracy of

the ephemeris of Mars, suggesting that one could use MoMo to correct those data when producing the new solution of the Mars trajectory.

The effect on the Doppler measurements is maximum when the TEC suddenly varies. This appends twice per day (see right panels of Fig. 10) during local sunrise and sunset (around 06:00 and 18:45 SLT). The maximum shift is of the order of ±1 mHz (corresponding to ±3.6 × 10⁻¹ mm s⁻¹ in velocity change) in UHF, ±2 × 10⁻¹ mHz (±1.5 × 10⁻² mm s⁻¹) in S-band and ±5 × 10⁻² mHz (±1.0 × 10⁻³ mm s⁻¹) in X-band when considering 20° elevation with respect to the asset.
<figure><img src="image_12.png" /></figure>

**Fig. 10.** Effect of the Mars' ionosphere on radio wave propagation. Effect on UHF (0.4 GHz, bottom), S-band (2 GHz, middle) and X-band (8 GHz, top) radio signals. Left: effect on phase delay in centimetres; Right: impact on instant Doppler shift in mHz (assuming $dt = 1$ s). The red curves correspond to an asset situated in the Northern Hemisphere (E0°, N10°) while the blue curves correspond to the one situated in the Southern Hemisphere (E0°, S10°). The bold, dashed-dots and dashed lines correspond to a LOS elevation of 90°, 55° and 20° respectively (see *e the elevation* on Fig. 9). The solar activity level is fixed at $F10.7P = 120$ sfu (moderate solar activity, see Fig. 1).

Figure 10 shows the seasonal variations of the Mars' ionosphere contribution to the Doppler measurements of an asset located at (E0°, N10°) and (E0°, S10°). The computation is systematically performed at 06:00 SLT for a signal LOS elevation angle of 20° above the asset. The solar activity level is maintained constant ($F10.7P = 120$ sfu), although two periods are distinguished to account for the variations in the thermospheric reservoir between summer and winter (see Table 1). This change is responsible for the periodic signal clearly shown on Figure 10. Modelling such a seasonal signal can be of great importance for radio-science experiments from a lander on Mars as their main scientific objective is the accurate determination of the seasonal variations of Mars rotation rate and spin axis orientation variations called nutations (e.g., Folkner et al., 2018). Indeed, the rotation/orientation variations also affect the Doppler measurements with a periodic signature of one Martian year. Therefore, a large ionospheric contribution to the Doppler shift could in turn introduce bias in the rotation/orientation parameter estimates. Our model predicts a maximum ionospheric contribution at 0.05 mHz level (i.e., $0.001 \text{ mm s}^{-1}$) in

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

**Fig. 11.** Ionospheric signature on Doppler shift during a Martian year. Data extracted for different Ls, for a 20° elevation with respect to the asset at 06:00 SLT. The dashed red and dotted bleu lines correspond to an asset in the Northern and Southern hemispheres respectively. The ionospheric signature on Doppler shift is given for UHF (0.4 GHz, bottom), S-band (2 GHz, middle) and X-band (8 GHz, top) radio signals.

the frequency band of both RISE and LaRa (upper panel of Fig. 11). This is about one order of magnitude below the estimated noise of those instruments, but is of same order of magnitude as the contribution of the liquid core in the Doppler ($\sim 10^{-3}-10^{-2} \text{ mm s}^{-1}$, Yseboodt et al., 2017). Consequently, the predictions made with MoMo suggests that the radio-science teams will have to either correct their data using our model for instance or adapt the mission programmatic to avoid operating during sunrise and sunset when the TEC rapidly varies.

## 4 Conclusion

We have developed a new empirical model of the Mars Total Electron Content (TEC) called MoMo (available at https://lara.oma.be/marsatmo/iono/momo.html). The model provides values for the vertical TEC (vTEC) for a given solar zenith angle, solar activity, solar longitude Ls in the two Mars hemispheres. The model has been validated with external radio occultation data for high SZA and with the IPIM (IRAP Plasmasphere Ionosphere Model) physical complex model for low SZA. One of the main motivations for this paper was to provide Mars ionosphere corrections for radio-science experiments. We showed that, even if the expected noise of the radio-science instruments is large compare to the Mars ionospheric contribution, it is recommended to correct for this contribution, as its seasonal variations are of the same periods as those of the geophysical parameters determined by the
experiments. It is advised to apply ionospheric corrections even during moderate solar activity level, above all for long term studies (i.e., more than one Mars sidereal year), such as those aiming at studying the deep interior of the planet. The model provides key radio-science quantities like phase delay and Doppler shift, which can be used to calibrate radio science data like those acquired by RISE on the ongoing InSight mission, by LaRa on the future mission ExoMars 2020 mission, or by historical lander missions (Viking, Pathfinder, MER). In particular, the Mars ephemeris calculations could certainly benefit from the calibration using MOMO of the S-band Viking landers ranging data.

**Acknowledgements.** This work has been financially supported by the Solar and Terrestrial Centre of Excellence (STCE, http://www.stce.be/) as well as by the Belgian PRODEX program managed by the European Space Agency (ESA) in collaboration with the Belgian Federal Science Policy Office (BELSPO). Also, support from the ESA-ESTEC Faculty is acknowledged. We are grateful for the providers of the tools and findings from SPICE or NAIF/PDS resources, and of the teams that provided the data used to make SPICE kernels. The ensemble of the pictures is produce using the Generic Mapping Tool (GMT, Wessel et al., 2013). We also thank ESA for providing the TEC data set from MARSIS (ftp://psa.esac.esa.int/) and the NGDC/NOAA for providing the F10.7 data (ftp://ftp.ngdc.noaa.gov). Finally the authors would like to highlight the Rob Lillis from UC Berkeley Space Sciences Laboratory and the anonymous person for their constructive review which clearly contributed to the improvement of the paper. The editor thanks Robert Lillis and an anonymous referee for their assistance in evaluating this paper.

## References

Acton C, Bachman N, Semenov B, Wright E. 2018. A look towards the future in the handling of space science mission geometry. *Planet Space Sci* **150**: 9–12. DOI: 10.1016/j.pss.2017.02.013.

Acton CH Jr. 1996. Ancillary data services of NASA's navigation and ancillary information facility. *Planet Space Sci* **44**(1): 65–70.

Acuña MH, Connerney JEP, Wasilewski P, Lin RP, Mitchell D, Anderson KA, Carlson CW, McFadden J, Rème H, Mazelle C, Vignes D. 2001. Magnetic field of Mars: Summary of results from the aerobraking and mapping orbits. *J Geophys Res Planet* **106** (*E10*): 23403–23417. DOI: 10.1029/2000JE001404.

Bergeot N, Tsagouri I, Bruyninx C, Legrand J, Chevalier JM, Defraigne P, Baire Q, Pottiaux E. 2013. The influence of space weather on ionospheric total electron content during the 23rd solar cycle. *J Space Weather Space Clim* **3**: A25. DOI: 10.1051/swsc/2013047.

Bilitza D, McKinnell LA, Reinisch B, Fuller-Rowell T, Bilitza D, McKinnell LA, Reinisch B, Fuller-Rowell T. 2011. The international reference ionosphere today and in the future. *J Geodesy* **85** (*12*): 909–920. DOI: 10.1007/s00190-010-0427-x.

Chapman S. 1931. The absorption and dissociative or ionizing effect of monochromatic radiation in an atmosphere on a rotating earth. *Proc Phys Soc* **43**(1): 26. DOI: 10.1088/0959-5309/43/1/305.

Covington AE. 1969. Solar Radio Emission at 10.7 cm, 1947–1968. *J R Astron Soc Can* **63**: 125–132.

Dehant V, Folkner W, Renotte E, Orban D, Asmar S, Balmino G, Barriot JP, Benoist J, Biancale R, Biele J, Budnik F. 2009. Lander radioscience for obtaining the rotation and orientation of Mars. *Planet Space Sci* **57**(8–9): 1050–1067. DOI: 10.1016/j.pss.2008.08.009.

Dehant V, Le Maistre S, Rivoldini A, Yseboodt M, Rosenblatt P, Van Hoolst T, Mitrovic M, Karatekin Ö, Marty JC, Chicarro A. 2011. Revealing Mars' deep interior: Future geodesy missions using radio links between landers, orbiters, and the Earth. *Planet Space Sci* **59**(10): 1069–1081. DOI: 10.1016/j.pss.2010.03.014.

Dehant V, Le Maistre S, Balanda RM, Bergeot N, Karatekin O, et al. 2019. The radio-science LaRa instrument onboard ExoMars 2020 to investigate the rotation and interior of Mars. *Planet. Space Sci.* Accepted for Publication.

Flynn CL, Vogt MF, Withers P, Andersson L, England S, Liu G. 2017. MAVEN observations of the effects of crustal magnetic fields on electron density and temperature in the Martian dayside ionosphere. *Geophys Res Lett* **44**(21): 10812–10821. DOI: 10.1002/2017GL075367.

Folkner WM, Dehant V, Le Maistre S, Yseboodt M, Rivoldini A, Van Hoolst T, Asmar SW, Golombek MP. 2018. The rotation and interior structure experiment on the InSight mission to Mars. *Space Sci Rev* **214**(5): 100. DOI: 10.1007/s11214-018-0530-5.

Grima C., Kofman W. 2008. MARSIS derived enhanced ionospheric calibration data. Technical Report No 21646/08/NL/NR, European Space Agency 3, pp. 1–10.

Ho C, Golshan N, Kliore A. 2002. Radio wave propagation handbook for communication on and around Mars, Technical Report 02-5, Jet Propulsion Laboratory, Pasadena, CA.

Klobuchar JA, Parkinson BW, Spilker JJ. 1996. Ionospheric effects in global positioning system: Theory and applications, American Institute of Aeronautics and Astronautics, Washington, DC. Google Scholar.

Kuchynka P, Folkner WM, Konopliv AS, Parker TJ, Park RS, Le Maistre S, Dehant V. 2014. New constraints on Mars rotation determined from radiometric tracking of the Opportunity Mars Exploration Rover. *Icarus* **229**: 340–347. DOI: 10.1016/j.icarus.2013.11.015.

Lillis RJ, Brain DA, England SL, Withers P, Fillingim MO, Safaeinili A. 2010. Total electron content in the Mars ionosphere: Temporal studies and dependence on solar EUV flux. *J Geophys Res Space Phys* **115**, A11314. DOI: 10.1029/2010JA015698.

Liu L, Chen Y. 2009. Statistical analysis of solar activity variations of total electron content derived at Jet Propulsion Laboratory from GPS observations. *J Geophys Res Space Phys* **114**, A10311. DOI: 10.1029/2009JA014533.

Liu L, Wan W, Ning B, Pirog OM, Kurkin VI. 2006. Solar activity variations of the ionospheric peak electron density. *J Geophys Res* **111**(8): A08304. DOI: 10.1029/2006JA011598.

Marchaudon A, Bleilly PL. 2015. A new interhemispheric 16-moment model of the plasmasphere-ionosphere system: IPIM. *J Geophys Res Space Phys* **120**(7): 5728–5745. DOI: 10.1002/2015JA021193.

Mendillo M, Marusiak AG, Withers P, Morgan D, Gurnett D. 2013. A new semiempirical model of the peak electron density of the Martian ionosphere. *Geophys Res Lett* **40**(20): 5361–5365. DOI: 10.1002/2013GL057631.

Mendillo M, Narvaez C, Trovato J, Withers P, Mayyasi M, Morgan D, Kopf A, Gurnett D, Němec F, Campbell B. 2018. Mars Initial Reference Ionosphere (MIRI) model: Updates and validations using MAVEN, MEX, and MRO data sets. *J Geophys Res Space Phys* **123**(7): 5674–5683. DOI: 10.1029/2018JA025263.
Morel L, Witasse O, Warnant R, Cerisier JC, Blelly PL, Lilensten J. 2004. Diagnostic of the dayside ionosphere of Mars using the Total Electron Content measurement by the NEIGE/Netlander experiment: An assessment study. *Planet Space Sci* **52**(7): 603–611.

Morgan DD, Witasse O, Nielsen E, Gurnett DA, Duru F, Kirchner DL. 2013. The processing of electron density profiles from the Mars Express MARSIS topside sounder. *Radio Sci* **48**(3): 197–207. DOI: 10.1002/rds.20023.

Mouginot J, Kofman W, Safaeinili A, Hérique A. 2008. Correction of the ionospheric distortion on the MARSIS surface sounding echoes. *Planet Space Sci* **56**(7): 917–926. DOI: 10.1016/j.pss.2008.01.010.

Pätzold M, Neubauer FM, Carone L, Hagermann A, Stanzel C, Häusler B, Remus S, Selle J, Hagl D, Hinson DP, Simpson RA. 2004. MaRS: Mars express orbiter radio science. In: *Mars Express: The Scientific Payload*, ESA Publications Division **1240**, pp. 141–163.

Peter K, Pätzold M, Molina-Cuberos G, Witasse O, González-Galindo F, Withers P, Bird MK, Häusler B, Hinson DP, Tellmann S, Tyler GL. 2014. The dayside ionospheres of Mars and Venus: Comparing a one-dimensional photochemical model with MaRS (Mars Express) and VeRa (Venus Express) observations. *Icarus* **233**: 66–82. DOI: 10.1016/j.icarus.2014.01.028.

Safaeinili A, Kofman W, Mouginot J, Gim Y, Herique A, Ivanov AB, Plaut JJ, Picardi G. 2007. Estimation of the total electron content of the Martian ionosphere using radar sounder surface echoes. *Geophys Res Lett* **34**(23), L23204. DOI: 10.1029/2007GL032154.

Sánchez-Cano B, Witasse O, Herraiz M, Radicella SM, Bauer J, Blelly PL, Rodríguez-Caderot G. 2012. Retrieval of ionospheric profiles from the Mars Express MARSIS experiment data and

comparison with radio occultation data. *Geosci Instrum Meth Data Syst* **1**(1): 77–84. DOI: 10.5194/gi-1-77-2012.

Sánchez-Cano B, Radicella SM, Herraiz M, Witasse O, Rodríguez-Caderot G. 2013. NeMars: An empirical model of the Martian dayside ionosphere based on Mars Express MARSIS data. *Icarus* **225**(1): 236–247. DOI: 10.1016/j.icarus.2013.03.021.

Sánchez-Cano B, Morgan DD, Witasse O, Radicella SM, Herraiz MF, Orosei R, Cartacci M, Cicchetti A, Noschese R, Kofman W, Grima C. 2015. Total electron content in the Martian atmosphere: A critical assessment of the Mars Express MARSIS data sets. *J Geophys Res Space Phys* **120**(3): 2166–2182. DOI: 10.1002/2014JA020630.

Sánchez-Cano B, Lester M, Witasse O, Blelly PL, Indurain M, Cartacci M, González-Galindo F, Vicente-Retortillo Á, Cicchetti A, Noschese R. 2018. Spatial, seasonal and solar cycle variations of the Martian total electron content (TEC): Is the TEC a good tracer for atmospheric cycles? *J Geophys Res Planet* **123**(7): 1746–1759. DOI: 10.1029/2018JE005626.

Smith FL III, Smith C. 1972. Numerical evaluation of Chapman's grazing incidence integral ch (X, y). *J Geophys Res* **77**(19): 3592–3597. DOI: 10.1029/JA077i019p03592.

Wessel P, Smith WH, Scharroo R, Luis J, Wobbe F. 2013. Generic mapping tools: improved version released. *Eos Trans Am Geophys Union* **94**(45): 409–410.

Yseboodt M, Dehant V, Péters MJ. 2017. Signatures of the Martian rotation parameters in the Doppler and range observables. *Planet Space Sci* **144**: 74–88. DOI: 10.1016/j.pss.2017.05.008.

Zuber MT, Lemoine FG, Smith DE, Konopliv AS, Smrekar SE, Asmar SW. 2007. Mars reconnaissance orbiter radio science gravity investigation. *J Geophys Res* **112**(E5), E05S07. DOI: 10.1029/2006je002833.

**Cite this article as:** Bergeot N, Witasse O, Le Maistre S, Blelly P-L, Kofman W, et al. 2019. MoMo: a new empirical model of the Mars ionospheric total electron content based on Mars Express MARSIS data. *J. Space Weather Space Clim.* **9**, A36.