The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

An observation of Jupiter's tidal response is anticipated for the on-going Juno spacecraft mission. We combine self-consistent, numerical models of Jupiter's equilibrium tidal response with observed Doppler shifts from the Juno gravity science experiment to test the sensitivity of the spacecraft to tides raised by the Galilean satellites and the Sun. The concentric Maclaurin spheroid (CMS) method finds the equilibrium shape and gravity field of a rotating, liquid planet with the tide raised by a satellite, expanded in Love numbers ($k_{nm}$). We present improvements to CMS theory that eliminate an unphysical center of mass offset and study in detail the convergence behavior of the CMS approach. We demonstrate that the dependence of $k_{nm}$ with orbital distance is important when considering the combined tidal response for Jupiter. Conversely, the details of the interior structure have a negligible influence on $k_{nm}$, for models that match the zonal harmonics $J_2$, $J_4$ and $J_6$, already measured to high precision by Juno. As the mission continues, improved coverage of Jupiter's gravity field at different phases of Io's orbit is expected to yield an observed value for the degree-2 Love number ($k_{22}$) and potentially select higher--degree $k_{nm}$. We present a test of the sensitivity of the Juno Doppler signal to the calculated $k_{nm}$, which suggests the detectability of $k_{33}$, $k_{42}$ and $k_{31}$, in addition to $k_{22}$. A mismatch of a robust Juno observation with the remarkably small range in calculated Io equilibrium $k_{22}=0.58976\pm0.0001$, would indicate a heretofore uncharacterized dynamic contribution to the tides.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

LaRa (Lander Radioscience) is an experiment on the ExoMars 2020 mission that uses the Doppler shift on the radio link due to the motion of the ExoMars platform tied to the surface of Mars with respect to the Earth ground stations (e.g. the deep space network stations of NASA), in order to precisely measure the relative velocity of the lander on Mars with respect to the Earth. The LaRa measurements shall improve the understanding of the structure and processes in the deep interior of Mars by obtaining the rotation and orientation of Mars with a better precision compared to the previous missions. In this paper, we provide the analysis done until now for the best realization of these objectives. We explain the geophysical observation that will be reached with LaRa (Length-of-day variations, precession, nutation, and possibly polar motion). We develop the experiment set up, which includes the ground stations on Earth (so-called ground segment). We describe the instrument, i.e. the transponder and its three antennas. We further detail the link budget and the expected noise level that will be reached. Finally, we detail the expected results, which encompasses the explanation of how we shall determine Mars' orientation parameters, and the way we shall deduce Mars' interior structure and Mars' atmosphere from them. Lastly, we explain briefly how we will be able to determine the Surface platform position.

The Laser Ranging Interferometer (LRI) instrument on the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has provided the first laser interferometric range measurements between remote spacecraft, separated by approximately 220 km. Autonomous controls that lock the laser frequency to a cavity reference and establish the 5 degree of freedom two-way laser link between remote spacecraft succeeded on the first attempt. Active beam pointing based on differential wavefront sensing compensates spacecraft attitude fluctuations. The LRI has operated continuously without breaks in phase tracking for more than 50 days, and has shown biased range measurements similar to the primary ranging instrument based on microwaves, but with much less noise at a level of $1\,{\rm nm}/\sqrt{\rm Hz}$ at Fourier frequencies above 100 mHz.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

In-spiraling supermassive black holes should emit gravitational waves, which would produce characteristic distortions in the time of arrival residuals from millisecond pulsars. Multiple national and regional consortia have constructed pulsar timing arrays by precise timing of different sets of millisecond pulsars. An essential aspect of precision timing is the transfer of the times of arrival to a (quasi-)inertial frame, conventionally the solar system barycenter. The barycenter is determined from the knowledge of the planetary masses and orbits, which has been refined over the past 50 years by multiple spacecraft. Within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), uncertainties on the solar system barycenter are emerging as an important element of the NANOGrav noise budget. We describe what is known about the solar system barycenter, touch upon how uncertainties in it affect gravitational wave studies with pulsar timing arrays, and consider future trends in spacecraft navigation.

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.

