Orbits and masses of Saturn's coorbital and F-ring shepherding satellites

We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm⁻³ are , , , and , respectively.     

Thermal infrared and optical observations of four near-Earth asteroids

We present thermal infrared photometry and spectrophotometry of four near-Earth asteroids (NEAs), namely (433) Eros, (66063) 1998 RO₁, (137032) 1998 UO₁, and (138258) 2000 GD₂, using the United Kingdom Infrared Telescope (UKIRT) in 2002. For two objects, i.e. (433) Eros and (137032) 1998 UO₁, quasi-simultaneous optical observations were also obtained, using the Jacobus Kapteyn Telescope (JKT). For (127032) 1998 UO₁, we obtain a rotation period P=3.0±0.1 h and an absolute visual magnitude HV=16.7±0.4. The Standard Thermal Model (STM), Fast Rotating Model (FRM) and near-Earth asteroid Thermal Model (NEATM) have been fitted to the IR fluxes to determine effective diameters Deff, geometric albedos pv, and beaming parameters η. The derived values are (433) Eros: Deff=23.3±3.5 km (at lightcurve maximum), pv=0.24±0.07, η=0.95±0.19; (66063) 1998 RO₁: , ; (137032) 1998 UO₁: Deff<1.13 km, pv>0.29; (138258) 2000 GD₂: Deff=0.27±0.04 km, , η=0.74±0.15. (66063) 1998 RO₁ is a binary asteroid from lightcurve characteristics [Pravec, P., and 56 colleagues, 2006. Icarus 181, 63-93] and we estimate the effective diameter of the primary (Dp) and secondary (Ds) components: and . The diameter and albedo of (138258) 2000 GD₂ are consistent with the trend of decreasing diameter for S- and Q-type asteroids found by Delbó et al. [Delbó, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. A possible trend of increasing beaming parameter with diameter for small (less than about 3 km) S- and Q-type asteroids is found.     

Detection of C2HD and the D/H ratio on Titan

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

Saturn’s wayward shepherds: the peregrinations of Prometheus and Pandora

Saturn’s narrow F ring is flanked by two nearby small satellites, Prometheus and Pandora, discovered in Voyager images taken in 1980 and 1981 (Synnott et al., 1983, Icarus 53, 156-158). Observations with the Hubble Space Telescope (HST) during the ring plane crossings (RPX) of 1995 led to the unexpected finding that Prometheus was ∼19° behind its predicted orbital longitude, based on the Synnott et al. (1983) Voyager ephemeris (Bosh and Rivkin, 1996 Science 272, 518-521; Nicholson et al., 1996, Science 272, 509-515). Whereas Pandora was at its predicted location in August 1995, McGhee (2000, Ph.D. thesis, Cornell University) found from the May and November 1995 RPX data that Pandora also deviates from the Synnott et al. (1983) Voyager ephemeris. Using archival HST data from 1994, previously unexamined RPX images, and a large series of targeted WFPC2 observations between 1996 and 2002, we have determined highly accurate sky-plane positions for Prometheus, Pandora, and nine other satellites found in our images. We compare the Prometheus and Pandora measurements to the predictions of substantially revised and improved ephemerides for the two satellites based on an extensive analysis of a large set of Voyager images (Murray et al., 2000, Bull. Am. Astron. Soc. 32, 1090; Evans, 2001 Ph.D. thesis, Queen Mary College). From December 1994 to December 2000, Prometheus’ orbital longitude lag was changing by −0.71° year⁻¹ relative to the new Voyager ephemeris. In contrast, Pandora is ahead of the revised Voyager prediction. From 1994 to 2000, its longitude offset changed by +0.44° year⁻¹, showing in addition an ∼585 day oscillatory component with amplitude Δλ_CR0 = 0.65 ± 0.07° whose phase matches the expected perturbation due to the nearby 3:2 corotation resonance with Mimas, modulated by the 71-year libration in the longitude of Mimas due to its 4:2 resonance with Tethys. We determine orbital elements for freely precessing equatorial orbits from fits to the 1994-2000 HST observations, from which we conclude that Prometheus’ semimajor axis was 0.31 km larger, and Pandora’s was 0.20 km smaller, than during the Voyager epoch. Subsequent observations in 2001-2002 reveal a new twist in the meanderings of these satellites: Prometheus’ mean motion changed suddenly by an additional −0.77° year⁻¹, equivalent to a further increase in semimajor axis of 0.33 km, at the same time that Pandora’s mean motion changed by +0.92° year⁻¹, corresponding to a change of −0.42 km in its semimajor axis. There is an apparent anticorrelation of the motions of these two moons seen in the 2001-2002 observations, as well as over the 20-year interval since the Voyager epoch. This suggests a common origin for their wanderings, perhaps through direct exchange of energy between the satellites as the result of resonances, possibly involving the F ring.     

The C/C isotopic ratio in Titan hydrocarbons from Cassini/CIRS infrared spectra

We have analyzed infrared spectra of Titan recorded by the Cassini Composite Infrared Spectrometer (CIRS) to measure the isotopic ratio ¹²C/¹³C in each of three chemical species in Titan's stratosphere: CH₄, C₂H₂ and C₂H₆. This is the first measurement of ¹²C/¹³C in any C₂ molecule on Titan, and the first measurement of ¹²CH₄/¹³CH₄ (non-deuterated) on Titan by remote sensing. Our spectra cover five widely-spaced latitudes, 65° S to 71° N and we have searched for both latitude variability of ¹²C/¹³C within a given species, and also for differences between the ¹²C/¹³C in the three gases. For CH₄ alone, we find (1-σ), essentially in agreement with the ¹²CH₄/¹³CH₄ measured by the Huygens Gas Chromatograph/Mass Spectrometer instrument (GCMS) [Niemann, H.B., and 17 colleagues, 2005. Nature 438, 779-784]: 82.3±1.0, and also with measured values in H¹³CN and ¹³CH₃D by CIRS at lower precision [Bézard, B., Nixon, C., Kleiner, I., Jennings, D., 2007. Icarus 191, 397-400; Vinatier, S., Bézard, B., Nixon, C., 2007. Icarus 191, 712-721]. For the C₂ species, we find in C₂H₂ and 89.8±7.3 in C₂H₆, a possible trend of increasingly value with molecular mass, although these values are both compatible with the Huygens GCMS value to within error bars. There are no convincing trends in latitude. Combining all fifteen measurements, we obtain a value of , also compatible with GCMS. Therefore, the evidence is mounting that ¹²C/¹³C is some 8% lower on Titan than on the Earth (88.9, inorganic standard), and lower than typical for the outer planets (88±7 [Sada, P.V., McCabe, G.H., Bjoraker, G.L., Jennings, D.E., Reuter, D.C., 1996. Astrophys. J. 472, 903-907]). There is no current model for this enrichment, and we discuss several mechanisms that may be at work.     

Parametric analysis of modeled ion escape from Mars

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

Astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn

We describe the astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn discovered in images obtained using the Imaging Science Subsystem (ISS) of the Cassini spacecraft. Included are details of 63 observations, of which 28 were obtained with Cassini's narrow-angle camera (NAC) and 35 using its wide-angle camera (WAC), covering an observation time-span of approximately 3 years. We estimate the diameter of Anthe to be ∼1.8 km. Orbit modeling based on a numerical integration of the full equations of motion fitted to the observations show that Anthe is in a first-order 11:10 mean motion resonance with Mimas. Two resonant arguments are librating: ?₁=11λ^′−10λ−? and ?₂=11λ^′−10λ−?^′−Ω^′+Ω, where λ, ? and Ω refer to the mean longitude, longitude of pericenter and longitude of ascending node of Mimas and Anthe, with the primed quantities corresponding to Anthe. These resonances cause periodic variations in the orbital elements. The semi-major axis varies by ±26 km over a 913-day period. Anthe is also close to a second-order eccentricity-type mean motion resonant relationship of the form 77:75 with Methone. Since Methone is also in a first-order resonance with Mimas [Spitale, J.N., Jacobson, R.A., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 692-710], an additional indirect perturbation exists between Methone and Anthe via Mimas. Neither effect is detectable in the orbit fitting and the short-term dynamical evolution of Anthe is dominated by the Mimas-Anthe resonances alone. The expected modulation effect from the Mimas-Tethys 4:2 inclination resonance is also insignificant over this time period. By including Cassini ISS observations of Mimas in the numerical integration fit, we estimate the GM of Mimas to be , consistent with Jacobson et al. [Jacobson, R.A., Spitale, J., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 711-713].     

Detection of CH3D on Titan

We report the detection of ¹³CH₃D in Titan's stratosphere from Cassini/CIRS infrared spectra near 8.7 μm. Fitting simultaneously the ν₆ bands of both ¹³CH₃D and ¹²CH₃D and the ν₄ band of CH₄, we derive a D/H ratio equal to and a ¹²C/¹³C ratio in deuterated methane of , consistent with that measured in normal methane.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

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

The orbit and size distribution of small Solar System objects orbiting the Sun interior to the Earth's orbit

We present the first observational measurement of the orbit and size distribution of small Solar System objects whose orbits are wholly interior to the Earth's (Inner Earth Objects, IEOs, with aphelion <0.983 AU). We show that we are able to model the detections of near-Earth objects (NEO) by the Catalina Sky Survey (CSS) using a detailed parameterization of the CSS survey cadence and detection efficiencies as implemented within the Jedicke et al. [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J.M., Bottke, W.F., 2003. Icarus 161, 17-33] survey simulator and utilizing the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] model of the NEO population's size and orbit distribution. We then show that the CSS detections of 4 IEOs are consistent with the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] IEO model. Observational selection effects for the IEOs discovered by the CSS were then determined using the survey simulator in order to calculate the corrected number and H distribution of the IEOs. The actual number of IEOs with H<18 (21) is 36±26 (530±240) and the slope of the H magnitude distribution (∝¹⁰αH) for the IEOs is . The slope is consistent with previous measurements for the NEO population of αNEO=0.35±0.02 [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] and αNEO=0.39±0.013 [Stuart, J.S., Binzel, R.P., 2004. Icarus 170, 295-311]. Based on the agreement between the predicted and observed IEO orbit and absolute magnitude distributions there is no indication of any non-gravitational effects (e.g. Yarkovsky, tidal disruption) affecting the known IEO population.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

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.     

