Cassini SAR, radiometry, scatterometry and altimetry observations of Titan’s dune fields

Large expanses of linear dunes cover Titan’s equatorial regions. As the Cassini mission continues, more dune fields are becoming unveiled and examined by the microwave radar in all its modes of operation (SAR, radiometry, scatterometry, altimetry) and with an increasing variety of observational geometries. In this paper, we report on Cassini’s radar instrument observations of the dune fields mapped through May 2009 and present our key findings in terms of Titan’s geology and climate. We estimate that dune fields cover ∼12.5% of Titan’s surface, which corresponds to an area of ∼10 million km², roughly the area of the United States. If dune sand-sized particles are mainly composed of solid organics as suggested by VIMS observations (Cassini Visual and Infrared Mapping Spectrometer) and atmospheric modeling and supported by radiometry data, dune fields are the largest known organic reservoir on Titan. Dune regions are, with the exception of the polar lakes and seas, the least reflective and most emissive features on this moon. Interestingly, we also find a latitudinal dependence in the dune field microwave properties: up to a latitude of ∼11°, dune fields tend to become less emissive and brighter as one moves northward. Above ∼11° this trend is reversed. The microwave signatures of the dune regions are thought to be primarily controlled by the interdune proportion (relative to that of the dune), roughness and degree of sand cover. In agreement with radiometry and scatterometry observations, SAR images suggest that the fraction of interdunes increases northward up to a latitude of ∼14°. In general, scattering from the subsurface (volume scattering and surface scattering from buried interfaces) makes interdunal regions brighter than the dunes. The observed latitudinal trend may therefore also be partially caused by a gradual thinning of the interdunal sand cover or surrounding sand sheets to the north, thus allowing wave penetration in the underlying substrate. Altimetry measurements over dunes have highlighted a region located in the Fensal dune field (∼5° latitude) where the icy bedrock of Titan is likely exposed within smooth interdune areas. The hemispherical assymetry of dune field properties may point to a general reduction in the availability of sediments and/or an increase in the ground humidity toward the north, which could be related to Titan’s asymmetric seasonal polar insolation. Alternatively, it may indicate that either the wind pattern or the topography is less favorable for dune formation in Titan’s northern tropics.     

Elastic ice shells of synchronous moons: Implications for cracks on Europa and non-synchronous rotation of Titan

A number of synchronous moons are thought to harbor water oceans beneath their outer ice shells. A subsurface ocean frictionally decouples the shell from the interior. This has led to proposals that a weak tidal or atmospheric torque might cause the shell to rotate differentially with respect to the synchronously rotating interior. Applications along these lines have been made to Europa and Titan. However, the shell is coupled to the ocean by an elastic torque. As a result of centrifugal and tidal forces, the ocean would assume an ellipsoidal shape with its long axis aligned toward the parent planet. Any displacement of the shell away from its equilibrium position would induce strains thereby increasing its elastic energy and giving rise to an elastic restoring torque. In the investigation reported on here, the elastic torque is compared with the tidal torque acting on Europa and the atmospheric torque acting on Titan.Regarding Europa, it is shown that the tidal torque is far too weak to produce stresses that could fracture the ice shell, thus refuting an idea that has been widely advocated. Instead, it is suggested that the cracks arise from time-dependent stresses due to non-hydrostatic gravity anomalies from tidally driven, episodic convection in the satellite’s interior.Two years of Cassini RADAR observations of Titan’s surface have been interpreted as implying an angular displacement of ∼0.24° relative to synchronous rotation. Compatibility of the amplitude and phase of the observed non-synchronous rotation with estimates of the atmospheric torque requires that Titan’s shell be decoupled from its interior. We find that the elastic torque balances the seasonal atmospheric torque at an angular displacement ?0.05°, effectively coupling the shell to the interior. Moreover, if Titan’s surface were spinning faster than synchronous, the tidal torque tending to restore synchronous rotation would almost certainly be larger than the atmospheric torque. There must either be a problem with the interpretation of the radar observations, or with our basic understanding of Titan’s atmosphere and/or interior.     

Titan’s cloud seasonal activity from winter to spring with Cassini/VIMS

