Three tenuous rings/arcs for three tiny moons

Using Cassini images, we examine the faint material along the orbits of Methone, Anthe and Pallene, three small moons that reside between the orbits of Mimas and Enceladus. A continuous ring of material covers the orbit of Pallene; it is visible at extremely high phase angles and appears to be localized vertically to within ±25 km of Pallene's inclined orbit. By contrast, the material associated with Anthe and Methone appears to lie in longitudinally confined arcs. The Methone arc extends over ∼10° in longitude around the satellite's position, while the Anthe arc reaches ∼20° in length. The extents of these arcs are consistent with their confinement by nearby corotation eccentricity resonances with Mimas. Anthe has even been observed to shift in longitude relative to its arc in the expected manner given the predicted librations of the moon.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

Aegaeon (Saturn LIII), a G-ring object

Aegaeon (Saturn LIII, S/2008 S1) is a small satellite of Saturn that orbits within a bright arc of material near the inner edge of Saturn’s G-ring. This object was observed in 21 images with Cassini’s Narrow-Angle Camera between June 15 (DOY 166), 2007 and February 20 (DOY 051), 2009. If Aegaeon has similar surface scattering properties as other nearby small saturnian satellites (Pallene, Methone and Anthe), then its diameter is approximately 500 m. Orbit models based on numerical integrations of the full equations of motion show that Aegaeon’s orbital motion is strongly influenced by multiple resonances with Mimas. In particular, like the G-ring arc it inhabits, Aegaeon is trapped in the 7:6 corotation eccentricity resonance with Mimas. Aegaeon, Anthe and Methone therefore form a distinctive class of objects in the Saturn system: small moons in corotation eccentricity resonances with Mimas associated with arcs of debris. Comparisons among these different ring-arc systems reveal that Aegaeon’s orbit is closer to the exact resonance than Anthe’s and Methone’s orbits are. This could indicate that Aegaeon has undergone significant orbital evolution via its interactions with the other objects in its arc, which would be consistent with the evidence that Aegaeon’s mass is much smaller relative to the total mass in its arc than Anthe’s and Methone’s masses are.     

Prometheus and Pandora: masses and orbital positions during the Cassini tour

Hubble Space Telescope (HST) images of Prometheus and Pandora show longitude discrepancies of about 20° with respect to the Voyager ephemerides, with an abrupt change in mean motion at the end of 2000 (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07). These discrepancies are anti-correlated and arise from chaotic interactions between the two moons, occurring at interval of 6.2 yr, when their apses are anti-aligned (Goldreich and Rappaport, 2003a, Icarus 162, 391-399). This behavior is attributed to the overlap of four 121:118 apse-type mean motion resonances (Goldreich and Rappaport, 2003b, Icarus 166, 320-327). We study the Prometheus-Pandora system using a Radau-type integrator taking into account Saturn's oblateness up to and including terms in J₆, plus the effects of the major satellites. We first confirm the chaotic behavior of Prometheus and Pandora. By fitting the numerical integrations to the HST data (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07), we derive the satellite masses. The resulting GM values (with their standard 3-σ errors) for Prometheus and Pandora are respectively and . Using the nominal shape of the two moons (Thomas, 1989, Icarus 77, 248-274), we derive Prometheus and Pandora's densities, 0.40^+0.03_−0.07 and 0.49^+0.05_−0.09 g cm⁻³, respectively. Our numerical fits also enable us to constrain the time of the latest apse anti-alignment in 2000. Finally, using our fit, we predict the orbital positions of the two satellites during the Cassini tour, and provide a lower limit of the uncertainties due to chaos. These uncertainties amount to about 0.2° in mean longitude at the arrival of the Cassini spacecraft in July 2004, and to about 3° in 2008, at the end of the nominal tour.     

Search for and limits on plume activity on Mimas, Tethys, and Dione with the Cassini Visual Infrared Mapping Spectrometer (VIMS)

