We describe the prediction, design, execution and calibration of stellar and solar occultation observations of Saturn's rings by the Visual and Infrared Mapping Spectrometer (VIMS) instrument on the Cassini spacecraft. Particular attention is paid to the technique developed for onboard acquisition of the stellar target and to the geometric and photometric calibration of the data. Examples of both stellar and solar occultation data are presented, highlighting several aspects of the data as well as the different occultation geometries encountered during Cassini's 13 year orbital tour. Complete catalogs of ring stellar and solar occultations observed by Cassini-VIMS are presented, as a guide to the standard data sets which have been delivered to the Planetary Data System's Ring Moon Systems Node.

Density waves in Saturn's rings are usually tightly wrapped spiral patterns generated by resonances with either Saturn's moons or structures inside the planet. However, between the Barnard and Bessel Gaps in the Cassini Division (i.e. between 120,240 and 120,300 km from Saturn's spin axis), there are density variations that appear to form an axisymmetric density wave consisting of concentric zones of varying densities that propagate radially through the rings. Axisymmetric waves cannot be generated directly by a satellite resonance, but instead appear to be excited by interference between a nearby satellite resonance and normal mode oscillations on the inner edge of the Barnard Gap. Similar axisymmetric waves may exist just interior to other resonantly confined edges that exhibit a large number of normal modes, including the Dawes ringlet in the outer C ring and the outermost part of the B ring.

Recent studies of stellar occultations observed by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft have demonstrated that multiple spiral wave structures in Saturn's rings are probably generated by normal-mode oscillations inside the planet. Wavelet-based analyses have been able to unambiguously determine both the number of spiral arms and the rotation rate of many of these patterns. However, there are many more planetary normal modes that should have resonances in the rings, implying that many normal modes do not have sufficiently large amplitudes to generate obvious ring waves. Fortunately, recent advances in wavelet analysis allow weaker wave signals to be uncovered by combining data from multiple occultations. These new analytical tools reveal that a pattern previously identified as a single spiral wave actually consists of two superimposed waves, one with 5 spiral arms rotating at 1593.6 degrees/day and one with 11 spiral arms rotating at 1450.5 degrees/day. Furthermore, a broad search for new waves revealed four previously unknown wave patterns with 6, 7, 8 and 9 spiral arms rotating around the planet at 1538.2 degrees/day, 1492.5 degrees/day, 1454.2 degrees/day and 1421.8 degrees/day, respectively. These six patterns provide precise frequencies for another six fundamental normal modes inside Saturn, yielding what is now a complete sequence of fundamental sectoral normal modes with azimuthal wavenumbers from 2 to 10. These frequencies should place strong constraints on Saturn's interior structure and rotation rate, while the relative amplitudes of these waves should help clarify how the corresponding normal modes are excited inside the planet.

We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multispectral summary images have been computed for each single VIMS cube in order to rapidly identify their scientific content and assess their quality. These summary images are made available to the community on a public website (vims.univ-nantes.fr). The global mapping work faced several challenges due to the strong absorbing and scattering effects of the atmosphere coupled to the changing observing conditions linked to the orbital tour of the Cassini mission. We determined a surface photometric function which accounts for variations in incidence, emergence and phase angles, and which is able to mitigate brightness variations linked to the viewing geometry of the flybys. The atmospheric contribution has been reduced using the subtraction of the methane absorption band wings, considered as proxies for atmospheric haze scattering. We present a new global three color composite map of band ratios (red: 1.59/1.27 {\mu}m; green: 2.03/1.27 {\mu}m; blue: 1.27/1.08 {\mu}m), which has also been empirically corrected from an airmass (the solar photon path length through the atmosphere) dependence. This map provides a detailed global color view of Titan's surface partially corrected from the atmosphere and gives a global insight of the spectral variability, with the equatorial dunes fields appearing in brownish tones, and several occurrences of bluish tones localized in areas such as Sinlap, Menvra and Selk craters. (...)

We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We leave the detailed comparison with the Imaging Science Subsystem (ISS) and the Composite Infrared Spectrometer (CIRS) data sets to further studies. We have computed global hyperspectral mosaics for each of the 127 targeted flybys of Titan to produce synthetic color maps emphasizing the main atmospheric features. The north pole appears fully covered by a huge cloud as soon as the first observations in 2004 and up to the equinox in 2009 (Le Mou\'elic et al. 2012). The northern skies then became progressively clearer, after the circulation turnover in 2009, revealing the underlying lakes and seas to the optical instruments up to 2017. The reverse situation is observed over the south pole, which was mostly clear of such a high obscuring cloud during the first years of the mission, but started to develop a polar cloud in 2012. This feature grew up month after month until the end of the mission in 2017, with a poleward latitudinal extent of 75$^\circ$S in 2013 up to 58$^\circ$S in April 2017. Thanks to the spectral capabilities of VIMS, we have detected HCN spectral signatures over the north pole in almost all flybys between 2004 and 2008. These HCN signatures started then to show up over the south pole in almost all flybys between 2012 and 2017, so perfectly matching the timing and spatial extent of the northern and southern polar atmospheric features.

On day 138 of 2010, the plume of dust and gas emerging from Enceladus' South Polar Terrain passed between the Sun and the Cassini spacecraft. This solar occultation enabled Cassini's Ultraviolet Imaging Spectrograph (UVIS) and the Visual and Infrared Mapping Spectrometer (VIMS) to obtain simultaneous measurements of the plume's gas and dust components along the same lines of sight. The UVIS measurements of the plume's gas content are described in Hansen et al. (2011, GRL 38:11202) , while this paper describes the VIMS data and the information they provide about the plume's particle content. Together, the VIMS and UVIS measurements reveal that the plume material above Baghdad and Damascus sulci has a dust-to-gas mass ratio that is roughly an order of magnitude higher than the material above Alexandria and Cairo sulci. Similar trends in the plume's dust-to-gas ratio are also found in data obtained when Cassini flew through the plume in 2009, during which time the Ion and Neutral Mass Spectrometer (INMS), Radio and Plasma Wave Science instrument (RPWS) and Cosmic Dust Analyzer (CDA) instruments made in-situ measurements of the plume's gas and dust densities (Dong et al. 2015 JGR 120:915-937). These and other previously-published systematic differences in the material erupting from different fissures likely reflect variations in subsurface conditions across Encealdus' South Polar Terrain.

Planet-planet scattering is the leading mechanism to explain the broad eccentricity distribution of observed giant exoplanets. Here we study the orbital stability of primordial giant planet moons in this scenario. We use N-body simulations including realistic oblateness and evolving spin evolution for the giant planets. We find that the vast majority (80-90% across all our simulations) of orbital parameter space for moons is destabilized. There is a strong radial dependence, as moons past 0.1 Hill radii are systematically removed. Closer-in moons on Galilean-moon-like orbits (< 0.04 Hill radii) have a good ( 20-40%) chance of survival. Destabilized moons may undergo a collision with the star or a planet, be ejected from the system, be captured by another planet, be ejected but still orbiting its free-floating host planet, or survive on heliocentric orbits as "planets." The survival rate of moons increases with the host planet mass but is independent of the planet's final (post-scattering) orbits. Based on our simulations we predict the existence of an abundant galactic population of free- floating (former) moons.

The viscous spreading of planetary rings is believed to be counteracted by satellite torques, either through an individual resonance or through overlapping resonances. For the A ring of Saturn, it has been commonly believed that the satellite Janus alone can prevent the ring from spreading via its 7:6 Lindblad resonance. We discuss this common misconception and show that, in reality, the A ring is confined by the contributions from the group of satellites Pan, Atlas, Prometheus, Pandora, Janus, Epimetheus, and Mimas, whose cumulative torques from various resonances gradually decrease the angular momentum flux transported outward through the ring via density and bending waves. We further argue that this decrease in angular momentum flux occurs through 'flux reversal'. Furthermore, we use the magnitude of the satellites' resonance torques to estimate the effective viscosity profile across the A ring, showing that it decreases with radius from ~50 cm2 s-1 to less than ~10 cm2 s-1. The gradual estimated decrease of the angular momentum flux and effective viscosity are roughly consistent with results obtained by balancing the shepherding torques from Pan and Daphnis with the viscous torque at the edges of the Encke and Keeler gaps, as well as the edge of the A ring. On the other hand, the Mimas 2:1 Lindblad resonance alone seems to be capable of confining the edge of the B ring, and contrary to the situation in the A ring, we show that the effective viscosity across the B ring is relatively constant at ~24-30 cm2 s-1.

We analyze several thousand Cassini ISS images in order to study the inner edge of the Keeler gap in Saturn's outer A ring. We find strong evidence for an m=32 perturbation with a mean amplitude of radial variation of 4.5 km. Phase analysis yields a pattern speed consistent with the mean motion of Prometheus, indicating that this pattern is generated by the 32:31 Inner Lindblad resonance with Prometheus. In addition, we find evidence of 18-lobed and 20-lobed patterns with amplitudes of ~1.5 km. These patterns, whose rotation rates correspond to resonance locations ~4 km interior to the gap edge, are believed to be normal modes. The former is probably related to the nearby 18:17 (m=18) resonance with Pandora. In addition to these resonant and normal mode patterns, we also observe multiple localized features that appear to move at the local keplerian rate and that persist for only a few months. One hypothesis is that different groups of ring particles at the inner edge of the gap may be reacting differently to the resonance with Prometheus, with local variations in the forced eccentricity and/or pericenter; an alternative hypothesis is the existence of several unseen objects embedded at or near the inner edge of the Keeler gap, similar to those suspected to exist at the outer edges of the A and B rings (Spitale and Porco 2009, 2010). In either case, observations of the ring edge at opposite ansae demonstrate that the localized features must be on eccentric orbits.

We report the orbital distribution of the Trans-Neptunian objects (TNOs) discovered during the High Ecliptic Latitude (HiLat) extension of the Canada-France Ecliptic Plane Survey (CFEPS), conducted from June 2006 to July 2009. The HiLat component was designed to address one of the shortcomings of ecliptic surveys (like CFEPS), their lack of sensitivity to high-inclination objects. We searched 701~deg$^2$ of sky ranging from 12$^\circ$ to 85$^\circ$ ecliptic latitude and discovered \lKBO TNOs, with inclinations between 15$^\circ$ to 104$^\circ$. This survey places a very strong constraint on the inclination distribution of the hot component of the classical Kuiper Belt, ruling out any possibility of a large intrinsic fraction of highly inclined orbits. Using the parameterization of \citet{2001AJ....121.2804B}, the HiLat sample combined with CFEPS imposes a width $14^\circ \le \sigma \le 15.5^\circ$, with a best match for $\sigma = 14.5^\circ$. HiLat discovered the first retrograde TNO, 2008~KV$_{42}$, with an almost polar orbit with inclination 104$^\circ$, and (418993), a scattering object with perihelion in the region of Saturn's influence, with $a \sim 400$~AU and $i = 68^\circ$.

Saturn's B ring is the most opaque ring in our solar system, but many of its fundamental parameters, including its total mass, are not well constrained. Spiral density waves generated by mean-motion resonances with Saturn's moons provide some of the best constraints on the rings' mass density, but detecting and quantifying such waves in the B ring has been challenging because of this ring's high opacity and abundant fine-scale structure. Using a wavelet-based analyses of 17 occultations of the star gamma Crucis observed by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft, we are able to examine five density waves in the B ring. Two of these waves are generated by the Janus 2:1 and Mimas 5:2 Inner Lindblad Resonances at 96,427 km and 101,311 km from Saturn's center, respectively. Both of these waves can be detected in individual occultation profiles, but the multi-profile wavelet analysis reveals unexpected variations in the pattern speed of the Janus 2:1 wave that might arise from the periodic changes in Janus' orbit. The other three wave signatures are associated with the Janus 3:2, Enceladus 3:1 and Pandora 3:2 Inner Lindblad Resonances at 115,959 km, 115,207 km and 108,546 km. These waves are not visible in individual profiles, but structures with the correct pattern speeds can be detected in appropriately phase-corrected average wavelets. Estimates of the ring's surface mass density derived from these five waves fall between 40 and 140 g/cm^2, even though the ring's optical depth in these regions ranges from 1.5 to almost 5. This suggests that the total mass of the B ring is most likely between one-third and two-thirds the mass of Saturn's moon Mimas.

The spectral position of the 3.6 micron continuum peak measured on Cassini-VIMS I/F spectra is used as a marker to infer the temperature of the regolith particles covering the surfaces of Saturn's icy satellites. This feature is characterizing the crystalline water ice spectrum which is the dominant compositional endmember of the satellites' surfaces. Laboratory measurements indicate that the position of the 3.6 micron peak of pure water ice is temperature-dependent, shifting towards shorter wavelengths when the sample is cooled, from about 3.65 micron at T=123 K to about 3.55 micron at T=88 K. A similar method was already applied to VIMS Saturn's rings mosaics to retrieve ring particles temperature (Filacchione et al., 2014). We report here about the daytime temperature variations observed on the icy satellites as derived from three different VIMS observation types. Temperature maps are built by mining the complete VIMS dataset collected in years 2004-2009 (pre-equinox) and in 2009-2012 (post equinox) by selecting pixels with max 150 km/pixel resolution. VIMS-derived temperature maps allow to identify thermal anomalies across the equatorial lens of Mimas and Tethys.

We present the Cassini Atlas Of Stellar Spectra (CAOSS), comprised of near-infrared low-resolution spectra of bright stars recovered from space-based observations by the Cassini spacecraft. The 65 stellar targets in the atlas are predominately M, K and S giants. However it also contains spectra of other bright nearby stars including carbon stars and main sequence stars from A to F. The spectra presented are free of all spectral contamination caused by the Earth's atmosphere, including the detrimental telluric molecular bands which put parts of the near-infrared spectrum out of reach of terrestrial observations. With a single instrument, a spectro-photometric dataset is recovered that spans the near-infrared from 0.8 to 5.1 microns with spectral resolution ranging from R=53.5 to R=325. Spectra have been calibrated into absolute flux units after careful characterisation of the instrumental spectral efficiency. Spectral energy distributions for most stars match closely with literature values. All final data products have been made available online.

The James Webb Space Telescope (JWST) will provide unprecedented opportunities to observe the rings and small satellites in our solar system, accomplishing three primary objectives: 1) discovering new rings and moons, 2) unprecedented spectroscopy, and 3) time-domain observations. We give details on these science objectives and describe requirements that JWST must fulfill in order to accomplish the science objectives