Since Saturn orbital insertion in July 2004, the Cassini orbiter has been observing Titan throughout most of the northern winter season (October 2002–August 2009) and the beginning of spring, allowing a detailed monitoring of Titan’s cloud coverage at high spatial resolution with close flybys on a monthly basis. This study reports on the analysis of all the near-infrared images of Titan’s clouds acquired by the Visual and Infrared Mapping Spectrometer (VIMS) during 67 targeted flybys of Titan between July 2004 and April 2010.The VIMS observations show numerous sporadic clouds at southern high and mid-latitudes, rare clouds in the equatorial region, and reveal a long-lived cloud cap above the north pole, ubiquitous poleward of 60°N. These observations allow us to follow the evolution of the cloud coverage during almost a 6-year period including the equinox, and greatly help to further constrain global circulation models (GCMs). After 4 years of regular outbursts observed by Cassini between 2004 and 2008, southern polar cloud activity started declining, and completely ceased 1 year before spring equinox. The extensive cloud system over the north pole, stable between 2004 and 2008, progressively fractionated and vanished as Titan entered into northern spring. At southern mid-latitudes, clouds were continuously observed throughout the VIMS observing period, even after equinox, in a latitude band between 30°S and 60°S. During the whole period of observation, only a dozen clouds were observed closer to the equator, though they were slightly more frequent as equinox approached.We also investigated the distribution of clouds with longitude. We found that southern polar clouds, before disappearing in mid-2008, were systematically concentrated in the leading hemisphere of Titan, in particular above and to the east of Ontario Lacus, the largest reservoir of hydrocarbons in the area. Clouds are also non-homogeneously distributed with longitude at southern mid-latitudes. The n = 2-mode wave pattern of the distribution, observed since 2003 by Earth-based telescopes and confirmed by our Cassini observations, may be attributed to Saturn’s tides.Although the latitudinal distribution of clouds is now relatively well reproduced and understood by the GCMs, the non-homogeneous longitudinal distributions and the evolution of the cloud coverage with seasons still need investigation. If the observation of a few single clouds at the tropics and at northern mid-latitudes late in winter and at the start of spring cannot be further interpreted for the moment, the obvious shutdown of the cloud activity at Titan’s poles provides clear signs of the onset of the general circulation turnover that is expected to accompany the beginning of Titan’s northern spring. According to our GCM, the persistence of clouds at certain latitudes rather suggests a ‘sudden’ shift in near future of the meteorology into the more illuminated hemisphere. Finally, the observed seasonal change in cloud activity occurred with a significant time lag that is not predicted by our model. This may be due to an overall methane humidity at Titan’s surface higher than previously expected.     

Wave constraints for Titan’s Jingpo Lacus and Kraken Mare from VIMS specular reflection lightcurves

Stephan et al. (Stephan, K. et al. [2010]. Geophys. Res. Lett. 37, 7104-+.) first saw the glint of sunlight specularly reflected off of Titan’s lakes. We develop a quantitative model for analyzing the photometric lightcurve generated during a flyby in which the specularly reflected light flux depends on the fraction of the solar specular footprint that is covered by liquid. We allow for surface waves that spread out the geographic specular intensity distribution. Applying the model to the VIMS T58 observations shows that the waves on Jingpo Lacus must have slopes of no greater than 0.15°, two orders of magnitude flatter than waves on Earth’s oceans. Combining the model with theoretical estimates of the intensity of the specular reflection allows a tighter constraint on the waves: ?0.05°. Residual specular signal while the specular point lies on land implies that either the land is wetted, the wave slope distribution is non-Gaussian, or that 5% of the land off the southwest edge of Jingpo Lacus is covered in puddles. Another specular sequence off of Kraken Mare acquired during Cassini’s T59 flyby shows rapid flux changes that the static model cannot reproduce. Points just 1 min apart vary in flux by more than a factor of two. The present dataset does not uniquely determine the mechanism causing these rapid changes. We suggest that changing wind conditions, kilometer-wavelength waves, or moving clouds could account for the variability. Future specular observations should be designed with a fast cadence, at least 6 points per minute, in order to differentiate between these hypotheses. Such new data will further constrain the nature of Titan’s lakes and their interactions with Titan’s atmosphere.     

A primordial origin for the atmospheric methane of Saturn’s moon Titan

The origin of Titan’s atmospheric methane is a key issue for understanding the origin of the saturnian satellite system. It has been proposed that serpentinization reactions in Titan’s interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan’s planetesimals before its formation. Here, we point out that serpentinization reactions in Titan’s interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan’s water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan’s interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite’s planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan’s interior can be up to ∼1300 times the current mass of atmospheric methane.     

