Surface temperatures and retention of H2O frost on Ganymede and Callisto

Surface temperatures and ice evaporation rates are calculated for Ganymede and Callisto as a function of latitude, time of day, and albedo. The model uses surface thermal properties determined by eclipse radiometry (Morrison and Cruikshank, 1973Icarus18 224–236) and albedos determined from photometrically decalibrated Voyager images. Daytime temperatures on Callisto are roughly 8°K warmer than those in Ganymede's cratered terrain and 11°K warmer than those in Ganymede's grooved terrain. Diurnal mean ice evaporation rates are high enough on both bodies that the surface material probably consists of a very low density lag deposit of primarily silicate dust overlying a denser regolith of silicates and ice. The difference in temperature between Ganymede and Callisto is not great enough to account for the lack of bright polar caps on Callisto. This lack seems instead to reflect a real deficiency in the amount of available H₂O frost relative to Ganymede. The temperature difference between Ganymede's grooved and cratered terrains also cannot account for the strong concentration of bright ray craters in grooved terrain. This concentration suggests instead that an internal geologic process has enriched the grooved terrain in ice relative to the cratered terrain.     

Voyager photometry of surface features on Ganymede and Callisto

The photometric properties of selected surface features on Ganymede and Callisto have been studied using Voyager images over phase angles from 10 to 124° taken with the clear filter (effective wave wavelength ∽0.5 μm). Normal reflectences on Ganymede average 0.35 for the cratered terrain and 0.44 for the grooved terrain. The value for the ubiquitous cratered terrain on Callistro is 0.18. The photometric properties of these regions are described closely by a simple scattering function of the form I = Af(α)μ₀/(μ + μ₀), where A is a constant, μ is the cosine of the emission angle, μ₀ is the cosine of the incidence angle, and f(α) is a function of the phase angle, α, only. For these terrains the shape of f(α) is qualitatively similar to that for the moon—generally concave upward. By contrast, bright craters on both satellites have f(α)'s which are concave downward. The scattering properties of these bright features are definitely not Lambertian, but are described approximately by the scattering law given above. The brightest craters on Callisto have reflectances which are only 10% lower than the brightest craters on Ganymede; both have closely similar scattering laws. We estimate that the brightest craters on Ganymede may reach normal reflectances of 0.7. Our phase functions yield phase integrals of q = 0.8 and 0.6 for Ganymede and Callisto, respectively.     

Early thermal profiles and lithospheric strength of Ganymede from extensional tectonic features

Estimates of the brittle lithosphere thickness derived from the width and spacing of extensional tectonic features, coupled with lithospheric strength envelopes (brittle and ductile yield stress versus depth) appropriate for ice, allow the quantitative determination of early thermal profiles and lithospheric strength and stability on Ganymede. Furrows and grooves indicate brittle lithospheric thicknesses of 5–10 and 2–5 km, respectively, assuming that their spacing is controlled by an extensional instability or that their width is controlled by the intersection depth of their bounding faults. Plots of the brittle and ductile yield stress versus depth for the icy lithosphere of Ganymede show a linear increase in brittle strength with depth to a maximum at the brittle-ductile transition, followed by an exponential decrease in ductile yield stress with depth. Because the depth to the brittle-ductile transition depends primarily on the thermal gradient, the thickness of the brittle lithosphere can be used to calculate early thermal profiles of 1.5–6 and 4.5–20°/km during the formation of the furrows and grooves, respectively. Lithospheric strength, the integral of the yield stress versus depth curve, varied from 30–125 GPa m when the furrows formed to 5–30 GPa m when the grooves formed, which correspond to maximum yield stresses of 6–11 and 2.5–6 MPa, respectively. These results indicate that the thermal gradient and lithospheric strength varied laterally by as much as a factor of 5 and that Ganymede cooled in a highly inhomogeneous manner with significant lateral thermal anomalies. Finally, this analysis provides a straightforward explanation for the stability of large remnants of cratered terrain such as Galileo Regio that had a low thermal gradient and strong lithosphere in contrast to small remnants of cratered terrain that were fractured and broken up by grooved terrain as a result of higher thermal gradients and weaker lithospheres.     

A geometric analysis of surface deformation: Implications for the tectonic evolution of Ganymede