High angular resolution images of IRC+10216 are presented in several near infrared wavelengths spanning more than 8 years. These maps have been reconstructed from interferometric observations obtained at both Keck and the VLT, and also from stellar occultations by the rings of Saturn observed with the Cassini spacecraft. The dynamic inner regions of the circumstellar environment are monitored over eight epochs ranging between January 2000 and July 2008. The system is shown to experience substantial evolution within this period including the fading of many previously reported persistent features, some of which had been identified as the stellar photosphere. These changes are discussed in context of existing models for the nature of the underlying star and the circumstellar environment. With access to these new images, we are able to report that none of the previously identified bright spots in fact contain the star, which is buried in its own dust and not directly visible in the near infrared.

We use solar occultations observed by the Visual and Infrared Mapping Spectrometer aboard the Cassini Spacecraft to extract the 1 to 5 micron transmission spectrum of Saturn, as if it were a transiting exoplanet. We detect absorption from methane, ethane, acetylene, aliphatic hydrocarbons, and possibly carbon monoxide with peak-to-peak features of up to 90 parts-per-million despite the presence of ammonia clouds. We also find that atmospheric refraction, as opposed to clouds or haze, determines the minimum altitude that could be probed during mid-transit. Self-consistent exoplanet atmosphere models show good agreement with Saturn's transmission spectrum but fail to reproduce a large absorption feature near 3.4 microns likely caused by gaseous ethane and a C-H stretching mode of an unknown aliphatic hydrocarbon. This large feature is located in one of the Spitzer Space Telescope bandpasses and could alter interpretations of transmission spectra if not properly modeled. The large signal in Saturn's transmission spectrum suggests that transmission spectroscopy of cold, long-period gaseous exoplanets should be possible with current and future observatories. Motivated by these results, we briefly consider the feasibility of a survey to search for and characterize cold exoplanets analogous to Jupiter and Saturn using a target-of-opportunity approach.

