Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

The presence of rings around a transiting planet can cause its radius to be overestimated and lead to an underestimation of its density if the mass is known. We employ a Bayesian framework to show that the anomalously low density ($\sim$0.09 g cm${^{-3}}$) of the transiting long-period planet HIP$\,$41378$\,f$ might be due to the presence of opaque circum-planetary rings. Given our adopted model priors and data from the K2 mission, we find the statistical evidence for the ringed planet scenario to be comparable to that of the planet-only scenario. The ringed planet solution suggests a larger planetary density of $\sim$1.23$\,$g$\,$cm$^{-3}$ similar to Uranus. The associated ring extends from 1.05 to 2.59 times the planetary radius and is inclined away from the sky-plane by $\sim$25$^\mathrm{o}$. Future high-precision transit observations of HIP$\,$41378$\,f$ would be necessary to confirm/dismiss the presence of planetary rings.

Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown. We perform non-ideal MHD numerical simulations using the adaptive mesh refinement code Ramses, of a collapsing, one solar mass, molecular core to study the disc formation and early, up to 100 kyr, evolution, paying great attention to the impact of numerical resolution and accretion scheme. We found that while the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle, the disc mass, and to a lower extent its size, heavily depend on the accretion scheme, which we found, is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetized. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetized but its immediate surrounding is on the contrary highly magnetized. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object, suggesting that the disc mass is eventually controlled by the small scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are (abridged).

The isotopic heterogeneity of the Solar System shown by meteorite analyses is more pronounced for its earliest objects, the Calcium-Aluminum-rich Inclusions (CAIs). This suggests that it was inherited from spatial variations in different stardust populations in the protosolar cloud. We model the formation of the solar protoplanetary disk following its collapse and find that the solid-weighted standard deviation of different nucleosynthetic contributions in the disk is reduced by one order of magnitude compared to the protosolar cloud, whose successive isotopic signatures are fossilized by CAIs. The enrichment of carbonaceous chondrites in r-process components, whose proportions are inferred to have diminished near the end of infall, is consistent with their formation at large heliocentric distances, where the early signatures would have been preferentially preserved after outward advection. We also argue that thermal processing had little effect on the (mass-independent) isotopic composition of bulk meteorites for refractory elements.

The short-lived radionuclide $^{26}$Al is widely used to determine the relative ages of chondrite components and timescales of physical and thermal events that attended the formation of the Solar System. However, an important assumption for using $^{26}$Al as a chronometer is its homogeneous distribution in the disk. Yet, the oldest components in chondrites, the Ca-Al-rich inclusions (CAIs), which are usually considered as time anchors for this chronometer, show evidence of $^{26}$Al/$^{27}$Al variations independent of radioactive decay. Since their formation epoch may have been contemporaneous with the collapse of the parent cloud that formed the disk, this suggests that $^{26}$Al was heteregeneously distributed in the cloud. We model the collapse of such an heterogeneous cloud, using two different $^{26}$Al distributions (monotonic and non-monotonic), and follow its re-distribution in the first condensates and bulk dust that populate the forming disk. We find that CAIs inherit the $^{26}$Al/$^{27}$Al ratio of the matter infalling at the time of their formation, so that variations of $^{26}$Al/$^{27}$Al among primordial CAIs can be accounted for, independently of radioactive decay. The prevalence of a canonical ratio among them and its necessity for the differentiation of the first planetesimals suggest a (monotonic) scenario where $^{26}$Al sharply rose relatively close to the center of the protosolar cloud and essentially remained at a high level outward (rather than decreased since). As the $ ^{26}$Al abundance would be relatively homogeneous after cessation of infall, this would warrant the use of the Al-Mg chronometer from the formation of "regular" CAIs onward, to chondrules and chondrite accretion.

Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne observations of the Martian moons and of the more modern theories proposed to account for the formation of other moon systems. The most recent scenario of formation of the Martian moons relies on a giant impact occurring at early Mars history and having also formed the so-called hemispheric crustal dichotomy. This scenario accounts for the current orbits of both moons unlike the scenario of capture of small size asteroids. It also predicts a composition of disk material as a mixture of Mars and impactor materials that is in agreement with remote sensing observations of both moon surfaces, which suggests a composition different from Mars. The composition of the Martian moons is however unclear, given the ambiguity on the interpretation of the remote sensing observations. The study of the formation of the Martian moon system has improved our understanding of moon formation of terrestrial planets: The giant collision scenario can have various outcomes and not only a big moon as for the Earth. This scenario finds a natural place in our current vision of the early solar system when conditions were favorable for giant collisions to occur. The next step in exploration of Martian moon is a sample return mission to test the giant collision scenario for their origin, and to provide tests of models of early solar system dynamics since Mars may retain material exchanged between the inner and outer solar system.

