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

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 uses 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.


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, . 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., 2009Dehant et al., , 2011Dehant et al., , 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 Planetary Space Weather ESA-Roscosmos ExoMars 2020 mission. The X-band measurements provided by these two missions have unprecedented accuracy (and low power level in the case of LaRa). In particular, 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 À1 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 À1 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 radioscience data quality. Pretty low when Mars elongation is abovẽ 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 & Blelly, 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 applications to radio-science experiments.

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.

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 photoionization 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 considered as the new proxy which better represents the solar EUV N. Bergeot et al.: J. Space Weather Space Clim. 2019, 9, A36 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 F10:7 ¼ F 10:7 The 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 F10:7P using the relation: where F10:7A is the 81-previous days mean of daily 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 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 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: 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 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 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 F10:7P neither latitude nor longitude. In our model, we take into consideration this dichotomy by analysing separately the Northern and the Southern hemispheres.

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 N. Bergeot et al.: J. Space Weather Space Clim. 2019, 9, A36 incidence integral noted ch(X p , 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): where X p 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: vTEC SZA; Hem; Ls; F10:7P where ch(X p , SZA) is the function given in equation (2), a 1 and a 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 a 1 and a 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 a 1 coefficients can be considered as constant with respect to the F10:7P (see Fig. 3 left panels), while the a 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: TECðSZA; Hem; Ls; F10:7PÞ where mean(a 1 ) is the mean of the a 1 and a 2 is replaced by a linear equation parameterised by b 1 and b 2 and defined as 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  Table 1. The model (subroutine sub_momo.f) and the F10:7P index are available at https://lara.oma.be/marsatmo/iono/momo.html.

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 ± 0.05 TECu for SZA greater than 50° (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 ± 0.04 TECu, while the differences in the South are 0.03 ± 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°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 ( 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 ± 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°. 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°will be necessary to validate MoMo and all the models available.

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.

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° SZA 110°) of the  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), 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.   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.

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°w hich 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)  expressed in m s À1 (mm s À1 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 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 À1 mm s À1 in velocity change) in UHF, ±2 Â 10 À1 mHz (±1.5 Â 10 À2 mm s À1 ) in S-band and ±5 Â 10 À2 mHz (±1.0 Â 10 À3 mm s À1 ) in X-band when considering 20°elevation with respect to the asset.   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 mm s À1 ) in 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 (~10 À3 -10 À2 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.

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 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). 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.
N. Bergeot et al.: J. Space Weather Space Clim. 2019, 9, A36 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.