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.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.

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}$.

In 1984, three telescopes in South America recorded an occultation of a star near Neptune. It was attributed to the existence of a partial ring or ring arc. The existence of ring arcs around Neptune was confirmed during subsequent years via other occultation experiments and by the Voyager 2 spacecraft. The Voyager observations established that the Neptunian arcs are concentrations of particles embedded within Neptune's narrow Adams ring, the outermost of six tenuous rings discovered by Voyager and discussed here. Four ring arcs were identified: the trailing arc Fraternit\'e, a double-component arc Egalit\'e, dubbed Egalit\'e 1 and 2, Libert\'e, and the leading arc Courage. The arcs varied in extent from $\sim$ 1$^\circ$ to $\sim$ 10$^\circ$, and together were confined to a longitude range of 40$^\circ$, with typical radial widths of $\sim$ 15 km and optical depth of order 0.1. The properties of the dusty component of Neptune's rings are also discussed in this chapter. Although the arcs should have been destroyed in a few months time through differential Keplerian motion, they appeared to persist at least throughout the Voyager era, and well beyond. However, observations from Earth (both with the Hubble Space Telescope and ground-based adaptive optics) show an evolution in the last three decades, with the disappearance of both Courage and Libert\'e by 2009. This chapter reviews the constraints provided on the mean motion of the arcs and Galatea, the satellite possibly responsible for the arc confinement. This in turn constrains the various models that have been proposed to explain the arc longevity.

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

We searched for dust or debris rings in the Pluto-Charon system before, during, and after the New Horizons encounter. Methodologies included searching for back-scattered light during the approach to Pluto (phase $\sim15^\circ$), in situ detection of impacting particles, a search for stellar occultations near the time of closest approach, and by forward-scattered light during departure (phase $\sim165^\circ$). A search using HST prior to the encounter also contributed to the results. No rings, debris, or dust features were observed, but our detection limits provide an improved picture of the environment throughout the Pluto-Charon system. Searches for rings in back-scattered light covered 35,000-250,000 km from the system barycenter, a zone that starts interior to the orbit of Styx, and extends to four times the orbital radius of Hydra. We obtained our firmest limits using the NH LORRI camera in the inner half of this region. Our limits on the normal $I/F$ of an unseen ring depends on the radial scale of the rings: $2\times10^{-8}$ ($3\sigma$) for 1500 km wide rings, $1\times10^{-8}$ for 6000 km rings, and $7\times10^{-9}$ for 12,000 km rings. Beyond $\sim100,000$ km from Pluto, HST observations limit normal $I/F$ to $\sim8\times10^{-8}$. Searches for dust from forward-scattered light extended from the surface of Pluto to the Pluto-Charon Hill sphere ($r_{\rm Hill}=6.4\times10^6$ km). No evidence for rings or dust was detected to normal $I/F$ limits of $\sim8.9\times10^{-7}$ on $\sim10^4$ km scales. Four occulation observations also probed the space interior to Hydra, but again no dust or debris was detected. Elsewhere in the solar system, small moons commonly share their orbits with faint dust rings. Our results suggest that small grains are quickly lost from the system due to solar radiation pressure, whereas larger particles are unstable due to perturbations by the known moons.

A review of the jovian ring system

Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected to have tenuous dusty rings, though they have so far evaded detection. More remotely, transit events observed around a star in 2007 may have revealed for the first time exoplanetary rings around a giant planet orbiting that star. So, evidence is building to show that rings are more common features in the universe than previously thought. Several interesting issues arising from the discovery (or suspicion) of new ring systems are described in this chapter.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

The New Horizons mission has provided resolved measurements of Pluto's moons Styx, Nix, Kerberos, and Hydra. All four are small, with equivalent spherical diameters of $\approx$40 km for Nix and Hydra and ~10 km for Styx and Kerberos. They are also highly elongated, with maximum to minimum axis ratios of $\approx$2. All four moons have high albedos ( $\approx$50-90 %) suggestive of a water-ice surface composition. Crater densities on Nix and Hydra imply surface ages $\gtrsim$ 4 Ga. The small moons rotate much faster than synchronous, with rotational poles clustered nearly orthogonal to the common pole directions of Pluto and Charon. These results reinforce the hypothesis that the small moons formed in the aftermath of a collision that produced the Pluto-Charon binary.

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.

