Draft version July 20, 2021 Typeset using LATEX preprint2 style in AASTeX62

Evidence that a novel type of satellite wake might exist in Saturn’s E ring M.M. Hedman1 and M. Young1 1Department of Physics, University of Idaho, Moscow Idaho 83843

ABSTRACT Saturn’s E ring consists of micron-sized particles launched from Enceladus by that moon’s geological activity. A variety of small-scale structures in the E-ring’s brightness have been attributed to tendrils of material recently launched from Enceladus. However, one of these features occurs at a location where Enceladus’ gravitational perturbations should concentrate background E-ring particles into structures known as satellite wakes. While satellite wakes have been observed previously in ring material drifting past other moons, these E-ring structures would be the first examples of wakes involving particles following horseshoe orbits near Enceladus’ orbit. The predicted intensity of these wake signatures are particularly sensitive to the fraction E-ring particles’ on orbits with low eccentricities and semi-major axes just outside of Enceladus’ orbit, and so detailed analyses of these and other small-scale E-ring features should place strong constraints on the orbital properties and evolution of E-ring particles.

1. INTRODUCTION cal models of the E-ring’s large-scale structure Saturn’s E ring is a broad and tenuous ring suggested that its broad extent arises because composed of very small (< 5µm) ice grains the orbital eccentricities and inclinations of the that encompasses the orbits of many of Sat- small particles launched from Enceladus can urn’s icy moons (see Kempf et al. 2018; Hed- be greatly increased by a combination of so- man et al. 2018, for recent reviews). This ring lar radiation pressure and Saturn’s offset mag- is brightest around the orbit of the small moon netic field, while on longer timescales the parti- Enceladus, and so it was long suspected that it cles are slowly eroded and transported outwards had some connection with that moon (Showal- by plasma drag (Hor´anyi et al. 1992; Hamil- ter et al. 1991; Hor´anyi et al. 1992; Hamilton ton 1993a,b; Hamilton and Burns 1994; Juh´asz and Burns 1994; Burns et al. 2001). This suspi- and Hor´anyi 2002, 2004; Juh´aszet al. 2007; cion was dramatically confirmed by the Cassini Hor´anyi et al. 2008). Cassini observations have arXiv:2105.05294v2 [astro-ph.EP] 19 Jul 2021 mission, which revealed that Enceladus was geo- confirmed many aspects of these models, but logically active, with a series of fissures around have also identified large-scale asymmetries in the moon’s south pole launching micron-sized the ring’s brightness which imply that particle particles from beneath the moon’s surface into transport within this ring is more complex than orbit around Saturn (see Schenk et al. 2018, and expected (Hedman et al. 2012; Ye et al. 2016). references therein). Fine-scale structures in the E ring close to In addition to revealing the E-ring’s origin, Enceladus provide important opportunities to the Cassini mission clarfied several aspects of further explore the connections between Ence- the ring’s structure and dynamics. Theoreti- ladus and the E ring. Thus far, published stud- 2

W1537013154 Mitchell et al. 2015 Simulations

Interior Tendril Leading Tendril

W1725971847 Mitchell et al. 2015 Simulations

Leading Tendril

Interior Tendril

Figure 1. Comparisons of E-ring images with numerical simulations, adapted from Mitchell et al.(2015). The left panels show high-pass filtered and brightness inverted versions of two images of the E ring around Enceladus. The middle and right panels show two different simulations of material recently ejected from Enceladus, with different assumptions about Enceladus’ activity. Two concentrations of material near Ence- ladus are marked as the Interior and Leading Tendrils. (Note the direction of orbital motion is counterclock- wise in the top image and clockwise in the bottom image). For both observations, the leading tendril extends closer to Enceladus in the images than it does in the simulations, so that it is more prominent between the two dotted horizontal lines. 3 ies of these features have focused on whether the gravitational perturbations from a moon they can be attributed to the activity of spe- produce organized epicyclic motions in nearby cific sources near Enceladus’ south pole (Kempf ring particles, which in turn produce localized et al. 2010; Mitchell et al. 2015). However, in variations in the surface density with a similar this paper we will take a closer look at a sub-set orientation to the features seen here (Murray of these features that might represent structures and Dermott 2001). in the background E ring rather than freshly Mitchell et al.(2015) argued that it was un- ejected material from Enceladus. likely for any E-ring features to represent wakes Figure1 shows high-pass filtered and brightness- because the spacing of the features was incon- inverted versions of two images of the E ring sistent with those expected for standard wakes, first analyzed by Mitchell et al.(2015), along and because it is unlikely such structures would with numerical simulations performed by those form within a ring where most of the particles authors that track particles of various sizes have substantial orbital eccentricity. However, launched from various locations around Ence- a tenuous ring like the E ring can contain a ladus’ south pole. The particles in the simu- novel class of satellite wakes that occur in mate- lations are organized into elongated structures rial following horseshoe orbits. These “rebound called “tendrils”. One of these tendrils extends wakes” occur at very different locations from interior from the moon and is designated as the other satellite wakes, and the most prominent interior tendril, while others are found at vari- of these wakes should occur very close to the ous locations in front of the moon, with the first extensions of the leading tendril seen in these of the “leading tendrils” being the most promi- two images. nent in the images. If we take a closer look In the rest of this paper, we argue that satel- at the first leading tendril, we can note that it lite wakes are a reasonable potential explana- extends closer to Enceladus in the images than tion for these specific features, and so this pos- it does in the simulations. This is most obvi- sibility merits further examination with addi- ous in the earlier image W1537013154, where a tional numerical simulations. In Section2 we prominent band can be seen between the two use Hill’s Equations to demonstrate that mate- dashed lines that is not present in the sim- rial on horseshoe orbits can produce wake-like ulations. For W1725971847 the difference is density variations. In Section3 we describe how more subtle, but it still appears that the in- we process the two images in order to determine tensity of the band remains constant across the the locations and brightnesses of the relevant dotted line, whereas the simulations predict a features. In Section4 we show that the bright sharp edge to this feature. In principle, these features extending beyond the tips of the lead- discrepancies could potentially be resolved by ing tendrils are consistent with the expected lo- simulating plume particles with a broader range cations of these wakes. We also show that there of physical parameters. However, it turns out could be enough particles on the appropriate that these features could also be density varia- orbits to produce the observed brightness varia- tions in pre-existing E-ring particles known as tions. Finally, in Section5 we discuss the poten- satellite wakes. tial implications of this feature being composed Satellite wakes have been observed in Saturn’s of particles on horseshoe orbits for the orbital main rings near the small moons Pan and Daph- distribution of E-ring particles. nis (Showalter et al. 1986; Horn et al. 1996; Tis- . careno et al. 2007). These wakes arise because 4

