Tidal effects in planetary systems are the main driver in the orbital migration of natural satellites. They result from physical processes occurring deep inside celestial bodies, whose effects are rarely observable from surface imaging. For giant planet systems, the tidal migration rate is determined by poorly understood dissipative processes in the planet, and standard theories suggest an orbital expansion rate inversely proportional to the power 11/2 in distance, implying little migration for outer moons such as Saturn's largest moon, Titan. Here, we use two independent measurements obtained with the Cassini spacecraft to measure Titan's orbital expansion rate. We find Titan migrates away from Saturn at 11.3 $\pm$ 2.0 cm/year, corresponding to a tidal quality factor of Saturn of Q $\simeq$ 100, and a migration timescale of roughly 10 Gyr. This rapid orbital expansion suggests Titan formed significantly closer to Saturn and has migrated outward to its current position. Our results for Titan and five other moons agree with the predictions of a resonance locking tidal theory, sustained by excitation of inertial waves inside the planet. The associated tidal expansion is only weakly sensitive to orbital distance, motivating a revision of the evolutionary history of Saturn's moon system. The resonance locking mechanism could operate in other systems such as stellar binaries and exoplanet systems, and it may allow for tidal dissipation to occur at larger orbital separations than previously believed.

Context. Bright stellar positions are now known with an uncertainty below 1 mas thanks to Gaia DR2. Between 2019-2020, the Galactic plane will be the background of Jupiter. The dense stellar background will lead to an increase in the number of occultations, while the Gaia DR2 catalogue will reduce the prediction uncertainties for the shadow path. Aims. We observed a stellar occultation by the Galilean moon Europa (J2) and propose a campaign for observing stellar occultations for all Galilean moons. Methods. During a predicted period of time, we measured the light flux of the occulted star and the object to determine the time when the flux dropped with respect to one or more reference stars, and the time that it rose again for each observational station. The chords obtained from these observations allowed us to determine apparent sizes, oblatness, and positions with kilometre accuracy. Results. We present results obtained from the first stellar occultation by the Galilean moon Europa observed on 2017 March 31. The apparent fitted ellipse presents an equivalent radius of 1561.2 $\pm$ 3.6 km and oblatenesses 0.0010 $\pm$ 0.0028. A very precise Europa position was determined with an uncertainty of 0.8 mas. We also present prospects for a campaign to observe the future events that will occur between 2019 and 2021 for all Galilean moons. Conclusions. Stellar occultation is a suitable technique for obtaining physical parameters and highly accurate positions of bright satellites close to their primary. A number of successful events can render the 3D shapes of the Galilean moons with high accuracy. We encourage the observational community (amateurs included) to observe the future predicted events.

During the thirteen years in orbit around Saturn before its final plunge, the Cassini spacecraft provided more than ten thousand astrometric measurements. Such large amounts of accurate data enable the search for extremely faint signals in the orbital motion of the saturnian moons. Among these, the detection of the dynamical feedback of the rotation of the inner moons of Saturn on their respective orbits becomes possible. Using all the currently available astrometric data associated with Atlas, Prometheus, Pandora, Janus and Epimetheus, we first provide a detailed analysis of the Cassini Imaging Science Subsystem (ISS) data, with special emphasis on their statistical behavior and sources of bias. Then, we give updated estimates of the moons' averaged densities and try to infer more details about their interior properties by estimating the physical librations for Prometheus, Pandora, Epimetheus and Janus from anomalies in their apsidal precession. Our results are compatible with a homogeneous interior for Janus and Epimetheus, within the uncertainty of the measurements. On the other hand, we found some inconsistency for Pandora and Prometheus, which might result from a dynamical mismodeling of Saturn's gravity field. Last, we show how the synergistic introduction of libration measurements directly derived from imaging should allow the moons' moments of inertia to be better constrained.

