In this work, we model the radial density profile of the outer A ring of Saturn observed with Cassini cameras (Tiscareno & Harris 2018): An axisymmetric diffusion model has been developed, accounting for the outward viscous flow of the ring material, and the counteracting inward drift caused by resonances of the planet's large outer moons. It has been generally accepted that the 7:6 resonance of Janus alone confines the outer A ring which, however, was disproved by Tajeddine et al. (2017). We show that the step-like density profile of the outer A ring is predominantly defined by the discrete first-order resonances of Janus, the 5:3 second-order resonance of Mimas, and the overlapping resonances of Prometheus and Pandora.

One of the most intriguing facets of Saturn's rings are the sharp edges of gaps in the rings where the surface density abruptly drops to zero. This is despite of the fact that the range over which a moon transfers angular momentum onto the ring material is much larger. Recent UVIS-scans of the edges of the Encke and Keeler gap show that this drop occurs over a range approximately equal to the rings' thickness. Borderies et al. (1982, 1989) show that this striking feature is likely related to the local reversal of the usually outward directed viscous transport of angular momentum in strongly perturbed regions. In this article we revise the Borderies et al. (1989) model using a granular flow model to define the shear and bulk viscosities and incorporate the angular momentum flux reversal effect into the axisymmetric diffusion model we developed for gaps in dense planetary rings (Gr\"atz et al. 2018). Finally, we apply our model to the Encke and Keeler division in order to estimate the shear and bulk viscosities in the vicinity of both gaps.

Among the great discoveries of the Cassini mission are the propeller-shaped structures created by small moonlets embedded in Saturn's dense rings. We analyze images of the sunlit side of Saturn's outer A ring, which show the propeller Bleriot with clearly visible partial propeller gaps. By determining radial brightness profiles at different azimuthal locations, we obtain the evolution of the gap minimum downstream of the moonlet. From the radial separation of the partial propeller gaps we estimate the Hill radius of Bleriot to be about 400m. Further, we fit the analytic solution from Sremcevic et al. (2002) describing the azimuthal evolution of the surface mass density in the propeller gap region to the azimuthal gap evolution obtained from Cassini images. From these fits, we estimate a kinematic shear viscosity in the range of 60 cm$^2$/s to 100 cm$^2$/s in Bleriot's ring region. These values are consistent with the parametrization given by Daisaka et al. (2001) and agree well with values estimated for the Encke gap edge (Tajeddine et al. 2017; Graetz et al. 2018).

The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring (Tiscareno et al. 2010, Seiler et al. 2017, Spahn et al. 2018) raises the question, how the propeller responds to the wandering of the central embedded moonlet. Here, we study numerically how the induced propeller is changing for a librating moonlet. It turns out, that the libration of the moonlet induces an asymmetry in the structural imprint of the propeller, where the asymmetry is depending on the moonlet's libration period and amplitude. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with 400 m Hill radius, which is located in the outer A ring. In this way, we are able to apply our findings to the largest found propeller structures -- such as Bl\'eriot -- which are expected to be of similar size. For Bl\'eriot, we can conclude that, supposed the moonlet is librating with the largest observed period of 11.1 years and an azimuthal amplitude of about 1845 km (Seiler et al. 2017, Spahn et al. 2018), a small assymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The excess motions of the other giant propellers -- such as Santos Dumont and Earhart -- have similar amplitudes as Bl\'eriot and thus might allow the observation of larger asymmetries due to their smaller azimuthal extent. This would permit to scan the whole gap structure for asymmetries. Although the librational model of the moonlet is a simplification, our results are a first step towards the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.

One of the biggest successes of the Cassini mission is the detection of small moons (moonlets) embedded in Saturn's rings which cause S-shaped density structures in their close vicinity, called propellers (Spahn and Sremcevic 2000; Tiscareno et al. 2006; Sremcevic et al. 2007). Here, we present isothermal hydrodynamic simulations of moonlet-induced propellers in Saturn's A ring which denote a further development of the original model (Spahn and Sremcevic 2000). We find excellent agreement between these new hydrodynamic and corresponding N-body simulations. Furthermore, the hydrodynamic simulations confirm the predicted scaling laws (Spahn and Sremcevic 2000) and the analytical solution for the density in the propeller gaps (Sremcevic et al. 2002). Finally, this mean field approach allows us to simulate the pattern of the giant propeller Bleriot, which is too large to be modeled by direct N-body simulations. Our results are compared to two stellar occultation observations by the Cassini Ultraviolet Imaging Spectrometer (UVIS), that intersect the propeller Bleriot. Best fits to the UVIS optical depth profiles are achieved for a Hill radius of 590 m, which implies a moonlet diameter of about 860 m. Furthermore, the model favours a kinematic shear viscosity of the surrounding ring material of $\nu_0 = 340$ cm^2/s, a dispersion velocity in the range of 0.3 cm/s $< c_0 <$ 1.5 cm/s, and a fairly high bulk viscosity $7 < \xi_0/\nu_0 < 17$. These large transport values might be overestimated by our isothermal ring model and should be reviewed by an extended model including thermal fluctuations.

The interplanetary space probe Cassini/Huygens reached Saturn in July 2004 after seven years of cruise phase. The Cosmic Dust Analyzer (CDA) measures the interplanetary, interstellar and planetary dust in our solar system since 1999 and provided unique discoveries. In 1999, CDA detected interstellar dust in the inner solar system followed by the detection of electrical charges of interplanetary dust grains during the cruise phase between Earth and Jupiter. The instrument determined the composition of interplanetary dust and the nanometre sized dust streams originating from Jupiter's moon Io. During the approach to Saturn in 2004, similar streams of submicron grains with speeds in the order of 100 km/s were detected from Saturn's inner and outer ring system and are released to the interplanetary magnetic field. Since 2004 CDA measured more than one million dust impacts characterizing the dust environment of Saturn. The instrument is one of three experiments which discovered the active ice geysers located at the south pole of Saturn's moon Enceladus in 2005. Later, a detailed compositional analysis of the water ice grains in Saturn's E ring system lead to the discovery of large reservoirs of liquid water (oceans) below the icy crust of Enceladus. Finally, the determination of the dust- magnetosphere interaction and the discovery of the extended E ring (at least twice as large as predicted) allowed the definition of a dynamical dust model of Saturn's E ring describing the observed properties. This paper summarizes the discoveries of a ten year story of success based on reliable measurements with the most advanced dust detector flown in space until today. This paper focuses on cruise results and findings achieved at Saturn with a focus on flux and density measurements.

In this paper, the dynamical analysis of the Jovian dust originating from the four Galilean moons is presented. High accuracy orbital integrations of dust particles are used to determine their dynamical evolution. A variety of forces are taken into account, including the Lorentz force, solar radiation pressure, Poynting-Robertson drag, solar gravity, the satellites' gravity, plasma drag, and gravitational effects due to non-sphericity of Jupiter. More than 20,000 dust particles from each source moon in the size range from 0.05 micron to 1 cm are simulated over 8,000 (Earth) years until each dust grain hits a sink (moons, Jupiter, or escape from the system). Configurations of dust number density in the Jovicentric equatorial inertial frame are calculated and shown. In a Jovicentric frame rotating with the Sun the dust distributions are found to be asymmetric. For certain small particle sizes, the dust population is displaced towards the Sun, while for certain larger sizes, the dust population is displaced away from the Sun. The average lifetime as a function of particle size for ejecta from each source moon is derived, and two sharp jumps in the average lifetime are analyzed. Transport of dust between the Galilean moons and to Jupiter is investigated. Most of the orbits for dust particles from Galilean moons are prograde, while, surprisingly, a small fraction of orbits are found to become retrograde mainly due to solar radiation pressure and Lorentz force. The distribution of orbital elements is also analyzed.

The reconstruction of the orbital evolution of the propeller structure Bl\'eriot orbiting in Saturn's A ring from recurrent observations in Cassini ISS images yielded a considerable offset motion from the expected Keplerian orbit (Tiscareno et al., 2010). This offset motion can be composed by three sinusoidal harmonics with amplitudes and periods of 1845, 152, 58 km and 11.1, 3.7 and 2.2 years, respectively (Srem\v{c}evi\'c et al., 2014). In this paper we present results from N-Body simulations, where we integrated the orbital evolution of a moonlet, which is placed at the radial position of Bl\'eriot under the gravitational action of the Saturnian satellites. Our simulations yield, that especially the gravitational interactions with Prometheus, Pandora and Mimas is forcing the moonlet to librate with the right frequencies, but the libration-amplitudes are far too small to explain the observations. Thus, further mechanisms are needed to explain the amplitudes of the forced librations - e.g. moonlet-ring interactions. Here, we develop a model, where the moonlet is allowed to be slightly displaced with respect to its created gaps, resulting in a breaking point-symmetry and in a repulsive force. As a result, the evolution of the moonlet's longitude can be described by a harmonic oscillator. In the presence of external forcings by the outer moons, the libration amplitude gets the more amplified, the more the forcing frequency gets close to the eigenfrequency of the disturbed propeller oscillator. Applying our model to Bl\'eriot, it is possible to reproduce a libration period of 13 years with an amplitude of about 2000 km.

Most theoretical investigations of dust charging processes in space have treated the current balance condition as independent of grain size. However, for small grains, as they are often observed in space environments, a dependence on grain size is expected due to secondary electron emission (SEE). Here, by the term "small" we mean a particle size comparable to the typical penetration depth for given primary electron energy. The results are relevant for the dynamics of small, charged dust particles emitted by the volcanic moon Io, which form the Jovian dust streams. We revise the theory of charging of small (sub-micron sized) micrometeoroids to take into account a high production of secondary electrons for small grains immersed in an isotropic flux of electrons. We apply our model to obtain an improved estimate for the charge of the dust streams leaving the Jovian system, detected by several spacecraft. For the Jupiter plasma environment we derive the surface potential of grains composed of NaCl (believed to be the major constituent of Jovian dust stream particles) or silicates. For small particles the potential depends on grain size and the secondary electron current induces a sensitivity to material properties. As a result of the small particle effect, the estimates for the charging times and for the fractional charge fluctuations of NaCl grains obtained using our general approach to SEE give results qualitatively different from the analogous estimates derived from the traditional approach to SEE. We find that for the charging environment considered in this paper field emission does not limit the charging of NaCl grains.

A kinky and clumpy ringlet shares orbit with the moon Pan in the center of the 320-km wide Encke gap in Saturn's rings (Porco et al., 2005). The ringlet is mainly composed of micron-sized particles (Showalter, 1991, Hedman et al., 2011), implying that these particles may be significantly perturbed by non-gravitational forces, which can limit their lifetimes. We establish a kinetic model considering the birth, evolution, and death of dust in the Encke central ringlet allowing to evaluate the ringlet optical depth. First, we investigate the generation of dust by micrometeorite impacts (the `impact-ejecta' process) on putative, yet undetected embedded moonlets. Taking into account the orbital evolution under the influence of the relevant perturbation forces, the dominant loss mechanisms are collisions with ring particles in the gap edges, the putative moonlets in the gap, or erosion by sputtering in Saturn's plasma environment. However, our results show that this impact-ejecta process alone can only sustain a ringlet of optical depth 3-4 orders of magnitude smaller than the observed values. Consequently, other processes must be taken into account. For example, mutual collisions among putative moonlets should produce dust at about the same rate, while further disruption of ejecta, as proposed by Dikarev (1999), should increase the total amount of particles. Furthermore, observations show an azimuthal asymmetry of the material in the Encke gap ringlets (Ferrari and Brahic 1997, M. Srem\v{c}evi\'c, private communication). We investigate the scenario proposed by Hedman et al. (2013) that for the Encke central ringlet, the observed asymmetry is mainly due to the combined action of plasma drag and Pan's gravity causing a focusing of dust in the region leading Pan's orbit (Hedman et al. 2013).

Methone and Anthe are two tiny moons (with diameter $<3$ km) in the inner part of Saturn's E ring. Both moons are embedded in an arc of dust particles. To understand the amount of micron-sized dust and their spatial distribution in these arcs, we model the source, dynamical evolution, and sinks of these dust in the arc. We assume hypervelocity impacts of micrometeoroids on the moons as source of these dust (Hedman et al., 2009), the so called impact-ejecta process (Krivov et al., 2003; Spahn et al., 2006). After ejecting and escaping from the moons, these micron-sized particles are subject to several perturbing forces, including gravitational perturbation from Mimas, oblateness of Saturn, Lorentz force, solar radiation pressure, and plasma drag. Particles can be either confined in the arcs due to corotational resonance with Mimas, as their source moons (Spitale et al., 2006; Cooper et al., 2008; Hedman et al., 2009), or pushed outward by plasma drag. Particle sinks are recollisions with the source moon, collision with other moons, or migrate out of the radial zone of interest. In addition to that, the upper limit of the particle lifetimes are controlled by plasma sputtering, which decreases the particle size in a rate of order 1$\mu$m radius every 100 years (Johnson et al., 2008). Our simulation results show that ejecta from both moons can form the arcs of maximal optical depths $\tau$ in the order of $10^{-8} - 10^{-6}$, although the absolute amount of dust have uncertainties. We also find out the longitudinal extension of the arcs in our simulation are consistent with observation and the theory. Smaller particles are more likely to escape the arc because of the stronger influence of plasma drag. On the other hand, large particles can stay in arcs for longer time and therefore are more likely to collide with the source moon.

Saturn's rings consist of a huge number of water ice particles, with a tiny addition of rocky material. They form a flat disk, as the result of an interplay of angular momentum conservation and the steady loss of energy in dissipative inter-particle collisions. For particles in the size range from a few centimeters to a few meters, a power-law distribution of radii, $\sim r^{-q}$ with $q \approx 3$, has been inferred; for larger sizes, the distribution has a steep cutoff. It has been suggested that this size distribution may arise from a balance between aggregation and fragmentation of ring particles, yet neither the power-law dependence nor the upper size cutoff have been established on theoretical grounds. Here we propose a model for the particle size distribution that quantitatively explains the observations. In accordance with data, our model predicts the exponent $q$ to be constrained to the interval $2.75 \le q \le 3.5$. Also an exponential cutoff for larger particle sizes establishes naturally with the cutoff-radius being set by the relative frequency of aggregating and disruptive collisions. This cutoff is much smaller than the typical scale of micro-structures seen in Saturn's rings.

We study the vertical extent of propeller structures in Saturn's rings. Our focus lies on the gap region of the propeller and on non-inclined propeller moonlets. In order to describe the vertical structure of propellers we extend the model of Spahn and Sremcevic (2000) to include the vertical direction. We find that the gravitational interaction of ring particles with the non-inclined moonlet does not induce considerable vertical excursions of ring particles, but causes a considerable thermal motion in the ring plane. We expect ring particle collisions to partly convert the lateral induced thermal motion into vertical excursions of ring particles. For the gap region of the propeller, we calculate gap averaged propeller heights on the order of 0.7 Hill radii, which is of the order of the moonlet radius. In our model the propeller height decreases exponentially until viscous heating and collisional cooling balance. We estimate Hill radii of 370m and 615m for the propellers Earhart and Bleriot. Our model predicts about 120km for the azimuthal extent of the Earhart propeller at Saturn's 2009 equinox, being consistent with values determined from Cassini images.

We propose a new mechanism which explains the existence of enormously sharp edges in the rings of Saturn. This mechanism is based on the synchronization phenomenon due to which the epicycle rotational phases of particles in the ring, under certain conditions, become synchronized with the phase of external satellite, e.g. with the phase of Mimas in the case of the outer B ring edge. This synchronization eliminates collisions between particles and suppress the diffusion induced by collisions by orders of magnitude. The minimum of the diffusion is reached at the center of the synchronization regime corresponding to the ratio 2:1 between the orbital frequency at the edge of B ring and the orbital frequency of Mimas. The synchronization theory gives the sharpness of the edge in few tens of meters that is in agreement with available observations.