We use the planetary ephemeris INPOP17b to constrain the mass of the graviton in the Newtonian limit. We also give an interpretation of this result for a specific case of fifth force framework. We find that the residuals for the Cassini spacecraft significantly (90\% C.L.) degrade for Compton wavelengths of the graviton smaller than $1.83\times 10^{13}$ km, corresponding to a graviton mass bigger than $6.76\times 10^{-23} eV/c^2$. This limit is comparable in magnitude to the one obtained by the LIGO-Virgo collaboration in the radiative regime. We also use this specific example to illustrate why constraints on alternative theories of gravity obtained from postfit residuals are generically overestimated.

Lunar laser ranging (LLR) data and Apollo seismic data analyses, revealed independent evidence for the presence of a fluid lunar core. However, the size of the lunar fluid core remained uncertain by $\pm55$ km (encompassing two contrasting 2011 Apollo seismic data analyses). Here we show that a new description of the lunar interior's dynamical model provides a determination of the radius and geometry of the lunar core-mantle boundary (CMB) from the LLR observations. We compare the present-day lunar core oblateness obtained from LLR analysis with the expected hydrostatic model values, over a range of previously expected CMB radii. The findings suggest a core oblateness ($f_c=(2.2\pm0.6)\times10^{-4}$) that satisfies the assumption of hydrostatic equilibrium over a tight range of lunar CMB radii ($\mathcal{R}_{CMB}=381\pm12$ km). Our estimates of a presently-relaxed lunar CMB translates to a core mass fraction in the range of $1.59-1.77\%$ with a present-day Free Core Nutation (FCN) within $(367\pm100)$ years.

Context. CM-like asteroids (Ch and Cgh classes) are a major population within the broader C-complex, encompassing about 10% of the mass of the main asteroid belt. Their internal structure has been predicted to be homogeneous, based on their compositional similarity as inferred from spectroscopy (Vernazza et al., 2016, AJ 152, 154) and numerical modeling of their early thermal evolution (Bland & Travis, 2017, Sci. Adv. 3, e1602514). Aims. Here we aim to test this hypothesis by deriving the density of the CM-like asteroid (41) Daphne from detailed modeling of its shape and the orbit of its small satellite. Methods. We observed Daphne and its satellite within our imaging survey with the Very Large Telescope extreme adaptive-optics SPHERE/ZIMPOL camera (ID 199.C-0074, PI P. Vernazza) and complemented this data set with earlier Keck/NIRC2 and VLT/NACO observations. We analyzed the dynamics of the satellite with our Genoid meta-heuristic algorithm. Combining our high-angular resolution images with optical lightcurves and stellar occultations, we determine the spin period, orientation, and 3-D shape, using our ADAM shape modeling algorithm. Results. The satellite orbits Daphne on an equatorial, quasi-circular, prograde orbit, like the satellites of many other large main-belt asteroids. The shape model of Daphne reveals several large flat areas that could be large impact craters. The mass determined from this orbit combined with the volume computed from the shape model implies a density for Daphne of 1.77+/-0.26 g/cm3 (3 {\sigma}). This density is consistent with a primordial CM-like homogeneous internal structure with some level of macroporosity (~17%). Conclusions. Based on our analysis of the density of Daphne and 75 other Ch/Cgh-type asteroids gathered from the literature, we conclude that the primordial internal structure of the CM parent bodies was homogeneous.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

Asteroid (16) Psyche, that for long was the largest M-type with no detection of hydration features in its spectrum, was recently discovered to have a weak 3 micron band and thus it eventually was added to the group of hydrated asteroids. Its relatively high density in combination with the high radar albedo, led to classify the asteroid as a metallic object, possibly a core of a differentiated body, remnant of "hit-and-run" collisions. The detection of hydration is, in principle, inconsistent with a pure metallic origin of this body. Here we consider the scenario that the hydration on its surface is exogenous and was delivered by hydrated impactors. We show that impacting asteroids that belong to families whose members have the 3 m band can deliver the hydrated material to Psyche. We developed a collisional model with which we test all the dark carbonaceous asteroid families, which contain hydrated members. We find that the major source of hydrated impactors is the family of Themis, with a total implanted mass on Psyche to be of the order of 10^14 kg. However, the hydrated fraction could be only a few per cent of the implanted mass, as the water content in carbonaceous chondrite meteorites, the best analogue for the Themis asteroid family, is typically a few per cent of their mass.

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.