The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).

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.

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.

Lunar laser ranging (LLR) has made major contributions to our understanding of the Moon's internal structure and the dynamics of the Earth-Moon system. Because of the recent improvements of the ground-based laser ranging facilities, the present LLR measurement accuracy is limited by the retro-reflectors currently on the lunar surface, which are arrays of small corner-cubes. Because of lunar librations, the surfaces of these arrays do not, in general, point directly at the Earth. This effect results in a spread of arrival times, because each cube that comprises the retroreflector is at a slightly different distance from the Earth, leading to the reduced ranging accuracy. Thus, a single, wide aperture corner-cube could have a clear advantage. In addition, after nearly four decades of successful operations the retro-reflectors arrays currently on the Moon started to show performance degradation; as a result, they yield still useful, but much weaker return signals. Thus, fresh and bright instruments on the lunar surface are needed to continue precision LLR measurements. We have developed a new retro-reflector design to enable advanced LLR operations. It is based on a single, hollow corner cube with a large aperture for which preliminary thermal, mechanical, and optical design and analysis have been performed. The new instrument will be able to reach an Earth-Moon range precision of 1-mm in a single pulse while being subjected to significant thermal variations present on the lunar surface, and will have low mass to allow robotic deployment. Here we report on our design results and instrument development effort.

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.

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

This note suggests nomenclature for dealing with the units of various astronomical quantities that are used with the relativistic time scales TT, TDB, TCB and TCG. It is suggested to avoid wordings like "TDB units" and "TT units" and avoid contrasting them to "SI units". The quantities intended for use with TCG, TCB, TT or TDB should be called "TCG-compatible", "TCB-compatible", "TT-compatible" or "TDB-compatible", respectively. The names of the units second and meter for numerical values of all these quantities should be used with out any adjectives. This suggestion comes from a special discussion forum created within IAU Commission 52 "Relativity in Fundamental Astronomy".

We present the basic elements and first results of an end-to-end simulation package whose purpose is to test the validity of the Space Interferometer Mission design. The fundamental simulation time step is one millisecond, with substructure at 1/8 ms, and the total duration of the simulation is five years. The end product of a given "wide-angle" astrometry run is an estimated grid star catalog over the entire sky with an accuracy of about 4 micro-arcseconds. SIMsim is divided into five separate modules that communicate via data pipes. The first generates the 'truth' data on the spacecraft structure and laser metrology. The second module generates uncorrupted fringes for the three SIM interferometers, based on the current spacecraft orientation, target stars' positions, etc. The third module reads out the CCD white light fringe data at specified times, corrupting that and the metrology data with appropriate errors. The data stream out of this module represents the basic data stream on the simulated spacecraft. The fourth module performs fringe-fitting tasks on this data, recovering the total path delay, and the fifth and final module inverts the entire metrology/delay dataset to ultimately determine the instantaneous path delay on a fiducial baseline fixed in space. (Pathlength feed forward is used every few milliseconds to re-position the interferometer to keep the fringes in the delay window.) The average of all such delays over an integration time (typically 30 s) is reported as one of several hundred thousand measured stellar delays over the five-year period, which are then inverted to produce the simulated catalog. Future plans include taking into account more sources of error from the SIM error budget and including narrow angle observations in the observing plan.

We present the enabling concept and technology for a dual spacecraft formation-flying optical interferometer, to be launched into a deep space orbit as Space Technology 3. The combiner spacecraft makes use of a nested cat's eye delay line configuration that minimizes wavefront distortion and stores 20 m of optical pathlength in a package of \sim 1.5 m length. A parabolic trajectory for the secondary collector spacecraft enables baselines of up to 200 m for a fixed 20 m stored delay and spacecraft separations of up to 1 km.