Analysis of Preferred Directions in Phase Space for Tidal Measurements at Europa

The NASA Jupiter Europa Orbiter mission requires a circular, near-polar orbit to measure Europa's Love numbers, geophysical coefficients which give insight into whether a liquid ocean exists. This type of orbit about planetary satellites is known to be unstable. The effects of Jupiter's tidal gravity are seen in changes in Europa's gravity field and surface deformation, which are sensed through doppler tracking over time and altimetry measurements respectively. These two measurement types separately determine the h and k Love numbers, a combination of which bounds how thick the ice shell of Europa is and whether liquid water is present. This work shows how the properties of an unstable periodic orbit about Europa generate preferred measurement directions in position and velocity phase space for the orbit determination process. We generate an error covariance over seven days for the orbiter state and science parameters using a periodic orbit and then disperse the orbit initial conditions in a Monte Carlo simulation according to this covariance. The dispersed orbits are shown to have a bias toward longer lifetimes and we discuss this as an effect of the stable and unstable manifolds of the periodic orbit. Using an epoch formulation of a square-root information filter, measurements aligned with the unstable manifold mapped back in time add more information to the orbit determination process than measurements aligned with the stable manifold. This corresponds to a contraction in the uncertainty of the estimate of the desired parameters, including the Love numbers. We demonstrate this mapping mathematically using a representation of the State Transition Matrix involving its eigenvectors and eigenvalues. Then using the properties of left and right eigenvectors, we show how measurements in the orbit determination process are mapped in time leading to a concentration of information at epoch. We present examples of measurements taken on different time schedules to show the effect of preferred phase space directions in the estimation process. Manifold coordinate decomposition is applied to the orbit initial conditions as well as measurement partials in the filter to show the alignment of each with the stable and unstable manifolds of the periodic orbit. The connection between orbit lifetime and regions of increased information density in phase space is made using the properties of these manifolds. Low altitude, near-polar periodic orbits with these characteristics are discussed along with the estimation results for the Love numbers, orbiter state, and orbit lifetime. Different measurement schedules and the resulting estimation performance are presented along with an analysis of information content for single measurements with respect to manifold alignment. These results allow more sensitive estimation of the tidal Love numbers and may allow measurements to be taken less frequently or compensate for corrupted data arcs. Other measurement types will be mapped in the same way using the State Transition matrix and have increased information density at epoch if aligned with the unstable manifold. In the same way, these results are applicable to planetary satellite orbiters about Enceladus or Dione since they share the governing equations of motion and properties of unstable periodic orbits.