For many years, lunar laser ranging (LLR) observations using a green wavelength have suffered an inhomogeneity problem both temporally and spatially. This paper reports on the implementation of a new infrared detection at the Grasse LLR station and describes how infrared telemetry improves this situation. Our first results show that infrared detection permits us to densify the observations and allows measurements during the new and the full Moon periods. The link budget improvement leads to homogeneous telemetric measurements on each lunar retro-reflector. Finally, a surprising result is obtained on the Lunokhod 2 array which attains the same efficiency as Lunokhod 1 with an infrared laser link, although those two targets exhibit a differential efficiency of six with a green laser link.

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.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

To explain the unusual distribution of Kuiper Belt objects, several authors have advocated the existence of a super-Earth planet in the outer solar system. It has recently been proposed that a 10 M$_{\oplus}$ object with an orbit of 700 AU semi major axis and 0.6 eccentricity can explain the observed distribution of Kuiper Belt objects around Sedna. Here we use the INPOP planetary ephemerides model as a sensor for testing for an additional body in the solar system. We test the possibility of adding the proposed planet without increasing the residuals of the planetary ephemerides, fitted over the whole INPOP planetary data sample. We demonstrate that the presence of such an object is not compatible with the most sensitive data set, the Cassini radio ranging data, if its true anomaly is in the intervals $[-130^\circ:-100^\circ]$ or $[-65^\circ : 85^\circ]$. Moreover, we find that the addition of this object can reduce the Cassini residuals, with a most probable position given by a true anomaly $v = {117.8^\circ}^{ + 11^\circ}_{ - 10^\circ} $.

Planetary ephemerides are a good tool for studying general relativity at the scale of our solar system. We present here new evaluations of advances of perihelia for Mercury and Saturn.

As part of a reorganization of the International Astronomical Union (IAU), Commission 4 (Ephemerides) went out of existence after the IAU General Assembly in August 2015. This paper presents brief discussions of some of the developments in fundamental astronomy that have influenced and been influenced by the work of Commission 4 over its 96-year history. The paper also presents notes about some of the publications of the national institutions that have played an essential role in the commission's mission. The contents of this paper were submitted for Commission 4's final report, to appear in IAU Transactions Vol. XXIX-A.

Based on the use of MESSENGER radiotracking data in the construction of new Mercury ephemerides (Verma et al. 2014) a new planetary ephemerides INPOP13b was built including Mercury improvements but also improvements on the Mars orbit and on the tie of INPOP planetary ephemerides to ICRF in general.

The IAU Commission 4 Working Group on Standardizing Access to Ephemerides recommends the use of the Spacecraft and Planet Kernel (SPK) format as a standard format for the position ephemerides of planets and other natural solar system bodies, and the use of the Planetary Constants Kernel (PCK) format for the orientation of these bodies. It further recommends that other supporting data be stored in a text PCK. These formats were developed for use by the SPICE Toolkit by the Navigation and Ancillary Information Facility of NASA's Jet Propulsion Laboratory (JPL). The CALCEPH library developed by the Institut de mecanique celeste de calcul des ephemerides (IMCCE) is also able to make use of these files. High accuracy ephemerides available in files conforming to the SPK and PCK formats include: the Development Ephemerides (DE) from JPL, Integrateur Numerique Planetaire de l'Observatoire de Paris (INPOP) from IMCCE, and the Ephemerides Planets and the Moon (EPM), developed by the Institute for Applied Astronomy (IAA). The bulk of this report is a description of the portion of PCK and SPK formats required for these ephemerides. New SPK and PCK data types, both called Type 20: Chebyshev (Velocity Only), have been added. Other changes to the specification are (i) a new object identification number for coordinate time ephemerides and (ii) a set of three new data types that use the TCB rather than the TDB time scale for the ephemerides, but are otherwise identical to their TDB versions.

