Cross-sectional profiles of Baltis Vallis channel on Venus: Reconstructions from Magellan SAR brightness data

Baltis Vallis is a 6800-km long canali-type channel on Venus. Canali have a unique combination of morphological characteristics: extraordinary length, a single main conduit, and a degree of similarity to terrestrial rivers. These characteristics have given rise to intensive discussions on whether the origin of canali is erosional or constructional. Cross-sectional profiles of such channels reveal the detailed morphology of the structure and enable us to distinguish between these two possible origins; however, canali are just several kilometers wide and are therefore too small for the construction of cross-sectional profiles from Magellan altimetry data. Instead, we propose a new method for reconstructing short-wavelength topography using brightness data from Synthetic Aperture Radar images. We apply Muhleman's backscattering function to the backscatter intensity calculated from the brightness of Magellan Full-Resolution SAR Map images. The estimated vertical error of this new method is less than 5 m for a distance of 1 km across the channel. We studied 120 sites along an approximately 6000 km extent of Baltis Vallis. The channel profiles reveal that in nearly 90% of these sites, the bottom surface of the channel is lower than the surrounding plains by 20-100 m. Clear levee structures and intra-channel ridges are recognized in about 30 and 25%, respectively, of the sites analyzed within Baltis Vallis. Most of the levees occur in the upper segment of Baltis Vallis, while intra-channel ridges are mostly confined to the region between 1500 and 3000 km downstream from the probable source. The average depth and width of the channel are 46 m (standard deviation: 16 m) and 2.2 km (standard deviation: 0.4 km), respectively, and the depth profile along the channel is highly undulatory. The groove-like morphology and paucity of levee structures indicates an erosional origin. Furthermore, the observed undulations in depth along the channel indicate that Baltis Vallis most likely formed by mechanical erosion. The observed morphological transition from levees to intra-channel ridges suggests that the channel-forming processes changed across an area located approximately 1500 km from the source. Carbonatite is the most likely candidate material for the low-viscosity fluid that formed Baltis Vallis.     

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.     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

Depth profiles of venusian sinuous rilles and valley networks

More than 200 venusian channels and valleys have been mapped based on analyses of Magellan SAR images. Sinuous rilles are the most abundant channels among six types of venusian channels, and they are widely distributed on Venus. Morphological characteristics of venusian sinuous rilles include sinuous narrowing reaches, source depressions, and length of several 10s to a few 100s of km. This type of channels is known to exist on the Moon and possibly on Mars. Valley networks on Venus often occur in the vicinity of or in connection to sinuous rilles. Cross-sectional morphologies of sinuous rilles and valley networks are of special importance in discussing their formation processes both qualitatively and quantitatively. We reconstructed cross-sectional profiles of 6 sinuous rilles and 2 valley networks using a new radar clinometric method. Reconstructed cross-sections revealed that floors of the channels and valleys are clearly lower than the surrounding plains. This finding implies that the sinuous rilles and the valley networks have erosional origins. Longitudinal depth profiles of the sinuous rilles show distinct decreasing trends toward the termini. Such decreasing trends of depths are qualitatively in agreement with theoretical models and laboratory experiments of thermal erosion. In order to verify this assertion quantitatively, we conduct simple 1-dimensional model calculations under the assumption that both channel-forming lavas and ground substrate are tholeiitic basalt. For initial lava thicknesses in the range 2-6 m, the model calculations yield good matches to the depth profiles. Estimated duration of lava effusion ranges from several months to a few years. These numerical results support thermal erosion of the sinuous rilles but do not necessarily exclude contributions from mechanical erosion processes. Valley networks seem to have formed under a strong structural control in comparison to sinuous rilles. The valleys vary widely in characteristics of the depth profile and flow directions relative to surface slopes. Therefore valley networks appear to have originated from diverse formation mechanisms.     

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.     

The frequency-area distribution of volcanic units on Venus: Implications for planetary resurfacing

