Determining Titan surface topography from Cassini SAR data

A technique, referred to as SARTopo, has been developed for obtaining surface height estimates with 10 km horizontal resolution and 75 m vertical resolution of the surface of Titan along each Cassini Synthetic Aperture Radar (SAR) swath. We describe the technique and present maps of the co-located data sets. A global map and regional maps of Xanadu and the northern hemisphere hydrocarbon lakes district are included in the results. A strength of the technique is that it provides topographic information co-located with SAR imagery. Having a topographic context vastly improves the interpretability of the SAR imagery and is essential for understanding Titan.SARTopo is capable of estimating surface heights for most of the SAR-imaged surface of Titan. Currently nearly 30% of the surface is within 100 km of a SARTopo height profile. Other competing techniques provide orders of magnitude less coverage.We validate the SARTopo technique through comparison with known geomorphological features such as mountain ranges and craters, and by comparison with co-located nadir altimetry, including a 3000 km strip that had been observed by SAR a month earlier. In this area, the SARTopo and nadir altimetry data sets are co-located tightly (within 5-10 km for one 500 km section), have similar resolution, and as expected agree closely in surface height. Furthermore the region contains prominent high spatial resolution topography, so it provides an excellent test of the resolution and precision of both techniques.     

Ground-based radar observations of Titan: 2000-2008

We have observed Titan with the Arecibo Observatory’s 12.6 cm wavelength radar system during the last eight oppositions of the Saturn system with sufficient sensitivity to characterize its scattering properties as a function of sub-Earth longitude. In a few sessions the Green Bank Telescope was used as the receiving instrument in a bistatic configuration to boost sub-radar track length and integration time. Radar echo spectra have been obtained for a total of 92 viewing geometries with sub-Earth locations scattered through all longitudes and at latitudes between 7.6°S and 26.3°S, close to the maximum southern excursion of the sub-Earth track. We find Titan to have globally average radar albedos at this wavelength of 0.161 in the opposite circular polarization sense as that transmitted (OC) and 0.074 in the same sense (SC), giving a polarization ratio SC/OC of 0.46. These values are intermediate between lower reflectivity rocky surfaces and higher reflectivity clean icy surfaces. The variations with longitude in general mirror the surface brightness variations seen through the infrared atmospheric windows. Xanadu Regio’s radar reflectivity and polarization ratio are higher than the global averages, and suggest that its composition is relatively cleaner water ice or, possibly, some other material with low propagation loss at radio wavelengths. For all echo spectra most of the power is in a broad diffuse component but with a specular component whose strength and narrowness is highly variable as a function of surface location. For all data we fit a sum of the standard Hagfors scattering law describing the specular component and an empirical diffuse radar scattering model to extract bulk parameters of the surface. Many areas exhibit very narrow specular reflections implying terrain that are quite flat on centimeter to meter scales over spans of tens to perhaps hundreds of kilometers. The proportion of spectra showing these narrow specular echoes has fallen significantly over the observational time span, indicating either a latitudinal effect related to terrain differences or changing surface conditions over the past several years. A few radar tracks, especially those from the 2008 session, overlap some high resolution Cassini RADAR imagery swaths to allow a direct comparison with terrain.     

Stereo topographic models of Mercury after three MESSENGER flybys

From photogrammetric analysis of stereo images of Mercury obtained during three MESSENGER flybys, we have produced three digital terrain models (DTMs) that have a grid spacing of 1 km and together cover 30% of the planet's surface. The terrain models provide a rich source of information on the morphology of Mercury's surface, including details of tectonic scarp systems as well as impact craters and basins. More than 400 craters larger than 15 km in diameter are included in the models. Additionally, the models provide important test cases for the analysis of stereo image data to be collected during MESSENGER's orbital mission phase. Small lateral offsets and differences in trends between stereo DTMs and laser altimeter profiles may be due to remaining errors in spacecraft position, instrument pointing, or Mercury coordinate knowledge. Such errors should be resolved during the orbital mission phase, when more joint analyses of data and detailed orbit modeling will be possible.     

Topographic modeling of Phoebe using Cassini images