Figure 2. Trajectories of particles on orbits with the same eccentricity and pericenter location as a satellite (located at the origin), computed using Hill’s Equations. The top panel shows the trajectories of particles on that pass by the moon, while the bottom panel shows the trajectories of particles that undergo horseshoe motion. Note that the initial semi-major axis spacing of the trajectories is different in the two panels in order to better visualize the locations of the relevant wake patterns.

2. THEORY OF SATELLITE WAKES IN MP . Let the coordinates ξ and η be aligned MATERIAL ON HORSESHOE ORBITS with the local radial and azimuthal directions, The basic theory of satellite wakes was first measured in units of the satellite’s Hill Radius: developed to explain periodic brightness and op- tical depth variations found in the vicinity of 1/3 the narrow Encke Gap in Saturn’s main rings RH = as(ms/3MP ) . (1) (Showalter et al. 1986), and the relevant physics can be illustrated using Hill’s Equations. Consider a reference frame whose origin is at Note that for Enceladus the Hill Radius is the location of a satellite of mass ms on an orbit 939.74 km. In this coordinate system, the equa- with semi-major axis as around a planet of mass tions of motion for a test particle are given by Hill’s Equations (Henon and Petit 1986): 5  1  ¨ 2 passing the moon, the particles execute horse- ξ − 2nη˙ = 3n ξ 1 − 2 2 3/2 (2) (ξ + η ) shoe motion due to the semi-major axis evolu- tion induced by the satellite’s gravitational per- 3n2η η¨ + 2nξ˙ = − (3) turbations. However, as shown in the bottom (ξ2 + η2)3/2 panel of Figure2, wake-like structures can oc- where n is the satellite’s orbital mean motion. cur in this material. In this case, particles ap- Note these equations describe the motion of the proach the moon on initially circular orbits from test particle relative to the satellite even if the the right, but turn around before reaching the satellite is on an eccentric orbit (Henon and Pe- satellite’s location at distances that depend on tit 1986). the initial semi-major axis separation. When Standard satellite wakes arise when we con- the initial semi-major axis difference between sider particles initially on orbits with fixed ξ the particle and the moon is large enough the interior or exterior to the moon and when the particle also acquires a significant eccentricity semi-major axis difference between the parti- as it turns around in front of the moon. The cles and the moon ∆a is sufficiently large that resulting radial motions again cause the parti- the particles will not undergo horseshoe motion cles to follow oscillating trajectories with wave- (see top panel of Figure2), Due to Keplerian lengths that are proportional to the semi-major shear, these particles all drift in longitude rela- axis separation from the moon, giving rise to tive to the satellite. As the particles drift past diagonal bands of enhanced density similar to the moon, the gravitational pull from the moon those seen for material passing by the moon. To causes all the particles to acquire finite eccen- distinguish these wakes from the classic satellite tricities and thus execute radial motions with wakes, we call the wakes in material that drifts an amplitude that steadily decreases with in- past the moon “passing wakes” and the wakes in creasing ∆a (Julian and Toomre 1966; Lin and the material undergoing horseshoe motion “re- Papaloizou 1979; Murray and Dermott 2001). bound wakes”. As these particles continue to move downstream While the locations of both passing and re- from the moon, their radial epicyclic motion bound wakes can be most easily illustrated with causes them to follow quasi-sinusoidal paths the trajectories of particles initially on orbits with wavelength 3π∆a (Murray and Dermott with a constant ξ, most of the particles in the 2001). These systematic trends in the wave- E ring are not on these sorts of orbits. Tra- lengths of the particle’s trajectories in turn jectories with constant ξ in the E ring corre- cause the distances between particles on nearby spond to particles with the same orbital eccen- semi-major axes to vary, yielding variations in tricity and pericenter as Enceladus (Henon and the local particle density. More specifically, Petit 1986), and the particles in the E ring have there are diagonal bands where the particle tra- a wide range of eccentricities that are driven jectories get very close together and so the den- by non-gravitational forces like solar radiation sity should be exceptionally high. The first of pressure (Hor´anyi et al. 1992; Hamilton 1993a; these wakes is most obvious between 5 and 25 Hamilton and Burns 1994; Juh´aszand Hor´anyi Hill Radii in azimuth and -2 to -4 Hill Radii in 2002, 2004; Juh´asz et al. 2007; Hor´anyi et al. radius, while the second wake starts at 35 Hill 2008; Hedman et al. 2012; Ye et al. 2016). Par- Radii and extends off the page. ticles with orbital eccentricities and pericenters Material orbiting on semi-major axes closer to that differ from those of Enceladus would fol- the satellite is not normally considered suitable low cycloid-like trajectories spanning a range of for forming satellite wakes because instead of 6