Latitudinal variations in Titan’s methane and haze from Cassini VIMS observations

We analyze observations taken with Cassini’s Visual and Infrared Mapping Spectrometer (VIMS), to determine the current methane and haze latitudinal distribution between 60°S and 40°N. The methane variation was measured primarily from its absorption band at 0.61 μm, which is optically thin enough to be sensitive to the methane abundance at 20-50 km altitude. Haze characteristics were determined from Titan’s 0.4-1.6 μm spectra, which sample Titan’s atmosphere from the surface to 200 km altitude. Radiative transfer models based on the haze properties and methane absorption profiles at the Huygens site reproduced the observed VIMS spectra and allowed us to retrieve latitude variations in the methane abundance and haze. We find the haze variations can be reproduced by varying only the density and single scattering albedo above 80 km altitude. There is an ambiguity between methane abundance and haze optical depth, because higher haze optical depth causes shallower methane bands; thus a family of solutions is allowed by the data. We find that haze variations alone, with a constant methane abundance, can reproduce the spatial variation in the methane bands if the haze density increases by 60% between 20°S and 10°S (roughly the sub-solar latitude) and single scattering absorption increases by 20% between 60°S and 40°N. On the other hand, a higher abundance of methane between 20 and 50 km in the summer hemisphere, as much as two times that of the winter hemisphere, is also possible, if the haze variations are minimized. The range of possible methane variations between 27°S and 19°N is consistent with condensation as a result of temperature variations of 0-1.5 K at 20-30 km. Our analysis indicates that the latitudinal variations in Titan’s visible to near-IR albedo, the north/south asymmetry (NSA), result primarily from variations in the thickness of the darker haze layer, detected by Huygens DISR, above 80 km altitude. If we assume little to no latitudinal methane variations we can reproduce the NSA wavelength signatures with the derived haze characteristics. We calculate the solar heating rate as a function of latitude and derive variations of ∼10-15% near the sub-solar latitude resulting from the NSA. Most of the latitudinal variations in the heating rate stem from changes in solar zenith angle rather than compositional variations.     

Distribution and interplay of geologic processes on Titan from Cassini radar data

The Cassini Titan Radar Mapper is providing an unprecedented view of Titan’s surface geology. Here we use Synthetic Aperture Radar (SAR) image swaths (Ta-T30) obtained from October 2004 to December 2007 to infer the geologic processes that have shaped Titan’s surface. These SAR swaths cover about 20% of the surface, at a spatial resolution ranging from ∼350 m to ∼2 km. The SAR data are distributed over a wide latitudinal and longitudinal range, enabling some conclusions to be drawn about the global distribution of processes. They reveal a geologically complex surface that has been modified by all the major geologic processes seen on Earth - volcanism, tectonism, impact cratering, and erosion and deposition by fluvial and aeolian activity. In this paper, we map geomorphological units from SAR data and analyze their areal distribution and relative ages of modification in order to infer the geologic evolution of Titan’s surface. We find that dunes and hummocky and mountainous terrains are more widespread than lakes, putative cryovolcanic features, mottled plains, and craters and crateriform structures that may be due to impact. Undifferentiated plains are the largest areal unit; their origin is uncertain. In terms of latitudinal distribution, dunes and hummocky and mountainous terrains are located mostly at low latitudes (less than 30°), with no dunes being present above 60°. Channels formed by fluvial activity are present at all latitudes, but lakes are at high latitudes only. Crateriform structures that may have been formed by impact appear to be uniformly distributed with latitude, but the well-preserved impact craters are all located at low latitudes, possibly indicating that more resurfacing has occurred at higher latitudes. Cryovolcanic features are not ubiquitous, and are mostly located between 30° and 60° north. We examine temporal relationships between units wherever possible, and conclude that aeolian and fluvial/pluvial/lacustrine processes are the most recent, while tectonic processes that led to the formation of mountains and Xanadu are likely the most ancient.     

Impact craters on Titan