We present a study of the behavior of Saturn's A ring outer edge, using images and occultation data obtained by the Cassini spacecraft over a period of 8 years from 2006 to 2014. More than 5000 images and 170 occultations of the A ring outer edge are analyzed. Our fits confirm the expected response to the Janus 7:6 Inner Lindblad resonance (ILR) between 2006 and 2010, when Janus was on the inner leg of its regular orbit swap with Epimetheus. During this period, the edge exhibits a regular 7-lobed pattern with an amplitude of 12.8 km and one minimum aligned with the orbital longitude of Janus, as has been found by previous investigators. However, between 2010 and 2014, the Janus/Epimetheus orbit swap moves the Janus 7:6 LR away from the A ring outer edge, and the 7-lobed pattern disappears. In addition to several smaller-amplitudes modes, indeed, we found a variety of pattern speeds with different azimuthal wave numbers, and many of them may arise from resonant cavities between the ILR and the ring edge; also we found some other signatures consistent with tesseral resonances that could be associated with inhomogeneities in Saturn's gravity field. Moreover, these signatures do not have a fixed pattern speed. We present an analysis of these data and suggest a possible dynamical model for the behavior of the A ring's outer edge after 2010.

The spatial distribution of the tropospheric methane on Titan was measured using near-infrared spectroscopy. Ground-based observations at 1.5$\mu{\rm m}$ (H-band) were performed during the same night using instruments with adaptive optics at both the W. M. Keck Observatory and at the Paranal Observatory on 17 July 2014 UT. The integral field observations with SINFONI on the VLT covered the entire H-band at moderate resolving power, $R=\lambda/\Delta\lambda\approx1,500$, while the Keck observations were performed with NIRSPAO near 1.55254$\mu{\rm m}$ at higher resolution, $R\approx25,000$. The moderate resolution observations are used for flux calibration and for the determination of model parameters that can be degenerate in the interpretation of high resolution spectra. Line-by-line calculations of CH$_4$ and CH$_3$D correlated $k$ distributions from the HITRAN 2012 database were used, which incorporate revised line assignments near 1.5$\mu{\rm m}$. We fit the surface albedo and aerosol distributions in the VLT SINFONI observations that cover the entire H-band window and used these quantities to constrain the models of the high-resolution Keck NIRSPAO spectra when retrieving the methane abundances. Cassini VIMS images of the polar regions, acquired on 20 July 2014 UT, are used to validate the assumption that the opacity of tropospheric aerosol is relatively uniform below 10 km. We retrieved methane abundances at latitudes between 42$^{\circ}$ S and 80$^{\circ}$ N. The tropospheric methane in the Southern mid-latitudes was enhanced by a factor of $\sim$10--40% over the nominal profile that was measured using the GCMS on Huygens. The Northern hemisphere had $\sim$90% of the nominal methane abundance up to polar latitudes (80$^{\circ}$N). These measurements suggest that a source of saturated polar air is equilibrating with dryer conditions at lower latitudes.