ξ values. For such particles, it is difficult to vi- demonstrate that the locations of the rebound sualize the variations in the particle density us- wake signatures are insensitive to the orbital ing individual particle trajectories. Instead one parameter distribution, and their intensity is needs to generate maps of particle density for proportional to the fraction of particles on low- each specified distribution of orbital elements by eccentricity orbits with the appropriate range of simulating a sufficiently large number of parti- semi-major axes. For these simulations we trace cles. Mitchell et al.(2015) did two such simula- the trajectories of 105 particles randomly drawn tions, one where all the particles had negligible from a population with the following distribu- eccentricity and another where the particles had tions of initial orbital parameters: eccentricities of 0.02. The wake signals were not obvious in the latter simulation, implying that • A uniform distribution of initial eccentric- only particles on nearly circular orbits can pro- ities between 0 and 0.04 (corresponding to duce wake signatures. However, the features initial radial oscillations with amplitudes seen in the images are actually rather subtle between 0 and 10 Hill Radii). variations in the E-ring’s brightness, being only • A uniform distribution of pericenter loca- about 5% brighter than their surroundings (see tions. Section 4.2). Hence it is worth examining how the locations and intensity of potential wake sig- • A Gaussian distribution of semi-major natures vary among particle populations with axes with a center 1.5 Hill Radii external finite eccentricities. to the moon’s orbit and a Gaussian width The spatial density of particles with a range of of 0.5, 1.0 or 2.0 Hill radii. orbital properties can be derived from the tra- jectories of individual particles computed using Note that the eccentricity and semi-major axis Hill’s Equations. However, for any given initial distributions considered here are not intended particle orbital parameter distribution only a fi- to be fully representative of the E-ring, but in- nite number of these tracks can be calculated, stead serve to illustrate how the locations and and this will give rise to artificial small-scale intensity of wake features depend on the orbital density variations that can obscure the desired properties of the ring particles. The eccentric- wake signatures. This “sampling noise” can al- ity distribution was chosen because the wake ways be reduced by increasing the number of signatures see in Figure2 are only likely to be simulated trajectories, but this will also increase prominent for particles on orbits with initial ra- the time needed to do the simulations. In prac- dial motions less than 1 Hill Radii, which cor- tice, we find that 105-106 trajectories need to be responds to eccentricities less than about 0.005. simulated in order for wake-like signatures with The above eccentricity distribution therefore en- amplitudes of 5% to be detectable. These sorts sures that the fraction of particles on such orbits of calculations are therefore an inefficient way is small but not zero. Furthermore, this dis- to explore a broad parameter space of potential tribution gives a median eccentricity similar to orbital parameter distributions. Hence we will that used by Mitchell et al.(2015). The pericen- instead consider a small number of situations ter distribution is uniform both for the sake of that illustrate important aspects of the signals simplicity and because the observed large-scale from satellite wakes. structure of the E ring does not show strong First, we will simulate particles with a limited asymmetries in pericenter locations (Hedman range of semi-major axes and eccentricities to et al. 2012; Ye et al. 2016). Finally, the semi- major axis distributions are chosen so that there 7

Figure 3. Simulated particle density maps for par- Figure 4. The same simulated particle density ticle populations with uniform eccentricity distri- maps as shown in Figure3, except they have been butions between 0 and 0.04, uniform pericenter lo- high-pass filtered by removing a smoothed version cation distributions and Gaussian semi-major axis of the maps with a smoothing length of 1.5 Hill distributions with a center at +1.5 Hill Radii and Radii. The brightness range corresponds to frac- the three different widths indicated. All panels tional brightness variations of ±10%. Here we can show a broad distribution of ring particles, but in see a brightness maximum due to the rebound wake the top panel a density enhancement can be seen at in all three simulations. Note the strength of the the expected location of the rebound wake around feature decreases as the width of the semi-major -1 Hill Radii in radius and 15 Hill Radii in longi- axis distribution increases. tude. Note that no passing wakes are visible in this simulation because of the assumed semi-major axis dius of the moon. The resulting tables of parti- distribution. cle locations were then binned together to esti- mate the relative density of particles as func- are some particles on orbits likely to produce re- tions of radius and longitude relative to the bound wakes in all three cases, but the fraction moon. The resulting density maps are shown in of particles on these orbits differs. Figure3. These maps are clearly dominated by 5 The trajectories for the 10 particles drawn radial density variations that mostly reflect the from each of these populations were evaluated initial eccentricity distribution of the particles. at fixed time steps using Hill’s equations, start- The simulations with larger spreads in initial ing from a fixed longitude over 70 Hill Radii semi-major axis also show a gradual decrease in in front of the moon. These trajectories were the peak particle density near the moon’s loca- terminated if the particles got within 1 Hill Ra- 8 tion. This is because a substantial fraction of These numbers are comparable to the fraction the particles are on horseshoe orbits that turn of particles in the simulations on initial orbits around at various distances from the moon. For with eccentricities below 0.005 and semi-major the simulation with the tightest semi-major axis axes between 1 and 1.5 Hill radii exterior to the distribution, there is a ⊂ shape in the locations moon. Thus it is reasonable to conclude that of peak particle density that corresponds to the particles on these particular orbits generate a expected track of circular particles with initial wake signature at the location predicted by Fig- semi-major axes around 1.5 Hill Radii outside ure2, while particles on orbits with higher ec- of the moon’s orbit. The lower branch of this centricities and different semi-major axes sim- track is narrower than the upper branch because ply contribute a smooth background in this re- some fraction of the particles either pass by the gion. This naturally results in the observed frac- moon or are lost when they get too close to the tional density (or brightness) variations due to moon. the wake being roughly equal to the fraction Turning our attention to potential signals of particles on suitable orbits in this region be- from rebound wakes, we can note that the top cause the density variations due to the wake are panel in Figure3 shows a bright diagonal streak of order unity for the particles with the appro- running between -1 and -3 Hill Radii in radius priate orbital parameters. and 10 and 20 Hill Radii in longitude. This While the above simulations suggest that the streak falls at the expected location of the first amplitude of the wake signatures have a rel- rebound wake shown in Figure2, so signatures atively straightforward relationship with the of rebound wakes can indeed appear in particle particles’ orbit element distribution, complica- populations with a range of eccentricities and tions arise when we consider particle popula- semi-major axes. Similar wake signatures are tions with broad eccentricity distributions like also present in the other two simulations, but the E ring. Hedman et al.(2020) found that the they are harder to see against the background radial brightness profile of the E ring could be trends. Hence Figure4 shows the fractional reasonably well matched if the particles in the E difference between each of these maps and ver- ring had the following eccentricity distribution. sions of the maps smoothed over a scale of 1.5 F0 Hill Radii. In all of these high-pass-filtered F(e) = [1 − exp(−e/ec)] (4) 1 + e2/e2 maps, there is a localized density enhancement o around -1 Hill Radii in radius and 12 Hill Radii with eo ' 0.17 and ec ' 0.01. Such an ec- in longitude. These density enhancements all centricity distribution has only about 0.5% of occur where the first rebound wake is expected the particles on orbits with eccentricities less to be the most intense, so it is reasonable to than 0.005. At first, this would appear to im- attribute them all to this wake. The location ply that there would not be sufficient particles of this density enhancement therefore appears to produce wake signals of intensities compara- to be insensitive to the particles’ orbit distribu- ble to those seen in the images and the simu- tion. lations shown in Figure4. However, particles The intensity of these wake signatures also with small eccentricities spend more time close clearly declines as the width of the semi-major to Enceladus’ orbit than particles with large ec- axis distribution increases. The fractional den- centricities, so the particle density near Ence- sity variations associated with this feature are ladus’ orbit will have a larger fraction of parti- 10%, 5% and 2.5% for these three simulations. cles on low-eccentricity orbits. The density of particles near the E-ring’s core with a specified 9