Five certain impact craters and 44 additional nearly certain and probable ones have been identified on the 22% of Titan’s surface imaged by Cassini’s high-resolution radar through December 2007. The certain craters have morphologies similar to impact craters on rocky planets, as well as two with radar bright, jagged rims. The less certain craters often appear to be eroded versions of the certain ones. Titan’s craters are modified by a variety of processes including fluvial erosion, mass wasting, burial by dunes and submergence in seas, but there is no compelling evidence of isostatic adjustments as on other icy moons, nor draping by thick atmospheric deposits. The paucity of craters implies that Titan’s surface is quite young, but the modeled age depends on which published crater production rate is assumed. Using the model of Artemieva and Lunine (2005) suggests that craters with diameters smaller than about 35 km are younger than 200 million years old, and larger craters are older. Craters are not distributed uniformly; Xanadu has a crater density 2-9 times greater than the rest of Titan, and the density on equatorial dune areas is much lower than average. There is a small excess of craters on the leading hemisphere, and craters are deficient in the north polar region compared to the rest of the world. The youthful age of Titan overall, and the various erosional states of its likely impact craters, demonstrate that dynamic processes have destroyed most of the early history of the moon, and that multiple processes continue to strongly modify its surface. The existence of 24 possible impact craters with diameters less than 20 km appears consistent with the Ivanov, Basilevsky and Neukum (1997) model of the effectiveness of Titan’s atmosphere in destroying most but not all small projectiles.     

Determining a tilt in Titan’s north-south albedo asymmetry from Cassini images

Analysis of Titan’s hemispheric brightness asymmetry from mapped Cassini images reveals an axis of symmetry that is tilted with respect to the rotational axis of the solid body. Twenty images taken from 2004 through 2007 show a mean axial offset of 3.8 ± 0.9° relative to the solid body’s pole, directed 79 ± 24° to the west of the sub-solar longitude. These values are consistent with recent measurements of an implied atmospheric spin axis determined from isothermal mapping by [Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A., 2008. Icarus 197, 549-555].     

Particle size and abundance of HC3N ice in Titan’s lower stratosphere at high northern latitudes

Up to now, there has been no corroboration from Cassini CIRS of the Voyager IRIS-discovery of cyanoacetylene (HC₃N) ice in Titan’s thermal infrared spectrum. We report the first compelling spectral evidence from CIRS for the ν₆ HC₃N ice feature at 506 cm⁻¹ at latitudes 62°N and 70°N, from which we derive particle sizes and column abundances in Titan’s lower stratosphere. We find mean particle radii of 3.0 μm and 2.3 μm for condensed HC₃N at 62°N and 70°N, respectively, and corresponding ice phase molecular column abundances in the range 1-10 × 10¹⁶ mol cm⁻². Only upper limits for cloud abundances can be established at latitudes of 85°N, 55°N, 30°N, 10°N, and 15°S. Under the assumption that cloud tops coincide with the uppermost levels at which HC₃N vapor saturates, we infer geometric thicknesses for the clouds equivalent to 10-20 km or so, with tops at 165 km and 150 km at 70°N and 62°N, respectively.     

Nitrile gas chemistry in Titan’s atmosphere

This work presents the first study of the gaseous products resulting from the partial dissociation of methane and nitrogen in the PAMPRE experimental setup simulating Titan’s atmospheric chemistry.Using cryogenic trapping, the gaseous products generated from the chemical reactions occurring in the reactor have been trapped. Analyses of these products by gas chromatography coupled to mass spectrometry have allowed the detection and identification of more than 30 reaction products. Most of them are identified as nitrile species, accompanied by aliphatic hydrocarbons and a few aromatics compounds. The observed species are in agreement with the data from the recent Cassini-Huygens mission as well as from other laboratory setups capable of dissociating nitrogen and methane. This work emphasizes the probable importance of nitrogen-bearing compounds in the chemistry taking place in Titan’s atmosphere.Furthermore, a quantification of mono-nitriles with saturated alkyl chains has been performed relatively to hydrogen cyanide and shows a power law dependence in their concentration. This dependence is consistent with the Cassini-INMS data and Titan’s photochemical models.An empirical relationship has been extracted from our experimental data: [C_xH_2x−1N] = 100x⁻⁵, where x is the number of carbon atoms in the nitrile molecule. This relationship can be directly used in order to foretell the concentration of heavier nitriles induced by chemistry in Titan’s atmosphere.     

Titan’s atmosphere: An optimal gas mixture for aerosol production?