It is important to clarify where and when rocky and icy planetesimals are formed in a viscously evolving disk. We wish to understand how local runaway pile-up of solids occurs inside or outside the snow line. We assume an icy pebble contains micron-sized silicate grains that are uniformly mixed with ice and are released during the ice sublimation. Using a local one-dimensional code, we solve the radial drift and the turbulent diffusion of solids and the water vapor, taking account of their sublimation/condensation around the snow line. We systematically investigate effects of back-reactions of the solids to gas on the radial drift and diffusion of solids, scale height evolution of the released silicate particles, and possible difference in effective viscous parameters between that for turbulent diffusion ($\alpha_{\rm tur}$) and that for the gas accretion rate onto the central star ($\alpha_{\rm acc}$). We study the dependence on the ratio of the solid mass flux to the gas ($F_{\rm p/g}$). We show that the favorable locations for the pile-up of silicate grains and icy pebbles are the regions in the proximity of the water snow line inside and outside it, respectively. We found that runaway pile-ups occur when both the back-reactions for radial drift and diffusion are included. In the case with only the back-reaction for the radial drift, no runaway pile-up is found except for extremely high pebble flux, while the condition of streaming instability can be satisfied for relatively large $F_{\rm p/g}$ as found in the past literatures. If the back-reactions for radial diffusion is considered, the runaway pile-up occurs for reasonable value of pebble flux. The runaway pile-up of silicate grains that would lead to formation of rocky planetesimals occurs for $\alpha_{\rm tur} \ll \alpha_{\rm acc}$, while the runaway pile-up of icy pebbles is favored for $\alpha_{\rm tur} \sim \alpha_{\rm acc}$.

Chondritic meteorites, the building blocks of terrestrial planets, are made of an out-of-equilibrium assemblage of solids formed at high and low temperatures, either in our Solar system or previous generations of stars. This was considered for decades to result from large scale transport processes in the Sun's isolated accretion disk. However, mounting evidences suggest that refractory inclusions in chondrites formed contemporaneously with the disk building. Here we numerically investigate, using a 1D model and several physical and chemical processes, the formation and transport of rocky materials during the collapse of the Sun's parent cloud and the consequent Solar Nebula assembling. The interplay between the cloud collapse, the dynamics of gas and dust, vaporization, recondensation and thermal processing of different species in the disk, results in a local mixing of solids with different thermal histories. Moreover, our results also explains the overabundance of refractory materials far from the Sun and their short formation timescales, during the first tens of kyr of the Sun, corresponding to class 0-I, opening new windows into the origin of the compositional diversity of chondrites.

Recent works have shown that Martian moons Phobos and Deimos may have accreted within a giant impact-generated disk whose composition is about an equal mixture of Martian material and impactor material. Just after the giant impact, the Martian surface is heated up to $\sim3000-6000$ K and the building blocks of moons, including volatile-rich vapor, are heated up to $\sim2000$ K. In this paper, we investigate the volatile loss from the building blocks of Phobos and Deimos by hydrodynamic escape of vapor and radiation pressure on condensed particles. We show that a non-negligible amount of volatiles ($> 10\%$ of the vapor with temperature $> 1000$ K via hydrodynamic escape, and moderately volatile dusts that condense at $\sim700-2000$ K via radiation pressure) could be removed just after the impact during their first signle orbit from their pericenters to apocenters. Our results indicate that bulk Phobos and Deimos are depleted in volatile elements. Together with future explorations such as JAXA's MMX (Martian Moons eXploration) mission, our results would be used to constrain the origin of Phobos and Deimos.

Rings are ubiquitous around giant planets in our Solar System. They evolve jointly with the nearby satellite system. They could form either during the giant planet formation process or much later, as a result of large scale dynamical instabilities either in the local satellite system, or at the planetary scale. We review here the main characteristics of rings in our solar system, and discuss their main evolution processes and possible origin. We also discuss the recent discovery of rings around small bodies.

Several tens of white dwarfs are known to host circumstellar discs of dusty debris, thought to arise from the tidal disruption of rocky bodies originating in the star's remnant planetary system. This paper investigates the evolution of such discs if they are very massive, as may be the case if their progenitor was a terrestrial planet, moon, or dwarf planet. Assuming the discs are physically thin and flat, like Saturn's rings, their evolution is governed by Poynting-Robertson drag or viscous spreading, where the disc's effective viscosity is due to self-gravity wakes. For discs with masses >10^26 g, located in the outer parts of the tidal disruption zone, viscous spreading dominates the evolution, and mass is transported both in- and outwards. When outwards-spreading material flows beyond the Roche limit, it coagulates into new (minor) planets in a process analogous to the ongoing formation of moonlets at the outer edge of Saturn's rings. The newly formed bodies migrate outwards by exchanging angular momentum with the disc and coalesce into larger objects through mutual collisions. Eventually, the disc's Roche-limit overflow recycles tens of percent of the original disc mass; most ends up in a single large body near 2:1 mean-motion resonance with the disc's outer edge. Hence, the recycling of a tidally disrupted super-Earth, for example, could yield an Earth-mass planet on a ~10-h orbit, located in the habitable zone for 2-to-10-Gyr-old white dwarfs. The recycling process also creates a population of smaller bodies just outside the Roche limit, which may explain the minor planets recently postulated to orbit WD 1145+017.

The origin of Phobos and Deimos in a giant impact generated disk is gaining larger attention. Although this scenario has been the subject of many studies, an evaluation of the chemical composition of the Mars' moons in this framework is missing. The chemical composition of Phobos and Deimos is unconstrained. The large uncertainty about the origin of the mid-infrared features, the lack of absorption bands in the visible and near-infrared spectra, and the effects of secondary processes on the moons' surface make the determination of their composition very difficult from remote sensing data. Simulations suggest a formation of a disk made of gas and melt with their composition linked to the nature of the impactor and Mars. Using thermodynamic equilibrium we investigate the composition of dust (condensates from gas) and solids (from a cooling melt) that result from different types of Mars impactors (Mars-, CI-, CV-, EH-, comet-like). Our calculations show a wide range of possible chemical compositions and noticeable differences between dust and solids depending on the considered impactors. Assuming Phobos and Deimos as result of the accretion and mixing of dust and solids, we find that the derived assemblage (dust rich in metallic-iron, sulphides and/or carbon, and quenched solids rich in silicates) can be compatible with the observations. The JAXA's MMX (Martian Moons eXploration) mission will investigate the physical and chemical properties of the Maroons, especially sampling from Phobos, before returning to Earth. Our results could be then used to disentangle the origin and chemical composition of the pristine body that hit Mars and suggest guidelines for helping in the analysis of the returned samples.