The areas of volcanic units on Venus have been measured on the 1:5000000 geological maps published by NASA/USGS. These data were used to obtain a frequency-area distribution. The cumulative frequency-area distribution of 1544 specific occurrence of units cover six orders of magnitude from the largest unit (30 × 10⁶ km²) to the smallest (20 km²). The probability distribution function has been calculated. The medium and large volcanic units correlate well with a power-law (fractal) relation for the dependence of frequency on area with a slope of −1.83. There are fewer small units than the expected values provided by the power-law relation. Our measurements cover 21.02% of the planetary surface, 3.59% of the study area was found to be tessera terrain and is excluded from this study of volcanism. The measurements were restricted to areas where geological maps have been published. The analysis was performed on two independent areas of the planet, with a complete coverage of published maps. In both areas the largest volcanic unit covers a significant portion of the surface (58.75% and 63.64%, respectively). For the total measured volcanic units (excluding tessera), these two largest units (that could correspond to the same unit or not) cover the 61.18% and they are stratigraphically superimposed on older volcanic units which cover 3.37% of the area. The remaining area (35.45%) is occupied by younger volcanic units stratigraphically superimposed on the large volcanic unit(s). These results are based on the independent mapping of a large number of geologists with different ideas about the geodynamical evolution of Venus and different criteria for geological mapping. Despite this fact, the presence of these very large units is incompatible with the equilibrium resurfacing models, because their generation at different ages would destroy the crater randomness. Our frequency-area distribution of the mapped volcanic units supports a catastrophic resurfacing due to the emplacement of the largest unit(s) followed by a decay of volcanism. Our data for the frequency-area distribution of volcanic units provide new support for catastrophic resurfacing models. It is difficult to make our observations compatible with equilibrium, steady-state resurfacing models.     

Constraints on Titan’s topography through fractal analysis of shorelines

Titan’s north polar hydrocarbon lakes offer a unique opportunity to indirectly characterize the statistical properties of Titan’s landscape. The complexity of a shoreline can be related to the complexity of the landscape it is embedded in through fractal theory. We mapped the shorelines of 290 of the north polar titanian lakes in the Cassini synthetic aperture radar dataset. Out of these, we used a subset of 190 lake shorelines for our analysis. The fractal dimensions of the shorelines were calculated via two methods: the divider/ruler method and the box-counting method, at length scales of (1-10) km and found to average 1.27 and 1.32, respectively. The inferred power-spectral exponent of Titan’s topography (β) from theoretical and empirical relations is found to be ?2, which is lower than the values obtained from the global topography of the Earth or Venus. Some of the shorelines exhibit multi-fractal behavior (different fractal dimensions at different scales), which we interpret to signify a transition from one set of dominant surface processes to another. We did not observe any spatial variation in the fractal dimension with latitude; however we do report significant spatial variation of the fractal dimension with longitude. A systematic difference between the dimensions of orthogonal sections of lake shorelines is also noted, which signifies possible anisotropy in Titan’s topography. The topographic information thus gleaned can be used to constrain landscape evolution modeling to infer the dominant surface processes that sculpt the landscape of Titan.     

Generation and emplacement of fine-grained ejecta in planetary impacts

We report here on a survey of distal fine-grained ejecta deposits on the Moon, Mars, and Venus. On all three planets, fine-grained ejecta form circular haloes that extend beyond the continuous ejecta and other types of distal deposits such as run-out lobes or ramparts. Using Earth-based radar images, we find that lunar fine-grained ejecta haloes represent meters-thick deposits with abrupt margins, and are depleted in rocks ?1 cm in diameter. Martian haloes show low nighttime thermal IR temperatures and thermal inertia, indicating the presence of fine particles estimated to range from ∼10 μm to 10 mm. Using the large sample sizes afforded by global datasets for Venus and Mars, and a complete nearside radar map for the Moon, we establish statistically robust scaling relationships between crater radius R and fine-grained ejecta run-out r^* for all three planets. On the Moon, r^* ∼ R^−0.18 for craters 5-640 km in diameter. For Venus, radar-dark haloes are larger than those on the Moon, but scale as r^* ∼ R^−0.49, consistent with ejecta entrainment in Venus’ dense atmosphere. On Mars, fine-ejecta haloes are larger than lunar haloes for a given crater size, indicating entrainment of ejecta by the atmosphere or vaporized subsurface volatiles, but scale as R^−0.13, similar to the ballistic lunar scaling. Ejecta suspension in vortices generated by passage of the ejecta curtain is predicted to result in ejecta run-out that scales with crater size as R^1/2, and the wind speeds so generated may be insufficient to transport particles at the larger end of the calculated range. The observed scaling and morphology of the low-temperature haloes leads us rather to favor winds generated by early-stage vapor plume expansion as the emplacement mechanism for low-temperature halo materials.     

Radar imagery of the southern Caloris region, Mercury