A visual nonalignment of the furrows and the circularity of impact craters are used to study surface deformation on Ganymede. The furrow system is examined to test the hypothesis that lateral motion has taken place between areas of dark terrain. Results show that while lateral motion cannot be ruled out, it may not be required to explain the geometry of the system. Initial nonconcentricity of the furrows or an early period of penetrative deformation shortly after furrow formation could also account for the present configuration. Centers of curvature of the furrows in Galileo and Marius Regiones are numerically determined and it is shown that if lateral movement did occur, it is not possible to determine the amount of displacement. The axial ratios of impact craters in the Uruk Sulcus region which separates Galileo and Marius Regiones are determined and show that large scale shear deformation has not occured in that area since bright terrain was emplaced. Deformation of impact craters within Galileo Regio suggests that Ganymede's lithosphere has behaved rigidly throughout most of the satellite's evolution. The shapes and orientations of impact craters in dark terrain around wedges of bright terrain are used to place an upper limit on the amount of extension associated with bright terrain formation.     

Viscosity of the lithosphere of Enceladus

High resolution Voyager II images of Enceladus reveal that some regions on its surface are highly cratered; the most heavily cratered surfaces probably date back to a period of heavy bombardment. The forms of many of the craters on Enceladus are similar to those of fresh lunar craters, but many of the craters are much shallower in depth, and the floors of some craters are bowed up. The flattering of craters and bowing up of the floors are indicative of viscous relaxation of the topography. Analysis of the forms of the flattened craters suggests that the viscosity at the top of the lithosphere, in the most heavily cratered regions, is between 10²⁴ and 10²⁵ P. The exact time scale for the collapse of the craters is not known, but probably was between 100 my and 4 gy. The flattened craters are located in distinct zones that are adjacent to zones, of similar age, where craters have not flattened. The zones where flattened craters occur possibly are regions in which the heat flow was (or is) higher than that in the adjacent terrains. Because the temperature at the top of the lithosphere of Enceladus would be less than or equal to that of Ganymede and Callisto, if it is covered by a thick regolith, and because the required viscosity on Enceladus is one to two orders of magnitude less than that for Ganymede and Callisto, it can be concluded that the lithospheric material on Enceladus is different from that of Ganymede and Callisto. Enceladus probably has a mixture of ammonia ice and water ice in the lithosphere, whereas the lithospheres of Ganymede and Callisto are composed primarily of water ice.     

The effect of target properties on crater morphology: Comparison of central peak craters on the Moon and Ganymede

Abstract— We examine the morphology of central peak craters on the Moon and Ganymede in order to investigate differences in the near‐surface properties of these bodies. We have extracted topographic profiles across craters on Ganymede using Galileo images, and use these data to compile scaling trends. Comparisons between lunar and Ganymede craters show that crater depth, wall slope and amount of central uplift are all affected by material properties. We observe no major differences between similar‐sized craters in the dark and bright terrain of Ganymede, suggesting that dark terrain does not contain enough silicate material to significantly increase the strength of the surface ice. Below crater diameters of ?12 km, central peak craters on Ganymede and simple craters on the Moon have similar rim heights, indicating comparable amounts of rim collapse. This suggests that the formation of central peaks at smaller crater diameters on Ganymede than the Moon is dominated by enhanced central floor uplift rather than rim collapse. Crater wall slope trends are similar on the Moon and Ganymede, indicating that there is a similar trend in material weakening with increasing crater size, and possibly that the mechanism of weakening during impact is analogous in icy and rocky targets. We have run a suite of numerical models to simulate the formation of central peak craters on Ganymede and the Moon. Our modeling shows that the same styles of strength model can be applied to ice and rock, and that the strength model parameters do not differ significantly between materials.     

Topographic domes on Ganymede: Ice vulcanism or isostatic upwarping

Voyager images of Ganymede show two broad, gently sloping dome-shaped features. They lie in grooved terrain and have diameters of roughly 250 km. The one observed at high resolution has a summit elevation 2–2.5 km above the surrounding plains, and appears to be surrounded by a field of secondary craters. Two formation processes are considered: water vulcanism triggered by a major impact, and isostatic upwarping of a crater formed in a thin crust. Numerical simulation of nonadiabatic water vulcanism indicates that the volume of the domes is inconsistent with eruption through a conduit created by complete penetration of the crust by an impact. It is consistent, however, with eruption through fractures created by an impact that excavates partly through a thin crust. Isostatic upwarp rates calculated as a function of effective crustal temperature indicate that upwarping could also create such a dome if the impact excavated to depths where the crust was sufficiently warm and mobile. Both models require that the density of the crust slightly exceed that of a liquid water mantle for a short period of time. Morphologic evidence suggests that both processes may have been important. If either of the proposed models is correct, the situation of the domes in grooved terrain implies that grooved terrain formation occurred prior to the thickening and stiffening of Ganymede's crust.     