Figure 5. Simulated particle density maps for par- Figure 6. The same simulated particle density ticle populations with the Hedman et al.(2020) ec- maps as shown in Figure5, except they have been centricity distributions, uniform pericenter location high-pass filtered by removing a smoothed version distributions and uniform semi-major axis between of the maps with a smoothing length of 4 Hill the values indicated. All panels show a broad dis- Radii. The brightness range corresponds to frac- tribution of ring particles. Note that the simulated tional brightness variations of ±10%. Here we brightness distributions for all three populations can see a brightness maximum due to the rebound are very similar, although a weak rebound wake wake around +12 Hill Radii in azimuth and -1 Hill signal can be seen in the top panel. radii in Radius in at least the top two simulations. Again, the brightness of this feature decreases as range of eccentricities can be estimated using the width of the semi-major axis distribution in- the following formula (Hedman et al. 2020): creases.

Z emax rF(e) in the images (provided a sufficiently high frac- d(r) = de. (5) 2p 2 2 tion of the particles on low-eccentricity orbits emin πa e − (r/a − 1) had the appropriate semi-major axes). This formula indicates that around Enceladus’ In order to demonstrate that wake signals can orbit roughly 10% of the particles have orbital still be reasonably prominent in particle pop- eccentricities less than 0.005, and a similar frac- ulations with broad eccentricity distributions, tion should be within 0.005 of Enceladus’ orbital we performed simulations tracing the trajecto- parameters. This would imply that even a ring ries of 106 particles randomly drawn from the with a broad eccentricity distribution could pro- a population with the following distributions of duce wake signals comparable to those observed initial orbital parameters: 10

• The eccentricity distribution given by narrowest semi-major axis distribution shows a Equation4. clear brightness enhancement around -1 Hill Ra- dius in radius and 12 Hill Radii in longitude, • A uniform distribution of pericenter loca- consistent with the expected location of a re- tions. bound wake. A weak signal around the same • A uniform distribution of semi-major axes location can also be seen in the middle panel. between +0.1, 0.5 or 0.75 Hill Radii and For the bottom simulation the signal is too weak 1.5 Hill Radii of Enceladus to be robustly detectable, but there may be a slight brightness enhancement near the noise As with the previous simulations, these trajec- level at the appropriate location. The peak tories were evaluated at fixed time steps using amplitudes of the wake signals in these simula- Hill’s equations, starting from a fixed longitude tions are around 6%, 3% and 1.5%, respectively, over 70 Hill Radii in front of the moon, and which is again comparable to the brightness these trajectories were terminated if the par- variations associated with the small-scale struc- ticles got within 1 Hill Radius of the moon. tures in the E-ring. These simulations therefore Again, it is important to note that these cal- demonstrate that detectable wake signals can culations are not meant to be a fully realistic occur in particle populations with broad eccen- simulation of the E ring. In particular, Hill’s tricity distributions. equations are probably not accurate equations of motion for the particles with high eccentric- 3. OBSERVED PROPERTIES OF THE ities at large distances from the moon’s orbit. BRIGHT STREAKS IN THE E RING Forutnately, since we are only concerned with In order to compare the locations and bright- the particle distribution near Enceladus’ orbit nesses of the bright streaks seen in Figure1 with here, and particles on large eccentricities only the expected properties of rebound wakes, we contribute to the relatively uniform background re-processed the data from both these images. density in this region, such inaccuracies should Table1 provides a brief summary of the ob- not significantly affect the intensity of the wake servation parameters for both of these images, signatures in these simulations. Also, as before, including the phase angle, resolution and the we deliberately chose semi-major axis distribu- inertial longitudes (measured relative to the as- tions that would contain sufficient particles on cending node of the rings on the J2000 equator suitable orbits to produce wake signatures. plane) of the Sun, Enceladus, and Enceladus’ Figure5 shows the resulting density maps orbital pericenter. We did examine several other for these simulations. As expected, the over- images of this region, including those analyzed all brightness trends in these simulations are by Mitchell et al.(2015), but found that those smoother than the ones shown in Figure3, images either did not capture a wide enough and the differences between the three simula- range of longitudes or had much lower signal- tions are more subtle. More importantly, a to-noise, and so they were not included in this wake signature is faintly visible in the simu- investigation. lation with the narrowest semi-major axis dis- The raw data for both images were calibrated tribution. This signal can be seen much more using the standard CISSCAL routines (Porco clearly in Figure6, which shows the same simu- et al. 2004; West et al. 2010), and were geometri- lations after high-pass filtered by subtracting a cally navigated based on the locations of known version of the map smoothed over 4 Hill Radii. stars in the field of view. The data from both In these filtered maps, the simulation with the images were then re-projected to produce maps 11