Cassini Visual Infrared Mapping Spectrometer (VIMS) observations of Mimas, Tethys, and Dione obtained during the nominal and extended missions at large solar phase angles were analyzed to search for plume activity. No forward scattered peaks in the solar phase curves of these satellites were detected. The upper limit on water vapor production for Mimas and Tethys is one order of magnitude less than the production for Enceladus. For Dione, the upper limit is two orders of magnitude less, suggesting this world is as inert as Rhea (Pitman, K.M., Buratti, B.J., Mosher, J.A., Bauer, J.M., Momary, T., Brown, R.H., Nicholson, P.D., Hedman, M.M. [2008]. Astrophys. J. Lett. 680, L65-L68). Although the plumes are best seen at ∼2.0 μm, Imaging Science Subsystem (ISS) Narrow Angle Camera images obtained at the same time as the VIMS data were also inspected for these features. None of the Cassini ISS images shows evidence for plumes. The absence of evidence for any Enceladus-like plumes on the medium-sized saturnian satellites cannot absolutely rule out current geologic activity. The activity may below our threshold of detection, or it may be occurring but not captured on the handful of observations at large solar phase angles obtained for each moon. Many VIMS and ISS images of Enceladus at large solar phase angles, for example, do not contain plumes, as the active “tiger stripes” in the south pole region are pointed away from the spacecraft at these times. The 7-year Cassini Solstice Mission is scheduled to gather additional measurements at large solar phase angles that are capable of revealing activity on the saturnian moons.     

The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

How much oxygen is too much? Constraining Saturn's ring atmosphere

The discovery of a molecular oxygen atmosphere around Saturn's rings has important implications for the electrodynamics of the ring system. Its existence was inferred from the Cassini in situ detection of molecular oxygen ions above the rings during Saturn Orbit Insertion in 2004. Molecular oxygen is difficult to observe remotely, and theoretical estimates have yielded only a lower limit (Nn?¹⁰¹³ cm⁻²) to the O₂ column density. Comparison with observations has previously concerned matching ion densities at spacecraft altitudes far larger than the scale height of the neutral atmosphere. This is further complicated by charged particle propagation effects in Saturn's offset magnetic field. In this study we adopt a complementary approach, by focusing on bulk atmospheric properties and using additional aspects of the Cassini observations to place an upper limit on the column density. We develop a simple analytic model of the molecular atmosphere and its photo-ionization and dissociation products, with Nn a free parameter. Heating of the neutrals by viscous stirring, cooling by collisions with the rings, and torquing by collisions with pickup ions are all included in the model. We limit the neutral scale height to h?3000 km using the INMS neutral density nondetection over the A ring. A first upper limit to the neutral column is derived by using our model to reassess O₂ production and loss rates. Two further limits are then obtained from Cassini observations: corotation of the observed ions with the planet implies that the height-integrated conductivity of the ring atmosphere is less than that of Saturn's ionosphere; and the nondetection of fluorescent atomic oxygen over the rings constrains the molecular column from which it is produced via photo-dissociation. These latter limits are independent of production and loss rates and are only weakly dependent on temperature. From the three independent methods described, we obtain similar limits: Nn?2×¹⁰¹⁵ cm⁻². The mean free path for collisions between neutrals thus cannot be very much smaller than the scale height.     

Disk-integrated bolometric Bond albedos and rotational light curves of saturnian satellites from Cassini Visual and Infrared Mapping Spectrometer

We present values from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) of four fundamental disk-integrated spectrophotometric properties (bolometric Bond albedo, solar phase curve, phase integral, and geometric albedo at 7-15 different wavelengths in the λ = 0.35-5.1 μm range) for five mid-sized saturnian icy satellites: Rhea, Dione, Tethys, Mimas, and Enceladus. These values, which include data from the period 2004-2008 and add to past VIMS phase curves, include opposition surge effects at down to fractions of a degree in solar phase angle for several moons and extend to over double the solar phase angle coverage of the Voyager mission. We also present new rotational light curves for Rhea and Dione at 7 near-infrared bands not previously available in ground-based or spacecraft studies. The bolometric Bond albedos we derive are as follows: 0.48 ± 0.09 (Rhea), 0.52 ± 0.08 (Dione), 0.61 ± 0.09 (Tethys), 0.67 ± 0.10 (Mimas), and 0.85 ± 0.11 (Enceladus). We also provide breakdowns of the major photometric quantities in both leading and trailing hemispheres. These refined parameters can be used to construct future bolometric Bond albedo maps that will contribute to surface composition identification studies, as well as models of volatile transport and sublimation. Through such applications, these data will help to determine the physical properties of surface particles, how the E-ring affects the inner saturnian moons, what is responsible for the dark albedo patterns seen on Tethys, and if these moons (e.g., Dione) are geologically active.     