The technique of mutual approximations accurately gives the central instant at the maximum apparent approximation of two moving natural satellites in the sky plane. This can be used in ephemeris fitting to infer the relative positions between satellites with high precision. Only the mutual phenomena -- occultations and eclipses -- may achieve better results. However, mutual phenomena only occur every six years in the case of Jupiter. Mutual approximations do not have this restriction and can be observed at any time along the year as long as the satellites are visible. In this work, we present 104 central instants determined from the observations of 66 mutual approximations between the Galilean moons carried out at different sites in Brazil and France during the period 2016--2018. For 28 events we have at least two independent observations. All telescopes were equipped with a narrow-band filter centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The telescope apertures ranged between 25--120 cm. For comparison, the precision of the positions obtained with classical CCD astrometry is about 100 mas, for mutual phenomena it can achieve 10 mas or less and the average internal precision obtained with mutual approximations was 11.3 mas. This new kind of simple, yet accurate observations can significantly improve the orbits and ephemeris of Galilean satellites and thus be very useful for the planning of future space missions aiming at the Jovian system.

Natural satellite systems present a large variety of orbital configurations in the solar system. While some are clearly the result of known processes, others still have largely unexplained eccentricity and inclination values. Iapetus has a still unexplained 3% orbital eccentricity and its orbital plane is tilted with respect to its local Laplace plane. On the other hand, astrometric measurements of saturnian moons have revealed high tidal migration rates, corresponding to a quality factor Q of Saturn of around 1600 for the mid-sized icy moons. We show how a past crossing of the 5:1 mean motion resonance between Titan and Iapetus may be a plausible scenario to explain Iapetus' orbit. We have carried out numerical simulations of the resonance crossing using an N-Body code as well as using averaged equations of motion. A large span of migration rates were explored for Titan and Iapetus was started on its local Laplace plane with a circular orbit. The resonance crossing can trigger a chaotic evolution of the eccentricity and the inclination of Iapetus. The outcome of the resonance is highly dependent on the migration rate. For a quality factor Q of over around 2000, the chaotic evolution of Iapetus in the resonance leads in most cases to its ejection, while simulations with a quality factor between 100 and 2000 show a departure from the resonance with post-resonant eccentricities spanning from 0 up to 15%, and free inclinations capable of reaching 11 degrees. Usually high inclinations come with high eccentricities but some simulations show elements compatible with Iapetus' current orbit. A quality factor between 100 and 2000 at the frequency of Titan would bring Titan and Iapetus into a 5:1 resonance, which would perturb Iapetus' eccentricity and inclination to values observed today. Such rapid tidal migration would have avoided Iapetus' ejection around 40 to 800 million years ago.

The JUICE mission will perform detailed measurements of the properties of the Galilean moons, with a nominal in-system science-mission duration of about 3.5 years. Using both the radio tracking data, and (Earth- and JUICE-based) optical astrometry, the dynamics of the Galilean moons will be measured to unprecedented accuracy. This will provide crucial input to the determination of the ephemerides and physical properties of the system, most notably the dissipation in Io and Jupiter. The data from Planetary Radio Interferometry and Doppler Experiment (PRIDE) will provide the lateral position of the spacecraft in the ICRF. In this article, we analyze the relative quantitative influence of the JUICE-PRIDE observables to the determination of the ephemerides of the Jovian system and the associated physical parameters. We perform a covariance analysis for a broad range of mission and system characteristics. We analyze the influence of VLBI data quality, observation planning, as well as the influence of JUICE orbit determination quality. This provides key input for the further development of the PRIDE observational planning and ground segment development. Our analysis indicates that the VLBI data are especially important for constraining the dynamics of Ganymede and Callisto perpendicular to their orbital planes. Also, the use of the VLBI data makes the uncertainty in the ephemerides less dependent on the error in the orbit determination of the JUICE spacecraft itself. Furthermore, we find that optical astrometry data of especially Io using the JANUS instrument will be crucial for stabilizing the solution of the normal equations. Knowledge of the dissipation in the Jupiter system cannot be improved using satellite dynamics obtained from JUICE data alone, the uncertainty in Io's dissipation obtained from our simulations is similar to the present level of uncertainty.

Accurate tracking of probes is one of the key points of space exploration. Range and Doppler techniques are the most commonly used. In this paper we analyze the impact of the transponder delay, $i.e.$ the processing time between reception and re-emission of a two-way tracking link at the satellite, on tracking observables and on spacecraft orbits. We show that this term, only partially accounted for in the standard formulation of computed space observables, can actually be relevant for future missions with high nominal tracking accuracies or for the re-processing of old missions. We present several applications of our formulation to Earth flybys, the NASA GRAIL and the ESA BepiColombo missions.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 {\mu}m/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m.