The James Webb Space Telescope (JWST) will provide unprecedented opportunities to observe the rings and small satellites in our solar system, accomplishing three primary objectives: 1) discovering new rings and moons, 2) unprecedented spectroscopy, and 3) time-domain observations. We give details on these science objectives and describe requirements that JWST must fulfill in order to accomplish the science objectives

Images obtained by the Cassini spacecraft between 2012 and 2015 reveal a periodic brightness variation in a region of Saturn's D ring that previously appeared to be rather featureless. Furthermore, the intensity and radial wavenumber of this pattern have decreased steadily with time since it was first observed. Based on analogies with similar structures elsewhere in the D ring, we propose that this structure was created by some event that disturbed the orbital motions of the ring particles, giving them finite orbital eccentricities and initially aligned pericenters. Differential orbital precession then transformed this structure into a spiral pattern in the ring's optical depth that became increasingly tightly wound over time. The observed trends in the pattern's radial wavenumber are roughly consistent with this basic model, and also indicate that the ring-disturbing event occurred in early December 2011. Similar events in 1979 may have generated the periodic patterns seen in this same region by the Voyager spacecraft. The 2011 event could have been caused by debris striking the rings, or by a disturbance in the planet's electromagnetic environment. The rapid reduction in the intensity of the brightness variations over the course of just a few years indicates that some process is either damping orbital eccentricities in this region or causing the orbital pericenters of particles with the same semi-major axis to become misaligned.

As the 7 December 2007 equinox of Uranus approached, ring and atmosphere observers produced a substantial collection of observations using the 10-m Keck telescope and the Hubble Space Telescope. Those spanning the period from 7 June 2007 through 9 September 2007 we used to identify and track cloud features, determine atmospheric motions, characterize cloud morphology and dynamics, and define changes in atmospheric band structure. We confirmed the existence of the suspected northern hemisphere prograde jet, locating its peak near 58 N, and extended wind speed measurements to 73 N. For 28 cloud features we obtained extremely high wind-speed accuracy through extended tracking times. The new results confirm a small N-S asymmetry in the zonal wind profile, and the lack of any change in the southern hemisphere between 1986 (near solstice) and 2007 (near equinox) suggests that the asymmetry may be permanent rather than seasonally reversing. In the 2007 images we found two prominent groups of discrete cloud features with very long lifetimes. The one near 30 S has departed from its previous oscillatory motion and started a significant northward drift, accompanied by substantial morphological changes. The complex of features near 30 N remained at a nearly fixed latitude, while exhibiting some characteristics of a dark spot accompanied by bright companion features. Smaller and less stable features were used to track cloud motions at other latitudes, some of which lasted over many planet rotations, though many could not be tracked beyond a single transit. A bright band has developed near 45 N, while the bright band near 45 S has begun to decline, both events in agreement with the idea that the asymmetric band structure of Uranus is a delayed response to solar forcing, but with a surprisingly short delay of only a few years.

Previous investigations of Saturn's outer D ring (73,200-74,000 km from Saturn's center) identified periodic brightness variations whose radial wavenumber increased linearly over time. This pattern was attributed to a vertical corrugation, and its temporal variability implied that some event --possibly an impact with interplanetary debris-- caused the ring to become tilted out the planet's equatorial plane in 1983. This work examines these patterns in greater detail using a more extensive set of Cassini images in order to obtain additional insights into the 1983 event. These additional data reveal that the D ring is not only corrugated, but also contains a time-variable periodic modulation in its optical depth that probably represents organized eccentric motions of the D-ring's particles. This second pattern suggests that whatever event tilted the rings also disturbed the radial or azimuthal velocities of the ring particles. Furthermore, the relative amplitudes of the two patterns indicate that the vertical motions induced by the 1983 event were 2.3+/-0.5 times larger than the corresponding in-plane motions. If these structures were indeed produced by an impact, material would need to strike the ring at a steep angle (more than 60 degrees from the ring plane) to produce such motions. Meanwhile, the corrugation wavelengths in the D ring are about 0.7% shorter than one would predict based on extrapolations from similar structures in the nearby C ring. This could indicate that the D-ring was tilted/disturbed about 60 days before the C ring. Such a timing difference could be explained if the material that struck the ring was derived from debris released when some object broke up near Saturn some months earlier. This debris would need to have a substantial initial velocity dispersion and have its orbital properties perturbed by solar tides prior to its collision with the rings.

