The brightening of Saturn’s F ring

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 4800 images of Saturn’s F ring taken over a 5-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 τ ≈ 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.     

The three-dimensional structure of ‘loop-like’ coronal mass ejections

We have measured the brightness and latitudinal extent of the depleted (cavity) region behind the leading edge of three coronal mass ejection (CME) events from the Skylab epoch. The events chosen are among those believed to be typical of non-impulsive, loop-like CMEs (Sime, MacQueen, and Hundhausen, 1984). The pre-event coronal brightness has been matched by a model corona, assuming both a background contribution and a contribution from a range of hypothetical streamer models, distinguished by differing longitudinal extent. Then, assuming that the cavities are voided regions in which the local electron density is negligible, we estimate their minimum line-of-sight extent and find them to be comparable to, or greater than, their measured latitudinal extent. As a result, we suggest this unambiguously demonstrates the three-dimensional nature of these — and likely, this class of — events.     

The origin and evolution of Saturn’s 2011–2012 stratospheric vortex

The planet-encircling springtime storm in Saturn’s troposphere (December 2010–July 2011, Fletcher, L.N. et al. [2011]. Science 332, 1413–1414; Sánchez-Lavega, A. et al. [2011]. Nature 475, 71–74; Fischer, G. et al. [2011]. Nature 475, 75–77) produced dramatic perturbations to stratospheric temperatures, winds and composition at mbar pressures that persisted long after the tropospheric disturbance had abated. Thermal infrared (IR) spectroscopy from the Cassini Composite Infrared Spectrometer (CIRS), supported by ground-based IR imaging from the VISIR instrument on the Very Large Telescope and the MIRSI instrument on NASA’s IRTF, is used to track the evolution of a large, hot stratospheric anticyclone between January 2011 and March 2012. The evolutionary sequence can be divided into three phases: (I) the formation and intensification of two distinct warm airmasses near 0.5 mbar between 25 and 35°N (B1 and B2) between January–April 2011, moving westward with different zonal velocities, B1 residing directly above the convective tropospheric storm head; (II) the merging of the warm airmasses to form the large single ‘stratospheric beacon’ near 40°N (B0) between April and June 2011, disassociated from the storm head and at a higher pressure (2 mbar) than the original beacons, a downward shift of 1.4 scale heights (approximately 85 km) post-merger; and (III) the mature phase characterised by slow cooling (0.11 ± 0.01 K/day) and longitudinal shrinkage of the anticyclone since July 2011. Peak temperatures of 221.6 ± 1.4 K at 2 mbar were measured on May 5th 2011 immediately after the merger, some 80 K warmer than the quiescent surroundings. From July 2011 to the time of writing, B0 remained as a long-lived stable stratospheric phenomenon at 2 mbar, moving west with a near-constant velocity of 2.70 ± 0.04 deg/day (?24.5 ± 0.4 m/s at 40°N relative to System III longitudes). No perturbations to visible clouds and hazes were detected during this period.With no direct tracers of motion in the stratosphere, we use thermal windshear calculations to estimate clockwise peripheral velocities of 200–400 m/s at 2 mbar around B0. The peripheral velocities of the two original airmasses were smaller (70–140 m/s). In August 2011, the size of the vortex as defined by the peripheral collar was 65° longitude (50,000 km in diameter) and 25° latitude. Stratospheric acetylene (C₂H₂) was uniformly enhanced by a factor of three within the vortex, whereas ethane (C₂H₆) remained unaffected. The passage of B0 generated a new band of warm stratospheric emission at 0.5 mbar at its northern edge, and there are hints of warm stratospheric structures associated with the beacons at higher altitudes (p < 0.1 mbar) than can be reliably observed by CIRS nadir spectroscopy. Analysis of the zonal windshear suggests that Rossby wave perturbations from the convective storm could have propagated vertically into the stratosphere at this point in Saturn’s seasonal cycle, one possible source of energy for the formation of these stratospheric anticyclones.     

Saturn’s F ring as seen by Cassini UVIS: Kinematics and statistics