Radar imagery from July 2005 Arecibo observations has provided new information on surface relief over the southern portion of Caloris Basin and the smooth plains to the south of the basin. A lobe of smooth plains has been identified in the Mariner-unimaged region southwest of Mozart Crater that coincides precisely with topographically down-bowed terrain seen in earlier Arecibo radar altimetry. A 105-km-diameter crater has been found at 193.6° W, 25.6° N that appears to be the largest crater in the Caloris basin floor.     

Shallow radar (SHARAD) sounding observations of the Medusae Fossae Formation, Mars

The SHARAD (shallow radar) sounding radar on the Mars Reconnaissance Orbiter detects subsurface reflections in the eastern and western parts of the Medusae Fossae Formation (MFF). The radar waves penetrate up to 580 m of the MFF and detect clear subsurface interfaces in two locations: west MFF between 150 and 155° E and east MFF between 209 and 213° E. Analysis of SHARAD radargrams suggests that the real part of the permittivity is ∼3.0, which falls within the range of permittivity values inferred from MARSIS data for thicker parts of the MFF. The SHARAD data cannot uniquely determine the composition of the MFF material, but the low permittivity implies that the upper few hundred meters of the MFF material has a high porosity. One possibility is that the MFF is comprised of low-density welded or interlocked pyroclastic deposits that are capable of sustaining the steep-sided yardangs and ridges seen in imagery. The SHARAD surface echo power across the MFF is low relative to typical martian plains, and completely disappears in parts of the east MFF that correspond to the radar-dark Stealth region. These areas are extremely rough at centimeter to meter scales, and the lack of echo power is most likely due to a combination of surface roughness and a low near-surface permittivity that reduces the echo strength from any locally flat regions. There is also no radar evidence for internal layering in any of the SHARAD data for the MFF, despite the fact that tens-of-meters scale layering is apparent in infrared and visible wavelength images of nearby areas. These interfaces may not be detected in SHARAD data if their permittivity contrasts are low, or if the layers are discontinuous. The lack of closely spaced internal radar reflectors suggests that the MFF is not an equatorial analog to the current martian polar deposits, which show clear evidence of multiple internal layers in SHARAD data.     

Regional fracture patterns around volcanoes: Possible evidence for volcanic spreading on Venus

Magellan data show that the surface of Venus is dominated by volcanic landforms including large flow fields and a wide range of volcanic edifices that occur in different magmatic and tectonic environments. This study presents the results from a comprehensive survey of volcano-rift interaction in the BAT region and its surroundings. We carried out structural mapping of examples where interaction between volcanoes and regional fractures results in a deflection of the fractures around the volcanic features and discuss the nature of the local volcano-related stress fields that might be responsible for the observed variations of the regional fracture systems. We propose that the deflection of the regional fractures around these venusian volcanoes might be related to volcanic spreading, a process recognized as of great importance in the tectonic evolution of volcanoes on Earth and Mars, but not previously described on Venus.     

Mercury: Radar images of the equatorial and midlatitude zones

Radar imaging results for Mercury's non-polar regions are presented. The dual-polarization, delay-Doppler images were obtained from several years of observations with the upgraded Arecibo S-band (λ12.6-cm) radar telescope. The images are dominated by radar-bright features associated with fresh impact craters. As was found from earlier Goldstone-VLA and pre-upgrade Arecibo imaging, three of the most prominent crater features are located in the Mariner-unimaged hemisphere. These are: “A,” an 85-km-diameter crater (348° W, 34° S) whose radar ray system may be the most spectacular in the Solar System; “B,” a 95-km-diameter crater (343° W, 58° N) with a very bright halo but less distinct ray system; and “C,” an irregular feature with bright ejecta and rays distributed asymmetrically about a 125-km source crater (246° W, 11° N). Due south of “C” lies a “ghost” feature (242° W, 27° S) that resembles “A” but is much fainter. An even fainter such feature is associated with Bartok Crater. These may be two of the best mercurian examples of large ejecta/ray systems observed in an intermediate state of degradation. Virtually all of the bright rayed craters in the Mariner 10 images show radar rays and/or bright rim rings, with radar rays being less common than optical rays. Radar-bright craters are particularly common in the H-7 quadrangle. Some diffuse radar albedo variations are seen that have no obvious association with impact ejecta. In particular, some smooth plains regions such as the circum-Caloris plains in Tir, Budh, and Sobkou Planitiae and the interiors of Tolstoj and “Skinakas” basins show high depolarized brightness relative to their surroundings, which is the reverse of the mare/highlands contrast seen in lunar radar images. Caloris Basin, on the other hand, appears dark and featureless in the images.     

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.     

Cassini RADAR constraint on Titan's winter polar precipitation