Tidal dissipation is the main driver of orbital evolution of natural satellites and a key point to understand the exoplanetary system configurations. Despite its importance, its quantification from observations still remains difficult for most objects of our own Solar system. In this work, we overview the method that has been used to determine, directly from observations, the tidal parameters, with emphasis on the Love number k2 and the tidal quality factor Q. Up-to-date values of these tidal parameters are summarized. Last, an assessment on the possible determination of the tidal ratio k2/Q of Uranus and Neptune is done. This may be particularly relevant for coming astrometric campaigns and future space missions focused on these systems.

Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.

Accurate positional measurements of planets and satellites are used to improve our knowledge of their dynamics and to infer the accuracy of planet and satellite ephemerides. In the framework of the FP7 ESPaCE project, we provide the positions of Mars, Phobos, and Deimos taken with the U.S. Naval Observatory 26-inch refractor during the 1971 opposition of the planet. These plates were measured with the digitizer of the Royal Observatory of Belgium and reduced through an optimal process that includes image, instrumental, and spherical corrections to provide the most accurate data. We compared the observed positions of the planet Mars and its satellites with the theoretical positions from INPOP10 and DE430 planetary ephemerides, and from NOE and MAR097 satellite ephemerides. The rms residuals in RA and Dec. of one position is less than 60 mas, or about 20 km at Mars. This accuracy is comparable to the most recent CCD observations. Moreover, it shows that astrometric data derived from photographic plates can compete with those of old spacecraft (Mariner 9, Viking 1 and 2).