Near infrared imaging spectroscopy of Venus with the Anglo-Australian Telescope

We have obtained full-disk spatially resolved spectra of the Venus nightside at near-infrared wavelengths during July 2007 using the Anglo-Australian Telescope and Infrared Imager and Spectrograph 2 (IRIS2). The data have been used to map the intensity and rotational temperature of the O₂(a¹Δ_g) airglow band at . The temperatures agree with those obtained in earlier IRIS2 observations and are significantly higher than expected from the Venus International Reference Atmosphere (VIRA) profile. We also report the detection of the corresponding ν=0-1O₂ airglow band at with a similar spatial distribution to the ν=0-0 band. Observations in the thermal window have been used to image surface topography using two different methods of cloud correction. We have also obtained images that can be used to study cloud motion.     

Titan's hydrodynamically escaping atmosphere: Escape rates and the structure of the exobase region

In Strobel [Strobel, D.F., 2008. Icarus, 193, 588-594] a mass loss rate from Titan's upper atmosphere, , was calculated for a single constituent, N₂ atmosphere by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating due to CH₄ absorption. It was estimated, but not proven, that the hydrodynamic mass loss is essentially CH₄ and H₂ escape. Here the individual conservation of momentum equations for the three major components of the upper atmosphere (N₂, CH₄, H₂) are solved in the low Mach number limit and compared with Cassini Ion Neutral Mass Spectrometer (INMS) measurements to demonstrate that light gases (CH₄, H₂) preferentially escape over the heavy gas (N₂). The lightest gas (H₂) escapes with a flux 99% of its limiting flux, whereas CH₄ is restricted to ?75% of its limiting flux because there is insufficient solar power to support escape at the limiting rate. The respective calculated H₂ and CH₄ escape rates are 9.2×¹⁰²⁷ and 1.7×¹⁰²⁷ s−1, for a total of . From the calculated densities, mean free paths of N₂, CH₄, H₂, and macroscopic length scales, an extended region above the classic exobase is inferred where frequent collisions are still occurring and thermal heat conduction can deliver power to lift the escaping gas out of the gravitational potential well. In this region rapid acceleration of CH₄ outflow occurs. With the thermal structure of Titan's thermosphere inferred from INMS data by Müller-Wodarg et al. [Müller-Wodarg, I.C.F., Yelle, R.V., Cui, J., Waite Jr., J.H., 2008. J. Geophys. Res. 113, doi:10.1029/2007JE003033. E10005], in combination with calculated temperature profiles that include sputter induced plasma heating at the exobase, it is concluded that on average that the integrated, globally average, orbit-averaged, plasma heating rate during the Cassini epoch does not exceed ().     