The evolving flow of Jupiter’s White Ovals and adjacent cyclones

We present results regarding the dynamical meteorology of Jupiter’s White Ovals at different points in their evolution. Starting from the era with three White Ovals FA, BC, and DE (Galileo), continuing to the post-merger epoch with only one Oval BA (Cassini), and finally to Oval BA’s current reddened state (New Horizons), we demonstrate that the dynamics of their flow have similarly evolved along with their appearance. In the Galileo epoch, Oval DE had an elliptical shape with peak zonal wind speeds of ∼90 m s⁻¹ in both its northern and southern peripheries. During the post-merger epoch, Oval BA’s shape was more triangular and less elliptical than Oval DE; in addition to widening in the north-south direction, its northern periphery was 20 m s⁻¹ slower, and its southern periphery was 20 m s⁻¹ faster than Oval DE’s flow during the Galileo era. Finally, in the New Horizons era, the reddened Oval BA had evolved back to a classical elliptical form. The northern periphery of Oval BA increased in speed by 20 m s⁻¹ from Cassini to New Horizons, ending up at a speed nearly identical to that of the northern periphery of Oval DE during Galileo. However, the peak speeds along the southern rim of the newly formed Oval BA were consistently faster than the corresponding speeds in Oval DE, and they increased still further between Cassini and New Horizons, ending up at ∼140-150 m s⁻¹. Relative vorticity maps of Oval BA reveal a cyclonic ring surrounding its outer periphery, similar to the ring present around the Great Red Spot. The cyclonic ring around Oval BA in 2007 appears to be moderately stronger than observed in 1997 and 2001, suggesting that this may be associated with the coloration of the vortex. The modest strengthening of the winds in Oval BA, the appearance of red aerosols, and the appearance of a turbulent, cyclonic feature to Oval BA’s northwest create a strong resemblance with the Great Red Spot from both a dynamical and morphological perspective.In addition to the White Ovals, we also measure the winds within two compact cyclonic regions, one in the Galileo data set and one in the Cassini data set. In the images, these cyclonic features appear turbulent and filamentary, but our wind field reveals that the flow manifests as a coherent high-speed collar surrounding relatively quiescent interiors. Our relative vorticity maps show that the vorticity likewise concentrates in a collar near the outermost periphery, unlike the White Ovals which have peak relative vorticity magnitudes near the center of the vortex. The cyclones contain several localized bright regions consistent with the characteristics of thunderstorms identified in other studies. Although less studied than their anticyclonic cousins, these cyclones may offer crucial insights into the planet’s cloud-level energetics and dynamical meteorology.     

Revisiting the thermal inertia of Iapetus: Clues to the thickness of the dark material

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m⁻² K⁻¹ s^−1/2 while the light material varies between 15 and 25 J m⁻² K⁻¹ s^−1/2. Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.     

Structure of self-gravity wakes in Saturn's A ring as measured by Cassini CIRS

The CIRS infrared spectrometer onboard the Cassini spacecraft has scanned Saturn's A ring azimuthally from several viewing angles since its orbit insertion in 2004. A quadrupolar asymmetry has been detected in this ring at spacecraft elevations ranging between 16° to 37°. Its fractional amplitude decreases from 22% to 8% from 20° to 37° elevations. The patterns observed in two almost complete azimuthal scans at elevations 20° and 36° strongly favor the self-gravity wakes as the origin of the asymmetry. The elliptical, infinite cylinder model of Hedman et al. [Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D., Buratti, B.J., Baines, K.H., Brown, R.H., Clark, R.N., 2007. Astron. J. 133, 2624-2629] can reproduce the CIRS observations well. Such wakes are found to have an average height-to-spacing ratio H/λ=0.1607±0.0002, a width-over-spacing W/λ=0.3833±0.0008. Gaps between wakes, which are filled with particles, have an optical depth τG=0.1231±0.0005. The wakes mean pitch angle ΦW is 70.70°±0.07°, relative to the radial direction. The comparison of ground-based visible data with CIRS observations constrains the A ring to be a monolayer. For a surface mass density of 40 g cm⁻² [Tiscarino, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34], the expected spacing of wakes is λ≈60 m. Their height and width would then be H≈10 m and W≈24 m, values that match the maximum size of particles in this ring as determined from ground-based stellar occultations [French, R.G., Nicholson, P.D., 2000. Icarus 145, 502-523].     

