Saturn’s zonal wind profile in 2004-2009 from Cassini ISS images and its long-term variability

Five years of Cassini Imaging Science Subsystem images, from 2004 to 2009, are analyzed in this work to retrieve global zonal wind profiles of Saturn’s northern and southern hemispheres in the methane absorbing bands at 890 and 727 nm and in their respective adjacent continuum wavelengths of 939 and 752 nm. A complete view of Saturn’s global circulation, including the equator, at two pressure levels, in the tropopause (60 mbar to 250 mbar with the MT filters) and in the upper troposphere (from ∼350 mbar to ∼500 mbar with the CB filter set), is presented. Both zonal wind profiles (available at the Supplementary Material Section), show the same structure but with significant differences in the peak of the eastward jets and the equatorial region, including a region of positive vertical shear symmetrically located around the equator between the 10° < |φ_c| < 25° where zonal velocities close to the tropopause are higher than at 500 mbar. A comparison of previously published zonal wind sets obtained by Voyager 1 and 2 (1980-1981), Hubble Space Telescope, and ground-based telescopes (1990-2004) with the present Cassini profiles (2004-2009) covering a full Saturn year shows that the shape of the zonal wind profile and intensity of the jets has remained almost unchanged except at the equator, despite the seasonal insolation cycle and the variability of Saturn’s emitted power. The major wind changes occurred at equatorial latitudes, perhaps following the Great White Spot eruption in 1990. It is not evident from our study if the seasonal insolation cycle and its associated ring shadowing influence the equatorial circulation at cloud level.     

Empirical models of pressure and density in Saturn's interior: Implications for the helium concentration, its depth dependence, and Saturn's precession rate

We present ‘empirical’ models (pressure vs. density) of Saturn's interior constrained by the gravitational coefficients J₂, J₄, and J₆ for different assumed rotation rates of the planet. The empirical pressure-density profile is interpreted in terms of a hydrogen and helium physical equation of state to deduce the hydrogen to helium ratio in Saturn and to constrain the depth dependence of helium and heavy element abundances. The planet's internal structure (pressure vs. density) and composition are found to be insensitive to the assumed rotation rate for periods between 10h:32m:35s and 10h:41m:35s. We find that helium is depleted in the upper envelope, while in the high pressure region (P?1 Mbar) either the helium abundance or the concentration of heavier elements is significantly enhanced. Taking the ratio of hydrogen to helium in Saturn to be solar, we find that the maximum mass of heavy elements in Saturn's interior ranges from ∼6 to 20 M_⊕. The empirical models of Saturn's interior yield a moment of inertia factor varying from 0.22271 to 0.22599 for rotation periods between 10h:32m:35s and 10h:41m:35s, respectively. A long-term precession rate of about ^0.754″ yr⁻¹ is found to be consistent with the derived moment of inertia values and assumed rotation rates over the entire range of investigated rotation rates. This suggests that the long-term precession period of Saturn is somewhat shorter than the generally assumed value of 1.77×¹⁰⁶ years inferred from modeling and observations.     

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.     

On the vertical wind shear of Saturn's Equatorial Jet at cloud level

In this work we analyze and compare the vertical cloud structure of Saturn's Equatorial Zone in two different epochs: the first one close to the Voyagers flybys (1979-1981) and the second one in 2004, when the Cassini spacecraft entered its orbit around the planet. Our goal is to retrieve the altitude of cloud features used as zonal wind tracers in both epochs. We reanalyze three different sets of photometrically calibrated published data: ground-based in 1979, Voyager 2 PPS and ISS observations in 1981, and we analyze a new set of Hubble Space Telescope images for 2004. For all situations we reproduced the observed reflectivity by means of a similar vertical model with three layers. The results indicate the presence of a changing tropospheric haze in 1979-1981 (Ptop∼100 mbar, τ∼10) and in 2004 (Ptop∼50 mbar, τ∼15) where the tracers are embedded. According to this model the Voyager 2 ISS images locate cloud tracers moving with zonal velocities of 455 to 465 (±2) m/s at a pressure level of 360 ± 140 mbar. For HST observations, our previous works had showed cloud tracers moving with zonal wind speeds of 280±10 m/s at a pressure level of about 50±10 mbar. All these values are calculated in the same region (3°±2° N). This speed difference, if interpreted as a vertical wind shear, requires a change of per scale height, two times greater than that estimated from temperature observations. We also perform an initial guess on Cassini ISS vertical sounding levels, retrieving values compatible with HST ones and Cassini CIRS derived vertical wind shear, but not with Voyager wind measurements. We conclude that the wind speed velocity differences measured between 1979-1981 and 2004 cannot be explained as a wind shear effect alone and demand dynamical processes.     