The geology of Dione

Dione is one of the more geologically complex of Saturn's satellites. Several geologic units have been identified including ancient heavily cratered terrain; two plains units: cratered plains and lightly cratered plains; lobate deposits; crater rim deposits; and bright wispy material. The only structural features observed are a series of troughs which cross portions of the surface and subtle northeast and northwest trending lineaments. The troughs are associated with volcanic deposits and are interpreted to be the vents through which material was erupted. Correlations exist between telescopically observed albedo patterns and the distribution of geologic units.     

Distribution of small channels on the Martian surface

The distribution of small channels on Mars has been mapped from Mariner 9 images, at the 1:5 000 000 scale, by the author. The small channels referred to here are small valleys ranging in width from the resolution limit of the Mariner 9 wide-angle images (～1 km) to about 10 km. The greatest density of small band occurs in dark cratered terrain. This dark zone forms a broad subequatorial band around the planet. The observed distribution may be the result of decreased small-channel visibility in bright areas due to obscuration by a high albedo dust or sediment mantle. Crater densities within two small-channel segments show crater size-frequency distributions consistent with those of the oldest of the heavily cratered plains units. Such crater densities coupled with the almost exclusive occurrence of small channels in old cratered terrain and the generally degraded appearance of small channels in the high-resolution images (～100 m) imply a major episode of small-channel formation early in Martian geologic history.     

Thermal segregation of water ice on the Galilean satellites

Consideration of the thermal sublimation of ice on the Galilean satellites suggests that dirty-ice surfaces are susceptible to a process of cold-trapping of water in local bright patches and its preferential removal from dark areas. The result may be very rapid (decade time scale) segregation on the surface into bright icy regions and regions covered by dark ice-free lag deposits. Ion sputtering and micrometeorite bombardment are probably insufficient to prevent this process at low latitudes on Ganymede and Callisto. Sputtering on Europa may prevent segregation, especially on the trailing side. Segregated regions must be mostly smaller than the kilometer resolution of the Voyager images, but larger than centimeter size.     

2.7- to 4.1-μm spectrophotometry of icy satellites of Saturn and Jupiter

We present spectrophotometry in the 27–41 μm spectral region for icy satellites of Saturn (Tethys, Dione, Rhea, Iapetus, and Hyperion) and Jupiter (Europa, Ganymede, and Callisto). The 3.6-μm reflectance peak characteristic of fine-grained water ice is observed prominently on the satellites of Saturn, faintly on the leading side of Europa, and not all on Ganymede, Callisto, or the dark side of Iapetus. The spectral reflectances of these icy satellites may be affected by their equilibrium surface temperatures and magnetospheric effects.     

Surface evolution of two-component stone/ice bodies in the Jupiter region

Observational and theoretical data converge on the conclusion that planetesimals in Jupiter's region of the solar nebula were initially composed predominantly of a mixture of roughly 39–70% H₂O ice by volume, and 30–61% dark stony material resembling carbonaceous chondrites. Recent observations emphasize a division of most asteroid and satellite surfaces in this region into two distinct groups: bright icy material and dark stony material. The present model accounts for these by two main processes: an impact-induced buildup of a dark stony regolith in the absence of surface thermal disturbance, and thermal-disturbance-induced eruption of “water magmas” that create icy surfaces. “Thermal disturbances” include tidal and radiative effects caused by nearness of a planet. A correlation of crater density and albedo, Ganymede's dark-ray craters, and other observed phenomena (listed in the summary) appear consistent with the model discussed here.     

Ring geometry on Ganymede and Callisto

Arguments have been made, based on geometry, for both an impact and an internal origin for the ancient, partially preserved, major furrow system of Ganymede. Zuber and Parmentier concluded that furrows are not concentric but could be impact related if multiringed structures on icy satellites are initially noncircular. We examine the geometry of the Valhalla ring structure on Callisto in order to assess the circularity of an unmodified ring system. Despite prominent local meandering, the only gross deviations from concentricity in the Valhalla system are found in the outer northeast quadrant of the system. Here, a number of ring segments intersect small circles about the center at angles up to 30°. The Ganymede furrow system was remapped to make use of improvements in coordinate control. The least-squares center of curvature (determined using natural weighting) for all furrows in Marius and Galileo Regio is 20.7 ± 1.1°S, 179.2 ± 1.1°W (2σ uncertainty). Furrows in Marius and Galileo Regio are reasonably concentric, are much more circular than previously estimated, and probably once covered at least an entire hemisphere of Ganymede. In addition, at least three other multiring systems of varying size are identified on Ganymede, indicative of a projectile population. Thus we find furrow geometry and occurrence are consistent with an impact origin. Deviations of some furrows from concentricity about the center of curvature, on the scale of those found at Valhalla, do exist. As in the case of Valhalla these variations are principally confined to outer regions of the structure and are interpreted as inherent properties of multiringed structures on icy satellites. The cause(s) of this may be in the ring formation mechanism itself, but are more likely due to variations in preexisting lithospheric mechanical properties. The perceived present nonalignment of the assumed originally concentric furrows has been used to argue for large-scale lateral motion of dark terrain blocks in Ganymede's crust, presumably in association with bright terrain formation. The overall alignment of furrows as well as the inherent scatter in centers of curvature for subregions of Galileo and Marius Regio do not support this hypothesis.     