Table 1. Parameters for images used in this study.

Image ID Filter Date Phase Angle Sub-solar Enceladus Enceladus Orbital Pixel Scale (degrees) Longitude (deg.) Longitude (deg.) Pericenter Long. (deg.) (km/pixel) W1537013154 Clear 2006-258T11:34:10 174.4 -169.8 -72.1 -81.6 141.9 W1752971847 Red 2013-200T23:41:28 169.7 -91.7 -12.3 44.9 72.5

Figure 7. Re-projected images of the E-ring structures around Enceladus. Note that both longitude and radius are measured in Hill Radii and the origin is located at Enceladus. Unlike Figure1, these images are neither spatially filtered nor brightness inverted (although they are independently stretched), and so regions of higher particle density appear brighter than the background. For this analysis we are particularly interested in the feature centered at a radius of -2 Hill Radii and a longitude of +20 Hill Radii. of the ring’s projected brightness as a function the observed locations of these bright bands are of radius and inertial longitude. Finally, we inconsistent with predicted locations of passing converted the relevant coordinates to radial and satellite wakes shown in Figure2. However, the azimuthal distances from Enceladus in units of first and most prominent rebound wake extends Hill Radii in order to facilitate comparisons with between 10 and 20 Hill Radii in azimuth and the predicted wake locations, between -1 and -3 Hill Radii in Radius, and so Figure7 shows the resulting maps of the E this feature occurs quite close to the observed ring in the vicinity of Enceladus on a com- location of a bright band in both images. mon scale. Note that these images are neither In order to more precisely compare the ob- brightness-inverted nor high-pass filtered, so re- served locations of these bright bands with the gions of higher density appear as bright streaks expected locations of the first rebound wake, in these images. The feature Mitchell et al. we used the following procedures to trace the (2015) identified as the interior tendril appears locations of the bright streaks in both maps. as a diagonal streak just below Enceledus, while First, we generated a high-pass filtered version the first leading tendril is the feature centered of each map by subtracting a version of the rel- around 2 Hill Radii interior to Enceladus’ orbit evant map smoothed over a spatial scale of 4000 and 20 Hill Radii ahead of the moon. Consis- tent with prior findings by Mitchell et al.(2015), 12 km, comparable to that used in Figure6. 1 As behind the wake’s location. Furthermore, the shown in the top panels of Figures8 and9, this fact that this feature’s position shifts between filtering removed large-scale brightness trends the two observations is inconsistent with the but preserved the signals from the bands. Next, expected behavior of satellite wakes. Hence, for each longitude in each map, we identified these structures are most likely tendrils consist- the radius of the brightness maximum within a ing of material recently launched from Ence- selected range of radii that contained the inte- ladus, whose precise positions can change de- rior and leading tendrils. The middle panels of pending on the locations of the active sources, Figures8 and9 show the resulting wake loca- the assumed launch velocity of the particles, tions overlaid on the high-pass-filtered images, and the observation geometry (Mitchell et al. demonstrating that this process properly traces 2015). the location of the structures. We also quan- At the same time, it is important to note that tified the brightness variations associated with there are no obvious brightness enhancements these structures in terms of the peak fractional at the predicted locations of the first passing brightness variation at each longitude for each wake in image W1537013154 (see Figure8), feature in both images. The peak fractional even though the observed tendril is well dis- brightness variation is simply the ratio of the placed from the expected wake position. The brightness of the peak in the high-pass-filtered lack of any observable wake signatures implies image to the brightness of the unfiltered image, that there are not sufficient particles on suitable and these numbers are plotted in Figure 10. orbits to form detectable interior passing wakes in the E ring (cf. Mitchell et al. 2015). 4. COMPARISONS BETWEEN OBSERVED However, if we instead consider the bright BRIGHT STREAKS AND PREDICTED streaks lying between 10 and 30 Hill Radii ahead PROPERTIES OF REBOUND WAKES of the moon, the situation is more complicated. 4.1. Locations of bright streaks overlap with The material more than 20 Hill Radii in front of the expected positions of rebound wakes Enceladus falls at the predicted locations of the The bottom panels of Figures8 and9 compare leading tendrils (Mitchell et al. 2015), and these the observed positions of the bright streaks in parts of the streaks extend beyond the expected the two maps with the predicted locations for locations of rebound wakes, so the distal parts both the passing and rebound wakes. of these streaks are almost certainly tendrils. If we first consider the bright streaks directly However, the portions of these bright bands that interior to Enceladus, it is clear that their lo- lie around 15-20 Hill Radii in front of Enceladus cations are inconsistent with those expected for and around 2 Hill Radii interior to Enceladus’ any satellite wakes. For image W1537013154, orbit are located almost exactly where the first the bright band lies well behind the expected rebound wake should occur. Furthermore, in location of the wake, and even extends slightly image W1537103154, there appears to be two behind Enceladus itself. For W1725971847 the overlapping streaks in this region (see Figure8), bright feature is closer to the expected location one situated between 10 and 20 Hill Radii in of the wake, but is still roughly 2 Hill Radii azimuth and around -1.5 Hill Radii in radius, while the other extends between 20 and 30 Hill 1 This smoothing scale is larger than the one used Radii in azimuth and between -3 and -4 Hill for the simulations in Figure4 because the background Radii in radius. The lower feature extends over brightness variations in the ring are on a larger scale due a range of radii similar to the tendril below the to the broader range of eccentricities. 13