The irregular satellites of the giant planets are believed to have been captured during the evolution of the solar system. Knowing their physical parameters, such as size, density, and albedo is important for constraining where they came from and how they were captured. The best way to obtain these parameters are observations in situ by spacecrafts or from stellar occultations by the objects. Both techniques demand that the orbits are well known. We aimed to obtain good astrometric positions of irregular satellites to improve their orbits and ephemeris. We identified and reduced observations of several irregular satellites from three databases containing more than 8000 images obtained between 1992 and 2014 at three sites (Observat\'orio do Pico dos Dias, Observatoire de Haute-Provence, and European Southern Observatory - La Silla). We used the software PRAIA (Platform for Reduction of Astronomical Images Automatically) to make the astrometric reduction of the CCD frames. The UCAC4 catalog represented the International Celestial Reference System in the reductions. Identification of the satellites in the frames was done through their ephemerides as determined from the SPICE/NAIF kernels. Some procedures were followed to overcome missing or incomplete information (coordinates, date), mostly for the older images. We managed to obtain more than 6000 positions for 18 irregular satellites: 12 of Jupiter, 4 of Saturn, 1 of Uranus (Sycorax), and 1 of Neptune (Nereid). For some satellites the number of obtained positions is more than 50\% of what was used in earlier orbital numerical integrations. Comparison of our positions with recent JPL ephemeris suggests there are systematic errors in the orbits for some of the irregular satellites. The most evident case was an error in the inclination of Carme.

Tidal dissipation, which is directly linked to internal structure, is one of the key physical mechanisms that drive systems evolution and govern their architecture. A robust evaluation of its amplitude is thus needed to predict evolution time for spins and orbits and their final states. The purpose of this paper is to refine recent model of the anelastic tidal dissipation in the central dense region of giant planets, commonly assumed to retain a large amount of heavy elements, which constitute an important source of dissipation. The previous paper evaluated the impact of the presence of the static fluid envelope on the tidal deformation of the core and on the associated anelastic tidal dissipation, through the tidal quality factor Qc. We examine here its impact on the corresponding effective anelastic tidal dissipation, through the effective tidal quality factor Qp. We show that the strength of this mechanism mainly depends on mass concentration. In the case of Jupiter- and Saturn-like planets, it can increase their effective tidal dissipation by, around, a factor 2.4 and 2 respectively. In particular, the range of the rheologies compatible with the observations is enlarged compared to the results issued from previous formulations. We derive here an improved expression of the tidal effective factor Qp in terms of the tidal dissipation factor of the core Qc, without assuming the commonly used assumptions. When applied to giant planets, the formulation obtained here allows a better match between the an elastic core's tidal dissipation of a two-layer model and the observations.

We present astrometric observations of the Saturnian satellites Mimas, Enceladus, Tethys, Dione and Rhea from Cassini Imaging Science Subsystem (ISS) narrow-angle camera (NAC) images. Image sequences were designed to observe mutual occultations between these satellites. The positions of satellite centres were estimated by fitting ellipsoidal shape models to the measured limbs of the imaged satellites. Spacecraft pointing corrections were computed using the UCAC2 star catalogue. We provide a total of 2303 astrometric observations, resulting in 976 pairs, the remainder consisting of observations of a single satellite. Mean residuals for the individual satellite positions relative to the SAT360 ephemeris were 4.3 km in the line direction and -2.4 km in the sample direction, with standard deviations of 5.6 and 7.0 km respectively, an order of magnitude improvement in precision compared to published HST observations. By considering inter-satellite separations, uncertainties in camera pointing and spacecraft positioning along with possible biases in the individual positions of the satellites can be largely eliminated, resulting in an order-of-magnitude increase in accuracy compared to that achievable using the individual satellite positions. We show how factors relating to the viewing geometry cause small biases in the individual positions of order 0.28 pixel to become systematic across the dataset as a whole and discuss options for reducing their effects . The reduced astrometric data are provided in the form of individual positions for each satellite, together with the measured positions of reference stars, in order to allow more flexibility in the processing of the observations, taking into account possible future advances in limb-fitting techniques as well as the future availability of more accurate star catalogues, such as those from the GAIA mission.

During the last decade, the precision in the tracking of spacecraft has constantly improved. The discovery of few astrometric anomalies, such as the Pioneer and Earth flyby anomalies, stimulated further analysis of the operative modeling currently adopted in Deep Space Navigation (DSN). Our study shows that some traditional approximations lead to neglect tiny terms that could have consequences in the orbit determination of a probe in specific configurations such as during an Earth flyby. Therefore, we suggest here a way to improve the light time calculation used for probe tracking.

In this paper, we focus on the possibility to test General Relativity in the Solar System with radioscience measurements. To this aim, we present a new software that simulates Range and Doppler signals directly from the space-time metric. This flexible approach allows one to perform simulations in General Relativity and in alternative metric theories of gravity. In a second step, a least-squares fit of the different initial conditions involved in the situation is performed in order to compare anomalous signals produced by a given alternative theory with the ones obtained in General Relativity. This software provides orders of magnitude and signatures stemming from hypothetical alternative theories of gravity on radioscience signals. As an application, we present some simulations done for the Cassini mission in Post-Einsteinian Gravity and in the context of MOND External Field Effect. We deduce constraints on the Post-Einsteinian parameters but find that the considered arc of the Cassini mission is not useful to constrain the MOND External Field Effect.

Earth-like planets have anelastic mantles, whereas giant planets may have anelastic cores. As for the fluid parts of a body, the tidal dissipation of such solid regions, gravitationally perturbed by a companion body, highly depends on its internal friction, and thus on its internal structure. Therefore, modelling this kind of interaction presents a high interest to provide constraints on planet interiors, whose properties are still quite uncertain. Here, we examine the equilibrium tide in the solid central region of a planet, taking into account the presence of a fluid envelope. We first present the equations governing the problem, and show how to obtain the different Love numbers that describe its deformation. We discuss how the quality factor Q depends on the rheological parameters, and the size of the core. Taking plausible values for the anelastic parameters, and examinig the frequency-dependence of the solid dissipation, we show how this mechanism may compete with the dissipation in fluid layers, when applied to Jupiter- and Saturn-like planets. We also discuss the case of the icy giants Uranus and Neptune.

GAIA will provide observations of some multiple asteroid and dwarf systems. These observations are a way to determine and improve the quantification of dynamical parameters, such as the masses and the gravity fields, in these multiple systems. Here we investigate this problem in the cases of Pluto's and Eugenia's system. We simulate observations reproducing an approximate planning of the GAIA observations for both systems, as well as the New Horizons observations of Pluto. We have developed a numerical model reproducing the specific behavior of multiple asteroid system around the Sun and fit it to the simulated observations using least-square method, giving the uncertainties on the fitted parameters. We found that GAIA will improve significantly the precision of Pluto's and Charon's mass, as well as Petit Prince's orbital elements and Eugenia's polar oblateness.

Earth-like planets have viscoelastic mantles, whereas giant planets may have viscoelastic cores. The tidal dissipation of such solid regions, gravitationally perturbed by a companion body, highly depends on their rheology and on the tidal frequency. Therefore, modelling tidal interactions presents a high interest to provide constraints on planets' properties and to understand their history and their evolution, in our Solar System or in exoplanetary systems. We examine the equilibrium tide in the anelastic parts of a planet whatever the rheology, taking into account the presence of a fluid envelope of constant density. We show how to obtain the different Love numbers that describe its tidal deformation. Thus, we discuss how the tidal dissipation in solid parts depends on the planet's internal structure and rheology. Finally, we show how the results may be implemented to describe the dynamical evolution of planetary systems. The first manifestation of the tide is to distort the shape of the planet adiabatically along the line of centers. Then, the response potential of the body to the tidal potential defines the complex Love numbers whose real part corresponds to the purely adiabatic elastic deformation, while its imaginary part accounts for dissipation. This dissipation is responsible for the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which we derive the characteristic evolution times. The rate at which the system evolves depends on the physical properties of tidal dissipation, and specifically on how the shear modulus varies with tidal frequency, on the radius and also the rheological properties of the solid core. The quantification of the tidal dissipation in solid cores of giant planets reveals a possible high dissipation which may compete with dissipation in fluid layers.

Tidal interactions between Saturn and its satellites play a crucial role in both the orbital migration of the satellites and the heating of their interiors. Therefore constraining the tidal dissipation of Saturn (here the ratio k2/Q) opens the door to the past evolution of the whole system. If Saturn's tidal ratio can be determined at different frequencies, it may also be possible to constrain the giant planet's interior structure, which is still uncertain. Here, we try to determine Saturn's tidal ratio through its current effect on the orbits of the main moons, using astrometric data spanning more than a century. We find an intense tidal dissipation (k2/Q= (2.3 \pm 0.7) \times 10-4), which is about ten times higher than the usual value estimated from theoretical arguments. As a consequence, eccentricity equilibrium for Enceladus can now account for the huge heat emitted from Enceladus' south pole. Moreover, the measured k2/Q is found to be poorly sensitive to the tidal frequency, on the short frequency interval considered. This suggests that Saturn's dissipation may not be controlled by turbulent friction in the fluid envelope as commonly believed. If correct, the large tidal expansion of the moon orbits due to this strong Saturnian dissipation would be inconsistent with the moon formations 4.5 Byr ago above the synchronous orbit in the Saturnian subnebulae. But it would be compatible with a new model of satellite formation in which the Saturnian satellites formed possibly over longer time scale at the outer edge of the main rings. In an attempt to take into account for possible significant torques exerted by the rings on Mimas, we fitted a constant rate da/dt on Mimas semi-major axis, also. We obtained an unexpected large acceleration related to a negative value of da/dt= -(15.7 \pm 4.4) \times 10-15 au/day.

The laws of gravitation have been tested for a long time with steadily improving precision, leading at some moment of time to paradigmatic evolutions. Pursuing this continual effort is of great importance for science. In this communication, we focus on Solar System tests of gravity and more precisely on possible tests that can be performed with radio science observations (Range and Doppler). After briefly reviewing the current tests of gravitation at Solar System scales, we give motivations to continue such experiments. In order to obtain signature and estimate the amplitude of anomalous signals that could show up in radio science observables because of modified gravitational laws, we developed a new software that simulates Range/Doppler signals. We present this new tool that simulates radio science observables directly from the space-time metric. We apply this tool to the Cassini mission during its cruise from Jupiter to Saturn and derive constraints on the parameters entering alternative theories of gravity beyond the standard Parametrized Post Newtonian theory.

The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis).