Cratering rates on the Galilean satellites

We exploit recent theoretical advances toward the origin and orbital evolution of comets and asteroids to obtain revised estimates for cratering rates in the jovian system. We find that most, probably more than 90%, of the craters on the Galilean satellites are caused by the impact of Jupiter-family comets (JFCs). These are comets with short periods, in generally low-inclination orbits, whose dynamics are dominated by Jupiter. Nearly isotropic comets (long period and Halley-type) contribute at the 1-10% level. Trojan asteroids might also be important at the 1-10% level; if they are important, they would be especially important for smaller craters. Main belt asteroids are currently unimportant, as each 20-km crater made on Ganymede implies the disruption of a 200-km diameter parental asteroid, a destruction rate far beyond the resources of today's asteroid belt. Twenty-kilometer diameter craters are made by kilometer-size impactors; such events occur on a Galilean satellite about once in a million years. The paucity of 20-km craters on Europa indicates that its surface is of order 10 Ma. Lightly cratered surfaces on Ganymede are nominally of order 0.5-1.0 Ga. The uncertainty in these estimates is about a factor of five. Callisto is old, probably more than 4 Ga. It is too heavily cratered to be accounted for by the current flux of JFCs. The lack of pronounced apex-antapex asymmetries on Ganymede may be compatible with crater equilibrium, but it is more easily understood as evidence for nonsynchronous rotation of an icy carapace.     

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.     

Geology of the Selk crater region on Titan from Cassini VIMS observations

Observations of Titan obtained by the Cassini Visual and Infrared Mapping Spectrometer (VIMS) have revealed Selk crater, a geologically young, bright-rimmed, impact crater located ∼800 km north-northwest of the Huygens landing site. The crater rim-crest diameter is ∼90 km; its floor diameter is ∼60 km. A central pit/peak, 20-30 km in diameter, is seen; the ratio of the size of this feature to the crater diameter is consistent with similarly sized craters on Ganymede and Callisto, all of which are dome craters. The VIMS data, unfortunately, are not of sufficient resolution to detect such a dome. The inner rim of Selk crater is fluted, probably by eolian erosion, while the outer flank and presumed ejecta blanket appear dissected by drainages (particularly to the east), likely the result of fluvial erosion. Terracing is observed on the northern and western walls of Selk crater within a 10-15 km wide terrace zone identified in VIMS data; the terrace zone is bright in SAR data, consistent with it being a rough surface. The terrace zone is slightly wider than those observed on Ganymede and Callisto and may reflect differences in thermal structure and/or composition of the lithosphere. The polygonal appearance of the crater likely results from two preexisting planes of weakness (oriented at azimuths of 21° and 122° east of north). A unit of generally bright terrain that exhibits similar infrared-color variation and contrast to Selk crater extends east-southeast from the crater several hundred kilometers. We informally refer to this terrain as the Selk “bench.” Both Selk and the bench are surrounded by the infrared-dark Belet dune field. Hypotheses for the genesis of the optically bright terrain of the bench include: wind shadowing in the lee of Selk crater preventing the encroachment of dunes, impact-induced cryovolcanism, flow of a fluidized-ejecta blanket (similar to the bright crater outflows observed on Venus), and erosion of a streamlined upland formed in the lee of Selk crater by fluid flow. Vestigial circular outlines in this feature just east of Selk’s ejecta blanket suggest that this might be a remnant of an ancient, cratered crust. Evidently the southern margin of the feature has sufficient relief to prevent the encroachment of dunes from the Belet dune field. We conclude that this feature either represents a relatively high-viscosity, fluidized-ejecta flow (a class intermediate to ejecta blankets and long venusian-style ejecta flows) or a streamlined upland remnant that formed downstream from the crater by erosive fluid flow from the west-northwest.     