We present a new orbital model of Saturn’s F ring core based on 93 occultations by the Cassini Ultraviolet Imaging Spectrograph (UVIS) and the Voyager radio and stellar occultations. We demonstrate that the core, despite its intrinsic variability, is well-described as an inclined, freely precessing ellipse. We find that post-fit residuals with a root-mean-square of 24 km are genuine, representing the well-known non-Keplerian features observed in the ring. Over the nearly 4 years of UVIS observations we find the residual variance to increase, coincident with the apse anti-alignment of Prometheus and F ring core in December 2009. This increase in dynamical F ring core temperature most likely reflects the ever-stronger perturbations by Prometheus. Our results are in good agreement with Earth-based and HST observations as well as Voyager imaging.Cassini UVIS stellar occultations resolve the F ring at unprecedented resolutions of a few meters and we identify the F ring core and inner and outer strands. We infer their normal optical depth and full width at half maximum (FWHM) and show that core and strands form distinct morphological groups. Typically, a strand is about ten times wider than the core (average FWHM is ～10 km) while having a ten times smaller optical depth. Unlike in pre-Cassini occultations the F ring core displays significant optical depth with in some cases >3. In many cases we find a narrow optically thick component (～ few km and τ > 0.5) embedded in the F ring core. Entertaining the possibility that this is the actual, “true” F ring core then UVIS results suggest that this “true” core is highly non-continuous. In addition, we report the detection of a previously unknown structure – dubbed the “secondary” as it visually resembles the F ring core. Its morphology is similar to that of the core in optical depth and FWHM and it displays individual opaque features. Despite its core-like appearance, we show that its kinematics is consistent with that of strands. We conclude that it is the most prominent strand seen to date. It represents a striking example of strand creation resulting in what could be called a morphological “small-scale” version of the F ring core. This extraordinary object should be one of the prime targets of future F ring studies.     

Solar Cycle Effects on the Dynamics of Jupiter’s and Saturn’s Magnetospheres

The giant planetary magnetospheres surrounding Jupiter and Saturn respond in quite different ways, compared to Earth, to changes in upstream solar wind conditions. Spacecraft have visited Jupiter and Saturn during both solar cycle minima and maxima. In this paper we explore the large-scale structure of the interplanetary magnetic field (IMF) upstream of Saturn and Jupiter as a function of solar cycle, deduced from solar wind observations by spacecraft and from models. We show the distributions of solar wind dynamic pressure and IMF azimuthal and meridional angles over the changing solar cycle conditions, detailing how they compare to Parker predictions and to our general understanding of expected heliospheric structure at 5 and 9 AU. We explore how Jupiter’s and Saturn’s magnetospheric dynamics respond to varying solar wind driving over a solar cycle under varying Mach number regimes, and consider how changing dayside coupling can have a direct effect on the nightside magnetospheric response. We also address how solar UV flux variability over a solar cycle influences the plasma and neutral tori in the inner magnetospheres of Jupiter and Saturn, and estimate the solar cycle effects on internally driven magnetospheric dynamics. We conclude by commenting on the effects of the solar cycle in the release of heavy ion plasma into the heliosphere, ultimately derived from the moons of Jupiter and Saturn.     

Small-scale variability in Saturn’s lower ionosphere

We perform and present a wavelet analysis on all 31 Cassini electron density profiles published to date (Nagy, A.F. et al. [2006]. J. Geophys. Res. 111 (A6), CiteID A06310; Kliore, A.J. et al. [2009]. J. Geophys. Res. 114 (A4), CiteID A04315). We detect several discrete scales of variability present in the observations. Small-scale variability (S < 700 km) is observed in almost all data sets at different latitudes, both at dawn and dusk conditions. The most typical scale of variability is 300 km with scales between 200 km and 450 km being commonly present in the vast majority of the profiles. A low latitude dawn/dusk asymmetry is noted in the prevalent scales with the spectrum peaking sharply at the 300 km scale at dusk conditions and being broader at dawn conditions. Compared to dawn conditions the dusk ionosphere also shows more significant variability at the 100 km scale. The 300 km vertical scale is also present in the few available profiles from the northern hemisphere. Early observations from 2005 show a dominant scale at 350 km whereas later in 2007–2008 the spectrum shifts to the shorter scales with the most prominent scale being 300 km. The performed wavelet analysis and the obtained results are independent of assumptions about the nature of the layers and do not require a definition for a “background” electron density profile.In the second part of the paper we present a gravity wave propagation/dissipation model for Saturn’s upper atmosphere and compare the wave properties to the characteristics of the observed electron density variability at different scales. The general features observed in the data are consistent with gravity waves being present in the lower ionosphere and causing layering of the ions and the electrons. The wave-driving mechanism provides a simultaneous explanation for several of the properties of the observed variability: (i) lack of variability in the electron density above the predicted region of wave dissipation; (ii) in most cases the peak amplitude of variability occurs within the altitude range for dissipation of gravity waves or below; (iii) shorter scales have smaller amplitudes than the longer scales; (iv) shorter scales are present at lower altitudes whereas longer scales persist to higher altitudes; and (v) several layers often form a system of equally spaced maxima and minima that can be traced over a large altitude range.     

