A double planet system or planet binary undergoes eclipses that modify the reflective light curve. In the time domain, the eclipse events are fast and weak. This would make their signal difficult to find and recognize in the phase light curve, even for small inclinations when eclipses happen frequently. However, due to the quasiperiodic nature of the phenomenon, the Fourier transform of the direct reflection signal consists of a double sum of sharp peaks. These peaks can be resolved for large close binaries and sufficiently long observation times with a star coronagraph. Eclipses modulate the phase curve, having an orbital period $2\pi/\omega$, with a contribution from the relative motion in the binary plane of a period $2\pi/\Omega$. This leads to a spectral structure with basis frequencies $\omega$ and $\Omega$. We aim to characterize these spectra. We studied the regime of short eclipses that occur when the planet radii are small compared to the planet separation. We derived formulas for the peak amplitudes applicable to homogeneous (Lambertian) planet binaries in circular orbit with small inclination. The effects of an eclipse and of double reflection appear as first- and second-order contributions (in planet radius over separation) in the reflection signal respectively. Small peaks appear as observable side bands in the spectrum. Identical structures around $m\Omega$ are characteristic of short-duration eclipses. Deceasing side bands could indicate double reflection between companions. Fourier analysis of the light curve of non-transiting planets can be used to find planets and their moons. Difficulties in interpreting the structures arise for small planet separation and when there are several moons in mean-motion resonance.

The Weak Equivalence Principle (WEP), stating that two bodies of different compositions and/or mass fall at the same rate in a gravitational field (universality of free fall), is at the very foundation of General Relativity. The MICROSCOPE mission aims to test its validity to a precision of $10^{-15}$, two orders of magnitude better than current on-ground tests, by using two masses of different compositions (titanium and platinum alloys) on a quasi-circular trajectory around the Earth. This is realised by measuring the accelerations inferred from the forces required to maintain the two masses exactly in the same orbit. Any significant difference between the measured accelerations, occurring at a defined frequency, would correspond to the detection of a violation of the WEP, or to the discovery of a tiny new type of force added to gravity. MICROSCOPE's first results show no hint for such a difference, expressed in terms of E\"otv\"os parameter $\delta(Ti,Pt)=[-1\pm{}9{\rm (stat)}\pm{}9{\rm (syst)}] \times{}10^{-15}$ (both 1$\sigma$ uncertainties) for a titanium and platinum pair of materials. This result was obtained on a session with 120 orbital revolutions representing 7\% of the current available data acquired during the whole mission. The quadratic combination of 1$\sigma$ uncertainties leads to a current limit on $\delta$ of about $1.3\times{}10^{-14}$.

At present, tracking data for planetary missions largely consists of radio observables: range-rate range and angular position. Future planetary missions may use Interplanetary Laser Ranging (ILR) as a tracking observable. Two-way ILR will provide range data that are about 2 orders of magnitude more accurate than radio-based range data. ILR does not produce Doppler data, however. In this article, we compare the relative strength of radio Doppler and laser range data for the retrieval of parameters of interest in planetary missions, to clarify and quantify the science case of ILR, with a focus on geodetic observables. We first provide an overview of the near-term attainable quality of ILR, in terms of both the realization of the observable and the models used to process the measurements. Subsequently, we analyze the sensitivity of radio-Doppler and laser-range measurements in representative mission scenarios. We use both an analytical approximation and numerical analyses of the relative sensitivity of ILR and radio Doppler observables for more general cases. We show that mm-precise range normal points are feasible for ILR, but mm-level accuracy and stability is unlikely to be attained, due to a combination of instrumental and model errors. ILR has the potential for superior performance in observing signatures in the data with a characteristic period of greater than 0.33-1.65 hours This indicates that Doppler tracking will typically remain the method of choice for gravity field determination and spacecraft orbit determination in planetary missions. Laser ranging data, however, are shown to have a significant advantage for the retrieval of rotational and tidal characteristics from landers. Similarly, laser ranging data will be superior for the construction of planetary ephemerides and the improvement of solar system tests of gravitation.

According to the Weak Equivalence Principle, all bodies should fall at the same rate in a gravitational field. The MICROSCOPE satellite, launched in April 2016, aims to test its validity at the $10^{-15}$ precision level, by measuring the force required to maintain two test masses (of titanium and platinum alloys) exactly in the same orbit. A non-vanishing result would correspond to a violation of the Equivalence Principle, or to the discovery of a new long-range force. Analysis of the first data gives $\delta\rm{(Ti,Pt)}= [-1 \pm 9 (\mathrm{stat}) \pm 9 (\mathrm{syst})] \times 10^{-15}$ (1$\sigma$ statistical uncertainty) for the titanium-platinum E\"otv\"os parameter characterizing the relative difference in their free-fall accelerations.

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.

Many physical parameters that can be estimated from space mission tracking data influence both the translational dynamics and proper time rates of observers. These different proper time rates cause a variability of the time transfer observable beyond that caused by their translational (and rotational) dynamics. With the near-future implementation of (interplanetary) transponder laser ranging, these effects will become increasingly important, requiring a re-evaluation of the common data analysis practice of using a priori time ephemerides, which is the goal of this paper. We develop a framework for the simultaneous estimation of the initial translational state and the initial proper time of an observer, with the goal of facilitating robust tracking data analysis from next-generation space missions carrying highly accurate clocks and tracking equipment. Using our approach, the influence of physical parameters on both translational and time dynamics are considered at the same level in the analysis, and mutual correlations between the signatures of the two are automatically identified. We perform a covariance analysis using our proposed method with simulated laser data from Earth-based stations to both a Mars and Mercury lander. Using four years of tracking data for the Mars lander simulations, we find a difference between our results using the simultaneous space-time dynamics estimation and the classical analysis technique (with an ta priori time ephemeris) of around 0.1 % in formal errors and correlation coefficients. For a Mercury lander this rises to around 1% for a 1-month mission and 10 % for a 4-year mission. By means of Monte Carlo simulation, we find that using an a priori time ephemeris of representative accuracy will result in estimation errors that are orders of magnitude above the formal error when processing highly accurate laser time transfer data.

The weak orbital-phase dependent reflection signal of an exoplanet contains information on the planet surface, such as the distribution of continents and oceans on terrestrial planets. This light curve is usually studied in the time domain, but because the signal from a stationary surface is (quasi)periodic, analysis of the Fourier series may provide an alternative, complementary approach. We study Fourier spectra from reflected light curves for geometrically simple configurations. Depending on its atmospheric properties, a rotating planet in the habitable zone could have circular polar ice caps. Tidally locked planets, on the other hand, may have symmetric circular oceans facing the star. These cases are interesting because the high-albedo contrast at the sharp edges of the ice-sheets and the glint from the host star in the ocean may produce recognizable light curves with orbital periodicity, which could also be interpreted in the Fourier domain. We derive a simple general expression for the Fourier coefficients of a quasiperiodic light curve in terms of the albedo map of a Lambertian planet surface. Analytic expressions for light curves and their spectra are calculated for idealized situations, and dependence of spectral peaks on the key parameters inclination, obliquity, and cap size is studied.