Detection of C2HD and the D/H ratio on Titan

We report here the first detection of mono-deuterated acetylene (acetylene-d1, C₂HD) in Titan's atmosphere from the presence of two of its emission bands at 678 and 519 cm⁻¹ as observed in CIRS spectral averages of nadir and limb observations taken between July 2004 and mid-2007. By using new laboratory spectra for this molecule, we were able to derive its abundance at different locations over Titan's disk. We find the C₂HD value () to be roughly constant with latitude from the South to about 45° N and then to increase slightly in the North, as is the case for C₂H₂. Fitting the 678 cm⁻¹ν₅ band simultaneously with the nearby C₂H₂ 729 cm⁻¹ν₅ band, allows us to infer a D/H ratio in acetylene on Titan with an average of the modal values of 2.09±0.45×¹⁰⁻⁴ from the nadir observations, the uncertainties being mainly due to the vertical profile used for the fit of the acetylene band. Although still subject to significant uncertainty, this D/H ratio appears to be significantly larger than the one derived in methane from the CH₃D band (upper limit of 1.5×¹⁰⁻⁴; Bézard, B., Nixon, C.A., Kleiner, I., Jennings, D.E., 2007. Icarus, 191, 397-400; Coustenis, A., Achterberg, R., Conrath, B., Jennings, D., Marten, A., Gautier, D., Bjoraker, G., Nixon, C., Romani, P., Carlson, R., Flasar, M., Samuelson, R.E., Teanby, N., Irwin, P., Bézard, B., Orton, G., Kunde, V., Abbas, M., Courtin, R., Fouchet, Th., Hubert, A., Lellouch, E., Mondellini, J., Taylor, F.W., Vinatier, S., 2007. Icarus 189, 35-62). From the analysis of limb data we infer D/H values of (at 54° S), (at 15° S), (at 54° N) and (at 80° N), which average to a mean value of 1.63±0.27×¹⁰⁻⁴.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.     

Interpretation of Solar Magnetic Field Strength Observations

This study based on longitudinal Zeeman effect magnetograms and spectral line scans investigates the dependence of solar surface magnetic fields on the spectral line used and the way the line is sampled to estimate the magnetic flux emerging above the solar atmosphere and penetrating to the corona from magnetograms of the Mt. Wilson 150-foot tower synoptic program (MWO). We have compared the synoptic program λ5250 Å line of Fe?i to the line of Fe?i at λ5233 Å since this latter line has a broad shape with a profile that is nearly linear over a large portion of its wings. The present study uses five pairs of sampling points on the λ5233 Å line. Line profile observations show that the determination of the field strength from the Stokes V parameter or from line bisectors in the circularly polarized line profiles lead to similar dependencies on the spectral sampling of the lines, with the bisector method being the less sensitive. We recommend adoption of the field determined with the line bisector method as the best estimate of the emergent photospheric flux and further recommend the use of a sampling point as close to the line core as is practical. The combination of the line profile measurements and the cross-correlation of fields measured simultaneously with λ5250 Å and λ5233 Å yields a formula for the scale factor δ ^?1 that multiplies the MWO synoptic magnetic fields. By using ρ as the center-to-limb angle (CLA), a fit to this scale factor is δ ^?1=4.15?2.82sin?²(ρ). Previously δ ^?1=4.5?2.5sin?²(ρ) had been used. The new calibration shows that magnetic fields measured by the MDI system on the SOHO spacecraft are equal to 0.619±0.018 times the true value at a center-to-limb position 30°. Berger and Lites (2003, Solar Phys. 213, 213) found this factor to be 0.64±0.013 based on a comparison using the Advanced Stokes Polarimeter.     