On the thermal history of Saturn’s satellites Titan and Enceladus

The thermal histories of two geologically active satellites of Saturn—Titan and Enceladus—are discussed. During the Cassini mission, it was found that there are both nitrogen-containing compounds—NH₃ and N₂-and CO₂ and CH₄ in the water plumes of Enceladus; at that, ammonia is the prevailing form. This may testify that during evolution, the material of the satellite was warmed up to T ∼ 500–600 K, when NH₃ (the form of nitrogen capable of being accreted) could only be partly converted into N₂. Contrary to Enceladus, the temperature inside Titan probably reached values higher than 800 K or even higher than 1000 K, since the process of the chemical dissociation of ammonia was completely finished on this satellite and its atmosphere contains only molecular nitrogen. While the internal heating of Titan up to high temperatures can be explained by its large mass, the heating source for Enceladus’ interior is far from evident. Such traditional heating sources as the energy of gravitational differentiation and the radiogenic heating due to shortliving ²⁶Al and ⁶⁰Fe could not be effective. The first one is because of the small size of Enceladus (RE ≈ 250 km), and the inefficiency of the second one is caused by the fact that the satellite was formed not earlier than 8–10 Myr after the formation of calcium and aluminum-enriched inclusions in carbonaceous chondrites (CAI), i.e., after ²⁶Al had completely decayed. In the present paper, we propose other heating mechanisms-the heat of long-living radioactive elements and tidal heat, which could provide the observed chemical composition of the water plumes of Enceladus rather than only the differentiation of its protomatter into the ironstone core and the ice mantle.     

Processes forming and sustaining Saturn’s proton radiation belts

Saturn’s proton radiation belts extend over the orbits of several moons that split this region of intense radiation into several distinct belts. Understanding their distribution requires to understand how their particles are created and evolve. High-energy protons are thought to be dominantly produced by cosmic ray albedo neutron decay (CRAND). The source of the lower energies and the role of other effects such as charge exchange with the gas originating from Enceladus is still an open question. There is also no certainty so far if the belts exist independently from each other and the rest of the magnetosphere or if and how particles are exchanged between these regions. We approach these problems by using measurements acquired by the MIMI/LEMMS instrument onboard the Cassini spacecraft. Protons in the range from 500 keV to 40 MeV are considered. Their intensities are averaged over 7 years of the mission and converted to phase space densities at constant first and second adiabatic invariant. We reproduce the resulting radial profiles with a numerical model that includes radial diffusion, losses from moons and interactions with gas, and a phenomenological source. Our results show that the dominating effects away from the moon sweeping corridors are diffusion and the source, while interactions with gas are secondary. Based on a GEANT4 simulation of the interaction of cosmic rays with Saturn’s rings, we conclude that secondary particles produced within the rings can only account for the high-energy part of the source. A comparison with the equivalent processes within Earth’s atmosphere shows that Saturn’s atmosphere can contribute to the production of the lower energies and might be even dominating at the higher energies. Other possibilities to supply the belts and exchange particles between them, as diffusion and injections from outside the belts, or stripping of ENAs, can be excluded.     

Diurnal variation of electron density in Saturn’s ionosphere: Model comparisons with Saturn Electrostatic Discharge (SED) observations

Using the Saturn Thermosphere Ionosphere Model (STIM), we present a study of the diurnal variation of electron density, with a focus on comparisons with peak electron densities (N_MAX) inferred from the low-frequency cutoff of radio emission due to lightning in the lower atmosphere, called Saturn Electrostatic Discharges (SEDs). It is demonstrated that photochemistry in Saturn’s ionosphere cannot reproduce the SED-inferred diurnal variation in N_MAX unless additional production and loss sources outside of the current best estimates are considered. Additional explanations of the SED-inferred diurnal variation of N_MAX are presented and analyzed, such as the possibility that the low-frequency cutoff seen in SEDs is due to the presence of sharp low-altitude layers of plasma, as frequently seen in radio occultation measurements. Finally, we outline the observational constraints that must be fulfilled by any candidate explanations of the SED-inferred diurnal variation of N_MAX.     