Figure 8. Comparisons of the observed positions Figure 9. Comparisons of the observed locations of features in image W1537013154 with those ex- of features in image W1752971847 with those ex- pected from passing and rebound wakes. The top pected from passing and rebound wakes. The top panel shows a high-pass filtered version of the re- panel shows a high-pass filtered version of the re- projected image, where the brightness variations projected image, where the brightness variations can be more clearly seen (brightness scale runs can be more clearly seen (brightness scale runs from from ±10%). The second panel shows the same ±10%). The second panel shows the same image image with the estimated locations of the bright- with the estimated locations of the brightness max- ness maxima marked as lime crosses and green dia- ima marked as orange crosses and red diamonds. monds. The bottom panel overlays these same posi- The bottom panel overlays these same positions tions with the trajectories predicted by Hill’s Equa- with the trajectories predicted by Hill’s Equations. tions. Note that the structure connected to Ence- Note that the structure connected to Enceladus is ladus is displaced several Hill Radii to the left of displaced a couple of Hill Radii to the left of the the first passing wake, the second brightness max- first passing wake, the second brightness maxima ima lie very close to the expected location of the lie very close to the expected location of the first first rebound wake, although it does appear to ex- rebound wake, although it does appear to extend tend beyond that wake. beyond that wake. 14 moon, and corresponds to the portion of this tions will prevent them from producing the ob- feature that was reproduced reasonably well as served small-scale density variations. We there- a tendril by Mitchell et al.(2015) (see Figure1). fore need to evaluate whether the brightness of By contrast, the upper part of this band was the these features is inconsistent with reasonable es- portion of this feature that those authors were timates of the orbital element distribution of the not able to reproduce in their simulations. Thus E ring. it is reasonable to consider the possibility that As shown in Figure 10, these features are actu- the lower feature is the actual leading tendril, ally quite subtle, representing variations in the while the upper feature is a rebound wake. The background E-ring density of only about 5%. structure in image W1752971847 cannot be so The fractional brightness variations of the rele- easily separated into two distinct components, vant feature in the images are therefore roughly but the fact that the portion of the band be- comparable to the fractional density variations tween 10 and 20 Hill Radii is at very similar shown in the simulations in Figures4 and6. locations in both images is consistent with the Thus the features in the wakes could potentially expected behavior of a satellite wake. be wake signatures so long as a large enough One further complication in interpreting this fraction of the particles in this part of the E feature is that the part of the feature closest ring have the requisite orbital elements. to Enceladus is not in exactly the same loca- As mentioned in Section2 above, the best con- tion in the two images. In the earlier image strained aspect of the particles’ orbital proper- W1537013154, the end of the band curves in- ties is the eccentricity distribution, and the sim- wards towards the planet and so falls interior to ulations discussed above demonstrate that this the expected wake location between 10 and 15 broad eccentricity distribution allows for wakes Hill Radii in front of the moon. By contrast in with the required intensity to form. Also, the the W1752971847 image the bright band curves distributions of pericenter locations and incli- outwards, and so falls exterior to the expected nations are also unlikely to strongly affect the wake location in this region and in fact extends intesnity of the wake signatures. While there all the way to Enceladus’ orbit. These varia- are longitudinal variations in the shape of the tions pose challenges to both the tendril and re- E-ring profile that could be due to asymmetries bound wake models, and will likely require more in the pericenter direction, these are relatively detailed simulations to explain. subtle and are tied to the direction of the Sun rather than the location of Enceladus (Hedman 4.2. The brightness of these features could be et al. 2012; Ye et al. 2016). Hence there are consistent with the E-ring’s orbital probably not strong asymmetries in the distri- element distribution bution of pericenter locations. In addition, even The other main argument made by Mitchell though the E-ring particles do have a range of et al.(2015) against the idea that these fea- orbital inclinations, this parameter is strongly tures could be satellite wakes is that such wakes correlated with the eccentricity and so it is rea- only form among particles that have nearly the sonable to expect that particles with low eccen- same eccentricity and pericenter as Enceladus, tricities also have low inclinations. Hence it is which is a small fraction of the particles in the probably reasonable to neglect orbital inclina- E ring. As discussed in Section2 above, par- tions for this preliminary investigation. ticles with epicylic motions relative to Ence- By contrast, the particles’ orbital semi-major ladus larger than around 1 Hill Radius are un- axis distribution certainly strongly affects the likely to form wakes because their radial mo- 15