Here we present the first quantitative study of the gas to solid particle conversion in a Radio Frequency dusty plasma experiment simulating the complex atmospheric reactivity on Titan.Analogs of Titan’s aerosols have been produced in different N₂-CH₄ gas mixtures. Using in situ mass spectrometry, it has been found that, by varying the initial methane concentration, aerosols could be produced in methane steady state concentrations similar to Titan’s atmospheric conditions. In our experiment, an initial ∼5% methane concentration is necessary to ensure a ∼1.5% methane steady state concentration in the plasma.The tholin mass production rate has been quantified as a function of the initial methane concentration. A maximum was found for a steady state CH₄ concentration in agreement with Titan’s atmospheric CH₄ concentrations. At this maximum, the tholin C/N ratio is about 1.45 and the carbon gas to solid conversion yield is about 35%.We have modeled the mass production rate by a parabolic function, highlighting two competitive chemical regimes controlling the tholin production efficiency: an efficient growth process which is proportional to the methane consumption, and an inhibiting process which opposes the growth process and dominates it for initial methane concentrations higher than ∼5%. To explain these two opposite effects, we propose two mechanisms: one involving HCN patterns in the tholins for the growth process, and one involving the increasing amount of atomic hydrogen in the plasma as well as the increase in aliphatic contributions in the tholins for the inhibiting process. This study highlights new routes for understanding the chemical growth of the organic aerosols in Titan’s atmosphere.     

A tropical haze band in Titan’s stratosphere

Inspection of near-infrared images from Cassini’s Imaging Science Subsystem and Visual and Infrared Mapping Spectrometer have revealed a new feature in Titan’s haze structure: a narrow band of increased scattering by haze south of the equator. The band seems to indicate a region of very limited mixing in the lower stratosphere, which causes haze particles to be trapped there. This could explain the sharp separation between the two hemispheres, known as the north-south asymmetry, seen in images. The separation of the two hemispheres can also be seen in the stratosphere above 150 km using infrared spectra measured by Cassini’s Composite Infrared Spectrometer. Titan’s behaviour in the lower tropical stratosphere is remarkably similar to that of the Earth’s tropical stratosphere, which hints at possible common dynamical processes.     

On Titan’s Xanadu region

A large, circular marking ∼1800 km across is seen in near-infrared images of Titan. The feature is centered at 10°S, 120°W on Titan, encompasses much of Titan’s western Xanadu region, and has an off-center, quasi-circular, inner margin about 700 km across, with lobate outer margins extending 200-500 km from the inner margin. On the feature’s southern flank is Tui Regio, an area that has very high reflectivity at 5 μm, and is hypothesized to exhibit geologically recent cryovolcanic flows (Barnes, J.W. et al. [2006]. Geophys. Res. Lett. 33), similar to flows seen in Hotei Regio, a cryovolcanic area whose morphology may be controlled by pre-existing, crustal fractures resulting from an ancient impact (Soderblom, L.A. et al. [2009]. Icarus, 204). The spectral reflectivity of the large, circular feature is quite different than that of its surroundings, making it compositionally distinct, and radar measurements of its topography, brightness temperature and volume scattering also suggest that the feature is quite distinct from its surroundings. These and several other lines of evidence, in addition to the feature’s morphology, suggest that it may occupy the site of an ancient impact.     

The photochemical products of benzene in Titan’s upper atmosphere

The Cassini spacecraft detected benzene high in Titan’s atmosphere as well as the presence of large mass positive and negative ions. Previous work has suggested that these large mass ions could be composed of fused-ring polycyclic aromatic hydrocarbon compounds. These fused-ring PAHs, such as naphthalene and anthracene, are usually the result of high temperature processes that may not occur in Titan’s thin, cold, upper thermosphere. Here we suggest that a different class of aromatic compounds, polyphenyls, may be a better explanation of the data. Polyphenyls can grow to be large polymeric structures and could condense to form the aerosols seen in Titan’s cloud and hazes. They have similar properties to fused-ring PAHs (for example, electron affinity, ionization potential) and could be the negative ion species seen in the CAPS instrument data from the Cassini spacecraft.     

Titan’s 5-μm window: observations with the Very Large Telescope