Saturn's F ring is subject to dynamic structural changes over short periods. Among the observed phenomena are diffuse extended bright clumps (ECs) ~ 3-40 degrees in longitudinal extent. These ECs appear, evolve, and disappear over a span of days to months. ECs have been seen by the two Voyager spacecraft, the Cassini orbiter, and various ground- and space-based telescopes. Showalter (2004, Icarus, 171, 356-371) analyzed all Voyager images of the F ring and found that there were 2-3 major and 20-40 minor ECs present in the ring at any given time. We expand upon these results by comparing the ECs seen by Voyager to those seen by Cassini in 2004-2010. We find that the number of minor ECs has stayed roughly constant and the ECs have similar distributions of angular width, absolute brightness, and semimajor axis. However, the common exceptionally bright ECs seen by Voyager are now exceedingly rare, with only two instances seen by Cassini in the six years, and they are now also much dimmer relative to the mean ring background. We hypothesize that these bright ECs are caused by the repeated impacts of small moonlets with the F ring core, and that these moonlets have decreased in number in the 25 years between missions.

We have explored the stability of the inner Uranian satellites using simulations based on the most recent observational data. We find that, across a wide range of mass assumptions, the system is unstable, resulting in the eventual crossing of orbits and probable subsequent collision of moons. Cupid and Belinda are usually the first satellites to cross orbits, and they do so on a time scale of 10^3-10^7 years. Cressida and Desdemona are generally the next pair to cross, on a time scale of 10^5-10^7 years. We show that the crossing times are highly sensitive to initial conditions and that Cupid's instability is related to its resonant interactions with Belinda. We also show that a previously discovered power law, which relates orbit crossing time to satellite mass, is valid across a wide range of masses. We generalize the power law to handle two unstable orbital pairs with overlapping lifetimes and show that it can be used to extend the time span of studies of orbital stability in a computationally efficient manner. Our results suggest that the current Uranian satellite system is in transition and that the moons will continue to collide and reaccrete for the foreseeable future.

Image photometry reveals that the F ring is approximately twice as bright during the Cassini tour as it was during the Voyager flybys of 1980 and 1981. It is also three times as wide and has a higher integrated optical depth. We have performed photometric measurements of more than 4,800 images of Saturn's F ring taken over a five-year period with Cassini's Narrow Angle Camera. We show that the ring is not optically thin in many observing geometries and apply a photometric model based on single-scattering in the presence of shadowing and obscuration, deriving a mean effective optical depth tau = 0.033. Stellar occultation data from Voyager PPS and Cassini VIMS validate both the optical depth and the width measurements. In contrast to this decades-scale change, the baseline properties of the F ring have not changed significantly from 2004 to 2009. However, we have investigated one major, bright feature that appeared in the ring in late 2006. This transient feature increased the ring's overall mean brightness by 84% and decayed with a half-life of 91 days.

D68 is a narrow ringlet located only 67,627 km (1.12 planetary radii) from Saturn's spin axis. Images of this ringlet obtained by the Cassini spacecraft reveal that this ringlet exhibits persistent longitudinal brightness variations and a substantial eccentricity (ae=25+/-1 km). By comparing observations made at different times, we confirm that the brightness variations revolve around the planet at approximately the local orbital rate (1751.6 degrees/day), and that the ringlet's pericenter precesses at 38.243+/-0.008 degrees/day, consistent with the expected apsidal precession rate at this location due to Saturn's higher-order gravitational harmonics. Surprisingly, we also find that the ringlet's semi-major axis appears to be decreasing with time at a rate of 2.4+/-0.4 km/year between 2005 and 2013. A closer look at these measurements, along with a consideration of earlier Voyager observations of this same ringlet, suggests that the mean radius of D68 moves back and forth, perhaps with a period of around 15 Earth years or about half a Saturn year. These observations could place important constraints on both the ringlet's local dynamical environment and the planet's gravitational field.