Phobos and Deimos are the two small Martian moons, orbiting almost on the equatorial plane of Mars. Recent works have shown that they can accrete within an impact-generated inner dense and outer light disk, and that the same impact potentially forms the Borealis basin, a large northern hemisphere basin on the current Mars. However, there is no a priori reason for the impact to take place close to the north pole (Borealis present location) nor to generate a debris disk in the equatorial plane of Mars (in which Phobos and Deimos orbit). In this paper, we investigate these remaining issues on the giant impact origin of the Martian moons. First, we show that the mass deficit created by the Borealis impact basin induces a global reorientation of the planet to realign its main moment of inertia with the rotation pole (True Polar Wander). This moves the location of the Borealis basin toward its current location. Next, using analytical arguments, we investigate the detailed dynamical evolution of the eccentric inclined disk from the equatorial plane of Mars that is formed by the Martian-moon-forming impact. We find that, as a result of precession of disk particles due to the Martian dynamical flattening $J_{2}$ term of its gravity field and particle-particle inelastic collisions, eccentricity and inclination are damped and an inner dense and outer light equatorial circular disk is eventually formed. Our results strengthen the giant impact origin of Phobos and Deimos that can finally be tested by a future sample return mission such as JAXA's Martian Moons eXploration (MMX) mission.

The ice giants Uranus and Neptune are the least understood class of planets in our solar system but the most frequently observed type of exoplanets. Presumed to have a small rocky core, a deep interior comprising ~70% heavy elements surrounded by a more dilute outer envelope of H2 and He, Uranus and Neptune are fundamentally different from the better-explored gas giants Jupiter and Saturn. Because of the lack of dedicated exploration missions, our knowledge of the composition and atmospheric processes of these distant worlds is primarily derived from remote sensing from Earth-based observatories and space telescopes. As a result, Uranus's and Neptune's physical and atmospheric properties remain poorly constrained and their roles in the evolution of the Solar System not well understood. Exploration of an ice giant system is therefore a high-priority science objective as these systems (including the magnetosphere, satellites, rings, atmosphere, and interior) challenge our understanding of planetary formation and evolution. Here we describe the main scientific goals to be addressed by a future in situ exploration of an ice giant. An atmospheric entry probe targeting the 10-bar level, about 5 scale heights beneath the tropopause, would yield insight into two broad themes: i) the formation history of the ice giants and, in a broader extent, that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. In addition, possible mission concepts and partnerships are presented, and a strawman ice-giant probe payload is described. An ice-giant atmospheric probe could represent a significant ESA contribution to a future NASA ice-giant flagship mission.

Phobos and Deimos are the two small moons of Mars. Recent works have shown that they can accrete within an impact-generated disk. However, the detailed structure and initial thermodynamic properties of the disk are poorly understood. In this paper, we perform high-resolution SPH simulations of the Martian moon-forming giant impact that can also form the Borealis basin. This giant impact heats up the disk material (around $\sim 2000$ K in temperature) with an entropy increase of $\sim 1500$ J K$^{-1}$ kg$^{-1}$. Thus, the disk material should be mostly molten, though a tiny fraction of disk material ($< 5\%$) would even experience vaporization. Typically, a piece of molten disk material is estimated to be meter sized due to the fragmentation regulated by their shear velocity and surface tension during the impact process. The disk materials initially have highly eccentric orbits ($e \sim 0.6-0.9$) and successive collisions between meter-sized fragments at high impact velocity ($\sim 3-5$ km s$^{-1}$) can grind them down to $\sim100 \mu$m-sized particles. On the other hand, a tiny amount of vaporized disk material condenses into $\sim 0.1 \mu$m-sized grains. Thus, the building blocks of the Martian moons are expected to be a mixture of these different sized particles from meter-sized down to $\sim 100 \mu$m-sized particles and $\sim 0.1 \mu$m-sized grains. Our simulations also suggest that the building blocks of Phobos and Deimos contain both impactor and Martian materials (at least 35%), most of which come from the Martian mantle (50-150 km in depth; at least 50%). Our results will give useful information for planning a future sample return mission to Martian moons, such as JAXA's MMX (Martian Moons eXploration) mission.