High-resolution Cassini stereo images of Saturn's moon Phoebe have been used to derive a regional digital terrain model (DTM) and an orthoimage mosaic of the surface. For DTM-control a network of 130 points measured in 14 images (70-390 m/pixel resolution) was established which was simultaneously used to determine the orientation of the spin-axis. The J2000 spin-axis was found at Dec=78.0°±0.1° and RA=356.6°±0.3°, substantially different from the former Voyager solution. The control points yield a mean figure radius of 107.2 km with RMS residuals of 6.2 km demonstrating the irregular shape of this body. The DTM was computed from densely spaced conjugate image points determined by methods of digital image correlation. It has a horizontal resolution of 1-2 km and vertical accuracies in the range 50-100 m. It is limited in coverage, but higher in resolution than the previously derived global shape model of Phoebe [Porco et al., 2005. Cassini imaging science: initial results on Phoebe and Iapetus. Science 307, 1237-1242] and allows us to study the morphology of the surface in more detail. There is evidence for unconsolidated material from a steep and smooth slope at the rim of a 100 km impact feature. There are several conically shaped craters on Phoebe, which may hint at highly porous and low compaction material on the surface.     

The topography of Ganymede's grooved terrain

Using the technique of photoclinometry, topographic profiles across areas of grooved terrain and several other features on Ganymede have been constructed. The grooved terrain examined consists of subparallel grooves spaced 3–10 km apart. Topographic amplitudes are typically 300–400 m, with a maximum of about 700 m. Slopes are very gentle and tend to be primarily concave upward. Very few major positive relief features exist on Ganymede. The most important of these is a broad, gently sloping dome-shaped feature 260 km in diameter and over 2 km high.     

Correlations between VIMS and RADAR data over the surface of Titan: Implications for Titan’s surface properties

We apply a multivariate statistical method to Titan data acquired by different instruments onboard the Cassini spacecraft. We have searched through Cassini/VIMS hyperspectral cubes, selecting those data with convenient viewing geometry and that overlap with Cassini/RADAR scatterometry footprints with a comparable spatial resolution. We look for correlations between the infrared and microwave ranges the two instruments cover. Where found, the normalized backscatter cross-section obtained from the scatterometer measurement, corrected for incidence angle, and the calibrated antenna temperature measured along with the scatterometry echoes, are combined with the infrared reflectances, with estimated errors, to produce an aggregate data set, that we process using a multivariate classification method to identify homogeneous taxonomic units in the multivariate space of the samples.In medium resolution data (from 20 to 100 km/pixel), sampling relatively large portions of the satellite’s surface, we find regional geophysical units matching both the major dark and bright features seen in the optical mosaic. Given the VIMS cubes and RADAR scatterometer passes considered in this work, the largest homogeneous type is associated with the dark equatorial basins, showing similar characteristics as each other on the basis of all the considered parameters.On the other hand, the major bright features seen in these data generally do not show the same characteristics as each other. Xanadu, the largest continental feature, is as bright as the other equatorial bright features, while showing the highest backscattering coefficient of the entire satellite. Tsegihi is very bright at 5 μm but it shows a low backscattering coefficient, so it could have a low roughness on a regional scale and/or a different composition. Another well-defined region, located southwest of Xanadu beyond the Tui Regio, seems to be detached from the surrounding terrains, being bright at 2.69, 2.78 and 5 μm but having a low radar brightness. In this way, other units can be found that show correlations or anti-correlations between the scatterometric response and the spectrophotometric behavior, not evident from the optical remote sensing data.     

Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

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

Radar Backscatter from Underdense Meteors and Diffusion Rates

Many meteoroids burn up between about 120 km and 70 km, deposit metals and dust and form ionized trails which are detected by radars. Model studies about the influence of neutral or positively charged background dust on the ambipolar diffusion indicate that significant smaller decay times should be observed for weak meteor echoes compared to strong meteor echoes which can affect the estimation of temperatures. The variation of meteor decay times in dependence on echo strength, height, and season was studied using radar observations at 69° N, 22° S, and 67° S. Significantly reduced decay times were found for weak echoes below about 88 km at low latitudes throughout the year, and at high latitudes with the exception of summer. In summer at high latitudes, decreasing decay times of weak and strong meteors are observed at altitudes below about 85 km during the appearance of noctilucent clouds. The impact of reduced decay times on the estimation of neutral temperatures from decay times is discussed.     

