Since the discovery of an ice particle plume erupting from the south polar terrain on Saturn's moon Enceladus, the geophysical mechanisms driving its activity have been the focus of substantial scientific research. The pattern and deposition rate of plume material on Enceladus' surface is of interest because it provides valuable information about the dynamics of the ice particle ejection as well as the surface erosion. Surface deposition maps derived from numerical plume simulations by Kempf et al. (2010) have been used by various researchers to interpret data obtained by various Cassini instruments. Here, an updated and detailed set of deposition maps is provided based on a deep-source plume model (Schmidt et al., 2008), for the eight ice-particle jets identified in Spitale and Porco (2007), the updated set of jets proposed in Porco et al. (2014), and a contrasting curtain-style plume proposed in Spitale et al. (2015). Methods for computing the surface deposition are detailed, and the structure of surface deposition patterns is shown to be consistent across changes in the production rate and size distribution of the plume. Maps are also provided of the surface deposition structure originating in each of the four Tiger Stripes. Finally, the differing approaches used in Porco et al. (2014) and Spitale et al. (2015) have given rise to a jets vs. curtains controversy regarding the emission structure of the Enceladus plume. Here we simulate each, leading to new insight that, over time, most emissions must be directed relatively orthogonal to the surface because jets "tilted" significantly away from orthogonal lead to surface deposition patterns inconsistent with surface images. Data for maps are available in HDF5 format for a variety of particle sizes at http://impact.colorado.edu/southworth_data.

A new model for the shape of the prominent eccentric ringlet in the gap exterior to Saturn's B-ring is developed based on Cassini imaging observations taken over about 8 years. Unlike previous treatments, the new model treats each edge of the ringlet separately. The Keplerian component of the model is consistent with results derived from Voyager observations, and $m=2$ modes forced by the nearby Mimas 2:1 Lindblad resonance are seen. Additionally, a free $m=2$ mode is seen on the outer edge of the ringlet. Significant irregular structure that cannot be described using normal-mode analysis is seen on the ringlet edges as well. Particularly on the inner edge, that structure remains coherent over multi-year intervals, moving at the local Keplerian rate. We interpret the irregular structure as the signature of embedded massive bodies. The long coherence time suggests the responsible bodies are concentrated near the edge of the ringlet. Long wake-like structures originate from two locations on the inner edge of the ringlet, revealing the locations of the two most massive embedded bodies in that region. As with the Voyager observations, the Cassini data sets showed no correlation between the width and the radius of the ringlet as would be expected for a self-gravitating configuration, except for a brief interval during late 2006, when the width-radius relation was similar to those seen in most other narrow eccentric ringlets in the Solar System.

A new symplectic N-body integrator is introduced, one designed to calculate the global 360 degree evolution of a self-gravitating planetary ring that is in orbit about an oblate planet. This freely-available code is called epi_int, and it is distinct from other such codes in its use of streamlines to calculate the effects of ring self-gravity. The great advantage of this approach is that the perturbing forces arise from smooth wires of ring matter rather than discreet particles, so there is very little gravitational scattering and so only a modest number of particles are needed to simulate, say, the scalloped edge of a resonantly confined ring or the propagation of spiral density waves. The code is applied to the outer edge of Saturn's B ring, and a comparison of Cassini measurements of the ring's forced response to simulations of Mimas' resonant perturbations reveals that the B ring's surface density at its outer edge is 195+-60 gm/cm^2 which, if the same everywhere across the ring would mean that the B ring's mass is about 90% of Mimas' mass. Cassini observations show that the B ring-edge has several free normal modes, which are long-lived disturbances of the ring-edge that are not driven by any known satellite resonances. Although the mechanism that excites or sustains these normal modes is unknown, we can plant such a disturbance at a simulated ring's edge, and find that these modes persist without any damping for more than ~10^5 orbits or ~100 yrs despite the simulated ring's viscosity of 100 cm^2/sec. These simulations also indicate that impulsive disturbances at a ring can excite long-lived normal modes, which suggests that an impact in the recent past by perhaps a cloud of cometary debris might have excited these disturbances which are quite common to many of Saturn's sharp-edged rings.

Voyager images and Cassini occultation data have previously shown that the behavior of the outer edge of Saturn's massive B ring is determined only in part by a static response to the 2:1 inner Lindblad resonance with Mimas. In Cassini images of this region, we find, in addition to the expected wavenumber-2 forced distortion, evidence for unforced self-excited wavenumber-3, wavenumber-2, and wavenumber-1 normal modes. These are the first observations to suggest substantial wave amplification in Saturn's broad rings. Moreover, the presence of these free modes strongly implicates viscous overstability as their underlying cause and, by inference, the cause for most if not all of the unforced structures throughout the high-mass-density B ring and in other high-mass-density regions in Saturn's rings. Analysis of each of the inferred waves reveals a consistent lower bound on the average surface mass density of ~ 44 g/cm^2 for the outer 250 km of the ring, though the true surface density could be as high as 100 g/cm^2 or higher. Interference between the forced and free wavenumber-2 modes yields a total wavenumber-2 pattern that varies in amplitude and orientation with a characteristic period of ~ 5.5 years. We also find localized disturbances, including 3.5-km-tall vertical structures, that provide circumstantial evidence for embedded massive bodies in the Mimas resonance zone. The presence of such bodies is supported by the presence of a shadow-casting moonlet ~ 0.3 km wide near the ring's edge.

We report the discovery of several large "propeller" moons in the outer part of Saturn's A ring, objects large enough to be followed over the 5-year duration of the Cassini mission. These are the first objects ever discovered that can be tracked as individual moons, but do not orbit in empty space. We infer sizes up to 1--2 km for the unseen moonlets at the center of the propeller-shaped structures, though many structural and photometric properties of propeller structures remain unclear. Finally, we demonstrate that some propellers undergo sustained non-keplerian orbit motion. (Note: This arXiv version of the paper contains supplementary tables that were left out of the ApJL version due to lack of space).

(Abridged) The following describes a model of a broad planetary ring whose sharp edge is confined by a satellite's m^th Lindblad resonance (LR). This model uses a streamline formalism to calculate the ring's internal forces, namely, ring gravity, pressure, viscosity, as well as a hypothetical drag force. The model calculates the streamlines' forced orbit elements and surface density throughout the perturbed ring. The model is then applied to the outer edge of Saturn's B ring, which is maintained by an m=2 inner LR with the satellite Mimas. Ring models are used to illustrate how a ring's perturbed state depends on the ring's physical properties: surface density, viscosity, dispersion velocity, and the hypothetical drag force. A comparison of models to the observed outer B ring suggests that the ring's surface density there is between 10 and 280 gm/cm^2. The ring's edge also indicates where the viscous torque counterbalances the perturbing satellite's gravitational torque on the ring. But an examination of seemingly conventional viscous B ring models shows that they all fail to balance these torques at the ring's edge. This is due ring self-gravity and the fact that a viscous ring tends to be nearly peri-aligned with the satellite, which reduces the satellite's torque on the ring and makes the ring's edge more difficult to maintain. Nonetheless, the following shows that a torque balance can still be achieved in a viscous B ring, but only in an extreme case where the ratio of the ring's bulk/shear viscosities satisfy ~10^4. However, if the dissipation of the ring's forced motions is instead dominated by a weak drag force, then the satellite can exert a much stronger torque that can counterbalance the ring's viscous torque.