Saturn's auroral/polar H3 infrared emission: II. A comparison with plasma flow models

We present a detailed analysis of the H⁺₃ intensity and velocity profiles crossing Saturn's auroral/polar region, as described by Stallard et al. [Stallard, T., Miller, S., Melin, H., Lystrup, M., Dougherty, M., Achilleos, N., 2007. Icarus 189, 1-13], with a view to understanding the magnetospheric processes with which they are connected. The data are not consistent with the theory that Saturn's main auroral oval is associated with corotation enforcement currents in the middle magnetosphere. This implies that the main auroral oval can be associated with the open-closed field line boundary [Cowley, S.W.H., Bunce, E.J., O'Rourke, J.M., 2004. J. Geophys. Res. 109. A05212]; a third model, by Sittler et al. [Sittler, E.C., Blanc, M.F., Richardson, J.D., 2006. J. Geophys. Res. 111. A06208] associates the main oval with centrifugal instabilities in the outer magnetosphere, but does not make predictions about ionospheric plasma flows with which we can compare our data. We do, however, tentatively identify emission at latitudes lower than the main auroral oval which may be associated with the corotation enforcement currents in the middle magnetosphere. We also find that at latitudes higher than the main auroral oval there is often a region of the ionosphere that is in rigid corotation with the planet. We suggest that this region corresponds to field lines embedded in the centre of the magnetotail which are shielded from the solar wind such that their rotation is controlled only by the neutral atmosphere.     

Long-term and large-scale viscous evolution of dense planetary rings

Planetary rings are common in the outer Solar System but their origin and long-term evolution is still a matter of debate. It is well known that viscous spreading is a major evolutionary process for rings, as it globally redistributes the disk’s mass and angular momentum, and can lead to the disk’s loosing mass by infall onto the planet or through the Roche limit. However, describing this process is highly dependent on the model used for the viscosity. In this paper we investigate the global and long-term viscous evolution of a circumplanetary disk. We have developed a simple 1D numerical code, but we use a physically realistic viscosity model derived from N-body simulations (Daisaka et al., 2001), and dependent on the disk’s local properties (surface mass density, particle size, distance to the planet). Particularly, we include the effects of gravitational instabilities (wakes) that importantly enhance the disk’s viscosity. This method allows to study the global evolution of the disk over the age of the Solar System.Common estimates of the disk’s spreading time-scales with constant viscosity significantly underestimate the rings’ lifetime. We show that, with a realistic viscosity model, an initially narrow ring undergoes two successive evolutionary stages: (1) a transient rapid spreading when the disk is self-gravitating, with the formation of a density peak inward and an outer region marginally gravitationally stable, and with an emptying time-scale proportional to (where M₀ is the disk’s initial mass), (2) an asymptotic regime where the spreading rate continuously slows down as larger parts of the disk become non-self-gravitating due to the decrease of the surface density, until the disk becomes completely non-self-gravitating. At this point its evolution dramatically slows down, with an emptying time-scale proportional to 1/M₀, which significantly increases the disk’s lifetime compared to the case with constant viscosity. We show also that the disk’s width scales like t^1/4 with the realistic viscosity model, while it scales like t^1/2 in the case of constant viscosity, resulting in much larger evolutionary time-scales in our model. We find however that the present shape of Saturn’s rings looks like a 100 million-years old disk in our simulations. Concerning Jupiter’s, Uranus’ and Neptune’s rings that are faint today, it is not likely that they were much more massive in the past and lost most of their mass due to viscous spreading alone.     

High-resolution spectroscopy of Saturn at 3 microns: CH4, CH3D, C2H2, C2H6, PH3, clouds, and haze

We report observation and analysis of a high-resolution 2.87-3.54 μm spectrum of the southern temperate region of Saturn obtained with NIRSPEC at Keck II. The spectrum reveals absorption and emission lines of five molecular species as well as spectral features of haze particles. The ν₂+ν₃ band of CH₃D is detected in absorption between 2.87 and 2.92 μm; and we derived from it a mixing ratio approximately consistent with the Infrared Space Observatory result. The ν₃ band of C₂H₂ also is detected in absorption between 2.95 and 3.05 μm; analysis indicates a sudden drop in the C₂H₂ mixing ratio at 15 mbar (130 km above the 1 bar level), probably due to condensation in the low stratosphere. The presence of the ν₃+ν₉+ν₁₁ band of C₂H₆ near 3.07 μm, first reported by Bjoraker et al. [Bjoraker, G.L., Larson, H.P., Fink, U., 1981. Astrophys. J. 248, 856-862], is confirmed, and a C₂H₆ condensation altitude of 10 mbar (140 km) in the low stratosphere is determined. We assign weak emission lines within the 3.3 μm band of CH₄ to the ν₇ band of C₂H₆, and derive a mixing ratio of 9±4×¹⁰⁻⁶ for this species. Most of the C₂H₆ 3.3 μm line emission arises in the altitude range 460-620 km (at ∼μbar pressure levels), much higher than the 160-370 km range where the 12 μm thermal molecular line emission of this species arises. At 2.87-2.90 μm the major absorber is tropospheric PH₃. The cloud level determined here and at 3.22-3.54 is 390-460 mbar (∼30 km), somewhat higher than found by Kim and Geballe [Kim, S.J., Geballe, T.R., 2005. Icarus 179, 449-458] from analysis of a low resolution spectrum. A broad absorption feature at 2.96 μm, which might be due to NH₃ ice particles in saturnian clouds, is also present. The effect of a haze layer at about 125 km (∼12 mbar level) on the 3.20-3.54 μm spectrum, which was not apparent in the low resolution spectrum, is clearly evident in the high resolution data, and the spectral properties of the haze particles suggest that they are composed of hydrocarbons.     

Modeling of global variations and ring shadowing in Saturn's ionosphere

A time-dependent one-dimensional model of Saturn's ionosphere has been developed as an intermediate step towards a fully coupled Saturn Thermosphere-Ionosphere Model (STIM). A global circulation model (GCM) of the thermosphere provides the latitude and local time dependent neutral atmosphere, from which a globally varying ionosphere is calculated. Four ion species are used (H⁺, H⁺₂, H⁺₃, and He⁺) with current cross-sections and reaction rates, and the SOLAR2000 model for the Sun's irradiance. Occultation data from the Voyager photopolarimeter system (PPS) are adapted to model the radial profile of the ultraviolet (UV) optical depth of the rings. Diurnal electron density peak values and heights are generated for all latitudes and two seasons under solar minimum and solar maximum conditions, both with and without shadowing from the rings. Saturn's lower ionosphere is shown to be in photochemical equilibrium, whereas diffusive processes are important in the topside. In agreement with previous 1-D models, the ionosphere is dominated by H⁺ and H⁺₃, with a peak electron density of ∼¹⁰⁴ electrons cm⁻³. At low- and mid-latitudes, H⁺ is the dominant ion, and the electron density exhibits a diurnal maximum during the mid-afternoon. At higher latitudes and shadowed latitudes (smaller ionizing fluxes), the diurnal maximum retreats towards noon, and the ratio of [H⁺]/[H⁺₃] decreases, with H⁺₃ becoming the dominant ion at altitudes near the peak (∼1200-1600 km) for noon-time hours. Shadowing from the rings leads to attenuation of solar flux, the magnitude and latitudinal structure of which is seasonal. During solstice, the season for the Cassini spacecraft's encounter with Saturn, attenuation has a maximum of two orders of magnitude, causing a reduction in modeled peak electron densities and total electron column contents by as much as a factor of three. Calculations are performed that explore the parameter space for charge-exchange reactions of H⁺ with vibrationally excited H₂, and for different influxes of H₂O, resulting in a maximum diurnal variation in electron density much weaker than the diurnal variations inferred from Voyager's Saturn Electrostatic Discharge (SED) measurements. Peak values of height-integrated Pedersen conductivities at high latitudes during solar maximum are modeled to be ∼42 mho in the summer hemisphere during solstice and ∼18 mho during equinox, indicating that even without ionization produced by auroral processes, magnetosphere-ionosphere coupling can be highly variable.     

Ion and neutral sources and sinks within Saturn's inner magnetosphere: Cassini results

Using ion-electron fluid parameters derived from Cassini Plasma Spectrometer (CAPS) observations within Saturn's inner magnetosphere as presented in Sittler et al. [2006a. Cassini observations of Saturn's inner plasmasphere: Saturn orbit insertion results. Planet. Space Sci., 54, 1197-1210], one can estimate the ion total flux tube content, N_IONL², for protons, H⁺, and water group ions, W⁺, as a function of radial distance or dipole L shell. In Sittler et al. [2005. Preliminary results on Saturn's inner plasmasphere as observed by Cassini: comparison with Voyager. Geophys. Res. Lett. 32(14), L14S04), it was shown that protons and water group ions dominated the plasmasphere composition. Using the ion-electron fluid parameters as boundary condition for each L shell traversed by the Cassini spacecraft, we self-consistently solve for the ambipolar electric field and the ion distribution along each of those field lines. Temperature anisotropies from Voyager plasma observations are used with (T_⊥/T_∥)_W⁺∼5 and (T_⊥/T_∥)_H⁺∼2. The radio and plasma wave science (RPWS) electron density observations from previous publications are used to indirectly confirm usage of the above temperature anisotropies for water group ions and protons. In the case of electrons we assume they are isotropic due to their short scattering time scales. When the above is done, our calculation show N_IONL² for H⁺ and W⁺ peaking near Dione's L shell with values similar to that found from Voyager plasma observations. We are able to show that water molecules are the dominant source of ions within Saturn's inner magnetosphere. We estimate the ion production rate S_ION∼10²⁷ ions/s as function of dipole L using N_H⁺, N_W⁺ and the time scale for ion loss due to radial transport τ_D and ion-electron recombination τ_REC. The ion production shows localized peaks near the L shells of Tethys, Dione and Rhea, but not Enceladus. We then estimate the neutral production rate, S_W, from our ion production rate, S_ION, and the time scale for loss of neutrals by ionization, τ_ION, and charge exchange, τ_CH. The estimated source rate for water molecules shows a pronounced peak near Enceladus’ L shell L∼4, with a value S_W∼2×10²⁸ mol/s.     

Titan impacts and escape

We report on hydrodynamic calculations of impacts of large (multi-kilometer) objects on Saturn’s moon Titan. We assess escape from Titan, and evaluate the hypothesis that escaping ejecta blackened the leading hemisphere of Iapetus and peppered the surface of Hyperion.We carried out two- and three-dimensional simulations of impactors ranging in size from 4 to 100 km diameter, impact velocities between 7 and 15 km s⁻¹, and impact angles from 0° to 75° from the vertical. We used the ZEUSMP2 hydrocode for the calculations. Simulations were made using three different geometries: three-dimensional Cartesian, two-dimensional axisymmetric spherical polar, and two-dimensional plane polar. Three-dimensional Cartesian geometry calculations were carried out over a limited domain (e.g. 240 km on a side for an impactor of size d_i = 10 km), and the results compared to ones with the same parameters done by Artemieva and Lunine (2005); in general the comparison was good. Being computationally less demanding, two-dimensional calculations were possible for much larger domains, covering global regions of the satellite (from 800 km below Titan’s surface to the exobase altitude 1700 km above the surface). Axisymmetric spherical polar calculations were carried out for vertical impacts. Two-dimensional plane-polar geometry calculations were made for both vertical and oblique impacts. In general, calculations among all three geometries gave consistent results.Our basic result is that the amount of escaping material is less than or approximately equal to the impactor mass even for the most favorable cases. Amounts of escaping material scaled most strongly as a function of velocity, with high-velocity impacts generating the largest amount, as expected. Dependence of the relative amount of escaping mass f_esc = m_esc/M_i on impactor diameter d_i was weak. Oblique impacts (impact angle θ_i > 45°) were more effective than vertical or near-vertical impacts; ratios of m_esc/M_i ∼ 1-2 were found in the simulations.     

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.     

Ion composition and chemistry in the coma of Comet 1P/Halley—A comparison between Giotto's Ion Mass Spectrometer and our ion-chemical network

In order to understand the cometary plasma environment it is important to track the closely linked chemical reactions that dominate ion evolution. We used a coupled MHD ion-chemistry model to analyze previously unpublished Giotto High Intensity Ion Mass Spectrometer (HIS-IMS) data. In this way we study the major species, but we also try to match some minor species like the CH_x and the NH_x groups. Crucial for this match is the model used for the electrons since they are important for ion-electron recombination. To further improve our results we included an enhanced density of supersonic electrons in the ion pile-up region which increases the local electron impact ionization. In this paper we discuss the results for the following important ions: C⁺, CH⁺, CH⁺₂, CH⁺₃, N⁺, NH⁺, NH⁺₂, NH⁺₃, NH⁺₄, O⁺, OH⁺, H₂O⁺, H₃O⁺, CO⁺, HCO⁺, H₃CO⁺, and CH₃OH⁺₂. We also address the inner shock which is very distinctive in our MHD model as well as in the IMS data. It is located just inside the contact surface at approximately 4550 km. Comparisons of the ion bulk flow directions and velocities from our MHD model with the data measured by the HIS-IMS give indication for a solar wind magnetic field direction different from the standard Parker angle at Halley's position. Our ion-chemical network model results are in a good agreement with the experimental data. In order to achieve the presented results we included an additional short lived inner source for the C⁺, CH⁺, and CH⁺₂ ions. Furthermore we performed our simulations with two different production rates to better match the measurements which is an indication for a change and/or an asymmetric pattern (e.g. jets) in the production rate during Giotto's fly-by at Halley's comet.     

Exchange of ejecta between Telesto and Calypso: Tadpoles, horseshoes, and passing orbits

We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn’s major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent’s escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction of launch.The slowest ejecta follow suborbital trajectories and re-impact their source moon in less than one day. Slightly faster debris barely escape their parent’s Hill sphere and are confined to tadpole orbits, librating about Tethys’ triangular Lagrange points L₄ (leading, near Telesto) or L₅ (trailing, near Calypso) with nearly the same orbital semi-major axis as Tethys, Telesto, and Calypso. These ejecta too eventually re-impact their source moon, but with a median lifetime of a few dozen years. Those which re-impact within the first 10 years or so have lifetimes near integer multiples of 348.6 days (half the tadpole period).Still faster debris with azimuthal velocity components ?10 m/s enter horseshoe orbits which enclose both L₄ and L₅ as well as L₃, but which avoid Tethys and its Hill sphere. These ejecta impact either Telesto or Calypso at comparable rates, with median lifetimes of several thousand years. However, they cannot reach Tethys itself; only the fastest ejecta, with azimuthal velocities ?40 m/s, achieve “passing orbits” which are able to encounter Tethys. Tethys accretes most of these ejecta within several years, but some 1% of them are scattered either inward to hit Enceladus or outward to strike Dione, over timescales on the order of a few hundred years.     

An analysis of auroral E region neutral winds based on incoherent scatter radar observations at Chatanika

Auroral E region neutral winds determined from incoherent scatter radar observations at Chatanika, AK, during geomagnetic disturbances (15 May 1974) are compared with detailed theoretical calculations of neutral velocities for these conditions. The theoretical velocities are obtained by numerically solving the ion and neutral momentum equations in the ion drag approximation, including coriolis and viscous forces, using observed electric fields and electron densities. Large vertical gradients are found in the calculated velocities for altitudes below about 130 km. As a consequence of this structure and fluctuations in the electron density profiles, the data analysis procedure of Brekke et al. (1973) for obtaining neutral winds from radar data is found to underestimate the wind speed by up to 40%, but it determines the direction and temporal structure reasonably well. Comparison of observed neutral velocities with calculated values shows that ion drag alone cannot account for the observations. An equation is derived to estimate the pressure gradients required to resolve the discrepancy between calculated and observed neutral winds. Accelerations due to these pressure gradients are of the same order as those due to ion drag, but at least an order of magnitude larger than those due to solar heating. Directions of the horizontal pressure gradients are consistent with expected locations of auroral heating. During geomagnetic disturbances, ion drag and auroral heating both appear to play important roles in the generation and modification of neutral winds.     

Observations of C2H6 and C2H2 in the stratosphere of Saturn

We have performed high-resolution spectral observations at mid-infrared wavelengths of C₂H₆ (12.16 μm), and C₂H₂ (13.45 μm) on Saturn. These emission features probe the stratosphere of the planet and provide information on the hydrocarbon photochemical processes taking place in that region of the atmosphere. The observations were performed using our cryogenic echelle spectrometer Celeste, in conjunction with the McMath-Pierce 1.5-m solar telescope in November and December 1994. We used Voyager IRIS CH₄ observations (7.67 μm) to derive a temperature profile on the saturnian atmosphere for the region of the stratosphere. This profile was then used in conjunction with height-dependent volume mixing ratios of each hydrocarbon to determine global abundances for ethane and acetylene. Our ground-based measurements indicate abundances of for C₂H₆ (1.0 mbar pressure level), and for C₂H₂ (1.6 mbar pressure level). We also derived new mixing ratios from the Voyager mid-latitude IRIS observations; 8.6±0.9×¹⁰⁻⁶ for C₂H₆ (0.1-3.0 mbar pressure level), and 1.6±0.2×¹⁰⁻⁷ for C₂H₂ (2.0 mbar pressure level).     

Effects of a large convective storm on Saturn's equatorial jet

Hubble Space Telescope observations revealed that Saturn's equatorial jet at the cloud level blows at ∼275 m s⁻¹ today, approximately half the ∼470 m s⁻¹ wind during the Voyager flybys in 1980-1981. Radiative transfer calculations estimate the clouds to be significantly higher today than in 1980. The higher clouds make it difficult to observationally isolate any true slowdown from the vertical wind shear because Voyager and Cassini observations show that the winds become slower with altitude. Here, we test the hypothesis that the large equatorial storm in 1990 called the Great White Spot (GWS) decelerated the equatorial jet. We first use order of magnitude estimates to show: (1) if the GWS triggers vertical momentum redistribution, a minor speed change in the troposphere can lead to a substantial stratospheric wind speed change; (2) storm-triggered turbulent mixing slows a prograde equatorial jet; and (3) a prograde equatorial jet inhibits turbulent mixing in latitude. To test whether a GWS-like large storm decelerates the equatorial jet, we perform numerical experiments using the Explicit Planetary Isentropic Coordinate (EPIC) atmosphere model. Our simulation results are consistent with our order of magnitude predictions. We show that the storm excites waves, and the waves transport westward momentum from the troposphere to the stratosphere and decelerate the equatorial jet by as much as ∼40 m s⁻¹ at the 10-mbar level. However, our results show that the storm's effect is too weak at the cloud levels to halve the jet's speed from ∼470 m s⁻¹. Our results suggest that a combination of higher clouds and a true slowdown is necessary to explain the apparent equatorial jet slowdown. We also analyze the effect of waves on the apparent cloud motions, and show that waves can influence cloud-tracking wind speed measurements.     

Vertical and meridional distribution of ethane, acetylene and propane in Saturn’s stratosphere from CIRS/Cassini limb observations

Measuring the spatial distribution of chemical compounds in Saturn’s stratosphere is critical to better understand the planet’s photochemistry and dynamics. Here we present an analysis of infrared spectra in the range 600-1400 cm⁻¹ acquired in limb geometry by the Cassini spacecraft between March 2005 and January 2008. We first determine the vertical temperature profiles from 3 to 0.01 hPa, at latitudes ranging from 70°N to 80°S. We infer a similar meridional temperature gradient at 1-2 hPa as in recent previous studies [Fletcher, L.N., Irwin, P.G.J., Teanby, N.A., Orton, G.S., Parrish, P.D., de Kok, R., Howett, C., Calcutt, S.B., Bowles, N., Taylor, F.W., 2007. Icarus 189, 457-478; Howett, C.J.A., Irwin, P.G.J., Teanby, N.A., Simon-Miller, A., Calcutt, S.B., Fletcher, L.N., de Kok, R., 2007. Icarus 190, 556-572]. We then retrieve the vertical profiles of C₂H₆ and C₂H₂ from 3 to 0.01 hPa and of C₃H₈ around 1 hPa. At 1 hPa, the meridional variation of C₂H₂ is found to follow the yearly averaged solar insolation, except for a strong equatorial mole fraction of 8×10^-7, nearly two times higher than expected. This enhancement in abundance can be explained by the descent of hydrocarbon-rich air, with a vertical wind speed at the equator of 0.25±0.1 mm/s at 1 hPa and 0.4±0.15 mm/s at 0.1 hPa. The ethane distribution is relatively uniform at 1 hPa, with only a moderate 25% increase from 35°S to 80°S. Propane is found to increase from north to south by a factor of 1.9, suggesting that its lifetime may be shorter than Saturn’s year at 1 hPa. At high altitudes (1 Pa), C₂H₂ and C₂H₆ abundances depart significantly from the photochemical model predictions of Moses and Greathouse [Moses, J.I., Greathouse, T.K., 2005. J. Geophys. Res. 110, 9007], except at high southern latitudes (62, 70 and 80°S) and near the equator. The observed abundances are found strongly depleted in the 20-40°S region and enhanced in the 20-30°N region, the latter coinciding with the ring’s shadow. We favor a dynamical explanation for these anomalies.     

The first detection of propane on Saturn

We report the first detection of propane, C₃H₈, in Saturn's stratosphere. Observations taken on September 8, 2002 UT at NASA's IRTF using TEXES, show multiple emission lines due to the 748 cm⁻¹ν₂₁ band of C₃H₈. Using a line-by-line radiative transfer code, we are able to fit the data by scaling the propane vertical mixing ratio profile from the photochemical model of Moses et al. [2000. Icarus 143, 244-298]. Multiplicative factors of 0.7 and 0.65 are required to fit the −20° and −80° planetocentric latitude spectra. The resultant profiles are characterized by a 5 mbar mixing ratio of 2.7±0.8×¹⁰⁻⁸ at −20° and at −80° latitude. These results suggest that the time scale for meridional circulation lies between the net photochemical lifetimes of C₂H₂ and C₃H₈, ≈30-600 years.     

Waves in Cassini UVIS stellar occultations: 2. The C ring

We performed a complete wavelet analysis of Saturn’s C ring on 62 stellar occultation profiles. These profiles were obtained by Cassini’s Ultraviolet Imaging Spectrograph High Speed Photometer. We used a WWZ wavelet power transform to analyze them. With a co-adding process, we found evidence of 40 wavelike structures, 18 of which are reported here for the first time. Seventeen of these appear to be propagating waves (wavelength changing systematically with distance from Saturn). The longest new wavetrain in the C ring is a 52-km-long wave in a plateau at 86,397 km. We produced a complete map of resonances with external satellites and possible structures rotating with Saturn’s rotation period up to the eighth order, allowing us to associate a previously observed wave with the Atlas 2:1 inner Lindblad resonance (ILR) and newly detected waves with the Mimas 6:2 ILR and the Pandora 4:2 ILR. We derived surface mass densities and mass extinction coefficients, finding σ = 0.22(±0.03) g cm⁻² for the Atlas 2:1 ILR, σ = 1.31(±0.20) g cm⁻² for the Mimas 6:2 ILR, and σ = 1.42(±0.21) g cm⁻² for the Pandora 4:2 ILR. We determined a range of mass extinction coefficients (κ = τ/σ) for the waves associated with resonances with κ = 0.13 (±0.03) to 0.28(±0.06) cm² g⁻¹, where τ is the optical depth. These values are higher than the reported values for the A ring (0.01-0.02 cm² g⁻¹) and the Cassini Division (0.07-0.12 cm² g⁻¹ from Colwell et al. (Colwell, J.E., Cooney, J.H., Esposito, L.W., Srem?evi?, M. [2009]. Icarus 200, 574-580)). We also note that the mass extinction coefficient is probably not constant across the C ring (in contrast to the A ring and the Cassini Division): it is systematically higher in the plateaus than elsewhere, suggesting smaller particles in the plateaus. We present the results of our analysis of these waves in the C ring and estimate the mass of the C ring to be between3.7(±0.9) × 10¹⁶ kg and 7.9(±2.0) × 10¹⁶ kg (equivalent to an icy satellite of radius between 28.0(±2.3) km and 36.2(±3.0) km with a density of 400 kg m⁻³, close to that of Pan or Atlas). Using the ring viscosity derived from the wave damping length, we also estimate the vertical thickness of the C ring between 1.9(±0.4) m and 5.6(±1.4) m, comparable to the vertical thickness of the Cassini Division.