We report on mid-resolution (R∼2000) spectroscopic observations of Titan, acquired in November 2000 with the Very Large Telescope and covering the range 4.75-5.07 μm. These observations provide a detailed characterization of the CO (1-0) vibrational band, clearly separating for the first time individual CO lines (P10 to P19 lines of ¹³CO). They indicate that the CO/N₂ mixing ratio in Titan’s troposphere is 32±10 ppm. Comparison with photochemical models indicates that CO is not in a steady state in Titan’s atmosphere. The observations confirm that Titan’s 5-μm continuum geometric albedo is ∼0.06, and further indicates a ∼20% albedo decrease over 4.98-5.07 μm. Nonzero flux is detected at the 0.01 geometric albedo level in the saturated core of the ¹²CO (1-0) band, at 4.75-4.85 μm, providing evidence for backscattering on the stratospheric haze. Finally, emission lines are detected at 4.75-4.835 μm, coinciding in position with lines from the CO(1-0) and/or CO(2-1) bands. Matching them by thermal emission would require Titan’s stratosphere to be much warmer (by ∼ 25 K at 0.1 mbar) than indicated by the methane 7.7-μm emission and the Voyager radio-occultation. We show instead that a nonthermal mechanism, namely solar-excited fluorescence, is a more plausible source for these emissions. Improved observations and laboratory measurements on the vibrational-translational relaxation of CO are needed for further interpretation of these emissions in terms of a CO stratospheric mixing ratio.     

Chemical reactions in the Titan’s troposphere during lightning

In the lower troposphere of the Titan the temperature is about 90 K, therefore the chemical production of compounds in the CH₄/N₂ atmosphere is extremely slow. However, atmospheric electricity could provide conditions at which chemical reactions are fast. This paper is based on the assumption that there are lightning discharges in the Titan’s lower atmosphere. The temporal temperature profile of a gas parcel after lightning was calculated at the conditions of 10 km above the Titan’s surface. Using this temperature profile, composition of the after-lightning atmosphere was simulated using a detailed chemical kinetic mechanism consisting of 1829 reactions of 185 species. The main reaction paths leading to the products were investigated. The main products of lighting discharges in the Titan’s atmosphere are H₂, HCN, C₂N₂, C₂H₂, C₂H₄, C₂H₆, NH₃ and H₂CN. The annual production of these compounds was estimated in the Titan’s atmosphere.     

A study of the propagation of electromagnetic waves in Titan’s atmosphere with the TLM numerical method

A numerical modeling of the electromagnetic characteristics of Titan’s atmosphere is carried out by means of the TLM numerical method, with the aim of calculating the Schumann resonant frequencies of Saturn’s satellite. The detection and measurement of these resonances by the Huygens probe, which will enter Titan’s atmosphere at the beginning of 2005, is expected to show the existence of electric activity with lightning discharges in the atmosphere of this satellite. As happens with the Schumann frequencies on Earth, losses associated with electric conductivity will make these frequencies lower than theoretically expected, the fundamental frequency being located between 11 and 15 Hz. This numerical study also shows that the strong losses associated to the high conductivity make it impossible for an electromagnetic wave with a frequency of 10 MHz or lower, generated near the surface, to reach the outer part of Titan’s atmosphere.     

Transient surface liquid in Titan’s polar regions from Cassini

Cassini RADAR images of Titan’s south polar region acquired during southern summer contain lake features which disappear between observations. These features show a tenfold increases in backscatter cross-section between images acquired one year apart, which is inconsistent with common scattering models without invoking temporal variability. The morphologic boundaries are transient, further supporting changes in lake level. These observations are consistent with the exposure of diffusely scattering lakebeds that were previously hidden by an attenuating liquid medium. We use a two-layer model to explain backscatter variations and estimate a drop in liquid depth of approximately 1-m-per-year. On larger scales, we observe shoreline recession between ISS and RADAR images of Ontario Lacus, the largest lake in Titan’s south polar region. The recession, occurring between June 2005 and July 2009, is inversely proportional to slopes estimated from altimetric profiles and the exponential decay of near-shore backscatter, consistent with a uniform reduction of 4 ± 1.3 m in lake depth.Of the potential explanations for observed surface changes, we favor evaporation and infiltration. The disappearance of dark features and the recession of Ontario’s shoreline represents volatile transport in an active methane-based hydrologic cycle. Observed loss rates are compared and shown to be consistent with available global circulation models. To date, no unambiguous changes in lake level have been observed between repeat images in the north polar region, although further investigation is warranted. These observations constrain volatile flux rates in Titan’s hydrologic system and demonstrate that the surface plays an active role in its evolution. Constraining these seasonal changes represents the first step toward our understanding of longer climate cycles that may determine liquid distribution on Titan over orbital time periods.