Thermal evolution of a differentiated Ganymede and implications for surface features

Thermal evolution models are presented for Ganymede, assuming a mostly differentiated initial state of a water ocean overlying a rock layer. The only heat sources are assumed to be primordial heat (provided by accretion) and the long-lived radiogenic heat sources in the rock component. As Ganymede cools, the ocean thins, and two ice layers develop, one above composed of ice I, and the other below composed of high-pressure polymorphs of ice. Subsolidus convection proceeds separately in each ice layer, its transport of heat calculated using a simple parameterized convection scheme and the most recent data on ice rheology. The model requires that the average entropy of the deep ice layer exceeds that of the ice I layer. If the residual ocean separating these layers becomes thin enough, then a Rayleigh-Taylor-like (“diapiric”) instability may ensue, driven by the greater entropy of the deeper ice and merging the two ice mantles into a single convective layer. This instability is not predicted by linear analysis but occurs for plausible finite amplitude perturbations associated with large Rayleigh number convection. The resulting warm ice diapirs may lead to a dramatic “heat pulse” at the surface and to fracturing of the lithosphere, and may be directly or indirectly responsible for resurfacing and grooved terrain formation on Ganymede. The timing of this event depends rather sensitively on poorly known rheological parameters, but could be consistent with chronologies deduced from estimated cratering rates. Irrespective of the occurrence or importance of the heat pulse, we find that lithospheric fracturing requires rapid stress loading (on a time scale ⪅10⁴ years). Such a time scale can be realized by warm ice diapirism, but not directly by gradual global expansion. In the absence of any quantitative and self-consistent model for the resurfacing of Ganymede by liquid water, we favor resurfacing by warm ice flows, which we demonstrate to be physically possible, a plausible consequence of our models, compatible with existing observations, and a hypothesis testable by Galileo. We discuss core formation as an alternative driver for resurfacing, and conclude that it is less attractive. We also consider anew the puzzle of why Callisto differs so greatly from Ganymede, offering several possible explanations. The models presented do not provide a compelling explanation for all aspects of Ganymedean geological evolution, since we have identified several potential problems, most notably the apparently extended period of grooved terrain formation (several hundred million years), which is difficult to reconcile with the heat pulse phenomenon.     

Global geological mapping of Ganymede

We have compiled a global geological map of Ganymede that represents the most recent understanding of the satellite based on Galileo mission results. This contribution builds on important previous accomplishments in the study of Ganymede utilizing Voyager data and incorporates the many new discoveries that were brought about by examination of Galileo data. We discuss the material properties of geological units defined utilizing a global mosaic of the surface with a nominal resolution of 1 km/pixel assembled by the USGS with the best available Voyager and Galileo regional coverage and high resolution imagery (100-200 m/pixel) of characteristic features and terrain types obtained by the Galileo spacecraft. We also use crater density measurements obtained from our mapping efforts to examine age relationships amongst the various defined units. These efforts have resulted in a more complete understanding of the major geological processes operating on Ganymede, especially the roles of cryovolcanic and tectonic processes in the formation of might materials. They have also clarified the characteristics of the geological units that comprise the satellite’s surface, the stratigraphic relationships of those geological units and structures, and the geological history inferred from those relationships. For instance, the characteristics and stratigraphic relationships of dark lineated material and reticulate material suggest they represent an intermediate stage between dark cratered material and light material units.     

Ganymede crater dimensions – Implications for central peak and central pit formation and development

The morphology of impact craters on the icy Galilean satellites differs from craters on rocky bodies. The differences are thought due to the relative weakness of ice and the possible presence of sub-surface water layers. Digital elevation models constructed from Galileo images were used to measure a range of dimensions of craters on the dark and bright terrains of Ganymede. Measurements were made from multiple profiles across each crater, so that natural variation in crater dimensions could be assessed and averaged scaling trends constructed. The additional depth, slope and volume information reported in this work has enabled study of central peak formation and development, and allowed a quantitative assessment of the various theories for central pit formation. We note a possible difference in the size-morphology progression between small craters on icy and silicate bodies, where central peaks occur in small craters before there is any slumping of the crater rim, which is the opposite to the observed sequence on the Moon. Conversely, our crater dimension analyses suggest that the size-morphology progression of large lunar craters from central peak to peak-ring is mirrored on Ganymede, but that the peak-ring is subsequently modified to a central pit morphology. Pit formation may occur via the collapse of surface material into a void left by the gradual release of impact-induced volatiles or the drainage of impact melt into sub-crater fractures.