The Encke Gap is a 320-km-wide opening in Saturn's outer A ring that contains the orbit of the small moon Pan and an array of dusty features composed of particles less than 100 microns across. In particular, there are three narrow ringlets in this region that are not longitudinally homogeneous, but instead contain series of bright clumps. Using images obtained by the Cassini spacecraft, we track the motions of these clumps and demonstrate that they do not follow the predicted trajectories of isolated ring particles moving under the influence of Saturn's and Pan's gravitational fields. We also examine the orbital properties of these ringlets by comparing images taken at different longitudes and times. We find evidence that the orbits of these particles have forced eccentricities induced by solar radiation pressure. In addition, the mean radial positions of the particles in these ringlets appear to vary with local co-rotating longitude, perhaps due to the combined action of drag forces, gravitational perturbations from Pan, and collisions among the ring particles. The dynamics of the dust within this gap therefore appears to be much more complex than previously appreciated.

Saturn's diffuse E ring consists of many tiny (micron and sub-micron) grains of water ice distributed between the orbits of Mimas and Titan. Various gravitational and non-gravitational forces perturb these particles' orbits, causing the ring's local particle density to vary noticeably with distance from the planet, height above the ring-plane, hour angle and time. Using remote-sensing data obtained by the Cassini spacecraft in 2005 and 2006, we investigate the E-ring's three-dimensional structure during a time when the Sun illuminated the rings from the south at high elevation angles (> 15 degrees). These observations show that the ring's vertical thickness grows with distance from Enceladus' orbit and its peak brightness density shifts from south to north of Saturn's equator plane with increasing distance from the planet. These data also reveal a localized depletion in particle density near Saturn's equatorial plane around Enceladus' semi-major axis. Finally, variations are detected in the radial brightness profile and the vertical thickness of the ring as a function of longitude relative to the Sun. Possible physical mechanisms and processes that may be responsible for some of these structures include solar radiation pressure, variations in the ambient plasma, and electromagnetic perturbations associated with Saturn's shadow.

Stellar occultations by Saturn's rings observed with the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft reveal that dusty features such as the F ring and the ringlets in the Encke and the Laplace Gaps have distinctive infrared transmission spectra. These spectra show a narrow optical depth minimum at wavelengths around 2.87 microns. This minimum is likely due to the Christiansen Effect, a reduction in the extinction of small particles when their (complex) refractive index is close to that of the surrounding medium. Simple Mie-scattering models demonstrate that the strength of this opacity dip is sensitive to the size distribution of particles between 1 and 100 microns across. Furthermore, the spatial resolution of the occultation data is sufficient to reveal variations in the transmission spectra within and among these rings. For example, in both the Encke Gap ringlets and F ring, the opacity dip weakens with increasing local optical depth, which is consistent with the larger particles being concentrated near the cores of these rings. The strength of the opacity dip varies most dramatically within the F ring; certain compact regions of enhanced optical depth lack an opacity dip and therefore appear to have a greatly reduced fraction of grains in the few-micron size range.Such spectrally-identifiable structures probably represent a subset of the compact optically-thick clumps observed by other Cassini instruments. These variations in the ring's particle size distribution can provide new insights into the processes of grain aggregation, disruption and transport within dusty rings. For example, the unusual spectral properties of the F-ring clumps could perhaps be ascribed to small grains adhering onto the surface of larger particles in regions of anomalously low velocity dispersion.

We present the first dynamical solution of the triple asteroid system (45) Eugenia and its two moons Petit-Prince (Diameter~7 km) and S/2004 (45) 1 (Diameter~5 km). The two moons orbit at 1165 and 610 km from the primary, describing an almost-circular orbit (e~6x10-3 and e~7x10-2 respectively). The system is quite different from the other known triple systems in the main belt since the inclinations of the moon orbits are sizeable (9 deg and 18 deg with respect to the equator of the primary respectively). No resonances, neither secular nor due to Lidov-Kozai mechanism, were detected in our dynamical solution, suggesting that these inclinations are not due to excitation modes between the primary and the moons. A 10-year evolution study shows that the orbits are slightly affected by perturbations from the Sun, and to a lesser extent by mutual interactions between the moons. The estimated J2 of the primary is three times lower than the theoretical one, calculated assuming the shape of the primary and an homogeneous interior, possibly suggesting the importance of other gravitational harmonics.