Based on the new developped planetary ephemerides INPOP13c, determinations of acceptable intervals of General Relativity violation in considering simultaneously the PPN parameters $\beta$, PPN $\gamma$, the flattening of the sun $J_{2}^\odot$ and time variation of the gravitational mass of the sun $\mu$ are obtained in using Monte Carlo simulation coupled with basic genetic algorithm. Possible time variations of the gravitational constant G are also deduced. Discussions are lead about the better choice of indicators for the goodness-of-fit for each run and limits consistent with general relativity are obtained simultaneously.

We analyze the impact of the installation of a lunar laser ranging device on the NTT 3.6m ESO telescope. With such an instrument, the scientific communities of fundamental physics and solar system formation will highly benefit of the only LLR station in the Southern Hemisphere. The quality of the NTT 3.6 meter telescope will also greatly complement the LLR 3.5 meter Apache Point telescope (3.5 m) instrument in the Northern Hemisphere (USA) which is the best instrument for tracking the Moon since 2006. Finally, we also consider the technical characteristics of such installation including the observational constraints.

The current knowledge of Mercury orbit has mainly been gained by direct radar ranging obtained from the 60s to 1998 and by five Mercury flybys made by Mariner 10 in the 70s, and MESSENGER made in 2008 and 2009. On March 18, 2011, MESSENGER became the first spacecraft to orbit Mercury. The radioscience observations acquired during the orbital phase of MESSENGER drastically improved our knowledge of the orbit of Mercury. An accurate MESSENGER orbit is obtained by fitting one-and-half years of tracking data using GINS orbit determination software. The systematic error in the Earth-Mercury geometric positions, also called range bias, obtained from GINS are then used to fit the INPOP dynamical modeling of the planet motions. An improved ephemeris of the planets is then obtained, INPOP13a, and used to perform general relativity tests of PPN-formalism. Our estimations of PPN parameters (beta and gamma?) are more stringent than previous results.

The INPOP ephemerides have known several improvements and evolutions since the first INPOP06 release (Fienga et al. 2008) in 2008. In 2010, anticipating the IAU 2012 resolutions, adjustement of the gravitational solar mass with a fixed astronomical unit (AU) has been for the first time implemented in INPOP10a (Fienga et al. 2011) together with improvements in the asteroid mass determinations. With the latest INPOP10e version, such advancements have been enhanced and studies about solar corona have also been investigated (Verma et al. 2012). The use of planetary ephemerides for several physical applications are presented here from electronic densities of solar slow and fast winds to asteroid mass determinations and tests of general relativity operated with INPOP10a. Perspectives will also be drawn especially related to the analysis of the Messenger spacecraft data for the planetary orbits and future computation of the time variations of the gravitational mass of the sun.

We review here the tests of fundamental physics based on the dynamics of solar system objects.

Compared to the previous INPOP versions, the INPOP10a planetary and lunar ephemeris has several improvements. For the planets of our solar system, no big change was brought in the dynamics but improvements were implemented in the fitting process, the data sets used in the fit and in the selection of fitted parameters. We report here the main characteristics of the planetary part of INPOP10a like the fit of the product of the Solar mass with the gravitational constant (GM$_{\odot}$) instead of the astronomical unit. Determinations of PPN parameters as well as adjustments of the Sun J2 and of asteroid masses are also presented. New advances of nodes and perihelia of planets were also estimated and are given here. As for INPOP08, INPOP10a provides to the user, positions and velocities of the planets, the moon, the rotation angles of the Earth and the Moon as well as TT-TDB chebychev polynomials at http://www.imcce.fr/inpop