The geology of Hotei Regio, Titan: Correlation of Cassini VIMS and RADAR

Joint Cassini VIMS and RADAR SAR data of ∼700-km-wide Hotei Regio reveal a rich collection of geological features that correlate between the two sets of images. The degree of correlation is greater than anywhere else seen on Titan. Central to Hotei Regio is a basin filled with cryovolcanic flows that are anomalously bright in VIMS data (in particular at 5 μm) and quite variable in roughness in SAR. The edges of the flows are dark in SAR data and appear to overrun a VIMS-bright substrate. SAR-stereo topography shows the flows to be viscous, 100-200 m thick. On its southern edge the basin is ringed by higher (∼1 km) mountainous terrain. The mountains show mixed texture in SAR data: some regions are extremely rough, exhibit low and spectrally neutral albedo in VIMS data and may be partly coated with darker hydrocarbons. Around the southern margin of Hotei Regio, the SAR image shows several large, dendritic, radar-bright channels that flow down from the mountainous terrain and terminate in dark blue patches, seen in VIMS images, whose infrared color is consistent with enrichment in water ice. The patches are in depressions that we interpret to be filled with fluvial deposits eroded and transported by liquid methane in the channels. In the VIMS images the dark blue patches are encased in a latticework of lighter bands that we suggest to demark a set of circumferential and radial fault systems bounding structural depressions. Conceivably the circular features are tectonic structures that are remnant from an ancient impact structure. We suggest that impact-generated structures may have simply served as zones of weakness; no direct causal connection, such as impact-induced volcanism, is implied. We also speculate that two large dark features lying on the northern margin of Hotei Regio could be calderas. In summary the preservation of such a broad suite of VIMS infrared color variations and the detailed correlation with features in the SAR image and SAR topography evidence a complex set of geological processes (pluvial, fluvial, tectonic, cryovolcanic, impact) that have likely remained active up to very recent geological time (<10⁴ year). That the cryovolcanic flows are excessively bright in the infrared, particularly at 5 μm, might signal ongoing geological activity. One study [Nelson, R.M., and 28 colleagues, 2009. Icarus 199, 429-441] reported significant 2-μm albedo changes in VIMS data for Hotei Arcus acquired between 2004 and 2006, that were interpreted as evidence for such activity. However in our review of that work, we do not agree that such evidence has yet been found.     

A low cost meteor observation system using radio forward scattering and the interferometry technique

We present a low cost meteor observation system based on the radio forward scattering and interferometry technique at Kochi University of Technology (KUT). The system can be a suitable model for low budget educational institutes that target practical learning of astronomical objects and upper atmospheric characteristics. The system methodology for the automatic counting of meteor echoes, filtering noise and detecting meteor echo directions is described. Detection of the meteor echo directions, which is the basic element for determining the meteor trajectories and the orbital parameters of parent comets, is based on a software system developed for analysis of phase differences detected by interferometry. Randomly selected observation samples measured by the radio interferometer are compared to simultaneous optical observations by video cameras to verify the system accuracy. Preliminary error analysis revealed that the system accuracy is directly related to the duration of observed meteor echoes. Eighty percent of meteor echo samples with durations longer than 3 s showed agreement in azimuth and elevation angles measurements to within a 10° error range, while meteor echo samples with shorter durations showed lower agreement levels probably due to the low system sampling resolution of 0.1 s. The reasonable agreement level of meteor echoes with duration longer than 3 s demonstrated the applicability of the system methodology. Accurate observation of shorter duration meteor echoes could possibly be achieved by improving the system resolution.     

Evaluating digital elevation models for glaciologic applications: An example from Nevado Coropuna, Peruvian Andes

