We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

The numerical modeling of co-existing circumplanetary disks/rings and satellites is particularly challenging because each part of the system requires a very different approach. Disks are generally well represented by a fluid-like dense medium, whose evolution can be calculated by a hydrocode. On the other hand, the orbital evolution of satellites is generally performed using N-body integrators. We have developed a new numerical model that combines a 1-dimensional hydrocode with the N-body integrator SyMBA. The disk evolves due to its viscosity, and resonant torques from satellites. The latter is applied to the satellites as an additional "kick" to their accelerations. The integrator also includes the ability to spawn new moonlets at the disk's outer edge if the latter expands beyond a material-dependent Roche limit, as well as the effects of tidal dissipation in the planet and/or the satellite on the satellite orbits. The resulting integrator allows one to accurately model the evolution of an inner circumplanetary disk, and the formation of satellites by accumulation of disk material, all within a single self-consistent framework. Potential applications include the formation of Earth's Moon, the evolution of the inner Saturn system, the martian and uranian moons, and compact exoplanetary systems.

Saturn's rings are rock-poor, containing 90%-95% ice by mass. As a group, Saturn's moons interior to and including Tethys are also about 90% ice. Tethys itself contains <6% rock by mass, in contrast to its similar-mass outer neighbor Dione, which contains >40% rock. Here we simulate the evolution of a massive primordial ice-rich ring and the production of satellites as ring material spreads beyond the Roche limit. We describe the Roche-interior ring with an analytic model, and use an N-body code to describe material beyond the Roche limit. We track the accretion and interactions of spawned satellites, including tidal interaction with the planet, assuming a tidal dissipation factor for Saturn of $Q~10^4$. We find that ring torques and capture of moons into mutual resonances produce a system of ice-rich inner moons that extends outward to approximately Tethys's orbit in $10^9$ years, even with relatively slow orbital expansion due to tides. The resulting mass and semimajor axis distribution of spawned moons resembles that of Mimas, Enceladus, and Tethys. We estimate the mass of rock delivered to the moons by external cometary impactors during a late heavy bombardment. We find that the inner moons receive a mass in rock comparable to their current total rock content, while Dione and Rhea receive an order-of-magnitude less rock than their current rock content. This suggests that external ontamination may have been the primary source of rock in the inner moons, and that Dione and Rhea formed from much more rock-rich source material. Reproducing the distribution of rock among the current inner moons is challenging, and appears to require large impactors stochasticity and/or the presence of some rock in the initial ring.

Impacts that leave the Earth-Moon system with a large excess in angular momentum have recently been advocated as a means of generating a protolunar disc with a composition that is nearly identical to that of the Earth's mantle. We here investigate the accretion of the Moon from discs generated by such "non-canonical" impacts, which are typically more compact than discs produced by canonical impacts and have a higher fraction of their mass initially located inside the Roche limit. Our model predicts a similar overall accretional history for both canonical and non-canonical discs, with the Moon forming in three consecutive steps over hundreds of years. However, we find that, to yield a lunar-mass Moon, the more compact non-canonical discs must initially be more massive than implied by prior estimates, and only a few of the discs produced by impact simulations to date appear to meet this condition. Non-canonical impacts require that capture of the Moon into the evection resonance with the Sun reduced the Earth-Moon angular momentum by a factor of 2 or more. We find that the Moon's semimajor axis at the end of its accretion is approximately $7R_\oplus$, which is comparable to the location of the evection resonance for a post-impact Earth with a 2.5 h rotation period in the absence of a disc. Thus, the dynamics of the Moon's assembly may directly affect its ability to be captured into the resonance.

We use a hybrid numerical approach to simulate the formation of the Moon from an impact-generated disk, consisting of a fluid model for the disk inside the Roche limit and an N-body code to describe accretion outside the Roche limit. As the inner disk spreads due to a thermally regulated viscosity, material is delivered across the Roche limit and accretes into moonlets that are added to the N-body simulation. Contrary to an accretion timescale of a few months obtained with prior pure N-body codes, here the final stage of the Moon's growth is controlled by the slow spreading of the inner disk, resulting in a total lunar accretion timescale of ~10^2 years. It has been proposed that the inner disk may compositionally equilibrate with the Earth through diffusive mixing, which offers a potential explanation for the identical oxygen isotope compositions of the Earth and Moon. However, the mass fraction of the final Moon that is derived from the inner disk is limited by resonant torques between the disk and exterior growing moons. For initial disks containing < 2.5 lunar masses (ML), we find that a final Moon with mass > 0.8ML contains < 60% material derived from the inner disk, with this material preferentially delivered to the Moon at the end of its accretion.