The hypothesis of a recent origin of Saturn's rings and its mid-sized moons is actively debated. It was suggested that a proto-Rhea and a proto-Dione might have collided recently, giving birth to the modern system of mid-sized moons. It is also suggested that the rapid viscous spreading of the debris may have implanted mass inside Saturn's Roche limit, giving birth to the modern Saturn's ring system. However, this scenario has been only investigated in very simplified way for the moment. This paper investigates it in detail to assess its plausibility by using $N$-body simulations and analytical arguments. When the debris disk is dominated by its largest remnant, $N$-body simulations show that the system quickly re-accrete into a single satellite without significant spreading. On the other hand, if the disk is composed of small particles, analytical arguments suggest that the disk experiences dynamical evolutions in three steps. The disk starts significantly excited after the impact and collisional damping dominates over the viscous spreading. After the system flattens, the system can become gravitationally unstable when particles are smaller than $\sim$ 100 m. However, the particles grow faster than spreading. Then, the system becomes gravitationally stable again and accretion continues at a slower pace, but spreading is inhibited. Therefore, the debris is expected to re-accrete into several large bodies. In conclusion, our results show that such a scenario may not form the today's ring system. In contrast, our results suggest that today's mid-sized moons are likely re-accreted from such a catastrophic event.

This chapter of the book Planetary Ring Systems addresses the origin of planetary rings, one of the least understood processes related to planet formation and evolution. Whereas rings seem ubiquitous around giant planets, their great diversity of mass, structure and composition is a challenge for any formation scenario. Recent advances in our understanding of ring and satellite formation and destruction suggest that these processes are interconnected, so that rings and satellites may be two aspects of the same geological system. However, no single theory seems able to explain the origin of the different planetary rings known in our Solar System, and it now seems evident that rings may result from a variety of processes like giant collisions, tidal stripping of comets or satellites, as well as planet formation itself. In order to build any theory of ring formation it is important to specify physical processes that affect the long-term evolution of rings, as well as to describe the different observations that any ring formation model should explain. This is the topic of section 2. In section 3, we focus our attention on Saturn's rings and their main properties, and then discuss the pros and cons of a series of ring formation models. We also discuss the link between rings and satellites. In section 4, we extend the discussion to the other giant planets (Jupiter, Uranus, and Neptune). Section 5 is devoted to new types of rings -- the recent discovery of rings orbiting small outer Solar System bodies (Centaurs), and the possible rings around extrasolar planets. In section 6, we conclude and try to identify critical observations and theoretical advances needed to better understand the origin of rings and their significance in the global evolution of planets.

The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulations of the tidal destruction of big differentiated objects ($M_{\rm body}=10^{21-23}$) that experience close encounters with Saturn or Uranus. We find that about $0.1-10$% of the mass of the passing body is gravitationally captured around the planet. However, these fragments are initially big chunks and have highly eccentric orbits around the planet. Then, we perform N-body simulations including the planet's oblateness, starting with data obtained from the SPH simulations. Our N-body simulations show that the chunks are tidally destroyed during their next several orbits. Their individual orbits then start to precess incoherently around the planet's equator, which enhances their encounter velocities on longer-term evolution, resulting in more destructive impacts. These collisions would damp their eccentricities resulting in a progressive collapse of the debris cloud into a thin equatorial and low-eccentricity ring. These high energy impacts are expected to be catastrophic enough to produce small particles. Our numerical results also show that the mass of formed rings is large enough to explain current rings including inner regular satellites around Saturn and Uranus. In the case of Uranus, a body can go deeper inside the planet's Roche limit resulting in a more efficient capture of rocky material compared to Saturn's case in which mostly ice is captured. Thus, our results can naturally explain the compositional difference between the rings of Saturn, Uranus and Neptune.

Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the tidal disruption of a differentiated object that experiences a close encounter with a giant planet could naturally form diverse ring-satellite systems around the Centaurs. During the close encounter, the icy mantle of the passing object is preferentially ripped off by the planet's tidal force and the debris is distributed mostly within the Roche limit of the largest remnant body. Assuming the existence of $20-50$wt% silicate core below the icy mantle, a disk of particles is formed when the objects pass within $0.4-0.8$ of the planet's Roche limit with the relative velocity at infinity $3-6$km s$^{-1}$ and 8h initial spin period of the body. The resultant ring mass is $0.1-10$% of the central object's mass. Such particle disks are expected to spread radially, and materials spreading beyond the Roche limit would accrete into satellite(s). Our numerical results suggest that ring formation would be a natural outcome of such extreme close encounters and Centaurs can naturally have such ring systems because they cross the orbits of the giant planets.

Planet traps are necessary to prevent forming planets from falling onto their host star by type I migration. Surface mass density and temperature gradient irregularities favor the apparition of traps and deserts. Such features are found at the dust sublimation lines and heat transition barriers. We study how planets may remain trapped or escape as they grow and as the disk evolves. We model the temporal viscous evolution of a protoplanetary disk by coupling its dynamics, thermodynamics, geometry and composition. The resulting mid-plane density and temperature profiles allow the modeling of the interactions of such an evolving disk with potential planets, even before the steady state is reached. We follow the viscous evolution of a MMSN and compute the Lindblad and corotation torques that such a disk would exert on potential planets of various masses located within the planetary formation region. We determine the position of planet traps and deserts in relationship with the sublimation lines, shadowed regions and heat transition barriers. Planets that are a few tens of Earth masses can be trapped at the sublimation lines until they reach a certain mass while planets more massive than 100ME can only be trapped permanently at the heat transition barriers. Coupling a bimodal planetary migration model with a self-consistent evolved disk, we were able to distinguish several potential planet populations after 5 million years of evolution: two populations of giant planets that could stay trapped around 5.5 and 9 au and possibly open gaps, some super-Earths trapped around 5 and 7.5 au and a population of close-in super-Earths trapped inside 1 au. The traps corresponding to the last group could help validating the in-situ formation scenarios of the observed close-in super-Earths.