Appearance of Saturn’s F ring azimuthal channels for the anti-alignment configuration between the ring and Prometheus

In this article we explore the aspect of the F ring with respect to the anti-alignment configuration between the ring and Prometheus. We focus our attention on the shape of the F ring’s azimuthal channels which were first reported by Porco et al. (Porco, C.C., Baker, E., Barbara, J., Beurle, K., Brahic, A., Burns, J.A., Charnoz, S., Cooper, N., Dawson, D.D., Del Genio, A.D., Denk, T., Dones, L., Dyudina, U., Evans, M.W., Giese, B., Grazier, K., Helfenstein, P., Ingersoll, A.P., Jacobson, R.A., Johnson, T.V., McEwen, A., Murray, C.D., Neukum, G., Owen, W.M., Perry, J., Roatsch, T., Spitale, J., Squyres, S., Thomas, P., Tiscareno, M., Turtle, E., Vasavada, A.R., Veverka, J., Wagner, R., West, R. [2005] Science, 307, 1226-1236) and numerically explored by Murray et al. (Murray, C.D., Chavez, C., Beurle, K., Cooper, N., Evans, M.W., Burns, J.A., Porco, C.C. [2005] Nature 437, 1326-1329) who found excellent agreement between Cassini’s ISS reprojected images and their numerical model via a direct comparison. We find that for anti-alignment the channels are wider and go deeper inside the ring material. From our numerical model we find a new feature, an island in the middle of the channel. This island is made up of the particles that have been perturbed the most by Prometheus and only appears when this satellite is close to apoapsis. In addition, plots of the anti-alignment configuration for different orbital stages of Prometheus are obtained and discussed here.     