We present an advance in the use of Cassini observations of stellar occultations by the rings of Saturn for stellar studies. Stewart et al. (2013) demonstrated the potential use of such observations for measuring stellar angular diameters. Here, we use these same observations, and tomographic imaging reconstruction techniques, to produce two dimensional images of complex stellar systems. We detail the determination of the basic observational reference frame. A technique for recovering model-independent brightness profiles for data from each occulting edge is discussed, along with the tomographic combination of these profiles to build an image of the source star. Finally we demonstrate the technique with recovered images of the {\alpha} Centauri binary system and the circumstellar environment of the evolved late-type giant star, Mira.

This work shows the dynamical instability that can happen to close-in satellites when planet oblateness is not accounted for in non-coplanar multiplanet systems. Simulations include two secularly interacting Jupiter-mass planets mutually inclined by 10 degrees, with the host planet either oblate or spherical. With a spherical host planet, moons within a critical planetocentric distance experience high inclinations and in some cases high eccentricities, while more distant moons orbit stably with low inclinations and eccentricities, as expected. These counter-intuitive dynamical phenomena disappear with an oblate host planet, in which case the moons' Laplace plane transitions from the host planet's equatorial plane to the host planet's precessing orbital plane as their semi-major axes increase, and all moons are dynamically stable with very mild changes in orbits. Direct perturbation from the perturbing planet has been investigated and ruled out as an explanation for the behavior of the innermost satellites, therefore leaving the central star's perturbation as the cause. Instability occurs while the nodal precession of the satellite and the central star (as seen from the host planet's frame) approaches the 1:1 secular resonance. In non-coplanar systems, around a non-oblate planet, the nodal precession of the moon becomes slow and comparable to that of the planet, giving rise to resonant configurations. The above effect needs to be taken into account in setting up numerical simulations.

We present an analysis of the VIMS solar occultations dataset, which allows us to extract vertically resolved information on the characteristics of Titan's atmosphere between 100-700 km with a characteristic vertical resolution of 10 km. After a series of data treatment procedures, 4 occultations out of 10 are retained. This sample covers different seasons and latitudes of Titan. The transmittances show clearly the evolution of the haze and detect the detached layer at 310 km in Sept. 2011 at mid-northern latitudes. Through the inversion of the transmission spectra with a line-by-line radiative transfer code we retrieve the vertical distribution of CH$_4$ and CO mixing ratio. The two methane bands at 1.4 and 1.7 {\mu}m are always in good agreement and yield an average stratospheric abundance of $1.28\pm0.08$%. This is significantly less than the value of 1.48% obtained by the GCMS/Huygens instrument. The analysis of the residual spectra after the inversion shows that there are additional absorptions which affect a great part of the VIMS wavelength range. We attribute many of these additional bands to gaseous ethane, whose near-infrared spectrum is not well modeled yet. Ethane contributes significantly to the strong absorption between 3.2-3.5 {\mu}m that was previously attributed only to C-H stretching bands from aerosols. Ethane bands may affect the surface windows too, especially at 2.7 {\mu}m. Other residual bands are generated by stretching modes of C-H, C-C and C-N bonds. In addition to the C-H stretch from aliphatic hydrocarbons at 3.4 {\mu}m, we detect a strong and narrow absorption at 3.28 {\mu}m which we tentatively attribute to the presence of PAHs in the stratosphere. C-C and C-N stretching bands are possibly present between 4.3-4.5 {\mu}m. Finally, we obtain the CO mixing ratio between 70-170 km. The average result of $46\pm16$ ppm is in good agreement with previous studies.

In a previous paper (Hedman and Nicholson 2013), we developed tools that allowed us to confirm that several of the waves in Saturn's rings were likely generated by resonances with fundamental sectoral normal modes inside Saturn itself. Here we use these same tools to examine eight additional waves that are probably generated by structures inside the planet. One of these waves appears to be generated by a resonance with a fundamental sectoral normal mode with azimuthal harmonic number m=10. If this attribution is correct, then the m=10 mode must have a larger amplitude than the modes with m=5-9, since the latter do not appear to generate strong waves. We also identify five waves with pattern speeds between 807 degrees/day and 834 degrees/day. Since these pattern speeds are close to the planet's rotation rate, they probably are due to persistent gravitational anomalies within the planet. These waves are all found in regions of enhanced optical depth known as plateaux, but surprisingly the surface mass densities they yield are comparable to the surface mass densities of the background C ring. Finally, one wave appears to be a one-armed spiral pattern whose rotation rate suggests it is generated by a resonance with a structure inside Saturn, but the nature of this perturbing structure remains unclear. Strangely, the resonant radius for this wave seems to be drifting inwards at an average rate of 0.8 km/year over the last thirty years, implying that the relevant planetary oscillation frequency has been steadily increasing.

Water-ice-poor, 5-$\mu$m-bright material on Saturn's moon Titan has previously been geomorphologically identified as evaporitic. Here we present a global distribution of the occurrences of the 5-$\mu$m-bright spectral unit, identified with Cassini's Visual Infrared Mapping Spectrometer (VIMS) and examined with RADAR when possible. We explore the possibility that each of these occurrences are evaporite deposits. The 5-$\mu$m-bright material covers 1\% of Titan's surface and is not limited to the poles (the only regions with extensive, long-lived surface liquid). We find the greatest areal concentration to be in the equatorial basins Tui Regio and Hotei Regio. Our interpretations, based on the correlation between 5-$\mu$m-bright material and lakebeds, imply that there was enough liquid present at some time to create the observed 5-$\mu$m-bright material. We address the climate implications surrounding a lack of evaporitic material at the south polar basins: if the south pole basins were filled at some point in the past, then where is the evaporite?

Image photometry reveals that the F ring is approximately twice as bright during the Cassini tour as it was during the Voyager flybys of 1980 and 1981. It is also three times as wide and has a higher integrated optical depth. We have performed photometric measurements of more than 4,800 images of Saturn's F ring taken over a five-year period with Cassini's Narrow Angle Camera. We show that the ring is not optically thin in many observing geometries and apply a photometric model based on single-scattering in the presence of shadowing and obscuration, deriving a mean effective optical depth tau = 0.033. Stellar occultation data from Voyager PPS and Cassini VIMS validate both the optical depth and the width measurements. In contrast to this decades-scale change, the baseline properties of the F ring have not changed significantly from 2004 to 2009. However, we have investigated one major, bright feature that appeared in the ring in late 2006. This transient feature increased the ring's overall mean brightness by 84% and decayed with a half-life of 91 days.

Certain regions of Saturn's rings exhibit periodic opacity variations with characteristic radial wavelengths of up to a few hundred meters that have been attributed to viscous overstabilities. The Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft observed two stellar occultations of the star gamma Crucis that had sufficient resolution to discern a sub-set of these periodic patterns in a portion of the A ring between 124,000 and 125,000 km from Saturn center. These data reveal that the wavelengths and intensities of the patterns vary systematically across this region, but that these parameters are not strictly determined by the ring's average optical depth. Furthermore, our observations indicate that these opacity variations have an azimuthal coherence scale of around 3000 km.

Radio occultations of Saturn's main rings by spacecraft suggest a power law particle size-distribution down to sizes of the order of 1 cm (Marouf et al., 1983), (Zebker et al., 1985). The lack of optical depth variations between ultraviolet and near-IR wavelengths indicate a lack of micron-sized particles. Between these two regimes, the particle-size distribution is largely unknown. A cutoff where the particle-size distribution turns over must exist, but the position and shape of it is not clear from existing studies. Using a series of solar occultations performed by the VIMS instrument on-board Cassini in the near-infrared, we are able to measure light forward scattered by particles in the A and C rings. With a model of diffraction by ring particles, and the previous radio work as a constraint on the slope of the particle size distribution, we estimate the minimum particle size using a truncated power-law size distribution. The C Ring shows a minimum particle size of $4.1^{+3.8}_{-1.3}$ mm, with an assumed power law index of q=3.1 and a maximum particle size of 10 m. The A Ring signal shows a similar level of scattered flux, but modeling is complicated by the presence of self-gravity wakes and higher optical depths. If q<3, our A Ring model requires a minimum particle size below one millimeter (< 0.34 mm for an assumed q=2.75, or $0.56^{+0.35}_{-0.16}$ mm for a steeper q=2.9) to be consistent with VIMS observations. These results might seem to contradict previous optical(Dones et al., 1993) and infrared (French and Nicholson, 2000) work, which implied that there were few particles in the A Ring smaller than 1 cm. But, because of the shallow power law, relatively little optical depth (between 0.03 and 0.16 in extinction, or 0.015 - 0.08 in absorption) is provided by these particles.

We present a photometric model of the rings of Saturn which includes the main rings and an F ring, inclined to the main rings, with a Gaussian vertical profile of optical depth. This model reproduces the asymmetry in brightness between the east and west ansae of the rings of Saturn that was observed by the Hubble Space Telescope (HST) within a few hours after the Earth ring-plane crossing (RPX) of 10 August 1995. The model shows that during this observation the inclined F ring unevenly blocked the east and west ansae of the main rings. The brightness asymmetry produced by the model is highly sensitive to the vertical thickness and radial optical depth of the F ring. The F ring model that best matches the observations has a vertical full width at half maximum of 13 +/- 7 km and an equivalent depth of 10 +/- 4 km. The model also reproduces the shape of the HST profiles of ring brightness vs. distance from Saturn, both before and after the time of ring-plane crossing. Smaller asymmetries observed before the RPX, when the Earth was on the dark side of the rings, cannot be explained by blocking of the main rings by the F ring or vice versa and are probably instead due to the intrinsic longitudinal variation exhibited by the F ring.

Saturn's C ring contains multiple spiral patterns that appear to be density waves driven by periodic gravitational perturbations. In other parts of Saturn's rings, such waves are generated by Lindblad resonances with Saturn's various moons, but most of the wave-like C-ring features are not situated near any strong resonance with any known moon. Using stellar occultation data obtained by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft, we investigate the origin of six unidentified C-ring waves located between 80,900 and 87,200 km from Saturn's center. By measuring differences in the waves' phases among the different occultations, we are able to determine both the number of arms in each spiral pattern and the speeds at which these patterns rotate around the planet. We find that all six of these waves have between 2 and 4 arms and pattern speeds between 1660 degrees/day and 1861 degrees/day. These speeds are too large to be attributed to any satellite resonance. Instead they are comparable to the predicted pattern speeds of waves generated by low-order normal-mode oscillations within the planet [Marley & Porco 1993, Icarus 106, 508]. The precise pattern speeds associated with these waves should therefore provide strong constraints on Saturn's internal structure. Furthermore, we identify multiple waves with the same number of arms and very similar pattern speeds, indicating that multiple m=3 and m=2 sectoral (l=m) modes may exist within the planet.

We present novel observations utilising the Cassini spacecraft to conduct an observing campaign for stellar astronomy from a vantage point in the outer solar system. By exploiting occultation events in which Mira passed behind the Saturnian ring plane as viewed by Cassini, parametric imaging data were recovered spanning the near-infrared. From this, spatial information at extremely high angular resolution was recovered enabling a study of the stellar atmospheric extension across a spectral bandpass spanning the 1 - 5 {\mu}m spectral region in the near-infrared. The resulting measurements of the angular diameter of Mira were found to be consistent with existing observations of its variation in size with wavelength. The present study illustrates the validity of the technique; more detailed exploration of the stellar physics obtained by this novel experiment will be the subject of forthcoming papers.

Over the last eight years, the Visual and Infrared Mapping Spectrometer (VIMS) aboard the Cassini orbiter has returned hyperspectral images in the 0.35-5.1 micron range of the icy satellites and rings of Saturn. These very different objects show significant variations in surface composition, roughness and regolith grain size as a result of their evolutionary histories, endogenic processes and interactions with exogenic particles. The distributions of surface water ice and chromophores, i.e. organic and non-icy materials, across the saturnian system, are traced using specific spectral indicators (spectral slopes and absorption band depths) obtained from rings mosaics and disk-integrated satellites observations by VIMS.

Saturn's main rings exhibit variations in both their opacity and spectral properties on a broad range of spatial scales, and the correlations between these parameters can provide insights into the processes that shape the composition and dynamics of the rings. The Visual and Infrared Mapping Spectrometer (VIMS) instrument onboard the Cassini Spacecraft has obtained spectra of the rings between 0.35 and 5.2 microns with sufficient spatial resolution to discern variations on scales below 200 km. These relatively high-resolution spectral data reveal that both the depths of the near-infrared water-ice absorption bands and the visible spectral slopes are often correlated with structural parameters such as the rings' optical depth. Using a simplified model for the ring-particles' regolith properties, we have begun to disentangle the trends due to changes in the gross composition of the ring particles from those that may be due to shifts in the texture of the ring particles' regolith. Consistent with previous studies, this analysis finds that the C ring and the Cassini Division possess enhanced concentrations of a contaminant that absorbs light over a broad range of wavelengths. On the other hand, a second contaminant that preferentially absorbs at short visible and near-ultraviolet wavelengths is found to be more evenly distributed throughout the rings. The optical activity of this short-wavelength absorber increases in the inner B ring inwards of 100,000 km from Saturn center, which may provide clues to the origin of this contaminant. The spectral variations identified as shifts in the regolith texture are in some places clearly correlated with the ring's optical depth, and in other locations they appear to be associated with the disturbances generated by strong mean-motion resonances with Saturn's various moons.

High-inclination circumplanetary orbits that are gravitationally perturbed by the central star can undergo Kozai oscillations---large-amplitude, coupled variations in the orbital eccentricity and inclination. We first study how this effect is modified by incorporating perturbations from the planetary oblateness. Tremaine et al. (2009) found that, for planets with obliquities > 68.875 degrees, orbits in the equilibrium local Laplace plane are unstable to eccentricity perturbations over a finite radial range, and execute large-amplitude chaotic oscillations in eccentricity and inclination. In the hope of making that treatment more easily understandable, we analyze the problem using orbital elements, confirming this threshold obliquity. Furthermore, we find that orbits inclined to the Laplace plane will be unstable over a broader radial range, and that such orbits can go unstable for obliquities less than 68.875 degrees. Finally, we analyze the added effects of radiation pressure, which are important for dust grains and provide a natural mechanism for particle semimajor axes to sweep via Poynting-Robertson drag through any unstable range. We find that generally the effect persists; however, the unstable radial range is shifted and small retrograde particles can avoid the instability altogether. We argue that this is occurs because radiation pressure modifies the equilibrium Laplace plane.

The transneptunian objects (TNOs) trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant-planet formation process. We perform detailed modelling of the resonant objects detected in the Canada-France Ecliptic Plane Survey (CFEPS) in order to provide population estimates and, for some resonances, constrain the complex internal orbital element distribution. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. This paper discusses the 3:2, 5:2, 2:1, 3:1, 5:1, 4:3, 5:3, 7:3, 5:4, and 7:4 mean-motion resonances, all of which had CFEPS detections, along with our upper limit on 1:1 Neptune Trojans (which is consistent with their small population estimated elsewhere). For the plutinos (TNOs in the 3:2 resonance) we refine the orbital element distribution given in Kavelaars et al. (2009) and show that steep H-magnitude distributions (N(H) proportional to 10aH, with a=0.8-0.9) are favoured in the range Hg=8-9, and confirm that this resonance does not share the inclination distribution of the classical Kuiper Belt. We give the first population estimate for the 5:2 resonance and find that, to within the uncertainties, the population is equal to that of the 3:2 (13,000 TNOs with Hg < 9.16), whereas the 2:1 population is smaller by a factor of 3-4 compared to the other two resonances. We also measure significant populations inhabiting the 4:3, 5:3, 7:3, 5:4, 7:4, 3:1, and 5:1 resonances, with Hg < 9.16 (D >100 km) populations in the thousands. We compare our intrinsic population and orbital-element distributions with several published models of resonant-TNO production; the most striking discrepancy is that resonances beyond the 2:1 are in reality more heavily populated than in published models.

In the last few years Cassini-VIMS, the Visible and Infared Mapping Spectrometer, returned to us a comprehensive view of the Saturn's icy satellites and rings. After having analyzed the satellites' spectral properties (Filacchione et al. (2007a)) and their distribution across the satellites' hemispheres (Filacchione et al. (2010)), we proceed in this paper to investigate the radial variability of icy satellites (principal and minor) and main rings average spectral properties. This analysis is done by using 2,264 disk-integrated observations of the satellites and a 12x700 pixels-wide rings radial mosaic acquired with a spatial resolution of about 125 km/pixel. The comparative analysis of these data allows us to retrieve the amount of both water ice and red contaminant materials distributed across Saturn's system and the typical surface regolith grain sizes. These measurements highlight very striking differences in the population here analyzed, which vary from the almost uncontaminated and water ice-rich surfaces of Enceladus and Calypso to the metal/organic-rich and red surfaces of Iapetus' leading hemisphere and Phoebe. Rings spectra appear more red than the icy satellites in the visible range but show more intense 1.5-2.0 micron band depths. The correlations among spectral slopes, band depths, visual albedo and phase permit us to cluster the saturnian population in different spectral classes which are detected not only among the principal satellites and rings but among co-orbital minor moons as well. Finally, we have applied Hapke's theory to retrieve the best spectral fits to Saturn's inner regular satellites using the same methodology applied previously for Rhea data discussed in Ciarniello et al. (2011).

We report the orbital distribution of the trans-neptunian objects (TNOs) discovered during the Canada-France Ecliptic Plane Survey, whose discovery phase ran from early 2003 until early 2007. The follow-up observations started just after the first discoveries and extended until late 2009. We obtained characterized observations of 321 sq.deg. of sky to depths in the range g ~ 23.5--24.4 AB mag. We provide a database of 169 TNOs with high-precision dynamical classification and known discovery efficiency. Using this database, we find that the classical belt is a complex region with sub-structures that go beyond the usual splitting of inner (interior to 3:2 mean-motion resonance [MMR]), outer (exterior to 2:1 MMR), and main (in between). The main classical belt (a=40--47 AU) needs to be modeled with at least three components: the `hot' component with a wide inclination distribution and two `cold' components (stirred and kernel) with much narrower inclination distributions. The hot component must have a significantly shallower absolute magnitude (Hg) distribution than the other two components. With 95% confidence, there are 8000+1800-1600 objects in the main belt with Hg <= 8.0, of which 50% are from the hot component, 40% from the stirred component and 10% from the kernel; the hot component's fraction drops rapidly with increasing Hg. Because of this, the apparent population fractions depend on the depth and ecliptic latitude of a trans-neptunian survey. The stirred and kernel components are limited to only a portion of the main belt, while we find that the hot component is consistent with a smooth extension throughout the inner, main and outer regions of the classical belt; the inner and outer belts are consistent with containing only hot-component objects. The Hg <= 8.0 TNO population estimates are 400 for the inner belt and 10,000 for the outer belt within a factor of two.

Stellar occultations by Saturn's rings observed with the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft reveal that dusty features such as the F ring and the ringlets in the Encke and the Laplace Gaps have distinctive infrared transmission spectra. These spectra show a narrow optical depth minimum at wavelengths around 2.87 microns. This minimum is likely due to the Christiansen Effect, a reduction in the extinction of small particles when their (complex) refractive index is close to that of the surrounding medium. Simple Mie-scattering models demonstrate that the strength of this opacity dip is sensitive to the size distribution of particles between 1 and 100 microns across. Furthermore, the spatial resolution of the occultation data is sufficient to reveal variations in the transmission spectra within and among these rings. For example, in both the Encke Gap ringlets and F ring, the opacity dip weakens with increasing local optical depth, which is consistent with the larger particles being concentrated near the cores of these rings. The strength of the opacity dip varies most dramatically within the F ring; certain compact regions of enhanced optical depth lack an opacity dip and therefore appear to have a greatly reduced fraction of grains in the few-micron size range.Such spectrally-identifiable structures probably represent a subset of the compact optically-thick clumps observed by other Cassini instruments. These variations in the ring's particle size distribution can provide new insights into the processes of grain aggregation, disruption and transport within dusty rings. For example, the unusual spectral properties of the F-ring clumps could perhaps be ascribed to small grains adhering onto the surface of larger particles in regions of anomalously low velocity dispersion.

The Cassini Division in Saturn's rings contains a series of eight named gaps, three of which contain dense ringlets. Observations of stellar occultations by the Visual and Infrared Mapping Spectrometer onboard the Cassini spacecraft have yielded ~40 accurate and precise measurements of the radial position of the edges of all of these gaps and ringlets. These data reveal suggestive patterns in the shapes of many of the gap edges: the outer edges of the 5 gaps without ringlets are circular to within 1 km, while the inner edges of 6 of the gaps are eccentric, with apsidal precession rates consistent with those expected for eccentric orbits near each edge. Intriguingly, the pattern speeds of these eccentric inner gap edges, together with that of the eccentric Huygens ringlet,form a series with a characteristic spacing of 0.06 degrees/day. The two gaps with non-eccentric inner edges lie near first-order Inner Lindblad Resonances (ILRs) with moons. One such edge is close to the 5:4 ILR with Prometheus. The other resonantly confined edge is the outer edge of the B ring, which lies near the 2:1 Mimas ILR. Detailed investigation of the B-ring-edge data confirm the presence of an m=2 perturbation on the B-ring edge, but also suggest that this pattern moves or librates relative to Mimas. The B-ring edge also has an m=1 component that rotates around the planet at a rate close to the expected apsidal precession rate. The pattern speeds of the eccentric edges in the Cassini Division can potentially be generated from various combinations of the pattern speeds of structures observed on the edge of the B ring. We therefore suggest that the locations of most of the gaps in the Cassini Division may be determined by resonances involving a combination of perturbations from Mimas and the massive edge of the B ring.

Clouds on Titan result from the condensation of methane and ethane and, as on other planets, are primarily structured by circulation of the atmosphere. At present, cloud activity mainly occurs in the southern (summer) hemisphere, arising near the pole and at mid-latitudes from cumulus updrafts triggered by surface heating and/or local methane sources, and at the north (winter) pole, resulting from the subsidence and condensation of ethane-rich air into the colder troposphere. General circulation models predict that this distribution should change with the seasons on a 15-year timescale, and that clouds should develop under certain circumstances at temperate latitudes (~40\degree) in the winter hemisphere. The models, however, have hitherto been poorly constrained and their long-term predictions have not yet been observationally verified. Here we report that the global spatial cloud coverage on Titan is in general agreement with the models, confirming that cloud activity is mainly controlled by the global circulation. The non-detection of clouds at latitude ~40\degree N and the persistence of the southern clouds while the southern summer is ending are, however, both contrary to predictions. This suggests that Titan's equator-to-pole thermal contrast is overestimated in the models and that its atmosphere responds to the seasonal forcing with a greater inertia than expected.

Titan is one of the primary scientific objectives of the NASA ESA ASI Cassini Huygens mission. Scattering by haze particles in Titan's atmosphere and numerous methane absorptions dramatically veil Titan's surface in the visible range, though it can be studied more easily in some narrow infrared windows. The Visual and Infrared Mapping Spectrometer (VIMS) instrument onboard the Cassini spacecraft successfully imaged its surface in the atmospheric windows, taking hyperspectral images in the range 0.4 5.2 ?m. On 26 October (TA flyby) and 13 December 2004 (TB flyby), the Cassini Huygens mission flew over Titan at an altitude lower than 1200 km at closest approach. We report here on the analysis of VIMS images of the Huygens landing site acquired at TA and TB, with a spatial resolution ranging from 16 to14.4 km/pixel. The pure atmospheric backscattering component is corrected by using both an empirical method and a first-order theoretical model. Both approaches provide consistent results. After the removal of scattering, ratio images reveal subtle surface heterogeneities. A particularly contrasted structure appears in ratio images involving the 1.59 and 2.03 ?m images north of the Huygens landing site. Although pure water ice cannot be the only component exposed at Titan's surface, this area is consistent with a local enrichment in exposed water ice and seems to be consistent with DISR/Huygens images and spectra interpretations. The images show also a morphological structure that can be interpreted as a 150 km diameter impact crater with a central peak.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r~118,800 km, 4.5 m at r~120,700 km) and the inner A Ring (10 to 15 m for r<127,000 km).

We describe a model that accounts for the complex morphology of spiral density waves raised in Saturn's rings by the co-orbital satellites, Janus and Epimetheus. Our model may be corroborated by future Cassini observations of these time-variable wave patterns.

We report the discovery of the first trans-neptunian object, designated 2004 XR190, with a nearly-cirular orbit beyond the 2:1 mean-motion resonance. Fitting an orbit to 23 astrometric observations spread out over 12 months yields an orbit of a=57.2\pm0.4, e=0.08\pm0.04, and i=46.6 deg. All viable orbits have perihelia distances q>49 AU. The very high orbital inclination of this extended scattered disk object might be explained by several models, but its existence again points to a large as-yet undiscovered population of transneptunian objects with large orbital perihelia and inclination.

We present the data acquisition strategy and characterization procedures for the Canada-France Ecliptic Plane Survey (CFEPS), a sub-component of the Canada-France-Hawaii Telescope Legacy Survey. The survey began in early 2003 and as of summer 2005 has covered 430 square degrees of sky within a few degrees of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a magnitude limit of $m_R$=23 -- 24 (depending on the image quality). To track as large a sample as possible and avoid introducing followup bias, we have developed a multi-epoch observing strategy that is spread over several years. We present the evolution of the uncertainties in ephemeris position and orbital elements as the objects progress through the epochs. We then present a small 10-object sample that was tracked in this manner as part of a preliminary survey starting a year before the main CFEPS project. We describe the CFEPS survey simulator, to be released in 2006, which allows theoretical models of the Kuiper Belt to be compared with the survey discoveries since CFEPS has a well-documented pointing history with characterized detection efficiencies as a function of magnitude and rate of motion on the sky. Using the pre-survey objects we illustrate the usage of the simulator in modeling the classical Kuiper Belt.

By telescopic tracking, we have established that the orbit of the trans-neptunian object (2000 CR$_{105}$) has a perihelion of $\simeq$44 AU, and is thus outside the domain controlled by strong gravitational close encounters with Neptune. Because this object is on a very large, eccentric orbit (with semimajor axis $a\simeq$216 AU and eccentricity $e\simeq$0.8) this object must have been placed on this orbit by a gravitational perturbation which is {\it not} direct gravitational scattering off of any of the giant planets (on their current orbits). The existence of this object may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced it, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of trans-neptunian objects in an `extended scattered disk' with perihelia above the previously-discussed 38 AU boundary.

We have conducted pencil-beam searches for outer solar system objects to a limiting magnitude of R ~ 26. Five new trans-neptunian objects were detected in these searches. Our combined data set provides an estimate of ~90 trans-neptunian objects per square degree brighter than ~ 25.9. This estimate is a factor of 3 above the expected number of objects based on an extrapolation of previous surveys with brighter limits, and appears consistent with the hypothesis of a single power-law luminosity function for the entire trans-neptunian region. Maximum likelihood fits to all self-consistent published surveys with published efficiency functions predicts a cumulative sky density Sigma(<R) obeying log10(Sigma) = 0.76(R-23.4) objects per square degree brighter than a given magnitude R.