Whereas it is generally accepted that calcium-aluminum-rich inclusions (CAIs) from chondritic meteorites formed in a hot environment in the solar protoplanetary disk, the conditions of their formation remain debated. Recent laboratory studies of CAIs have provided new kind of data: their size distributions. We show that size distributions of CAIs measured in laboratory from sections of carbonaceous chondrites have a power law size distribution with cumulative size exponent between -1.7 and -1.9, which translates into cumulative size exponent between -2.5 and -2.8 after correction for sectioning. To explain these observations, numerical simulations were run to explore the growth of CAIs from micrometer to centimeter sizes, in a hot and turbulent protoplanetary disk through the competition of coagulation and fragmentation. We show that the size distributions obtained in growth simulations are in agreement with CAIs size distributions in meteorites. We explain the CAI sharp cut-off of their size distribution at centimeter sizes as the direct result from the famous fragmentation barrier, provided that CAI fragment for impact velocities larger than 10 m/s. The growth/destruction timescales of millimeter- and centimeter-sized CAIs is inversely proportional to the local dust/gas ratio and is about 10 years at 1300 K and up to 104 years at 1670K. This implies that the most refractory CAIs are expected to be smaller in size owing to their long growth timescale compared to less refractory CAIs. Conversely, the least refractory CAIs could have been recycled many times during the CAI production era which may have profound consequences for their radiometric age.

Aims. We track the time evolution of planet traps and snowlines in a viscously evolving protoplanetary disk using an opacity table that accounts for the composition of the dust material. Methods. We coupled a dynamical and thermodynamical disk model with a temperature-dependent opacity table (that accounts for the sublimation of the main dust components) to investigate the formation and evolution of snowlines and planet traps during the first million years of disk evolution. Results. Starting from a minimum mass solar nebula (MMSN), we find that the disk mid-plane temperature profile shows several plateaux (0.1-1 AU wide) at the different sublimation temperatures of the species that make up the dust. For water ice, the correspond- ing plateau can be larger than 1 AU, which means that this is a snow "region" rather than a snow "line". As a consequence, the surface density of solids in the snow region may increase gradually, not abruptly. Several planet traps and desert regions appear naturally as a result of abrupt local changes in the temperature and density profiles over the disk lifetime. These structures are mostly located at the edges of the temperature plateaux (surrounding the dust sublimation lines) and at the heat-transition barrier where the disk stellar heating and viscous heating are of the same magnitude (around 10 AU after 1 Myr). Conclusions. Several traps are identified: although a few appear to be transient, most of them slowly migrate along with the heat- transition barrier or the dust sublimation lines. These planet traps may temporarily favor the growth of planetary cores.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Violent stochastic collisional events have been invoked as a possible explanation for some debris discs displaying pronounced asymmetries or having a great luminosity excess. So far, no thorough modelling of the consequences of such events has been carried out, mainly because of the extreme numerical challenge of coupling the dynamical and collisional evolution of dust. We perform the first fully self-consistent modelling of the aftermath of massive breakups in debris discs. We follow the collisional and dynamical evolution of dust released after the breakup of a Ceres-sized body at 6 AU from its central star. We investigate the duration, magnitude and spatial structure of the signature left by such a violent event, as well as its observational detectability. We use the recently developed LIDT-DD code (Kral et al., 2013), which handles the coupled collisional and dynamical evolution of debris discs. The main focus is placed on the complex interplay between destructive collisions, Keplerian dynamics and radiation pressure forces. We use the GRaTer package to estimate the system's luminosity at different wavelengths. The breakup of a Ceres-sized body at 6 AU creates an asymmetric dust disc that is homogenized, by the coupled action of collisions and dynamics, on a timescale of a few $10^5$ years. The luminosity excess in the breakup's aftermath should be detectable by mid-IR photometry, from a 30 pc distance, over a period of $\sim 10^6$ years that exceeds the duration of the asymmetric phase of the disc (a few $10^5$ years). As for the asymmetric structures, we derive synthetic images for the SPHERE/VLT and MIRI/JWST instruments, showing that they should be clearly visible and resolved from a 10 pc distance. Images at 1.6$\mu$m (marginally), 11.4 and 15.5$\mu$m would show the inner disc structures while 23$\mu$m images would display the outer disc asymmetries.

Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as $r^{-1}$, or aspect ratios about 0.03 to 0.23. Accretion rates are also recurrently found around $10^{-8} - 10^{-6}~\mathrm{M_{\odot}~yr^{-1}}$ for disks already evolved (Isella et al., 2009, Andrews et al., 2009, 2010). Several models have been developed in order to recover these properties. However, most of them usually simplify the disk geometry if not its mid-plane temperature. This has major consequences for modeling the disk evolution over million years and consequently planet migration. In the present paper, we develop a viscous evolution hydrodynamical numerical code that determines simultaneously the disk photosphere geometry and the mid-plane temperature. We then compare our results of long-term simulations with similar simulations of disks with a constrained geometry along the Chiang & Goldreich (1997) prescription (dlnH/dlnr = 9/7). We find that the constrained geometry models provide a good approximation of the disk surface density evolution. However, they differ significantly regarding the temperature time evolution. In addition, we find that shadowed regions naturally appear at the transition between viscously dominated and radiation dominated regions that falls in the region of planetary formation. We show that $\chi$ (photosphere height to pressure scale height ratio) cannot be considered as a constant, consistently with Watanabe et al. (2008). Comparisons with observations show that all disk naturally evolve toward a shallow surface density disk ($\Sigma \propto r^{-1}$). The mass flux across the disk stabilizes in about 1 million year typically.