Turbulent viscosity and lifetime of Saturn's rings

The viscosity (the angular momentum flux) in the disk of mutually gravitating particles of Saturn's rings is investigated. The hydrodynamic theory of the gravitational Jeans-type instability of small gravity perturbations (e.g., those produced by spontaneous disturbances) of the disk is developed. It is suggested that in such a system the hydrodynamic turbulence may arise as a result of the instability. The turbulence is related to stochastic motions of “fluid” elements. The objective of this paper is to show that in the Jeans-unstable Saturnian ring disk the turbulent viscosity may exceed the ordinary microscopic viscosity substantially. The main result of local N-body simulations of planetary rings by Daisaka et al. (2001. Viscosity in a dense planetary ring with self-gravitating particles. Icarus 154, 296-312) is explained: in the presence of the gravitationally unstable density waves, the effective turbulent viscosity ν_eff is given as ν_eff=CG²Σ²/Ω³, where G, Σ, and Ω are the gravitational constant, the surface mass density of a ring, and the angular velocity, respectively, and the nondimensional correction factor C≈10. We argue that both Saturn's main rings and their irregular of the order of 100 m or even less fine-scale structure (being recurrently created and destroyed on the time scale of an order of Keplerian period ) are not likely much younger than the solar system.     

Chaotic motions of prometheus and pandora