Parametric analysis of modeled ion escape from Mars

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors - that of the EUV flux R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , where ?(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O⁺, we obtain the values, 1.73 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ?(i) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.     

Impact cratering records of the mid-sized, icy saturnian satellites

Resolution of Voyager 1 and 2 images of the mid-sized, icy saturnian satellites was generally not much better than 1 km per line pair, except for a few, isolated higher resolution images. Therefore, analyses of impact crater distributions were generally limited to diameters (D) of tens of kilometers. Even with the limitation, however, these analyses demonstrated that studying impact crater distributions could expand understanding of the geology of the saturnian satellites and impact cratering in the outer Solar System. Thus to gain further insight into Saturn’s mid-sized satellites and impact cratering in the outer Solar System, we have compiled cratering records of these satellites using higher resolution CassiniISS images. Images from Cassini of the satellites range in resolution from tens m/pixel to hundreds m/pixel. These high-resolution images provide a look at the impact cratering records of these satellites never seen before, expanding the observable craters down to diameters of hundreds of meters. The diameters and locations of all observable craters are recorded for regions of Mimas, Tethys, Dione, Rhea, Iapetus, and Phoebe. These impact crater data are then analyzed and compared using cumulative, differential and relative (R) size-frequency distributions. Results indicate that the heavily cratered terrains on Rhea and Iapetus have similar distributions implying one common impactor population bombarded these two satellites. The distributions for Mimas and Dione, however, are different from Rhea and Iapetus, but are similar to one another, possibly implying another impactor population common to those two satellites. The difference between these two populations is a relative increase of craters with diameters between 10 and 30 km and a relative deficiency of craters with diameters between 30 and 80 km for Mimas and Dione compared with Rhea and Iapetus. This may support the result from Voyager images of two distinct impactor populations. One population was suggested to have a greater number of large impactors, most likely heliocentric comets (Saturn Population I in the Voyager literature), and the other a relative deficiency of large impactors and a greater number of small impactors, most likely planetocentric debris (Saturn Population II). Meanwhile, Tethys’ impact crater size-frequency distribution, which has some similarity to the distributions of Mimas, Dione, Rhea, and Iapetus, may be transitional between the two populations. Furthermore, when the impact crater distributions from these older cratered terrains are compared to younger ones like Dione’s smooth plains, the distributions have some similarities and differences. Therefore, it is uncertain whether the size-frequency distribution of the impactor population(s) changed over time. Finally, we find that Phoebe has a unique impact crater distribution. Phoebe appears to be lacking craters in a narrow diameter range around 1 km. The explanation for this confined “dip” at D = 1 km is not yet clear, but may have something to do with the interaction of Saturn’s irregular satellites or the capture of Phoebe.     

A procedure to analyze nonlinear density waves in Saturn's rings using several occultation profiles

Cassini radio science experiments have provided multiple occultation optical depth profiles of Saturn's rings that can be used in combination to analyze density waves. This paper establishes an accurate procedure of inversion of the wave profiles to reconstruct the wave kinematic parameters as a function of semi-major axis, in the nonlinear regime. This procedure is established using simulated data in the presence of realistic noise perturbations, to control the reconstruction error. It is then applied to the Mimas 5:3 density wave. There are two important concepts at the basis of this procedure. The first one is that it uses the nonlinear representation of density waves, and the second one is that it relies on a combination of optical depth profiles instead of just one profile. A related method to analyze density waves was devised by Longaretti and Borderies [Longaretti, P.-Y., Borderies, N., 1986. Icarus 67, 211-223] to study the nonlinear density wave associated with the Mimas 5:3 resonance, but the single photopolarimetric profile provided limited constraints. Other studies of density waves analyzing Cassini data [Colwell, J.E., Esposito, L.W., 2007. Bull. Am. Astron. Soc. 39, 461; Tiscareno, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34] are based on the linear theory and find inconsistent results from profile to profile. Multiple cuts of the rings are helpful in a fundamental way to ensure the accuracy of the procedure by forcing consistency among the various optical depth profiles. By way of illustration we have applied our procedure to the Mimas 5:3 density wave. We were able to recover precisely the kinematic parameters from the radio experiment occultation data in most of the propagation region; a preliminary analysis of the pressure-corrected dispersion allowed us to determine new but still uncertain values for the opacity (K?0.02 cm²/g) and velocity dispersion of (c₀?0.6 cm/s) in the wave region. Our procedure constitutes the first step in our planned analysis of the density waves of Saturn's rings. It is very accurate and efficient in the far-wave region. However, improvements are required within the first wavelength. The ways in which this method can be used to establish diagnostics of ring physics are outlined.     

Low-speed impacts between rubble piles modeled as collections of polyhedra, 2

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.     

Jovian chromophore characteristics from multispectral HST images

The chromophores responsible for coloring the jovian atmosphere are embedded within Jupiter’s vertical aerosol structure. Sunlight propagates through this vertical distribution of aerosol particles, whose colors are defined by ?₀(λ), and we remotely observe the culmination of the radiative transfer as I/F(λ). In this study, we employed a radiative transfer code to retrieve ?₀(λ) for particles in Jupiter’s tropospheric haze at seven wavelengths in the near-UV and visible regimes. The data consisted of images of the 2008 passage of Oval BA to the south of the Great Red Spot obtained by the Wide Field Planetary Camera 2 on-board the Hubble Space Telescope. We present derived particle colors for locations that were selected from 14 weather regions, which spanned a large range of observed colors. All ?₀(λ) curves were absorbing in the blue, and ?₀(λ) increased monotonically to approximately unity as wavelength increased. We found accurate fits to all ?₀(λ) curves using an empirically derived functional form: ?₀(λ) = 1 − A exp(−Bλ). The best-fit parameters for the mean ?₀(λ) curve were A = 25.4 and B = 0.0149 for λ in units of nm. We performed a principal component analysis (PCA) on our ?₀(λ) results and found that one or two independent chromophores were sufficient to produce the variations in ?₀(λ). A PCA of I/F(λ) for the same jovian locations resulted in principal components (PCs) with roughly the same variances as the ?₀(λ) PCA, but they did not result in a one-to-one mapping of PC amplitudes between the ?₀(λ) PCA and I/F(λ) PCA. We suggest that statistical analyses performed on I/F(λ) image cubes have limited applicability to the characterization of chromophores in the jovian atmosphere due to the sensitivity of I/F(λ) to horizontal variations in the vertical aerosol distribution.