In most current debris disc models, the dynamical and the collisional evolutions are studied separately, with N-body and statistical codes, respectively, because of stringent computational constraints. We present here LIDT-DD, the first code able to mix both approaches in a fully self-consistent way. Our aim is for it to be generic enough so as to be applied to any astrophysical cases where we expect dynamics and collisions to be deeply interlocked with one another: planets in discs, violent massive breakups, destabilized planetesimal belts, exozodiacal discs, etc. The code takes its basic architecture from the LIDT3D algorithm developed by Charnoz et al.(2012) for protoplanetary discs, but has been strongly modified and updated in order to handle the very constraining specificities of debris discs physics: high-velocity fragmenting collisions, radiation-pressure affected orbits, absence of gas, etc. In LIDT-DD, grains of a given size at a given location in a disc are grouped into "super-particles", whose orbits are evolved with an N-body code and whose mutual collisions are individually tracked and treated using a particle-in-a-box prescription. To cope with the wide range of possible dynamics, tracers are sorted and regrouped into dynamical families depending on their orbits. The code retrieves the classical features known for debris discs, such as the particle size distributions in unperturbed discs, the outer radial density profiles (slope in -1.5) outside narrow collisionally active rings, and the depletion of small grains in "dynamically cold" discs. The potential of the new code is illustrated with the test case of the violent breakup of a massive planetesimal within a debris disc. The main potential future applications of the code are planet/disc interactions, and more generally any configurations where dynamics and collisions are expected to be intricately connected.

The first attempt at developing a fully self-consistent code coupling dynamics and collisions to study debris discs (Kral, Thebault, Charnoz, 2013) is presented. So far, these two crucial mechanisms were studied separately, with N-body and statistical collisional codes respectively, because of stringent computational constraints. In particular, incorporating collisional effects (especially destructive collisions) into an N-body scheme was deemed an impossible task because of the exponential increase of particles it would imply. We present here a new model named LIDT-DD which is able to follow over long timescales the coupled evolution of dynamics (including radiation pressure) and collisions in a self-consistent way.

When a planetary tidal disk -like Saturn's rings- spreads beyond the Roche radius (inside which planetary tides prevent aggregation), satellites form and migrate away. Here, we show that most regular satellites in the solar system probably formed in this way. According to our analytical model, when the spreading is slow, a retinue of satellites appear with masses increasing with distance to the Roche radius, in excellent agreement with Saturn's, Uranus', and Neptune's satellite systems. This suggests that Uranus and Neptune used to have massive rings that disappeared to give birth to most of their regular satellites. When the spreading is fast, only one large satellite forms, as was the case for Pluto and Earth. This conceptually bridges the gap between terrestrial and giant planet systems.

Dust is a major component of protoplanetary and debris disks as it is the main observable signature of planetary formation. However, since dust dynamics is size-dependent (because of gas-drag or radiation pressure) any attempt to understand the full dynamical evolution of circumstellar dusty-disks that neglect the coupling of collisional evolution with dynamical evolution is thwarted because of the feedback between these two processes. Here, a new hybrid lagrangian/eulerian code is presented that overcomes some of these difficulties. The particles representing "dust-clouds" are tracked individually in a lagrangian way. This system is then mapped on an eulerian spatial grid, inside the cells of which the local collisional evolutions are computed. Finally, the system is remapped back in a collection of discrete lagrangian particles keeping constant their number. An application example on dust growth in a turbulent protoplanetary disk at 1 AU is presented. First the growth of dust is considered in the absence of a dead-zone and the vertical distribution of dust is self-consistently computed. It is found that the mass is rapidly dominated by particles about a fraction of millimeter in size. Then the same case with an embedded dead-zone is investigated and It is found that coagulation is much more efficient and produces, in a short time scale, 1cm-10cm dust pebbles that dominate the mass. These pebbles may then be accumulated into embryos sized objects inside large-scale turbulent structures as shown recently (see e.g. Johansen et al., 2007).

Tidal interactions between Saturn and its satellites play a crucial role in both the orbital migration of the satellites and the heating of their interiors. Therefore constraining the tidal dissipation of Saturn (here the ratio k2/Q) opens the door to the past evolution of the whole system. If Saturn's tidal ratio can be determined at different frequencies, it may also be possible to constrain the giant planet's interior structure, which is still uncertain. Here, we try to determine Saturn's tidal ratio through its current effect on the orbits of the main moons, using astrometric data spanning more than a century. We find an intense tidal dissipation (k2/Q= (2.3 \pm 0.7) \times 10-4), which is about ten times higher than the usual value estimated from theoretical arguments. As a consequence, eccentricity equilibrium for Enceladus can now account for the huge heat emitted from Enceladus' south pole. Moreover, the measured k2/Q is found to be poorly sensitive to the tidal frequency, on the short frequency interval considered. This suggests that Saturn's dissipation may not be controlled by turbulent friction in the fluid envelope as commonly believed. If correct, the large tidal expansion of the moon orbits due to this strong Saturnian dissipation would be inconsistent with the moon formations 4.5 Byr ago above the synchronous orbit in the Saturnian subnebulae. But it would be compatible with a new model of satellite formation in which the Saturnian satellites formed possibly over longer time scale at the outer edge of the main rings. In an attempt to take into account for possible significant torques exerted by the rings on Mimas, we fitted a constant rate da/dt on Mimas semi-major axis, also. We obtained an unexpected large acceleration related to a negative value of da/dt= -(15.7 \pm 4.4) \times 10-15 au/day.

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).