This paper evaluates the suitability of readily available elevation data derived from recent sensors – the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and the Shuttle Radar Topography Mission (SRTM) – for glaciological applications. The study area is Nevado Coropuna (6426 m), situated in Cordillera Ampato of Southern Peru. The glaciated area was 82.6 km² in 1962, based on aerial photography. We estimate the glacier area to be ca. 60.8 km² in 2000, based on analysis of the ASTER L1B scene.We used two 1:50,000 topographic maps constructed from 1955 aerial photography to create a digital elevation model with 30 m resolution, which we used as a reference dataset. Of the various interpolation techniques examined, the TOPOGRID algorithm was found to be superior to other techniques, and yielded a DEM with a vertical accuracy of ± 14.7 m. The 1955 DEM was compared to the SRTM DEM (2000) and ASTER DEM (2001) on a cell-by-cell basis. Steps included: validating the DEM's against field GPS survey points on rock areas; visualization techniques such as shaded relief and contour maps; quantifying errors (bias) in each DEM; correlating vertical differences between various DEM's with topographic characteristics (elevation, slope and aspect) and subtracting DEM elevations on a cell-by-cell basis.The RMS error of the SRTM DEM with respect to GPS points on non-glaciated areas was 23 m. The ASTER DEM had a RMS error of 61 m with respect to GPS points and displayed 200–300 m horizontal offsets and elevation ‘spikes’ on the glaciated area when compared to the DEM from topographic data.Cell-by-cell comparison of SRTM and ASTER-derived elevations with topographic data showed ablation at the toes of the glaciers (− 25 m to − 75 m surface lowering) and an apparent thickening at the summits. The mean altitude difference on glaciated area (SRTM minus topographic DEM) was − 5 m, pointing towards a lowering of the glacier surface during the period 1955–2000. Spurious values on the glacier surface in the ASTER DEM affected the analysis and thus prevented us from quantifying the glacier changes based on the ASTER data.     

Analysis of a cryolava flow-like feature on Titan