The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis).

Hundred meter sized objects have been identified by the Cassini spacecraft in Saturn's A ring through the so-called "propeller" features they create in the ring. These moonlets should migrate, due to their gravitational interaction with the ring ; in fact, some orbital variation have been detected. The standard theory of type I migration of planets in protoplanetary disks can't be applied to the ring system, as it is pressureless. Thus, we compute the differential torque felt by a moonlet embedded in a two-dimensional disk of solid particles, with flat surface density profile, both analytically and numerically. We find that the corresponding migration rate is too small to explain the observed variations of the propeller's orbit in Saturn's A-ring. However, local density fluctuations (due to gravity wakes in the marginally gravitationally stable A-ring) may exert a stochastic torque on a moonlet. Our simulations show that this torque can be large enough to account for the observations, depending on the parameters of the rings. We find that on time scales of several years the migration of propellers is likely to be dominated by stochastic effects (while the former, non-stochastic migration dominates after ~ 10^{4-5} years). In that case, the migration rates provided by observations so far suggests that the surface density of the A ring should be of the order of 700 kg/m^2. The age of the propellers shouldn't exceed 1 to 100 million years, depending on the dominant migration regime.

We investigate the long-term and large-scale viscous evolution of dense planetary rings using a simple 1D numerical code. We use a physically realistic viscosity model derived from N-body simulations (Daisaka et al., 2001), and dependent on the disk's local properties (surface mass density, particle size, distance to the planet). Particularly, we include the effects of gravitational instabilities (wakes) that importantly enhance the disk's viscosity. We show that common estimates of the disk's spreading time-scales with constant viscosity significantly underestimate the rings' lifetime. With a realistic viscosity model, an initially narrow ring undergoes two successive evolutionary stages: (1) a transient rapid spreading when the disk is self-gravitating, with the formation of a density peak inward and an outer region marginally gravitationally stable, and with an emptying time-scale proportional to 1/M_0^2 (where M_0 is the disk's initial mass) (2) an asymptotic regime where the spreading rate continuously slows down as larger parts of the disk become not-self-gravitating due to the decrease of the surface density, until the disk becomes completely not-self-gravitating. At this point its evolution dramatically slows down, with an emptying time-scale proportional to 1/M_0, which significantly increases the disk's lifetime compared to the case with constant viscosity. We show also that the disk's width scales like t^{1/4} with the realistic viscosity model, while it scales like t^{1/2} in the case of constant viscosity, resulting in much larger evolutionary time-scales in our model. We find however that the present shape of Saturn's rings looks like a 100 million-years old disk in our simulations. Concerning Jupiter's, Uranus' and Neptune's rings that are faint today, it is not likely that they were much more massive in the past and lost most of their mass due to viscous spreading alone.

The origin of Saturn\' s massive ring system is still unknown. Two popular scenarios - the tidal splitting of passing comets and the collisional destruction of a satellite - rely on a high cometary flux in the past. In the present paper we attempt to quantify the cometary flux during the Late Heavy Bombardment (LHB) to assess the likelihood of both scenarios. Our analysis relies on the so-called Nice model of the origin of the LHB (Tsiganis et al., 2005; Morbidelli et al., 2005; Gomes et al., 2005) and on the size distribution of the primordial trans-Neptunian planetesimals constrained in Charnoz & Morbidelli (2007). We find that the cometary flux on Saturn during the LHB was so high that both scenarios for the formation of Saturn rings are viable in principle. However, a more detailed study shows that the comet tidal disruption scenario implies that all four giant planets should have comparable ring systems whereas the destroyed satellite scenario would work only for Saturn, and perhaps Jupiter. This is because in Saturn\'s system, the synchronous orbit is interior to the Roche Limit, which is a necessary condition for maintaining a satellite in the Roche zone up to the time of the LHB. We also discuss the apparent elimination of silicates from the ring parent body implied by the purity of the ice in Saturn \' s rings. The LHB has also strong implications for the survival of the Saturnian satellites: all satellites smaller than Mimas would have been destroyed during the LHB, whereas Enceladus would have had from 40% to 70% chance of survival depending on the disruption model. In conclusion, these results suggest that the LHB is the sweet moment for the formation of a massive ring system around Saturn.