We report integrated orbital fits for the inner regular moons of Neptune based on the most complete astrometric data set to date, with observations from Earth-based telescopes, Voyager 2, and the Hubble Space Telescope covering 1981-2016. We summarize the results in terms of state vectors, mean orbital elements, and orbital uncertainties. The estimated masses of the two innermost moons, Naiad and Thalassa, are $GM_{Naiad}$= 0.0080 $\pm$ 0.0043 $km^3 s^{-2}$ and $GM_{Thalassa}$=0.0236 $\pm$ 0.0064 $km^3 s^{-2}$, corresponding to densities of 0.80 $\pm$ 0.48 $g cm^{-3}$ and 1.23 $\pm$ 0.43 $g cm^{-3}$, respectively. Our analysis shows that Naiad and Thalassa are locked in an unusual type of orbital resonance. The resonant argument 73 $\dot{\lambda}_{Thalassa}$-69 $\dot{\lambda}_{Naiad}$-4 $\dot{\Omega}_{Naiad}$ $\approx$ 0 librates around 180 deg with an average amplitude of ~66 deg and a period of ~1.9 years for the nominal set of masses. This is the first fourth-order resonance discovered between the moons of the outer planets. More high precision astrometry is needed to better constrain the masses of Naiad and Thalassa, and consequently, the amplitude and the period of libration. We also report on a 13:11 near-resonance of Hippocamp and Proteus, which may lead to a mass estimate of Proteus provided that there are future observations of Hippocamp. Our fit yielded a value for Neptune's oblateness coefficient of $J_2$=3409.1$\pm$2.9 $\times 10^{-6}$.

We report the discovery of 12 new satellites of Jupiter, giving Jupiter 79 known satellites. The new finds are between 23rd-24th mag in the r-band and 1-3 km in diameter assuming dark albedos. Nine of the discoveries are in the distant retrograde satellite groupings. Two of the new satellites are in the closer Himalia prograde group near 28 degrees in inclination. S/2016 J2, nicknamed Valetudo, has an orbit unlike any other known outer satellite and is the most distant prograde satellite around any planet at 0.36 Hill radii. Numerical simulations show S/2016 J2 is very stable, with average and range of i=34.2+-3 deg, e=0.216+-0.125, and a=18.9+-0.7 million km over 100 Myrs. Our stability simulations show a S/2016 J2 like orbit would be stable out to a=21.8 million km or 0.41 Hill radii, but no further, unlike more distant and eccentric retrograde satellites. S/2016 J2's large semi-major axis means it significantly overlaps the orbits of the distant retrogrades. A prograde-retrograde moon-moon collision between outer satellites of Jupiter has likely happened over the age of the solar system.

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.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

Planetary ephemerides have been developed and improved over centuries. They are a fundamental tool for understanding solar system dynamics, and essential for planetary and small body mass determinations, occultation predictions, high-precision tests of general relativity, pulsar timing, and interplanetary spacecraft navigation. This paper presents recent results from a continuing program of high-precision astrometric very-long-baseline interferometry (VLBI) observations of the Cassini spacecraft orbiting Saturn, using the Very Long Baseline Array (VLBA). We have previously shown that VLBA measurements can be combined with spacecraft orbit determinations from Doppler and range tracking and VLBI links to the inertial extragalactic reference frame (ICRF) to provide the most accurate barycentric positions currently available for Saturn. Here we report an additional five years of VLBA observations along with improved phase reference source positions, resulting in an improvement in residuals with respect to the Jet Propulsion Laboratory's dynamical ephemeris.

The MOdified Newtonian Dynamics (MOND) is an attempt to modify the gravitation theory to solve the Dark Matter problem. This phenomenology is very successful at the galactic level. The main effect produced by MOND in the Solar System is called the External Field Effect parametrized by the parameter $Q_2$. We have used 9 years of Cassini range and Doppler measurements to constrain $Q_2$. Our estimate of this parameter based on Cassini data is given by $Q_2=(3 \pm 3)\times 10^{-27} \ \rm{s^{-2}}$ which shows no deviation from General Relativity and excludes a large part of the relativistic MOND theories. This limit can also be interpreted as a limit on a external tidal potential acting on the Solar System coming from the internal mass of our galaxy (including Dark Matter) or from a new hypothetical body.

In this communication, the current tests of gravitation available at Solar System scales are recalled. These tests rely mainly on two frameworks: the PPN framework and the search for a fifth force. Some motivations are given to look for deviations from General Relativity in other frameworks than the two extensively considered. A recent analysis of Cassini data in a MOND framework is presented. Furthermore, possibilities to constrain Standard Model Extension parameters using Solar System data are developed.

The planetary ephemeris is an essential tool for interplanetary spacecraft navigation, studies of solar system dynamics (including, for example, barycenter corrections for pulsar timing ephemeredes), the prediction of occultations, and tests of general relativity. We are carrying out a series of astrometric VLBI observations of the Cassini spacecraft currently in orbit around Saturn, using the Very Long Baseline Array (VLBA). These observations provide positions for the center of mass of Saturn in the International Celestial Reference Frame (ICRF) with accuracies ~0.3 milli-arcsecond (1.5 nrad), or about 2 km at the average distance of Saturn. This paper reports results from eight observing epochs between 2006 October and 2009 April. These data are combined with two VLBA observations by other investigators in 2004 and a Cassini-based gravitational deflection measurement by Fomalont et al. in 2009 to constrain a new ephemeris (DE 422). The DE 422 post-fit residuals for Saturn with respect to the VLBA data are generally 0.2 mas, but additional observations are needed to improve the positions of all of our phase reference sources to this level. Over time we expect to be able to improve the accuracy of all three coordinates in the Saturn ephemeris (latitude, longitude, and range) by a factor of at least three. This will represent a significant improvement not just in the Saturn ephemeris but also in the link between the inner and outer solar system ephemeredes and in the link to the inertial ICRF.