This paper reports on the analysis of the highest spatial resolution hyperspectral images acquired by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft during its prime mission. A bright area matches a flow-like feature coming out of a caldera-like feature observed in Synthetic Aperture Radar (SAR) data recorded by the Cassini radar experiment [Lopes et al., 2007. Cryovolcanic features on Titan's surface as revealed by the Cassini Titan Radar Mapper. Icarus 186, 395-412, doi:10.1016/j.icarus.2006.09.006]. In this SAR image, the flow extends about 160 km east of the caldera. The contrast in brightness between the flow and the surroundings progressively vanishes, suggesting alteration or evolution of the composition of the cryolava during the lifetime of the eruptions. Dunes seem to cover part of this flow on its eastern end. We analyze the different terrains using the Spectral Mixing Analysis (SMA) approach of the Multiple-Endmember Linear Unmixing Model (MELSUM, Combe et al., 2008). The study area can be fully modeled by using only two types of terrains. Then, the VIMS spectra are compared with laboratory spectra of known materials in the relevant atmospheric windows (from 1 to 2.78 μm). We considered simple molecules that could be produced during cryovolcanic events, including H₂O, CO₂ (using two different grain sizes), CH₄ and NH₃. We find that the mean spectrum of the cryoflow-like feature is not consistent with pure water ice. It can be best fitted by linear combinations of spectra of the candidate materials, showing that its composition is compatible with a mixture of H₂O, CH₄ and CO_2.     

High-resolution lunar radar map at 70-cm wavelength

New radar observations of the Moon in 1981–1984 were made using the 430 MHz (70 cm wavelength) radar at the Arecibo observatory, Puerto Rico. The new observations have produced a high resolution lunar radar map with radar cell-sizes near 2–5 km. This new resolution is a three-fold improvement over the previous mapping done in the late 1960's. Since the Arecibo radar antenna beam is only ten arc-minutes (about one-third of the width of the lunar disk), this new map is a mosaic of some eighteen observations. A radarmetric control between the various pieces of the mosaic was obtained via a beam-swing, limb-to-limb calibration.When the limb-to-limb calibration was combined with the mosaic, there were significant radar scattering differences across the maria. Eastern Mare Tranquillitatis and western Oceanus Procellarum have weaker echoes than other maria, while the central portion of Mare Serenitatis and northern Mare Imbrium have stronger echoes. There is a radar scattering difference across the southern terra as areas nearer Mare Orientale have stronger echoes than areas further from Mare Orientale.     

Radarclinometry of the sand seas of Africa’s Namibia and Saturn’s moon Titan

Radarclinometry is a powerful technique for estimating heights of landforms in synthetic aperture radar (SAR) images of planetary surfaces. In particular, it has been used to estimate heights of dunes in the sand seas of Saturn’s moon Titan (Lorenz, R.D., and 39 colleagues [2006]. Science 312, 724-727). In this work, we verify the technique by comparing dune heights derived from radarclinometry to known topography of dune fields in the Namib sand sea of western Africa. We compared results from three different image grid spacings, and found that 350 m/pixel (the same spacing at which the Cassini RADAR data was processed) is sufficient to determine dune height for dunes of similar morphometry to those of the Namib sand sea. At this grid spacing, height estimates derived from radarclinometry are largely representative of, though may underestimate by as much as 30%, or overestimate by as much as 40%, true dune height. Applying the technique to three regions on Titan, we estimate dune heights of 45-180 m, and dune spacings of 2.3-3.3 km. Obtaining accurate heights of Titan’s dunes will help to constrain the total organic inventory on Titan.     

A close look at Saturn's rings with Cassini VIMS

Soon after the Cassini-Huygens spacecraft entered orbit about Saturn on 1 July 2004, its Visual and Infrared Mapping Spectrometer obtained two continuous spectral scans across the rings, covering the wavelength range 0.35-5.1 μm, at a spatial resolution of 15-25 km. The first scan covers the outer C and inner B rings, while the second covers the Cassini Division and the entire A ring. Comparisons of the VIMS radial reflectance profile at 1.08 μm with similar profiles at a wavelength of 0.45 μm assembled from Voyager images show very little change in ring structure over the intervening 24 years, with the exception of a few features already known to be noncircular. A model for single-scattering by a classical, many-particle-thick slab of material with normal optical depths derived from the Voyager photopolarimeter stellar occultation is found to provide an excellent fit to the observed VIMS reflectance profiles for the C ring and Cassini Division, and an acceptable fit for the inner B ring. The A ring deviates significantly from such a model, consistent with previous suggestions that this region may be closer to a monolayer. An additional complication here is the azimuthally-variable average optical depth associated with “self-gravity wakes” in this region and the fact that much of the A ring may be a mixture of almost opaque wakes and relatively transparent interwake zones. Consistently with previous studies, we find that the near-infrared spectra of all main ring regions are dominated by water ice, with a typical regolith grain radius of 5-20 μm, while the steep decrease in visual reflectance shortward of 0.6 μm is suggestive of an organic contaminant, perhaps tholin-like. Although no materials other than H₂O ice have been identified with any certainty in the VIMS spectra of the rings, significant radial variations are seen in the strength of the water-ice absorption bands. Across the boundary between the C and B rings, over a radial range of ∼7000 km, the near-IR band depths strengthen considerably. A very similar pattern is seen across the outer half of the Cassini Division and into the inner A ring, accompanied by a steepening of the red slope in the visible spectrum shortward of 0.55 μm. We attribute these trends—as well as smaller-scale variations associated with strong density waves in the A ring—to differing grain sizes in the tholin-contaminated icy regolith that covers the surfaces of the decimeter-to-meter sized ring particles. On the largest scale, the spectral variations seen by VIMS suggest that the rings may be divided into two larger ‘ring complexes,’ with similar internal variations in structure, optical depth, particle size, regolith texture and composition. The inner complex comprises the C and B rings, while the outer comprises the Cassini Division and A ring.     

Cloud formation along mountain ridges on Titan

Cassini radar passes have shown a number of mountain ranges on Titan. Radar data covering approximately one quarter of Titan's surface places mountains in primarily equatorial regions with the mean height of about 900 m. The flow of air over topographic features can both trigger and enhance cloud formation. Orographically induced clouds near terrestrial mountain ranges include shallow wave clouds produced from upslope flow as well as precipitating stratus and cumulus type clouds; mountains can provide the perturbations needed to trigger convective clouds. The Titan regional atmospheric modeling system (TRAMS) has been used to explore a number of convective cloud properties and is now used to report on clouds formed when a mountain peak is placed within the model domain. Using a range of heights and surface winds compatible with Cassini/Huygens data, constraints can be placed on the scenarios in which clouds can be expected to form. Given sufficiently humid conditions (at least 50% humidity), convection is triggered. For drier environments similar to the Huygens landing site, short-lived, optically thin clouds form from air rising upslope. Precipitation is also seen in the cases of the convective clouds, which could have implications for the eroded appearance of Titan's mountains.     

A Proposed lunar orbiting gravity gradiometer experiment

Analysis of the gravity gradiometer developed by R. L. Forward and C. C. Bell at the Hughes Research Laboratories suggest than an accuracy, in the range 0.1 to 0.5 EU can be expected in a lunar orbiter application. This accuracy will allow gradient anomalies associated with mascons to be mapped with 1% accuracy and should reveal a great deal of new information about the lunar gravity field.The proposed experiment calls for putting such a gradiometer into a closely circular polar orbit at an average height of about 30 km above the lunar surface. This orbit allows the entire lunar surface to be covered in fourteen days, the gradiometer to be checked twice per revolution and results in successive passes above the lunar surface being spaced at about the resolution limit of about 30 km set both by the satellite altitude and instrumental integration time. Doppler tracking will be employed and the spacecraft will carry an electromagnetic altimeter. Gradient and altitude data from the far side of the Moon can be stored for replay when communication is re-established.     

A harmonic analysis of lunar topography

A global lunar topographic map has been derived from existing Earth-based and orbital observations supplemented in areas without data by a linear autocovariance predictor. Of 2592 bins, each 5° square, 1380 (64.7% by area) contain at least one measurement. A spherical harmonic analysis to degree 12 yields a mean radius of 1737.53 ± 0.03 km (formal standard error) and an offset of the center of figure of 1.98±0.06 km toward (19±2)°S, (194±1)°E. A Bouguer gravity map, derived from a 12-degree free-air gravity model and the present topography data, is presented for an elevation of 100 km above the mean surface. It is confirmed that the low-degree gravity harmonics are determined primarily by surface height variations and only secondarily by lateral density variations.     

Multiple radio echoes in the solar corona

Echo-type solar radio bursts are associated with preceding short-lived bursts in double events. The peculiar and rather rare decameter echoes are observed with a UTR-2 antenna. The initial narrow-band burst is followed, some 5 to 10 s later, by an echo-like burst at the same frequency. The observational data obtained for decameter echo evens are, on the whole, consistent with the model of a pulsed source emitting radio signals at the plasma-frequency harmonic, which is placed in a non-uniform corona and rotates together with the Sun.Intensity-time profiles of 25 MHz echo bursts are of an unusual shape, featuring an extended leading edge and an abrupt decay at the trailing edge and also showing some fine structures in the form of an additional, weakly pronounced maximum or a step at the final stage of the burst. Time parameters characterizing the profiles are evaluated. The step is delayed with respect to the main pulse at about two times longer than the principal echo maximum. At the same time, the time delays depend essentially on the heliolongitude of the active region and achieve their maximum values at the meridian. The step height does not exceed 0.5 of the echo maximum. At this level, the echo-decay time almost coincides with the initial burst duration but is about 1.7 times less than the echo-rise time. The feature at the echo tail can be interpreted as a result of a repeated reflection of the burst from the source region. The causes and conditions for the formation of multiple echoes are discussed. The extended leading edge of the echo permits us to assume a quasi-radial fibrous structure of the corona, capable of back-scattering the incident radiation.     

Cryovolcanic features on Titan's surface as revealed by the Cassini Titan Radar Mapper

The Cassini Titan Radar Mapper obtained Synthetic Aperture Radar images of Titan's surface during four fly-bys during the mission's first year. These images show that Titan's surface is very complex geologically, showing evidence of major planetary geologic processes, including cryovolcanism. This paper discusses the variety of cryovolcanic features identified from SAR images, their possible origin, and their geologic context. The features which we identify as cryovolcanic in origin include a large (180 km diameter) volcanic construct (dome or shield), several extensive flows, and three calderas which appear to be the source of flows. The composition of the cryomagma on Titan is still unknown, but constraints on rheological properties can be estimated using flow thickness. Rheological properties of one flow were estimated and appear inconsistent with ammonia-water slurries, and possibly more consistent with ammonia-water-methanol slurries. The extent of cryovolcanism on Titan is still not known, as only a small fraction of the surface has been imaged at sufficient resolution. Energetic considerations suggest that cryovolcanism may have been a dominant process in the resurfacing of Titan.