Hot temperature minerals have been detected in a large number of comets and were also identified in the samples of Comet Wild 2 that were returned by the Stardust mission. Meanwhile, observations of the distribution of hot minerals in young stellar systems suggest that these materials were produced in the inner part of the primordial nebula and have been transported outward in the formation zone of comets. We investigate the possibility that photophoresis provides a viable mechanism to transport high-temperature materials from the inner solar system to the regions in which the comets were forming. We use a grid of time-dependent disk models of the solar nebula to quantify the distance range at which hot minerals can be transported from the inner part of the disk toward its outer regions as a function of their size and density. The particles considered here are in the form of aggregates that presumably were assembled from hot mineral individual grains ranging down to submicron sizes and formed by condensation within the hottest portion of the solar nebula. Our particle-transport model includes the photophoresis, radiation pressure, and gas drag. Depending on the postulated disk parameters and the density of particles, 10-2 to 10-1 m aggregates can reach heliocentric distances up to 35 AU in the primordial nebula over very short timescales (no more than a few hundred thousand years). 10-3 m particles follow the same trajectory as the larger ones but their maximum migration distance does not exceed 26 AU and is reached at later epochs in the disks. On the other hand, 10-5 to 10-4 m aggregates are continuously pushed outward during the evolution of the solar nebula. Our simulations suggest that irrespective of the employed solar nebula model, photophoresis is a mechanism that can explain the presence of hot temperature minerals in the formation region of comets.

In order to understand how the chemical and isotopic compositions of dust grains in a gaseous turbulent protoplanetary disk are altered during their journey in the disk, it is important to determine their individual trajectories. We study here the dust-diffusive transport using lagrangian numerical simulations using the the popular "turbulent diffusion" formalism. However it is naturally expressed in an Eulerian form, which does not allow the trajectories of individual particles to be studied. We present a simple stochastic and physically justified procedure for modeling turbulent diffusion in a Lagrangian form that overcomes these difficulties. We show that a net diffusive flux F of the dust appears and that it is proportional to the gas density ({\rho}) gradient and the dust diffusion coefficient Dd: (F=Dd/{\rho}\timesgrad({\rho})). It induces an inward transport of dust in the disk's midplane, while favoring outward transport in the disk's upper layers. We present tests and applications comparing dust diffusion in the midplane and upper layers as well as sample trajectories of particles with different sizes. We also discuss potential applications for cosmochemistry and SPH codes.

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 and long-term evolution of Saturn's rings is still an unsolved problem in modern planetary science. In this chapter we review the current state of our knowledge on this long-standing question for the main rings (A, Cassini Division, B, C), the F Ring, and the diffuse rings (E and G). During the Voyager era, models of evolutionary processes affecting the rings on long time scales (erosion, viscous spreading, accretion, ballistic transport, etc.) had suggested that Saturn's rings are not older than 100 My. In addition, Saturn's large system of diffuse rings has been thought to be the result of material loss from one or more of Saturn's satellites. In the Cassini era, high spatial and spectral resolution data have allowed progress to be made on some of these questions. Discoveries such as the ''propellers'' in the A ring, the shape of ring-embedded moonlets, the clumps in the F Ring, and Enceladus' plume provide new constraints on evolutionary processes in Saturn's rings. At the same time, advances in numerical simulations over the last 20 years have opened the way to realistic models of the rings's fine scale structure, and progress in our understanding of the formation of the Solar System provides a better-defined historical context in which to understand ring formation. All these elements have important implications for the origin and long-term evolution of Saturn's rings. They strengthen the idea that Saturn's rings are very dynamical and rapidly evolving, while new arguments suggest that the rings could be older than previously believed, provided that they are regularly renewed. Key evolutionary processes, timescales and possible scenarios for the rings's origin are reviewed in the light of t

For its beautiful rings, active atmosphere and mysterious magnetic field, Saturn is a fascinating planet. It also holds some of the keys to understanding the formation of our Solar System and the evolution of giant planets in general. While the exploration by the Cassini-Huygens mission has led to great advances in our understanding of the planet and its moons, it has left us with puzzling questions: What is the bulk composition of the planet? Does it have a helium core? Is it enriched in noble gases like Jupiter? What powers and controls its gigantic storms? We have learned that we can measure an outer magnetic field that is filtered from its non-axisymmetric components, but what is Saturn's inner magnetic field? What are the rings made of and when were they formed? These questions are crucial in several ways: a detailed comparison of the compositions of Jupiter and Saturn is necessary to understand processes at work during the formation of these two planets and of the Solar System. This calls for the continued exploration of the second largest planet in our Solar System, with a variety of means including remote observations and space missions. Measurements of gravity and magnetic fields very close to the planet's cloud tops would be extremely valuable. Very high spatial resolution images of the rings would provide details on their structure and the material that form them. Last but not least, one or several probes sent into the atmosphere of the planet would provide the critical measurements that would allow a detailed comparison with the same measurements at Jupiter. [abridged abstract]