Figure 10. Plots showing the peak fractional brightness of the bright streaks in the two images as a function of radial displacement from Enceladus’ orbit. Note the color code matches that used in Figures8 and9. intensity of these sorts of wake signatures. Only stantial fraction of the E-ring particles on low- particles orbiting between 1 and 1.5 Hill Radii eccentricity orbits are on horseshoe orbits. On exterior to Enceladus can produce the required the other hand, the same simulations in Fig- rebound wake, so we need a substantial frac- ures4 and6 that have a strong rebound wake tion of the particles on low-eccentricity orbits also tend to have a ⊂-shaped feature that traces to have their semi-major axes in this range. At the average trajectory of the same particles re- the moment, there are no measurements of the sponsible for the wake. This ⊂-shaped feature semi-major axis distribution of particles close is not clearly visible in the images, which could to Enceladus’ orbit with high enough fidelity to place an upper limit on the fraction of parti- reliably estimate the fraction of particles that cles on orbits that could produce a rebound would be on these wake-producing orbits. wake. However, initial investigations show that In fact, the relative intensities of small-scale the relative strength of the wake signature and E-ring structures like these wakes probably pro- the ⊂-shaped feature depends on the details of vide the best opportunities to constrain the de- both the semi-major axis and eccentricity dis- tailed shape of the E-ring particles’ orbital pa- tributions. A broad range of orbital parameter rameter distribution. For example, the lack of distributions will therefore probably need to be obvious signatures of passing wakes indicates investigated to ascertain whether particle popu- that there are comparatively few particles on lations capable of produced detectable rebound low-eccentricity orbits with semi-major axes be- wakes are consistent with other aspects of the tween 2.5 and 4 Hill radii from Enceladus’ or- E-ring’s structure. bit, which could support the idea that a sub- 16