We present the first dynamical solution of the triple asteroid system (45) Eugenia and its two moons Petit-Prince (Diameter~7 km) and S/2004 (45) 1 (Diameter~5 km). The two moons orbit at 1165 and 610 km from the primary, describing an almost-circular orbit (e~6x10-3 and e~7x10-2 respectively). The system is quite different from the other known triple systems in the main belt since the inclinations of the moon orbits are sizeable (9 deg and 18 deg with respect to the equator of the primary respectively). No resonances, neither secular nor due to Lidov-Kozai mechanism, were detected in our dynamical solution, suggesting that these inclinations are not due to excitation modes between the primary and the moons. A 10-year evolution study shows that the orbits are slightly affected by perturbations from the Sun, and to a lesser extent by mutual interactions between the moons. The estimated J2 of the primary is three times lower than the theoretical one, calculated assuming the shape of the primary and an homogeneous interior, possibly suggesting the importance of other gravitational harmonics.

Any model of tides is based on a specific hypothesis of how lagging depends on the tidal-flexure frequency. For example, Gerstenkorn (1955), MacDonald (1964), and Kaula (1964) assumed constancy of the geometric lag angle, while Singer (1968) and Mignard (1979, 1980) asserted constancy of the time lag. Thus, each of these two models was based on a certain law of scaling of the geometric lag. The actual dependence of the geometric lag on the frequency is more complicated and is determined by the rheology of the planet. Besides, each particular functional form of this dependence will unambiguously fix the appropriate form of the frequency dependence of the tidal quality factor, Q. Since at present we know the shape of the dependence of Q upon the frequency, we can reverse our line of reasoning and single out the appropriate actual frequency-dependence of the angular lag. This dependence turns out to be different from those employed hitherto, and it entails considerable alterations in the time scales of the tide-generated dynamical evolution. Phobos' fall on Mars is an example we consider.