Photographic observations of Saturn’s moons at the MAO NAS of Ukraine in 1961–1990

A catalog of 1385 astrometric positions of Saturn’s moons S2–S9 has been compiled with Tycho-2 as a reference frame from photographic observations obtained at the Main Astronomical Observatory, National Academy of Sciences of Ukraine, in 1961–1990. Astronegatives have been digitized with an Epson Expression 10000XL commercial scanner in 16-bit grayscale with a resolution of 1200 dpi. Reduction has been performed in the LINUX-MIDAS-ROMAFOT software supplemented with additional modules. The internal positional accuracy of the reduction is 0.09…0.23′′ for both coordinates and 0.27…0.37^m for the photographic magnitudes of the Tycho-2 catalog. The calculated topocentric positions of the moons are compared online with the IMCCE ephemeris data (DE405 + TASS1,7). Moon-minus-moon differential coordinates are found for most of the moons and compared with theoretical data (http://lnfm1.sai.msu.ru/neb/nss/nssephmr.htm).     

Particle motions in Maxwell’s ring dynamical systems

The focus of this contribution is an effort to review and report the main results obtained so far, concerning the periodic motions of a small body in the combined gravitational field created by a regular ν-gon arrangement of ν big bodies with equal masses, where ν > 7, and another central primary with different mass. Various types of planar periodic motions are presented and networks of characteristic curves of families are depicted, in order to show their distribution in the space of the initial conditions, as well as the evolution of their members that are also examined under the variation of the parameters of the system. Furthermore, the regions of the allowed three-dimensional motions, as well as their variation, are illustrated by means of the zero-velocity surfaces. All this new material is added to the already existing data, and completes thus the profile of the dynamical behavior of the system.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.     

General circulation and transport in Saturn’s upper troposphere and stratosphere

We present a model for the general circulation and dynamical transport in Saturn’s upper troposphere and stratosphere and derive the effective advective circulation and eddy transport coefficients required for use in two-dimensional (latitude–altitude) photochemistry–transport models. A three-dimensional Outer-Planet General Circulation Model (OPGCM) is used to generate the transport data. We find that the OPGCM adequately captures the global-scale, pole-to pole temperature contrast, but overestimates mid- and high-latitude temperatures in the summer hemisphere by ～5 K. In addition, the model reproduces the local temperature minimum seen at the equator in Cassini Composite Infrared Spectrometer (CIRS) 0.1-mbar data but not the local maximum in 1-mbar temperatures, suggesting that it is capturing the phase of Saturn’s Semiannual Oscillation associated with a temperature minimum at the equator but not the opposite phase. The meridional circulation at low latitudes is found to be dominated by a seasonally reversing Hadley circulation, characterized by upwelling near the equator, cross-equatorial flow from summer to winter hemisphere, and strong subsidence centered near 25° latitude in the winter hemisphere. The cross-equatorial flow induces an asymmetry in which the equatorial jet is found to be stronger in the winter than in the summer stratosphere. The location of the subsidence near 25°N for L_s ～ 310° coincides with local maxima in acetylene, diacetylene, and methylacetylene mixing ratios measured by Cassini/CIRS (Guerlet, S., Fouchet, T., Bézard, B., Moses, J.I., Fletcher, L.N., Simon-Miller, A.A., Flasar, F.M. [2010]. Icarus 209, 682–695). This result supports the suggestion by Guerlet et al. (2010) that the hydrocarbon abundances are enhanced at this latitude by pronounced downward transport of hydrocarbon-rich air from above. The lateral eddy diffusion coefficient is found to typically be ～10⁵–10⁶ m² s^?1 at mid-latitudes, implying meridional eddy transport time scales of order 100–1000 years.     

Constraints on Saturn’s tropospheric general circulation from Cassini ISS images

An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of the mean zonal wind, and a population of winds comparable in size to that available for the much higher contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all latitudes sampled. Upper troposphere zonal winds derived from methane band images are ～10 m s^?1 weaker than cloud level winds in the cores of eastward jets and ～5 m s^?1 stronger on either side of the jet core, i.e., eastward jets appear to broaden with increasing altitude. In westward jet regions winds are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increasing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet regions (5.2 × 10^?5 m² s^?3) and smaller in westward jet regions (1.6 × 10^?5 m² s^?3) than the global mean value (4.1 × 10^?5 m² s^?3). Overall the results are consistent with theories that suggest that the jets and the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies due to instabilities of the large-scale flow near and/or below the cloud level.