tracking a lander on the surface of another planet with radio signals is an efficient way to observe its rotation ( dehant et al . , 2009 , 2011 ) .
this can be done by measuring the doppler shift of the radio signal between the lander and a large antenna on earth , like the ones of the deep space network ( dsn ) .
we consider here two - way x - band signal : the ground station transmits the radio signal to mars , which is coherently ( i.e.  without any phase shift ) sent back to earth by the transponder onboard the lander .
five martian landers have already been tracked from earth in order to get the martian rotation .
the first landers were the two viking landers in 1976 - 1982 ( see for example borderies et al .
, 1980 , yoder and standish , 1997 ) .
21 years later , the pathfinder lander stayed 3 months on mars surface in 1997 ( folkner et al .
, 1997a and 1997b ) .
more recently the two mars exploration rovers were tracked from earth when they were stationary ( le maistre , 2013 , kuchynka et al . , 2014 ) :
spirit was stuck at the end of its life in 2009 , and opportunity stopped moving during the winter for energy saving during four months in 2012 .    two martian missions including a geodesy experiment will be launched in the coming years .
the first one is the nasa mission insight ( interior exploration using seismic investigations , geodesy and heat transport ) to be launched in 2018 .
the insight spacecraft will have a small deep space transponder on board : the rise ( rotation and interior structure experiment , folkner et al . , 2012 ) transponder .
the second mission is the esa mission exomars2020 that will include the lara ( lander radioscience , dehant et al . , 2009 )
transponder on the surface platform .
this spacecraft is to be launched in 2020 . +
the martian orientation parameters ( mop ) are the nutations , the precession , the length - of - day variations and the polar motion .
the precession is the long - term drift of the rotation axis in space .
nutations are seasonal periodic motions of this axis observed from space , while the polar motion is the seasonal motion of the rotation axis in a frame tied to the planet .
the rotation rate variations , also called length - of - day variations , are the periodic variations of the martian diurnal rotation .
the periods are also the harmonics of the annual period . since these rotational motions depend on the interior structure of mars , in particular the core dimension , density and state , and on the dynamics of the atmosphere and the ice caps and the co@xmath0 sublimation / condensation process , accurate measurements of these angles are very useful because additional constraints on the geophysical models could be provided .    in this paper
we derive expressions of the signature of the rotation parameters in the doppler and range observables .
these expressions are useful in order to optimize the configuration parameters and to anticipate the observation times that maximize the signatures and the geophysical information .
this is of course not only helpful for the insight and exomars mission planning but also for other geodesy experiments on other planets .
+ the paper is organized as follows : in section 2 , we present the different angles and variables that describe the geometry of the doppler observable .
section 3 gives a summary of the models for the different mop . in section 4 ,
we give the analytical expressions of the signature of these mop in the doppler observable and their numerical values , while the section 5 focuses on the signatures in the range observable .
the instantaneous doppler frequency can be interpreted as the space probe radial velocity with respect to the tracking station . as a first order approximation , it is the projection on the line - of - sight of the difference of the velocity vector between the emitter and the receiver .
it is very sensitive to the configuration geometry , therefore we present in this section the different variables that are important for the signatures in the doppler observable .     with respect to mars as a function of time .
the sun - earth - probe ( sep ) angle is plotted on the upper graph in red with some arbitrary scaling , the black parts of the curve ( and the corresponding pink boxes ) being the points where the sep angle is smaller than @xmath1 , corresponding to radio data with a large noise due to the solar plasma .
the gray boxes correspond to time intervals where the earth declination is smaller than 10@xmath2 .
the insight and exomars nominal mission lifetime are also shown with the orange and red horizontal lines .
, width=680 ]    the exomars and insight missions have similar orbital configuration at launch and when arriving on mars . for both the distance between mars and the earth increases after the martian landing , being maximal a few months after the landing ( see fig .  [ fig_dist ] , upper graph ) .
the sun - earth - probe ( sep ) angle or elongation is small when the earth - mars distance is maximal , this is the conjunction time .
the pink boxes on fig .
[ fig_dist ] show where the sep angle is smaller than @xmath1 . at that time , the solar plasma perturbs the radio signal largely , and the noise on the radioscience measurements is much larger than the nominal noise . +
the earth declination @xmath3 with respect to mars is the angle between the mars - earth direction and the equatorial plane of mars .
because of the martian obliquity and the earth revolution in the ecliptic plane , the earth declination has variations due to the orbital motions of the two planets and due to the mars pole direction .
it varies between -25@xmath2 and 25@xmath2 ( see fig .
[ fig_dist ] , lower graph ) .
when the earth - mars distance is minimal , corresponding to the opposition , the earth - mars direction reverses its rotation : for a few weeks , this direction goes from a prograde to a retrograde motion .
this is due to the faster revolution velocity of the earth than for mars . at this time , in the earth s sky , mars is having an apparent retrograde motion with respect to the usual planetary rotation . at the same time
, the earth declination shows some irregularities because of this retrograde motion .
the maxima in the earth declination do not correspond to the distance minima because the direction of the martian spin axis has also to be taken into account .
the maxima have a periodicity longer than the earth - mars synodic period ( 780 days ) . on the time scale of a spacecraft mission ( about 2 years )
, the mop can be safely neglected for the computation of the earth declination , the largest mop effect being the martian precession that changes the earth declination by less of @xmath4 .
the gray boxes on fig .
[ fig_dist ] for the earth declination correspond to time intervals where the earth declination is smaller than 10@xmath2 . in section 4
we will show that these intervals correspond to a lower sensitivity to mars spin axis motion in space .
another useful equatorial coordinates that will appear in the mop signatures ( section 4 ) is the earth hour angle @xmath5 , the angle along the celestial equator from the mars local meridian to the hour circle passing through earth , in the retrograde direction .
if @xmath6 is the earth right ascension , @xmath7 the sidereal mars rotation angle ( the angle from the vernal point to the prime meridian ) and @xmath8 the lander longitude ( positive east from the prime meridian ) , the earth hour angle is related to the other angles in the equatorial plane by the following equation @xmath9 the earth hour angle has a diurnal variation ( period of 24h 37min for mars ) , increasing every day non - linearly from 0 to 360@xmath2 , and depends on the lander position .
every day the earth rises and crosses the lander sky . like the earth declination @xmath3
, the earth elevation @xmath10 has variations related to the orbital motion of the earth and mars .
if @xmath11 is the lander latitude , the maximal elevation @xmath12 of the earth in the lander sky varies for each day between @xmath13 and the angle corresponding to the minimum between @xmath14 and @xmath15 .
the rise instrument have two antennas designed to observe the earth within a @xmath1 to @xmath16 elevation range , with one antenna pointed nearly to the east and the other pointed nearly to the west .
lara antennas nominal range shall be @xmath17 elevation , for all azimuthal directions ( @xmath18 ) , with a maximum gain direction at @xmath19 of elevation .
the rotation matrix from the mars body - fixed reference frame to the inertial frame based on mars mean orbit of j2000 is ( folkner et al . , 1997a , konopliv et al . , 2006 , le maistre et al . , 2012 )
@xmath20 if @xmath21 is the lander position vector in the body - fixed reference frame , the position of the lander in the inertial frame @xmath22 is : @xmath23 @xmath24 , @xmath25 and @xmath26 are the classical elementary rotation matrices around the x , y , z axes respectively . as for the earth ,
five angles called here the martian orientation parameters ( mop ) are used to characterize the matrix @xmath27 : @xmath28 , @xmath29 , @xmath7 , @xmath30 and @xmath31 .
@xmath30 and @xmath31 are the crust - fixed coordinates of mars spin axis related to the polar motion .
@xmath7 is the martian sidereal rotation angle .
its variations @xmath32 are called the length - of - day ( lod ) variations . @xmath33 and @xmath34 are the precession - nutation matrices depending on the obliquity @xmath28 between the instantaneous mars equator and the mean orbital plane at j2000 and on the longitude of the node of the mars equator @xmath29 .
nutation in obliquity and longitude @xmath35 and @xmath36 are the periodic variations of @xmath28 and @xmath29 while the precession is @xmath37 the secular trend of @xmath29 .
we can write @xmath38 @xmath39 is the diurnal rotation rate .
the mop are small angles .
+ nutation is the periodic motion of the spin pole in space .
the nutations for a rigid planet can be computed from celestial mechanics and ephemerides .
a rigid nutation model is usually constructed as a sum of different frequencies with an amplitude and a phase for the nutation in obliquity ( @xmath35 ) and in longitude ( @xmath36 ) . in the following sections , we use the term `` nutation in longitude '' @xmath36 to mean the sum of the components at different frequencies of the nutation in longitude , and similarly the term `` nutation in obliquity '' @xmath35 as the sum of the main frequencies .
for mars , the largest nutations are the annual and semiannual nutations .
their maximal amplitude are about 550 and 1130 mas , corresponding to 9 and 19 m on mars surface .
we use roosbeek ( 1999 ) rigid mars nutations model .
however mars is not a rigid planet .
therefore the nutation amplitudes depend on the interior structure and parameters , like the state ( liquid or solid ) of the core or its moment of inertia .
we can model the non - rigid nutation amplitudes by applying a transfer function to the rigid nutations ( dehant et al .
, 2000 , van hoolst et al . , 2000a , 2000b ) .
this transfer function is based on the model developed for the earth by sasao et al .
( 1980 ) and includes the liquid core effects and other deformation effects .
for example , if the core is liquid , then an additional free mode , the free core nutation ( fcn ) amplifies the nutation response if the fcn period is close to one of the nutation periods .
current interior models estimate the fcn period to be between about -240 and -285 days .
therefore the fcn can mostly enhance the small retrograde terannual ( 229 days ) or the large prograde semiannual nutation ( 343.5 days ) .
considering the 3 sigma interval rather than the 1 sigma interval will make the fcn range larger than @xmath40 days and the fcn frequency can be even closer to the terannual frequency , increasing the amplification . in this study
, we use a transfer function and amplifications computed from up - to - date interior models ( panning et al . , 2016 ) .
the mantle mineralogy and the temperature profile are unique for all the models and the planet is assumed to be in hydrostatic equilibrium .
the variable quantities in the interior models are the size of the core and the density and thickness of the crust .
the models agree with the observed moment of inertia at one sigma .
the frequencies of the nutations are not affected by the transfer function , therefore the modeling we use is a trigonometric expansion with known periods but unknown amplitudes .
the amplifications are different for each nutation period and for the nutation in obliquity and in longitude ( see for example figure 10 of le maistre et al . , 2012 ) .
depending on the size of the core , the maximal change of the nutations due to the liquid core is between @xmath41  mas and @xmath42  mas ( corresponding to about @xmath43  cm and @xmath44  cm on mars surface ) for the nutation in obliquity and between @xmath42  mas and @xmath44  mas or @xmath45  cm and @xmath46  cm for the nutation in longitude .
the lower values correspond to a small core of @xmath47  km while the upper values correspond to a large core of about @xmath48  km .
an important objective of the geodesy experiments is to measure the nutations and particularly the transfer function because it is good tool to constraint the interior structure of mars and the core dimension and density in particular .
+ the precession is the secular motion of the spin axis about the normal to the orbit plane .
the spin axis completes one rotation in a period of about 170 300 years .
the precession is due to the gravitational pull of the sun on the oblate figure and depends on the moment of inertia about the spin axis . the present uncertainty on the precession rate ( konopliv et al . ,
2016 ) is 2.1 mas / y , corresponding to a period change of about 47 years .
since the precession rate @xmath37 is proportional to @xmath49 where @xmath50 is mars moment of inertia , the estimation of the precession rate gives a constraint on the interior mass distribution .
+ the rotation rate variations or length - of - day ( lod ) variations @xmath32 are modeled by a trigonometric series with the annual , semiannual , terannual and quaterannual periodicities .
the annual nutation has an amplitude of about @xmath51  mas or @xmath52  m on mars surface while the semiannual nutation is about @xmath53  mas or @xmath54  m ( konopliv et al . , 2016 ) .
in the rotation matrix around the z - axis ( @xmath7 angle , see eq .
[ eq_phi ] ) , above the lod variations , there is an additional term depending on the nutation @xmath55 , adding a nutation signature in the doppler observable .
the lod variations are of interest for global climate models , to constrain the atmosphere and ice caps dynamics .
+ each polar motion ( pm ) coordinate @xmath30 and @xmath31 is modeled by a trigonometric series ( annual , semiannual , terannual ) and a chandler wobble term .
the amplitudes of the seasonal terms of the polar motion are expected to be between 0 and 15 mas ( defraigne et al .
, 2000 , van den acker et al . , 2002 ) .
the period of the chandler wobble for mars is expected to be around 200 or 220 days .
its amplitude has not yet been measured but is expected to be somewhere between 10 and 100 mas .
the polar motion has both climatological and geodetic information ; particularly the chandler wobble period depends mainly on the dynamical flattening of the planet and it provides information on the planet s elasticity and on inelastic behavior .
& max value & analytical expression & max value + rigid nutations in long .
@xmath36 & 29.4 m & @xmath56 & 12.3 m + rigid nutations in obliq . @xmath35 & 10.3 m & @xmath57 & 10.3 m + & 0.4/0.7 m & & 0.4 m + & 0.2/0.4 m & & 0.2 m + precession : & 0.8 m & @xmath58 & 0.8 m + ( if @xmath59 2.1 mas / y ) & & @xmath60^{1/2}$ ] & + lod @xmath61 & 11.8 m & @xmath62 & 11.8 m + polar motion @xmath63 ( @xmath64 ) & 1.2 m & @xmath65 & 1.2 m + ( if chandler amp . @xmath66 m ) & & @xmath67 & +    the maximal value of each mop on a long time interval is given in the second column of table [ tab1 ] , assuming the nominal model .
the amplitude of the chandler wobble is assumed to be around @xmath68 meters .
for the liquid core effect in the nutations , the first numbers ( 0.4 m for a large core and 0.2 m for a small core ) are the components of the nutation in obliquity , while the second ( 0.7 and 0.4 m ) are the components of the nutation in longitude .
the precession rate is very large , we consider here only the effect of the present uncertainty on the precession rate at j2000 ( about 2.1 mas / y ) .
the displacement due to this uncertainty after about 20 years after j2000 is about 0.8 m.    each mop changes of the lander position .
if the position change vector is the difference between the lander position with and without the mop and using the rotation matrix @xmath27 ( see equation [ eq_r ] ) , we compute the norm of this position change vector .
the analytical expressions and the maximal value of the norm of this position change vector is given in the last columns of table [ tab1 ] .
since we use first order approximation , all the displacements are proportional to the mop .
some of the lander displacements are modulated by the diurnal rotation .
we see that the lander displacement due to the lod variations is larger when the lander is closer to the equator .
the maximal value of the displacement due to the rigid nutation is independent of the lander position ( because the square roots in table 1 are always equal or smaller than @xmath69 ) , but close to the pole , the displacement due to the nutations does not have any diurnal variations .
there is a reduction of a factor 2.3 from the maximal value of the nutation in longitude @xmath36 to its change in the lander motion ( from @xmath70 m to @xmath71 m ) because of the obliquity of mars ( @xmath72 ) induces a projection in the lander displacement and because of the rotation angle component ( - @xmath73 ) that removes a part of the signal ( see equation [ eq_phi ] ) . for the other
mop , the maximal displacement due to the mop is the same as the mop maximal value .
depending on the size of the core , the maximal change of the lander position due to the liquid core is between 19 cm ( small core ) and 42 cm ( large core ) .
the analytical expression is a combination of the transfer function and the change of lander position due to the rigid nutations in obliquity and in longitude , its expression is lengthly and therefore has not been included in table 1 .
the nutations in longitude and in obliquity together changes the lander by about 12.6 m.
in this section we give the analytical expressions of the mop signature in the doppler lander - earth observable (= direct - to - earth or dte observable ) .
we consider an instantaneous doppler observable ( with an infinite speed of light and no integration time ) .
the doppler measurements are highly sensitive to the position variations along the line - of - sight and not sensitive to position variations in the plane perpendicular to the line - of - sight .
+ the doppler measurements between the earth and a spacecraft on mars varies between about -17 and 17 km / s .
this corresponds to the value of the relative velocity of the two planets around the sun .
its main modulations come from the relative orbital motions of the earth and mars .
other modulations come from the variation of the orientation and the rotation of mars and the earth .
the diurnal rotations of mars and of earth impress a sinusoidal modulation upon the doppler signal , with a beat frequency of about 36 days .
this period of 36 days corresponds to the mars - earth synodic period of rotation given that the rotation period of earth and mars are very close , 23h56min versus 24h37min ( le maistre 2013 ) .
the doppler observable is positive when the lander is going away from the earth and is negative when getting closer .
+ the doppler observable @xmath74 can be modeled as the projection of the velocity difference @xmath75 on the line - of - sight @xmath76 , with @xmath77 being the angle between @xmath75 and the earth - lander line - of - sight .
the signature of a parameter in the doppler observable is the observable estimated taking into account a parameter minus the observable without this parameter .
to the first order in the mop , it can be expressed as ( yseboodt et al . , 2003 ) :
@xmath78 where @xmath79 is the change in @xmath80 due to the mop and @xmath81 is the change in the line - of - sight direction due to the mop . the first term in eq .
( [ eq_signatu ] ) is the `` geometrical effect '' or change of the line - of - sight direction .
this term expresses the change in direction of the emitter - receiver line induced by the rotation parameter multiplied by the amplitude of the velocity . the second term of eq .
( [ eq_signatu ] ) is the contribution of the velocity change induced by the rotation .
both terms are proportional to mars radius and to the mop .
the line - of - sight displacement because of the mop @xmath81 decreases if the emitter - receiver distance increases , therefore the geometrical effect is proportional to the relative velocity over the distance @xmath82 . the velocity effect ,
i.e.  the change of the lander velocity @xmath79 , is independent of the emitter - receiver distance but is proportional to the diurnal rotation rate @xmath39 .
if the relative velocity over the distance is smaller than the planet diurnal rotation rate , then the velocity effect dominates the signature in the doppler observable .
this happens when the earth - planet distance is large . since this is the case for mars and because mars diurnal rotation rate @xmath39 is large , therefore this velocity effect is the one that is modeled in this section while the line - of - sight displacement is neglected .
if the emitter - receiver distance is small , like for example if the receiver is on board an orbiter , then the dominant effect in the doppler signal is the change in the line - of - sight .
+ the signatures of the mop parameters in the observable have already been investigated in le maistre et al .
( 2012 ) for the doppler observable , konopliv et al .
( 2006 ) for the lod signature in the doppler observable and folkner et al .
( 1997a , eq .
1 ) for the range observable . previously , curkendall and mcreynolds ( 1969 ) gave first order analytical expressions for an orbiter - earth doppler and range observables and its link with the equatorial coordinates of the spacecraft while yseboodt et al .
( 2003 ) studied the mop signatures in the orbiter - earth doppler link in the frame of a geophysical mission including a network of landers on mars surface .
kawano et al . ( 1999 ) studied the mop signatures in the same beam interferometry ( sbi ) observable for two landers ( or more ) on mars surface . +    , the nutations in longitude @xmath36 , the liquid core effect in the nutations ( large core , free core nutation period of 242 days ) , a precession rate @xmath37 difference of 2 mas / y , the lod variations and the polar motion . on the left part of the graph ,
the insight lander is in western elysium planitia ( @xmath83n , @xmath84e ) and on the right part , the exomars lander is in oxia planum ( @xmath85n , @xmath86e ) . in the nutation plots ,
the dotted gray line shows the earth declination @xmath3 .
the gray boxes correspond to time intervals where the earth declination is smaller than @xmath87 while the pink boxes show where the sep angle is smaller than @xmath1 .
the red and purple curves are the surface displacements due to the nutation in longitude and nutation in obliquity , respectively .
they are scaled .
, width=755,height=604 ]    the signatures of the different mop in the doppler observable are plotted on fig .  [ fig_sensi ] .
since we use first order approximation , the mop signatures are proportional to the mop . therefore the mop periodicities ( annual , semiannual , chandler period ... ) can also be seen in the signatures plot .
the expected numerical values are summarized in section [ secamp ] .
if @xmath35 is the nutation in obliquity and @xmath36 the nutation in longitude , their signature in the lander - earth doppler observable is @xmath88 where @xmath89 is mars radius , @xmath39 the diurnal rotation rate of mars , @xmath7 the diurnal rotation angle of mars , @xmath3 the earth declination with respect to mars , @xmath28 mars obliquity , @xmath5 the earth hour angle and @xmath6 the right ascension .
the mathematical derivation of these expressions is given in appendix [ ap1 ] .
+ since the signature is proportional to the nutation , the annual , semiannual and terannual modulations in the @xmath35 and @xmath36 nutations also appear in the signature plots ( see also the red and purple lines in fig .
[ fig_sensi ] ) .
the time derivatives of the nutations @xmath90 and @xmath91 give second order effects in the signature since their period ( mostly annual and semiannual ) is much larger than the diurnal rotation rate @xmath39 .
the nutation signature is proportional to the distance between the lander and the spin axis @xmath92 , therefore the nutation signature is maximal at the equator and null at the pole .
the nutation signature in the doppler observable must be a function of the direction of its spin axis with respect to the earth direction , eqs .
( [ eq_signeps ] ) and ( [ eq_signpsi ] ) show that it is proportional to the sinus of the earth declination .
a tilt of the rotation axis because of the nutation tilts the diurnal velocity vector and the resulting velocity difference is , to the first order , along the z ( polar ) axis in a body - fixed reference frame .
the projection of this velocity vector on the line - of - sight is null if the earth in the equatorial plane of mars ( @xmath93 ) .
therefore the doppler observable sensitivity to nutations is low when mars spin axis is nearly perpendicular to the mars - earth direction ( see the gray line and gray boxes in fig .
[ fig_sensi ] ) . replacing @xmath94 by @xmath3 leads to a maximal error of 3 % because the earth declination is always smaller than mars obliquity ( @xmath95 ) .
the factor @xmath94 in the nutation signature varies between @xmath96 and @xmath97 , therefore decreasing the signatures . in eqs .
( [ eq_signeps ] ) and ( [ eq_signpsi ] ) , there is a diurnal periodicity in the signature through mars rotation angle @xmath7 or through the local earth hour angle @xmath5 .
the link between the angles in the equatorial plane is given in eq .
[ eq_he ] .
this diurnal variation is shown on the plot by the blue regions between the maximal and minimal curve . during one day
, the signature of @xmath36 and @xmath35 can not be maximal at the same time , there is a time shift of 6 hours between the maxima of the two nutations .
the earth elevation in the lander s sky does not appear in the nutation signature .
however there is an indirect link between both quantities , because the earth elevation depends on the chosen timing of observation .
for example if observations are chosen every early morning , the earth elevation will be always be small , and from day to day the value of the earth hour angle will almost be the same .
the liquid core signature in the doppler observable is the difference between the doppler observable with and without the nutation transfer function , each part being the sum of the @xmath35 and @xmath36 contributions ( eqs .
[ eq_signeps ] and [ eq_signpsi ] ) .
the liquid core signature is much smaller than the rigid nutations signature because the liquid core amplification is also smaller .
the signature of the non - rigid part of the nutations has more modulations than the other mop signatures because sometimes there is an amplification of the nutation in longitude due to the liquid core while later , it is the nutation in obliquity that is amplified . therefore , the liquid core signature is a combination of both signatures .
if plotted on the same figure , the rigid and non - rigid nutation signatures will be very close and indistinguishable on fig .
the dependence on the earth declination ( see the gray line and gray boxes ) and the lander latitude is the same as the rigid nutation signature . over an interval of a few years , the position of the maxima and of the roots moves if the free core nutation ( fcn ) period changes ( see fig .  [ fig_liquid ] ) .
( scaled ) .
, width=755,height=340 ]    this is because if the fcn period is close to the terannual period ( @xmath98 days ) , the terannual nutation will be largely amplified while if the fcn period is close to the semiannual period ( @xmath99 days ) , the terannual nutation amplification will be small but the semiannual nutation will be larger .
+ again , the diurnal variation is shown on the plot by the blue area . therefore depending on the lander position , each day
there is a different observation time that maximizes the liquid core signature .
12 hours after the time when the liquid core signature is maximum , the liquid core signature is minimum .
6 hours later , the liquid core signature is null .
having a large negative signature gives a geophysical information as important as a large positive signature .
the time where the liquid core signature is maximal during the day depends on the fcn value .
the time difference is less than an hour for the fcn range investigated here .      using a method similar to appendix [ ap1 ]
, we find the expression of the precession rate @xmath100 signature in the doppler observable . if the product @xmath101 is assumed to be small , the precession signature is @xmath102 the signature increases with time but remains relatively small on an interval like the mission duration , which is short with respect to the precession period of about @xmath103 years .
the precession signature has a diurnal variation via the earth hour angle @xmath5 .
+ the first term in eq .
[ eq_signpsip ] is the dominant one close to the equator because @xmath104 .
it depends on @xmath105 while the second term depends on @xmath106 .
the maximal signature is almost constant with time ( see figure [ fig_sensi ] ) , there is a long period variation with a small amplitude due to @xmath107 in the first term of eq .
( [ eq_signpsip ] ) and due to @xmath108 in the second term .    the second term ( but not the first one ) appears in the signature of the nutation in longitude @xmath36 .
there is a large difference for the analytical expressions of @xmath36 and @xmath109 because the nutation in longitude appears twice in the rotation matrix @xmath27 ( in @xmath110 and in mars rotation angle rotation @xmath111 ) while the precession appears once in @xmath110 .
the signature of martian rotation angle variations @xmath32 ( or lod variations ) in the doppler observable is @xmath112 there is a proportionality between the lod variations and its signature , therefore the main frequencies are annual , semiannual and quaterannual .
similarly as for the nutation , the lod signature is proportional to the distance from the lander to the spin axis @xmath92 , therefore the signature is maximal at the equator and null at the pole .
the lander longitude disappears from the signature .
however there is an indirect effect because the earth hour angle @xmath5 depends on the longitude .
the signature has again a diurnal periodicity via the earth hour angle .
the lod signature is very little affected by the variations of the earth declination because the @xmath113 stays between @xmath114 and @xmath69 .
+ the diurnal velocity vector is in the equatorial plane of the body - fixed reference frame and is perpendicular to the lander position vector @xmath115 .
the change of the velocity vector due to the lod variations is also in the equatorial plane and to the first order perpendicular to the diurnal velocity , it is in the direction of the lander position vector @xmath116 .
the scalar product between the velocity difference and the line - of - sight ( los ) is maximized when the two vectors are in the same direction , therefore when the los is `` aligned '' with the mars - lander vector .
each day , the lod signature is maximized when the mars - earth direction crosses the local meridian of the lander ( when the local earth hour angle is null ) or equivalently when the earth culminates in the lander sky .
typically this happens near the local noon .
when the lander is on the equator ( @xmath117 ) , the maximal earth elevation is @xmath118 .
therefore the maximal lod signature in the doppler observable is proportional to @xmath119 and on a long time interval , the signature increases as the maximal earth elevation increases or equivalently as the earth declination decreases . when the lander is not on the equator , the maximal lod signature in the doppler observable is again when the local earth hour angle is null , corresponding to the maximal earth elevation during that day .
however the lod signature is not anymore proportional to @xmath119 . @xmath120
the signature in the doppler observable of the equatorial plane coordinates of the lander @xmath121 and @xmath122 in the body - fixed reference frame is @xmath123 there is always a diurnal periodicity . on a daily basis ,
the @xmath121 signature is maximized when @xmath124 or @xmath125 while the @xmath122 signature is maximized when @xmath126 or @xmath127 .
if we express this signature as a function of the lander latitude @xmath11 , the lander longitude disappears from the equation @xmath128    it is known that the doppler measurements are not very sensitive to the @xmath129-coordinate of the lander ( along the polar axis ) , and therefore solving for this coordinate is difficult ( le maistre , 2016 ) . a change in the @xmath129-coordinate only does not change the lander diurnal velocity . therefore its signature is very small because this is a second order contribution : the direction of the line - of - sight is changed a little bit because of the lander displacement in the polar direction .
there is also a small second order velocity change .
however if we want to know the third component of the lander position , we can impose that the lander has to stay on mars surface and there is a relationship between the three coordinates . this additional constraint related to
the topography model helps to solve for the @xmath129-coordinate ( le maistre , 2016 ) .    using the same equations
, we can also evaluate the signature of a tidal displacement . a lander displacement along the x and y axis due to the tides
has a signature in the doppler observable given by eq .
[ eq_xysign ] .
the signature of the polar motion @xmath30 and @xmath31 in the doppler observable is @xmath130 to the first order in the polar motion parameters , the first two coordinates of the position of the spin pole are @xmath131 , explaining the sign difference when comparing this equation to eq .
[ eq_xysign ] .
the diurnal periodicity in the signal comes from the earth hour angle @xmath5 . as for the lod variations ,
the signature is not very sensitive to the earth declination . during one day
, the maxima and minima of the @xmath30 and @xmath31 signatures happen at different time , separated by 6 hours .
the polar motion signature is a function of @xmath132 , therefore the signature increases with the lander latitude , being maximal at the poles .
the reason is that the shift of the lander position due to the polar motion is @xmath133 .
since the diurnal velocity vector is in the equatorial plane , because of the scalar product of the doppler observable , the part of the lander displacement that is important is also in the equatorial plane .
this lander displacement in the equatorial plane is proportional to @xmath134 .
equivalently to the sensitivity to the lander @xmath129-coordinate , the doppler observable is not sensitive to the @xmath129 displacement due to the polar motion .
the maximal value of the mop signature in the doppler observable are given using the insight and exomars nominal missions configuration in table [ tablsi ] .
the nominal mission time coverage is taken into account : 2 earth years for the insight lander starting end of 2019 , 1 earth year for the exomars lander starting beginning of 2021 .
the insight landing site is expected to be close to the equator in western elysium planitia ( @xmath83n , @xmath84e ) while we pick the oxia planum location ( @xmath85n , @xmath86e ) for the exomars lander .
the typical noise on doppler data is @xmath135  mm / s for a doppler counting interval of count interval of 60 seconds .
@xmath136    the temporal evolution of the mop signature in the doppler observable is plotted on fig .
[ fig_sensi ] . changing the landers position from one landing site to another position further from the equator
does not affect a lot the value of the signatures ( this is barely noticeable on fig .  2 ,
see the left and the right parts ) , except for the polar motion signature .    from tab .
[ tablsi ] , we see that the maximal signatures are smaller for exomars than for insight , except for the polar motion signature , mostly because the exomars nominal mission does not cover one full martian cycle and because the chosen exomars location is farther from the equator than the insight landing site .
the signature of the precession rate are very close because of a small decrease due to a higher latitude compensated by a small increase because the beginning of the exomars observation time is later with respect to j2000 .
[ fig_sensilat ] shows that the maximal value of these signatures depends on the lander latitude .
some rotation parameters like the nutation , the precession and the lod variations vary with @xmath137 with a larger signature at the equator while the polar motion varies with @xmath132 with a larger signature at the pole .     and
@xmath36 , the liquid core effect through nutations , the lod variations and the polar motion ( pm ) in the doppler observable as a function of the lander latitude .
two different liquid core sizes have been considered .
the two vertical lines show the landers latitude while the horizontal line shows the typical noise on doppler measurements ( @xmath135  mm / s).,width=502,height=302 ]     + the maximal signature of the nutations in the doppler signal is up to @xmath138  mm / s , coming from the signature of the nutations in longitude @xmath36 .
the nutation in obliquity signature is smaller ( up to @xmath139  mm / s ) , see fig .
[ fig_sensilat ] .
there is a difference of 6 hours between the maxima in the nutation in obliquity and in longitude , therefore they can not be maximal at the same time .
the maximum of the nutation signature is smaller than the product of the maximal value of the component of equations ( [ eq_signeps ] ) and ( [ eq_signpsi ] ) ( i.e.  @xmath140 and the maximum of @xmath141 or @xmath142 ) because the declination and the nutation do not have their maximum at the same time .
+ the liquid core signature in the doppler observable is much smaller than the nutation signature and largely depend on the characteristics of the core because of the resonant effects in the transfer function .
the maximal signature of the liquid core in the doppler observable is plotted as a function of the fcn period on fig .
[ fig_sensinut ] .
the signature is larger ( up to 0.01  mm / s ) when the fcn is close to the terannual frequency .
this can also be seen on fig .
[ fig_liquid ] .        a shift on the precession equal to the present uncertainty on the precession rate ( 2.1 mas / y , corresponding to a period shift of @xmath143 years , konopliv et al . ,
2016 ) changes the doppler observable by about @xmath144  mm / s after @xmath145 years . + whatever the lander latitude below 80@xmath2 , the signature of the lod is well above the other signatures .
it culminates at about @xmath146  mm / s for a lander located close to the equator .
this explains why this rotation parameter has already been measured using the previous martian landers .
+ assuming a chandler wobble amplitude of @xmath68 meters , the maximal signature of the polar motion for a lander latitude of @xmath147 is 0.024  mm / s and @xmath148  mm / s for a lander latitude of @xmath149 .
the signature increases as the latitude increase and an equatorial lander is not sensitive to the polar motion .
if the lander is close to the pole , the maximal signature is close to 0.1  mm / s .
+ a tidal displacement of 1 cm on mars surface has a very small signature in the doppler observable of about @xmath150  mm / s for an equatorial lander .
+ the factor of proportionality between a mop ( expressed in surface displacement on the planet s surface ) and its signature in the doppler observable is close to @xmath39 for the lod variations , the x and y - coordinates , the tidal displacement and the polar motion .
this factor is @xmath151 for the nutation in obliquity and @xmath152 for the nutation in longitude . for mars ,
the rotation rate is large ( @xmath153 rad / s ) .
for example if a @xmath69 meter shift is applied to the x or y - coordinates , the doppler observable changes by maximum 0.07  mm / s .
however if the z - coordinate of the lander changes by 10 meters , the signature in the doppler is negligible ( less than 0.001  mm / s ) .
+ the expressions given before are first order expressions linearized with respect to the mop .
with respect to the full expression , it stays close up to 95 % of the signal .
we also compared these analytical expressions to outputs from the gins software ( godsie par intgrations numriques simultanes ) , a numerical software able to process planetary geodesy data , and the numerical difference between the full expression and our analytical expression is small .
we can apply the expressions of the rotational variations signature in the lander - earth direct radio link to other bodies of the solar system and estimate the order of magnitude of the signature , see tab .
[ tab_ss ] .
the lander has to be equipped with a transponder is on the planetary surface and the communication is direct - to - earth ( dte ) .
the previous signature expressions can not be applied to doppler observable between a lander and an orbiter , or between the orbiter and the earth .
the largest rotational variation of bodies in spin - orbit resonance is usually the annual forced libration .
this libration amplitude depends on the interior structure and properties and if there is a liquid layer or not .
additional forced librations with other frequencies also exist .
the obliquity is another parameter important to observe , because if the body is in the cassini state , the equilibrium obliquity depends on the body moment of inertia .
.amplitude of the libration at the orbital period and its maximal signature in the doppler observable between the earth and a lander on the planetary body .
the lander is assumed to be on the equator which maximized the libration displacement .
when the libration amplitude is not yet measured , we use estimation coming from a theoretical model . [
cols="^,^,<,^",options="header " , ]     mercury is in a 3:2 spin orbit resonance and orbits around the sun with a period of 88 days .
the main libration is the annual libration , with an amplitude of 38.5 as ( margot et al . , 2012 ) .
using eq .
( [ eq_utsign ] ) , the maximal signature in the doppler observable of this libration is 0.56  mm / s .
the earth declination seen from mercury varies between -11@xmath2 and 11@xmath2 .
mercury has an obliquity of 2.04 arcmin ( margot et al . , 2012 ) .
using eq .
( [ eq_signeps ] ) , the maximal signature of the obliquity in the doppler signal is about 0.33  mm / s , smaller than the libration signature . for mercury , like for mars ,
the change of the planet velocity due to the rotation parameter is also the dominant term in the doppler observable signature .
however since the distance to earth is usually smaller for mercury than for mars , the geometrical effect for mercury is not as small as for mars .
+ since we evaluate only the change of the lander velocity because of the rotation parameter , the previous equations are valid if the lander is far from the receiving antenna on earth .
if the lander - receiver distance is smaller , the change of the line - of - sight direction must be taken into account and increase the doppler observable signature . because this line - of - sight change is not modeled here and
because the geometry of the earth - moon configuration is very different from the earth - planet configuration due to the synchronous rotation , we do not evaluate the signatures for a lunar lander .
+ if the lander in on phobos surface , the signature of the libration amplitude will be very large ( about 100  mm / s ) , because phobos rotation rate and its libration amplitude are both large .
+ for the icy moons around jupiter and saturn , the velocity effect is larger than the geometry effect .
each libration signature in the doppler observable is proportional to the moon radius , to the diurnal frequency and to the moon libration amplitude .
the variation of the earth declination does not largely affect the libration signature .
for enceladus , the estimated signature is large ( about 28  mm / s ) because the measured libration amplitude is large ( @xmath154 , thomas et al . , 2016 ) .
for the other moons , the libration amplitude has not yet been measured but we use order of magnitude from theoretical estimations . the signature in the doppler observable varies between 0.07 and 2.2
mm / s , see tab .
[ tab_ss ] .
these signatures are larger than the typical noise on doppler data ( @xmath135  mm / s ) . + a tidal displacement of 60 cm on europa surface has a signature in the lander - earth doppler observable of about @xmath155  mm / s for an equatorial lander .
, the nutations in longitude @xmath36 , the liquid core effect in the nutations for a large core , a precession rate @xmath37 difference of 2 mas / y , the lod variations and the polar motion ( @xmath30 and @xmath31 ) . on the left part of the graph ,
the insight lander is in western elysium planitia ( @xmath83n , @xmath84e ) and on the right part , the exomars lander is in oxia planum ( @xmath85n , @xmath86e ) . in the nutation plots ,
the gray line shows the earth declination @xmath3 .
the red and purple curves are the surface displacements due to the nutation in longitude and nutation in obliquity , respectively .
, width=755,height=680 ]    in this section , we give the expressions of the mop and lander position signature in the range observable .
the instantaneous range observable is equivalent to the distance between the lander and the antenna on earth at a given time .
the mop signature in the range observable @xmath156 is , to the first order in the rotation parameter , the vector change of position due to the mop @xmath157 projected on the earth - mars line - of - sight @xmath158 . @xmath159
this expression is independent from the earth - mars distance .
the variables that appear in the following expressions are the position of the earth with respect to mars ( i.e.  the earth declination @xmath3 , the right ascension @xmath6 and the hour angle @xmath5 ) and the lander position .    because the doppler measurement can be seen as the time derivative of the range observable ,
the parameter signature in the doppler observable is the time derivative of the parameter signature in the range observable .
therefore the global shapes of the signature of these two observables are very similar , compare figs .
[ fig_sensi ] and [ fig_sensiran ] .
the range signatures may have additional shift or long - term variations that disappear in the doppler observable signature . again
the mop signature in the range observable has a diurnal variation and is proportional to the mop .
the signature of the nutation in longitude @xmath36 and nutation in obliquity @xmath35 in the range observable is @xmath160 the first term of these expressions changes slowly and is null if the lander is close to the equator while the second term has a diurnal variation and is dominant .
close to the equator , the dependence of the range signature to the declination is the same as the doppler signature : when the earth declination is null , the nutation signature in the range observable is null too .
the liquid core signature in the range observable is the difference between the range observable with and without the nutation transfer function , each part being the sum of the contributions of the nutation in longitude and in obliquity ( eqs .
[ eq_signraneps ] and [ eq_signranpsi ] ) .
the signature of the precession in the range observable is @xmath161 the precession signature increases with time but on an interval like the mission duration which is short with respect to the precession period of about @xmath103 years , this increase is relatively small .
the diurnal variations dominates the signature .
the second and third terms are equivalent to the signature of the nutation in longitude @xmath36 .
the additional term here with respect to the @xmath36 signature in the range observable ( the first one ) depends on @xmath105 while the others depend on @xmath106 .
the largest term is the first one ( @xmath162 ) , then the third one .
the second one ( the constant one ) is negligible .
the signature of the lod variations in the range observable is @xmath163 when the earth culminates in the lander sky ( @xmath164 ) , the lod signature is null .
the timing of the range maxima are therefore 6 hours later than the doppler maxima .
the signature of the lander coordinates in the range observable is @xmath165 the signatures of the x and y - coordinates have a diurnal variation .
the signature of the z - coordinate in the range observable , contrary to the doppler signature , is not null and has no diurnal variation . because of the @xmath108 , it has long term variations
. therefore some range measurements when the earth declination is not null can be useful to constrain the z - coordinate of a lander .
the signature of the polar motion in the range observable is @xmath166 the first part of these expressions ( @xmath167 or @xmath168 ) is large if the lander is close to the equator ( up to 0.3 m ) while the second term has a diurnal variation and is much smaller for an equatorial lander .
the maximal values of the mop signatures in the range observable are given in tab .
[ tablsiran ] .
they can also be viewed on fig .
[ fig_sensiran ] .
@xmath169    as for the doppler signature , the parameter that has the largest signature in the range observable is the lod variations with a signature of about 11 meters .
then the nutation has a maximal signature of about 3.7 to 5.5 meters .
the nutation signature in the range observable is smaller than the nutation amplitude and the lander displacement given in table 1 because of the reduction proportional to the sinus of the earth declination .
the liquid core signature is much smaller , about 0.17 meter .
the polar motion signature is larger than the liquid core signature ( up to 0.7 meter ) . if the precession shift by 2.1 mas / y , the precession signature in the range observable is max 0.74  m around 2020 .      using eq .
[ eq_utransign ] and similarly as the doppler observable study , we can numerically evaluate the signature in the range observable of an libration at the orbital period for a lander on another planet or moon of the solar system .
the typical noise on range data between the earth and a lander on another body is about a few meters .
the maximal value of the rotation parameter signature is simply the maximal displacement on the body surface due to the libration .
if we use the order of magnitude of the libration amplitude given in tab .
[ tab_ss ] , the maximal signature in the range observable for an equatorial lander is given by the second column .
it varies from a few meters up to 528 meters for enceladus .
the signature for mercury is also large , 455 meters .
the maximal signature of mercury obliquity in the range observable is about 270 m , smaller than the libration signature .
all the signatures are modulated by the diurnal rotation and by variations at at long period due to the configuration geometry .
the same beam interferometry ( sbi ) technique , also known as the inverse vlbi , is based on the simultaneous tracking of two martian landers with one large antenna on earth , for example a deep space network antenna , and measuring the delay between the two radio signals .
the advantage of the sbi technique is that the two signals cross almost the same media at the same time , and therefore taking the difference cancels a large part of the common errors , like the earth atmospheric and ionospheric errors or the plasma delays .
the sbi observable corresponds to the difference of the two ranges between the antenna and each lander .
therefore the rotation parameter signature in the sbi observable is the difference between the rotation signature in the range observable of the two landers .
this technique has never been used up to now it is expected to have an error budget up to 0.7 cm for martian landers ( gregnanin et al . , 2014 ) .
the mop signature in the sbi observable depends on the latitude and longitude of each lander .
we can replace the local angle @xmath5 by @xmath170 in the previous expressions and substitute the latitude and longitude of the lander by their rectangular coordinate @xmath171 .
for example , the lod signature in the sbi observable is @xmath172 where @xmath173 and @xmath174 are the difference of x and y - coordinates of the two landers .
the @xmath175 contribution does not appear in this signature .
the nutation signature in the sbi observable is @xmath176 we see that the parts proportional to @xmath173 and @xmath174 are affected by a diurnal modulation and have a negligible contribution if the earth declination is small .
the @xmath175 contribution has long term variations .
these expressions are similar to the sbi sensitivity given in kawano et al .
( 1999 , eq .  28 ) .
their variables describing the earth position with respect to mars are not equivalent to ours ( their model is simpler ) but it is possible to see the same modulations ( the diurnal and the long term ) by comparing the two methods .
the numerical values of the mop signatures in the sbi observable depend of course of the separation vector between the two landers . but
a general conclusion is that the signature is larger if the distance between the two landers is larger .
using the analytical expressions of the mars orientation parameters ( mop ) signatures in the doppler observable , we can anticipate the moments where the mop signatures are maximal , on a daily basis and on longer interval and we can predict if a mop can be observed or not . the mop signature in the doppler and the range observable is proportional to the mop signal .
therefore , the periodicities of the mop ( mostly annual , semiannual , terannual , the chandler period of about 200 days for the polar motion etc . )
can also be directly seen in the signatures . on a long time interval , for the mission operations planning ,
observation times should be privileged when the mop is large in order to get a larger signature . + additionally , some signatures like the nutations depend on the earth declination .
the signature is larger if the earth is outside the equatorial plane of mars .
the best time to observe the nutations is when the earth declination is large and when the sun - earth - probe angle is large because the noise on the doppler measurement due to the plasma will be smaller at that time .
this corresponds to the beginning ( 2018.9 - 2019.25 ) and end ( 2020.3 - 2021 ) of the insight mission .
for the exomars mission , this corresponds to approximately 2021 - 2021.1 and 2021.5 - 2021.8 .
+ in every signature , there is also a diurnal modulation because the main motion affecting the lander is the diurnal rotation of mars .
therefore every day , each mop signature has a maximum and a minimum signature .
these observation time should be favored .
however , the daily maxima in the signatures of the nutations in obliquity and in longitude are separated by 6 hours , therefore it is impossible to maximize both at the same time .
the daily maxima in the length - of - day variations signature in the doppler observable happen when the earth is at its highest point in the lander sky .
+ defining the best timing windows is useful for the mission planning of the geodesy experiments in the future insight and exomars missions .
we will use these results when testing different strategies for adjustments of the rotational parameters .
these strategies will be tested by performing numerical simulations for assessing the precision on the determination of the mop determination .
we thank a. rivoldini for the computation of the transfer functions and rose - marie baland and sebastien le maistre for the helpful comments .
this work was financially supported by the belgian prodex program managed by the european space agency in collaboration with the belgian federal science policy office .
borderies , n. , g. balmino , l. castel , and b. moynot , study of mars dynamics from lander tracking data analysis , _ the moon and the planets , 22 _ , 191 - 200 , 1980 .
curkendall , d. w. , and s. r. mcreynolds , a simplified approach for determining the information content of radio tracking data , _ j. spacecraft , 6_(5 ) , 520 - 525 , 1969 .
defraigne , p. , o. de viron , v. dehant , t. van hoolst , and f. hourdin , mars rotation variations induced by atmosphere and ice caps , _ j. geophys .
res . , 105_(e10 ) , 24563 - 24570 , 2000 .
dehant , v. , t. van hoolst , and p. defraigne , comparison between the nutations of the planet mars and the nutations of the earth , _ survey geophys .
, 21_(1 ) , 89 - 110 , 2000 .
dehant , v. and 27 co - authors , lander radioscience for obtaining the rotation and orientation of mars , _ planetary and space science _ 57,10501067 , 2009 .
dehant , v. , maistre , s.l . ,
rivoldini , a. , yseboodt , m. , rosenblatt , p. , hoolst , t.v . ,
mitrovic , m. , karatekin , o. , marty , j. , chicarro , a. , revealing mars deep interior : future geodesy missions using radiolinks between landers , orbiters , and the earth .
_ planetary and space science _ 59 , 10691081 , 2011 .
folkner , w. , r. kahn , r. preston , c. yoder , e. standish , j. williams , c. edwards , r. hellings , t. eubanks , and b. bills , mars dynamics from earth - based tracking of the mars pathfinder lander , _ j. geophys .
res . , 102_(e2 ) , 4057 - 4064 , 1997a .
folkner , w. , c. yoder , d. yuan , e. standish , and r. preston , interior structure and seasonal mass redistribution of mars from radio tracking of mars pathfinder , _ sci . , 278 _ , 1749 - 1751 , 1997b .
folkner , w.m . ,
asmar , s.w . ,
dehant , v. , warwick , r.w . , the rotation and interior structure experiment ( rise ) for the insight mission to mars . in : lunar and planetary institute
science conference abstracts , lunar and planetary institute technical report , vol .
43 , march , p. 1721
, 2012 .
gregnanin , m. , m. yseboodt , v. dehant , l. iess and t. van hoolst , estimation of mars geophysical information through same beam interferometry , european planetary science congress ( epsc ) , vol .
9 , epsc2014 - 395 , 2014 .
kawano , n. , m. hosokawa , h. hanada , and m. imae , inverse vlbi method for planetodesy , _ j. geod .
_ , 45(3 ) , 181203 , 1999 .
konopliv , a. , yoder , c. , standish , e.m . , yuan dah - ning , sjogren , w.l . , a global solution for the mars static and seasonal gravity , mars orientation , phobos and deimos masses , and mars ephemeris .
_ icarus _ 182 ( 1),2350 , 2006 .
konopliv , a.s .
park , w.m .
folkner , an improved jpl mars gravity field and orientation from mars orbiter and lander tracking data , _ icarus _ , 274 , 2016 .
kuchynka , p. , folkner , w.m . ,
konopliv , a.s . , park , r.s .
, le maistre s. , dehant , v. , new constraints on mars rotation determined from radiometric tracking of the opportunity mars exploration rover .
_ icarus _ 229 , 340347 , 2014 .
le maistre , s. , rosenblatt , p. , rivoldini , a. , dehant , v. , marty , j - c .
, karatekin , . , lander radio science experiment with a direct link between mars and the earth . _
space sci .
_ 68 , 2012 .
le maistre , s. , the rotation of mars and phobos from earth - based radio - tracking observations of a lander , phd dissertation , universit catholique de louvain , belgium , available on dial.uclouvain.be , 2013 .
le maistre , s. , insight coordinates determination from direct - to - earth radio - tracking and mars topography model , _ planet .
space sci .
_ 121 , 2016 .
margot , j .-
, s. j. peale , s. c. solomon , s. a. hauck ii , f. d. ghigo , r. f. jurgens , m. yseboodt , j. d. giorgini , s. padovan , and d. b. campbell , mercury s moment of inertia from spin and gravity data , _ j. geophys .
_ , 117 , e00l09 , doi:10.1029/2012je004161 , 2012 .    panning , m.p . , and 33 co - authors , planned products of the mars structure service for the insight mission to mars , submitted to space science reviews , 2016 .
roosbeek , f. , analytical developments of rigid mars nutation and tide generating potential series , _ celestial mechanics , 75 _ , 287 - 300 , 1999 .
sasao , t. , s. okubo , and m. saito , a simple theory on dynamical effects of a stratified fluid core upon nutational motion of the earth , proceedings of the agu symposium , no .
78 on _ nutation and the earth s rotation _ , fedorov ep et al .
( eds ) , kiew , d. reidel , norwell , mass . , 165 - 183 , 1980 .
thomas , p.c .
, r. tajeddine , m.s .
tiscareno , j.a .
burns , j. joseph , t.j .
loredo , p. helfenstein , c. porco , enceladus s measured physical libration requires a global subsurface ocean , _ icarus , 264 _ , 2016 .
van den acker , e. , t. van hoolst , o. de viron , p. defraigne , f. forget , f. hourdin , and v. dehant , influence of the winds and of the co@xmath0 mass exchange between the atmosphere and the polar ice caps on mars rotation ,
_ j. geophys .
, 107_(7 ) , doi:10.1029/2000je001539 , 2002 .
van hoolst , t. , v. dehant , and p. defraigne , chandler wobble and free core nutation for mars , _ planet .
space sci . , 48_(12 - 14 ) , 1145 - 1151 , 2000a .
van hoolst , t. , v. dehant , and p. defraigne , sensitivity of the free core nutation and the chandler wobble to changes in the interior structure of mars , _ phys .
earth planet .
inter . , 117 _ , 397 - 405 , 2000b .
yoder , c. f. , and e. m. standish , martian precession and rotation from viking lander range data , _ j. geophys .
res . , 102_(e2 ) , 4065 - 4080 , 1997 .
yseboodt , m. , j .-
barriot , and v. dehant , analytical modeling of the doppler tracking between a lander and a mars orbiter in terms of rotational dynamics , _ j. geophys .
res . , 108 _ , 2003 .
assuming there is no lod variations and no polar motion , we describe here the method to get the expression of the nutation signature in the doppler observable .
however this development can be applied similarly to each mop .
+ if the vector @xmath177 is the fixed lander position in the body fixed reference frame , the lander position in the inertial frame is @xmath178 .
its velocity in the inertial frame is : @xmath179 the terms proportional to @xmath100 and @xmath180 have been neglected here because they are much smaller than @xmath181 ( @xmath182 the diurnal rotation @xmath39 ) .    to the first order in the nutations in obliquity @xmath183 , the lander velocity in a non - rotating frame based on the mars
mean equator of j2000 ( mme2000 ) is @xmath184 the lander velocity @xmath185 is mostly a velocity in the xy plane .
but the lander velocity induced by a nutation in obliquity , in the non - rotating frame based on the mars mean equator of 2000 , is mainly along the z axis .
+ similarly , the main component of the signature of the nutations in longitude @xmath36 in the lander velocity is also the z component , because the dependence of the velocity is only in the z coordinate . @xmath186
this z component ( to the first order in @xmath36 ) is not as straightforward as for @xmath35 because there are @xmath36 dependent terms both in the products of 4 rotation matrices and in the @xmath187 matrix , and there is cancellation between some of these terms .
+ in the same inertial frame , the tracking station - lander direction is aproximated by the earth - mars direction .
@xmath6 is the earth right ascension .
the mop signature in the doppler observable is the scalar product between the velocity difference @xmath189 = ( @xmath190 or @xmath191 ) and the los vector : @xmath192