The Cassini RADAR instrument made a dedicated cloud backscatter observation near Titan's north pole, presently in winter darkness, to constrain the precipitation of material onto the surface. The detection limit is ∼5 orders of magnitude above that expected in methane rainstorms, and rules out ‘drizzle’ of more than , placing constraints on the winter accumulation of material on Titan's surface during polar winter.     

Serpentinization and the inorganic synthesis of H2 in planetary surfaces

The near-surface inorganic synthesis of molecular hydrogen (H₂) is a fundamental process relevant to the origins and to the sustenance of early life on Earth and potentially other planets. Hydrogen production through the decomposition of water is thought to be a principal reaction that occurs during hydrothermal alteration of olivine, an iron-magnesium silicate abundant near planetary surfaces. We demonstrate that copious amounts of H₂ are produced only when the olivine undergoing alteration (serpentinization) contains 1 to 50 mol% iron over a variety of planetary surface P-T conditions. This suggests that extrasolar Earth-like planets that are hosted by a star with iron contents up to two times the solar value could support life provided they are hydrothermally active and fall within the habitable zone around the star.     

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.     

A method to infer past surface mass balance and topography from internal layers in martian polar layered deposits

Internal layers in ice masses can be detected with ice-penetrating radar. In a flowing ice mass, each horizon represents a past surface that has been subsequently buried by accumulation, and strained by ice flow. These layers retain information about relative spatial patterns of accumulation and ablation (mass balance). Internal layers are necessary to accurately infer mass-balance patterns because the ice-surface shape only weakly reflects spatial variations in mass balance. Additional rate-controlling information, such as the layer age, the ice temperature, or the ice-grain sizes and ice-crystal fabric, can be used to infer the absolute rate of mass balance. To infer mass balance from the shapes of internal layers, we solve an inverse problem. The solution to the inverse problem is the best set or sets of unknown boundary conditions or initial conditions that, when used in our calculation of ice-surface elevation and internal-layer shape, generate appropriate predictions of observations that are available. We also show that internal layers can be used to infer martian paleo-surface topography from a past era of ice flow, even though the topography may have been largely altered by subsequent erosion. We have successfully inferred accumulation rates and surface topography from internal layers in Antarctica. Using synthetic data, we demonstrate the ability of this method to solve the corresponding inverse problem to infer accumulation and ablation rates, as well as the surface topography, for martian ice. If past ice flow has affected the shapes of martian internal layers, this method is necessary to infer the spatial pattern and rate of mass balance.     

Anomalous radar backscatter from Titan’s surface?

Since Cassini arrived at Saturn in 2004, its moon Titan has been thoroughly mapped by the RADAR instrument at 2-cm wavelength, in both active and passive modes. Some regions on Titan, including Xanadu and various bright hummocky bright terrains, contain surfaces that are among the most radar-bright encountered in the Solar System. This high brightness has been generally attributed to volume scattering processes in the inhomogeneous, low-loss medium expected for a cold, icy satellite surface. We can test this assumption now that the emissivity has been obtained from the concurrent radiometric measurements for nearly all the surface, with unprecedented accuracy (Janssen et al., and the Cassini RADAR Team [2009]. Icarus 200, 222-239). Kirchhoff’s law of thermal radiation relates the radar and radiometric properties in a way that has never been fully exploited. In this paper we examine here how this law may be applied in this case to better understand the nature of Titan’s radar-bright regions. We develop a quantitative model that, when compared to the observational data, allows us to conclude that either the reflective characteristics of the putative volume scattering subsurface must be highly constrained, or, more likely, organized structure on or in the surface is present that enhances the backscatter.     

Radar-bright channels on Titan

During Cassini’s T44 flyby of Titan (May 28, 2008), the Cassini SAR (synthetic aperture radar) revealed sinuous channels in the Southwest of Xanadu. These channels feature very large radar cross-sections, up to 5 dB, whereas the angle of incidence was relatively high, ∼20°. This backscatter is larger than allowed by the coherent backscatter model considered to explain the unusual reflective and polarization properties of the icy satellites and only a few radar scattering mechanisms can be responsible for such high radar returns. The presence of rounded (icy) pebbles with size larger than the radar wavelength (2.18 cm) is proposed to explain the large radar cross-sections measured in these units. The radar-bright channels are thus interpreted as riverbeds, where debris, likely shaped and transported by fluvial activity, have been deposited. Similar debris were observed in the landing site of the Huygens probe. This work may point the way to an explanation for the enhanced brightness of other fluvial regions of Titan.