We present here a new solution for the astronomical computation of the orbital motion of the Earth spanning from 0 to -250 Myr. The main improvement with respect to the previous numerical solution La2004 (Laskar et al. 2004) is an improved adjustment of the parameters and initial conditions through a fit over 1 Myr to a special version of the high accurate numerical ephemeris INPOP08 (Fienga et al. 2009). The precession equations have also been entirely revised and are no longer averaged over the orbital motion of the Earth and Moon. This new orbital solution is now valid over more than 50 Myr in the past or in the future with proper phases of the eccentricity variations. Due to chaotic behavior, the precision of the solution decreases rapidly beyond this time span, and we discuss the behavior of various solutions beyond 50 Myr. For paleoclimate calibrations, we provide several different solutions that are all compatible with the most precise planetary ephemeris. We have thus reached the time where geological data are now required to discriminate among planetary orbital solutions beyond 50 Myr.

The Planetary and Lunar ephemerides INPOP10a version has several improvements in the fitting process, the data sets used in the fit and in the general features of the solution. No big change was brought in the dynamics. As a consequence of these changes, the extrapolation capabilities of INPOP10a are improved compared to INPOP08, especially for the Earth, Mars, Mercury and Saturn orbits. As for INPOP08, INPOP10a provides to the user, positions, velocities of the planets and the moon and TT-TDB chebychev polynomials at \url{http://www.imcce.fr/inpop}.

We assess the ability of a solid ring to model a global perturbation induced by several thousands of main-belt asteroids. The ring is first studied in an analytical framework that provides an estimate of all the ring's parameters excepting mass. In the second part, numerically estimated perturbations on the Earth-Mars, Earth-Venus, and Earth-Mercury distances induced by various subsets of the main-belt population are compared with perturbations induced by a ring. To account for large uncertainties in the asteroid masses, we obtain results from Monte Carlo experiments based on asteroid masses randomly generated according to available data and the statistical asteroid model. The radius of the ring is analytically estimated at 2.8 AU. A systematic comparison of the ring with subsets of the main belt shows that, after removing the 300 most perturbing asteroids, the total main-belt perturbation of the Earth-Mars distance reaches on average 246 m on the 1969-2010 time interval. A ring with appropriate mass is able to reduce this effect to 38 m. We show that, by removing from the main belt ~240 asteroids that are not necessarily the most perturbing ones, the corresponding total perturbation reaches on average 472 m, but the ring is able to reduce it down to a few meters, thus accounting for more than 99% of the total effect.

The latest version of the planetary ephemerides developed at the Paris Observatory and at the Besancon Observatory is presented here. INPOP08 is a 4-dimension ephemeris since it provides to users positions and velocities of planets and the relation between TT and TDB. Investigations leading to improve the modeling of asteroids are described as well as the new sets of observations used for the fit of INPOP08. New observations provided by the European Space Agency (ESA) deduced from the tracking of the Mars Express (MEX) and Venus Express (VEX) missions are presented as well as the normal point deduced from the Cassini mission. We show the huge impact brought by these observations in the fit of INPOP08, especially in terms of Venus, Saturn and Earth-Moon barycenter orbits.

We summarise the scientific and technological aspects of the SAGAS (Search for Anomalous Gravitation using Atomic Sensors) project, submitted to ESA in June 2007 in response to the Cosmic Vision 2015-2025 call for proposals. The proposed mission aims at flying highly sensitive atomic sensors (optical clock, cold atom accelerometer, optical link) on a Solar System escape trajectory in the 2020 to 2030 time-frame. SAGAS has numerous science objectives in fundamental physics and Solar System science, for example numerous tests of general relativity and the exploration of the Kuiper belt. The combination of highly sensitive atomic sensors and of the laser link well adapted for large distances will allow measurements with unprecedented accuracy and on scales never reached before. We present the proposed mission in some detail, with particular emphasis on the science goals and associated measurements.