Construction of a theory of orbits about a precessing oblate planet, in terms of osculating elements defined in a frame of the equator of date, was started in Efroimsky and Goldreich (2004) and Efroimsky (2005, 2006). We now combine that analytical machinery with numerics. The resulting semianalytical theory is then applied to Deimos over long time scales. In parallel, we carry out a purely numerical integration in an inertial Cartesian frame. The results agree to within a small margin, for over 10 Myr, demonstrating the applicability of our semianalytical model over long timescales. This will enable us to employ it at the further steps of the project, enriching the model with the tides, the pull of the Sun, and the planet's triaxiality. Another goal of our work was to check if the equinoctial precession predicted for a rigid Mars could have been sufficient to repel the orbits away from the equator. We show that for low initial inclinations, the orbit inclination reckoned from the precessing equator of date is subject only to small variations. This is an extension, to non-uniform precession given by the Colombo model, of an old result obtained by Goldreich (1965) for the case of uniform precession and a low initial inclination. However, near-polar initial inclinations may exhibit considerable variations for up to +/- 10 deg in magnitude. Nevertheless, the analysis confirms that an oblate planet can, indeed, afford large variations of the equinoctial precession over hundreds of millions of years, without repelling its near-equatorial satellites away from the equator of date: the satellite inclination oscillates but does not show a secular increase. Nor does it show secular decrease, a fact that is relevant to the discussion of the possibility of high-inclination capture of Phobos and Deimos.

In this paper we investigate the possibility of measuring the post-Newtonian general relativistic gravitomagnetic Lense-Thirring effect in the Jovian system of its Galilean satellites Io, Europa, Ganymede and Callisto in view of recent developments in processing and modelling their optical observations spanning a large time interval (125 years). The present day best observations have an accuracy between several kilometers to few tens of kilometers, which is just the order of magnitude of the Lense-Thirring shifts of the orbits of the Galilean satellites over almost a century. From a comparison between analytical development and numerical integration it turns out that, unfortunately, most of the secular component of the gravitomagnetic signature is removed in the process of fitting the initial conditions. Indeed, an estimation of the magnitude of the Lense-Thirring effect in the ephemerides residuals is given; the resulting residuals have a maximum magnitude of 20 meters only (over 125 years).