5. DISCUSSION AND CONCLUSIONS Hedman et al. 2013). By symmetry, these par- The above analysis indicates that while there ticles would re-encounter Enceladus when their is no evidence for passing satellite wakes in the semi-major axes is about 2 Hill Radii exterior E ring, rebound wakes could contribute to the to Enceladus’ orbit. This would naturally cause bright features observed 10-20 Hill Radii in front many of these particles to approach Enceladus of Enceladus. Even though the distal ends of along semi-major axes close to those needed these features are most likely tendrils of mate- to generate rebound wakes. Of course, this rial recently launched from Enceladus (Mitchell process could also deliver particles onto orbits et al. 2015), the bright bands around 15 Hill that could produce passing wakes exterior to Radii in front of Enceladus are in the right lo- Enceladus’ orbit. Quantitatively evaluating the cations to be due to pre-existing E-ring particles expected strength of these different classes of following horseshoe orbits. Furthermore, these wakes willrequire sophisticated numerical simu- features are sufficiently subtle that they could lations including a variety of non-gravitational be generated by the small fraction of E-ring forces. particles that have suitable near-circular orbits. A larger suite of numerical simulations of ma- Properly quantifying and/or limiting the con- terial launched from Enceladus would also clar- tributions from rebound wakes to these features ify whether this material can ever produce fea- will require additional numerical simulations of tures as close to Enceladus as observed in these plume and E-ring material, as well as more de- images. The fact that the Mitchell et al.(2015) tailed analysis of the available data. simulations did not produce tendrils extending While such detailed investigations are beyond as close to Enceladus as the observations sug- the scope of this initial study, we can at least gests that, if these features are not rebound point out a potential mechanism that would wakes, then their morphology will place strong allow a large fraction of the particles on low- constraints on the particle launch velocity dis- eccentricity orbits to approach Enceladus on tribution from Enceladus. At the same time, trajectories that could produce the rebound simulations of various populations of E-ring ma- wake signature. The tendrils found in the terial are needed to determine which orbital Mitchell et al.(2015) simulations correspond to element distributions could produce structures material with semi-major axes between those of comparable in brightness to those observed, horseshoe and passing orbits, or around 2 Hill while at the same time avoiding producing other Radii interior and exterior to Enceladus’ orbit. structures in the rings (like signals from the After being launched from Enceladus, these par- passing wakes or ⊂-shaped patterns) that are ticles are expected to migrate outward due to too intense to be consistent with the data. Since plasma drag and other non-gravitational pro- only particles with a narrow range of orbital pa- cesses (Burns 1976; Gruen 1984; Hor´anyi et al. rameters can produce each of these structures, 1992; Juh´aszand Hor´anyi 2002, 2004; Juh´asz determining the fraction of particles with these et al. 2007; Hor´anyi et al. 2008). This outward properties should place strong constraints on migration will cause material initially launched the evolution rates of the particles’ semi-major Saturn-wards from Enceladus to approach and axes and eccentricities. then move past Enceladus’ orbit, so while these particles would first drift ahead of the moon, they should eventually reverse direction and drift back towards the moon (Murray 1994; 17

ACKNOWLEDGEMENTS NNX15AQ67G. The authors wish to thank C. This work was supported in part by NASA Mitchell and C. Porco for their helpful discus- Cassini Data Analysis Program grant number sions.

REFERENCES Burns, J. A., 1976. Elementary derivation of the Henon, M., Petit, J. M., 1986. Series Expansions perturbation equations of celestial mechanics. for Encounter-Type Solutions of Hill’s Problem. American Journal of Physics 44, 944–949. Celestial Mechanics 38 (1), 67–100. Burns, J. A., Hamilton, D. P., Showalter, M. R., Hor´anyi, M., Burns, J. A., Hamilton, D. P., 1992. 2001. Dusty rings and circumplanetary dust: The dynamics of Saturn’s E ring particles. Observations and simple physics. In: Gr¨un, E., Icarus 97, 248–259. Gustafson, B., Dermott, S., Fecthig, H. (Eds.), Hor´anyi, M., Juh´asz,A., Morfill, G. E., 2008. Interplanetary Dust, pp. 641–725. Springer. Large-scale structure of Saturn’s E-ring. Geophys. Res. Lett. 35, 4203. Gruen, E., 1984. Dynamics of Charged Small Horn, L. J., Showalter, M. R., Russell, C. T., Particles. In: Brahic, A. (Ed.), Planetary 1996. Detection and Behavior of Pan Wakes in Rings, pp. 535. Saturn’s A Ring. Icarus 124 (2), 663–676. Hamilton, D. P., 1993a. −Motion of dust in a Juh´asz,A., Hor´anyi, M., 2002. Saturn’s E ring: A planetary magnetosphere - Orbit-averaged dynamical approach. Journal of Geophysical equations for oblateness, electromagnetic, and Research (Space Physics) 107, 1066. radiation forces with application to Saturn’s E Juh´asz,A., Hor´anyi, M., 2004. Seasonal ring. Icarus 101, 244–264. variations in Saturn’s E-ring. Hamilton, D. P., 1993b. Erratum: Motion of dust Geophys. Res. Lett. 31, L19703. in a planetary magnetosphere: Orbit-averaged Juh´asz,A., Hor´anyi, M., Morfill, G. E., 2007. equations for oblateness, electromagnetic, and Signatures of Enceladus in Saturn’s E ring. radiation forces with application to Saturn’s Geophys. Res. Lett. 34, L09104. E-ring. Icarus 103, 161. Julian, W. H., Toomre, A., 1966. Hamilton, D. P., Burns, J. A., 1994. Origin of Non-Axisymmetric Responses of Differentially Saturn’s E ring: Self-sustained, naturally. Rotating Disks of Stars. ApJ 146, 810. Science 264, 550–553. Kempf, S., Beckmann, U., Schmidt, J., 2010. How Hedman, M. M., Burns, J. A., Hamilton, D. P., the Enceladus dust plume feeds Saturn’s E ring. Showalter, M. R., 2012. The three-dimensional Icarus 206, 446–457. structure of Saturn’s E ring. Icarus 217 (1), Kempf, S., Hor´anyi, M., Hsu, H. W., Hill, T. W., Juh´asz,A., Smith, H. T., 2018. Saturn’s 322–338. Diffuse E Ring and Its Connection with Hedman, M. M., Burns, J. A., Hamilton, D. P., Enceladus, pp. 195. Showalter, M. R., 2013. Of horseshoes and Lin, D. N. C., Papaloizou, J., 1979. Tidal torques heliotropes: Dynamics of dust in the Encke on accretion discs in binary systems with Gap. Icarus 223 (1), 252–276. extreme mass ratios. MNRAS 186, 799–812. Hedman, M. M., Helfenstein, P., Chancia, R. O., Mitchell, C. J., Porco, C. C., Weiss, J. W., 2015. Thomas, P., Roussos, E., Paranicas, C., Tracking the Geysers of Enceladus into Saturn’s Verbiscer, A. J., 2020. Photometric Analyses of E Ring. AJ 149 (5), 156. Saturn’s Small Moons: Aegaeon, Methone, and Murray, C. D., 1994. Dynamical Effects of Drag Pallene Are Dark; Helene and Calypso Are in the Circular Restricted Three-Body Problem. Bright. AJ 159 (4), 129. I. Location and Stability of the Lagrangian Hedman, M. M., Postberg, F., Hamilton, D. P., Equilibrium Points. Icarus 112 (2), 465–484. Renner, S., Hsu, H. W., 2018. Dusty Rings, pp. Murray, C. D., Dermott, S. F., 2001. Solar 308–337. System Dynamics. Cambridge University Press. 18

Porco, C. C., West, R. A., Squyres, S., McEwen, Tiscareno, M. S., Burns, J. A., Nicholson, P. D., A., Thomas, P., Murray, C. D., Del Genio, A., Hedman, M. M., Porco, C. C., 2007. Cassini Ingersoll, A. P., Johnson, T. V., Neukum, G., imaging of Saturn’s rings. II. A wavelet Veverka, J., Dones, L., Brahic, A., Burns, J. A., technique for analysis of density waves and Haemmerle, V., Knowles, B., Dawson, D., other radial structure in the rings. Icarus 189 Roatsch, T., Beurle, K., Owen, W., 2004. (1), 14–34. Cassini imaging science: Instrument characteristics and anticipated scientific West, R., Knowles, B., Birath, E., Charnoz, S., di investigations at Saturn. SSRv 115, 363–497. Nino, D., Hedman, M., Helfenstein, P., McEwen, Schenk, P. M., Clark, R. N., Howett, C. J. A., A., Perry, J., Porco, C., Salmon, J., Throop, H., Verbiscer, A. J., Waite, J. H., 2018. Enceladus Wilson, D., 2010. In-flight calibration of the and the Icy Moons of Saturn. Cassini imaging science sub-system cameras. Showalter, M. R., Cuzzi, J. N., Larson, S. M., Planet. Space Sci. 58, 1475–1488. 1991. Structure and particle properties of Ye, S. Y., Gurnett, D. A., Kurth, W. S., 2016. Saturn’s E Ring. Icarus 94, 451–473. In-situ measurements of Saturn’s dusty rings Showalter, M. R., Cuzzi, J. N., Marouf, E. A., based on dust impact signals detected by Esposito, L. W., 1986. Satellite “Wakes” and Cassini RPWS. Icarus 279, 51–61. the Orbit of the Encke Gap Moonlet. In: BAAS, Volume 18, pp. 839.