Recent HST images of the saturnian satellites Prometheus and Pandora show that their longitudes deviate from predictions of ephemerides based on Voyager images. Currently Prometheus is lagging and Pandora leading these predictions by somewhat more than 20°. We show that these discrepancies are fully accounted for by gravitational interactions between the two satellites. These peak every 24.8 days at conjunctions and excite chaotic perturbations. The Lyapunov exponent for the Prometheus-Pandora system is of order 0.3 year⁻¹ for satellite masses based on a nominal density of 0.63 g cm⁻³. Interactions are strongest when the orbits come closest together. This happens at intervals of 6.2 years when their apses are antialigned. In this context, we note the sudden changes of opposite signs in the mean motions of Prometheus and Pandora at the end of 2000 occurred around the time their apsidal lines were antialigned.     

The Titan N/N and C/C isotopic ratios in HCN from Cassini/CIRS

We report the detection of H¹³CN and HC¹⁵N in mid-infrared spectra recorded by the Composite Infrared Spectrometer (CIRS) aboard Cassini, along with the determination of the ¹²C/¹³C and ¹⁴N/¹⁵N isotopic ratios. We analyzed two sets of limb spectra recorded near 13-15° S (Tb flyby) and 83° N (T4 flyby) at 0.5 cm⁻¹ resolution. The spectral range 1210-1310 cm⁻¹ was used to retrieve the temperature profile in the range 145-490 km at 13° S and 165-300 km at 83° N. These two temperature profiles were then incorporated in the atmospheric model to retrieve the abundance profile of H¹²C¹⁴N, H¹³CN and HC¹⁵N from their bands at 713, 706 and 711 cm⁻¹, respectively. The HCN abundance profile was retrieved in the range 90-460 km at 15° S and 165-305 km at 83° N. There is no evidence for vertical variations of the isotopic ratios. Constraining the isotopic abundance profiles to be proportional to the HCN one, we find at 15° S, and at 83° N, two values that are statistically consistent. A combination of these results yields a ¹²C/¹³C value equal to 75±12. This global result, as well as the 15° S one, envelop the value in Titan's methane (82.3±1) [Niemann, H.B., and 17 colleagues, 2005. Nature 438, 779-784] measured at 10° S and is slightly lower than the terrestrial inorganic standard value (89). The ¹⁴N/¹⁵N isotopic ratio is found equal to at 15° S and at 83° N. Combining the two values yields ¹⁴N/¹⁵N = 56 ± 8, which corresponds to an enrichment in ¹⁵N of about 4.9 compared with the terrestrial ratio. These results agree with the values obtained from previous ground-based millimeter observations [Hidayat, T., Marten, A., Bézard, B., Gautier, D., Owen, T., Matthews, H.E., Paubert, G., 1997. Icarus 126, 170-182; Marten, A., Hidayat, T., Biraud, Y., Moreno, R., 2002. Icarus 158, 532-544]. The ¹⁵N/¹⁴N ratio found in HCN is ∼3 times higher than in N₂ [Niemann, H.B., and 17 colleagues, 2005. Nature 438, 779-784], which implies a large fractionation process in the HCN photochemistry.