Here we show that Iapetus can serve to discriminate between satellite formation models. Its accretion history can be understood in terms of a two-component gaseous subnebula, with a relatively dense inner region, and an extended tail out to the location of the irregular satellites, as in the SEMM model of Mosqueira and Estrada (2003a,b). Following giant planet formation, planetesimals in the feeding zone of Jupiter and Saturn become dynamically excited, and undergo a collisional cascade. Ablation and capture of planetesimal fragments crossing the gaseous circumplanetary disks delivers enough collisional rubble to account for the mass budgets of the regular satellites of Jupiter and Saturn. This process can result in rock/ice fractionation provided the make up of the population of disk crossers is non-homogeneous, thus offering a natural explanation for the marked compositional differences between outer solar nebula objects and those that accreted in the subnebulae of the giant planets. Consequently, our model leads to an enhancement of the ice content of Iapetus, and to a lesser degree those of Ganymede, Titan and Callisto, and accounts for the (non-stochastic) compositions of these large, low-porosity outer regular satellites of Jupiter and Saturn. (abridged)

In this paper, we characterize the morphology of the disk-integrated phase functions of satellites and rings around the giant planets of our Solar System. We find that the shape of the phase function is accurately represented by a logarithmic model (Bobrov, 1970, in Surfaces and Interiors of Planets and Satellites, Academic, edited by A. Dollfus). For practical purposes, we also parametrize the phase curves by a linear-exponential model (Kaasalainen et al., 2001, Journal of Quantitative Spectroscopy and Radiative Transfer, 70, 529-543) and a simple linear-by-parts model (Lumme and Irvine, 1976, Astronomical Journal, 81, 865-893), which provides three morphological parameters : the amplitude A and the Half-Width at Half-Maximum (HWHM) of the opposition surge, and the slope S of the linear part of the phase function at larger phase angles. Our analysis demonstrates that all of these morphological parameters are correlated with the single scattering albedos of the surfaces. By taking more accurately into consideration the finite angular size of the Sun, we find that the Galilean, Saturnian, Uranian and Neptunian satellites have similar HWHMs (0.5 degrees), whereas they have a wide range of amplitudes A. The Moon has the largest HWHM (2 degrees). We interpret that as a consequence of the solar size bias, via the finite size of the Sun which varies dramatically from the Earth to Neptune. By applying a new method that attempts to morphologically deconvolve the phase function to the solar angular size, we find that icy and young surfaces, with active resurfacing, have the smallest values of A and HWHM, whereas dark objects (and perhaps older surfaces) such as the Moon, Nereid and Saturn C ring have the largest A and HWHM.

The Cassini cameras have captured the opposition effect in Saturn rings with a high radial resolution at phase angles down to 0.01o in the entire main ring system. We derive phase functions from 0.01 to 25 degrees of phase angle in the Wide-Angle Camera clear filters with a central wavelength of 0.611microns and phase functions from 0.001 to 25 degrees of phase angle in the Narrow-Angle and Wide-Angle Cameras color filters (from the blue 0.451 microns to the near infrared 0.752 microns). We characterize the morphology of the phase functions of different features in the main rings. We find that the shape of the phase function is accurately represented by a logarithmic model (Bobrov 1970, in Surfaces and Interiors of Planets and Satellites, Academic, edited by A. Dollfus). For practical purposes, we also parametrize the phase curves by a simple linear-by-part model (Lumme and Irvine 1976, Astronomical Journal, 81, p865), which provides three morphological parameters : the amplitude and the Half-Width at Half-Maximum (HWHM) of the surge, and the slope S of the linear-part of the phase function at larger phase angles. Our analysis demonstrates that all of these parameters show trends with the optical depth of the rings. These trends imply that the optical depth is a key-element determining the physical properties which act on the opposition effect. Wavelength variations of the morphological parameters of the surge show important trends with the optical depth in the green filter (0.568 microns), which implies that grain size effects are maximum in this wavelength.

Since 2004, observations of Saturn's F ring have revealed that the ring's core is surrounded by structures with radial scales of hundreds of kilometers, called "spirals" and "jets". Gravitational scattering by nearby moons was suggested as a potential production mechanism; however, it remained doubtful because a population of Prometheus-mass moons is needed and, obviously, such a population does not exist in the F ring region. We investigate here another mechanism: dissipative physical collisions of kilometer-size moonlets (or clumps) with the F-ring core. We show that it is a viable and efficient mechanism for producing spirals and jets, provided that massive moonlets are embedded in the F-ring core and that they are impacted by loose clumps orbiting in the F ring region, which could be consistent with recent data from ISS, VIMS and UVIS. We show also that coefficients of restitution as low as ~0.1 are needed to reproduce the radial extent of spirals and jets, suggesting that collisions are very dissipative in the F ring region. In conclusion, spirals and jets would be the direct manifestation the ongoing collisional activity of the F ring region.

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.

Observations of the November 1995 Sun crossing of the Saturn's ring-plane made with the 3.6m CFH telescope, using the UHAO adaptive optics system, are presented here. We report the detection of four arcs located in the vicinity of the F ring. They can be seen one day later in HST images. The combination of both data sets gives accurate determinations of their orbits. Semi-major axes range from 140020 km to 140080 km, with a mean of 140060 +- 60 km. This is about 150 km smaller than previous estimates of the F ring radius from Voyager 1 and 2 data, but close to the orbit of another arc observed at the same epoch in HST images.