Tectonic framework of grooved terrain on ganymede

An analysis of the tectonic features (grooves and prominent troughs) comprising grooved terrain on Ganymede shows a global tectonic framework. Domains of grooves, bordered by deep structural boundaries that controlled their extent and orientation, developed along preferred orientations. Some of these directions are much more prominent than others and appear to coincide with at least two great circles with pole locations at 152° long. 55° lat. and 65° long. and 50° lat. The great circle hypothesis is in good agreement with several statistical tests of the data set. Several mechanisms for the emplacement of grooved areas along great circles are discussed. A plausible mechanism is a pattern of equatorial-symmetric rising and descending convection cells. Because the studied areas represent only a small portion of Ganymede, interpretations of global tectonic patterns and models of formation must remain tentative.     

Thermal-orbital evolution of Io and Europa

Coupled thermal-orbital evolution models of Europa and Io are presented. It is assumed that Io, Europa, and Ganymede evolve in the Laplace resonance and that tidal dissipation of orbital energy is an internal heat source for both Io and Europa. While dissipation in Io occurs in the mantle as in the mantle dissipation model of Segatz et al. (1988, Icarus 75, 187), two models for Europa are considered. In the first model dissipation occurs in the silicate mantle while in the second model dissipation occurs in the ice shell. In the latter model, ice shell melting and variations of the shell thickness above an ocean are explicitly included. The rheology of both the ice and the rock is cast in terms of a viscoelastic Maxwell rheology with viscosity and shear modulus depending on the average temperature of the dissipating layer. Heat transfer by convection is calculated using a parameterization for strongly temperature-dependent viscosity convection. Both models are consistent with the present orbital elements of Io, Europa, and Ganymede. It is shown that there may be phases of quasi-steady evolution with large or small dissipation rates (in comparison with radiogenic heating), phases with runaway heating or cooling and oscillatory phases during which the eccentricity and the tidal heating rate will oscillate. Europa's ice thickness varies between roughly 3 and 70 km (dissipation in the silicate layer) or 10 and 60 km (dissipation in the ice layer), suggesting that Europa's ocean existed for geological timescales. The variation in ice thickness, including both convective and purely conductive phases, may be reflected in the formation of different geological surface features on Europa. Both models suggest that at present Europa's ice thickness is several tens of km thick and is increasing, while the eccentricity decreases, implying that the satellites evolve out of resonance. Including lithospheric growth in the models makes it impossible to match the high heat flux constraint for Io. Other heat transfer processes than conduction through the lithosphere must be important for the present Io.     

On convection in ice I shells of outer Solar System bodies, with detailed application to Callisto

It has been argued that the dominant non-Newtonian creep mechanisms of water ice make the ice shell above Callisto's ocean, and by inference all radiogenically heated ice I shells in the outer Solar System, stable against solid-state convective overturn. Conductive heat transport and internal melting (oceans) are therefore predicted to be, or have been, widespread among midsize and larger icy satellites and Kuiper Belt objects. Alternatively, at low stresses (where non-Newtonian viscosities can be arbitrarily large), convective instabilities may arise in the diffusional creep regime for arbitrarily small temperature perturbations. For Callisto, ice viscosities are low enough that convection is expected over most of geologic time above the internal liquid layer for plausible ice grain sizes (?a few mm); the alternative for early Callisto, a conducting shell over a very deep ocean (>450 km), is not compatible with Callisto's present partially differentiated state. Moreover, if convection is occurring today, the stagnant lid would be quite thick (∼100 km) and compatible with the lack of active geology. Nevertheless, Callisto's steady-state heat flow may have fallen below the convective minimum for its ice I shell late in Solar System history. In this case convection ends, the ice shell melts back at its base, and the internal ocean widens considerably. The presence of such an ocean, of order 200 km thick, is compatible with Callisto's moment-of-inertia, but its formation would have caused an ∼0.25% radial expansion. The tectonic effects of such a late, slow expansion are not observed, so convection likely persists in Callisto, possibly subcritically. Ganymede, due to its greater size, rock fraction and full differentiation, has a substantially higher heat flow than Callisto and has not reached this tectonic end state. Titan, if differentiated, and Triton should be more similar to Ganymede in this regard. Pluto, like Callisto, may be near the tipping point for convective shutdown, but uncertainties in its size and rock fraction prevent a more definitive assessment.     

Oceans in the icy Galilean satellites of Jupiter?

Equilibrium models of heat transfer by heat conduction and thermal convection show that the three satellites of Jupiter—Europa, Ganymede, and Callisto—may have internal oceans underneath ice shells tens of kilometers to more than a hundred kilometers thick. A wide range of rheology and heat transfer parameter values and present-day heat production rates have been considered. The rheology was cast in terms of a reference viscosity ν₀ calculated at the melting temperature and the rate of change A of viscosity with inverse homologous temperature. The temperature dependence of the thermal conductivity k of ice I has been taken into account by calculating the average conductivity along the temperature profile. Heating rates are based on a chondritic radiogenic heating rate of 4.5 pW kg⁻¹ but have been varied around this value over a wide range. The phase diagrams of H₂O (ice I) and H₂O + 5 wt% NH₃ ice have been considered. The ice I models are worst-case scenarios for the existence of a subsurface liquid water ocean because ice I has the highest possible melting temperature and the highest thermal conductivity of candidate ices and the assumption of equilibrium ignores the contribution to ice shell heating from deep interior cooling. In the context of ice I models, we find that Europa is the satellite most likely to have a subsurface liquid ocean. Even with radiogenic heating alone the ocean is tens of kilometers thick in the nominal model. If tidal heating is invoked, the ocean will be much thicker and the ice shell will be a few tens of kilometers thick. Ganymede and Callisto have frozen their oceans in the nominal ice I models, but since these models represent the worst-case scenario, it is conceivable that these satellites also have oceans at the present time. The most important factor working against the existence of subsurface oceans is contamination of the outer ice shell by rock. Rock increases the density and the pressure gradient and shifts the triple point of ice I to shallower depths where the temperature is likely to be lower then the triple point temperature. According to present knowledge of ice phase diagrams, ammonia produces one of the largest reductions of the melting temperature. If we assume a bulk concentration of 5 wt% ammonia we find that all the satellites have substantial oceans. For a model of Europa heated only by radiogenic decay, the ice shell will be a few tens of kilometers thinner than in the ice I case. The underlying rock mantle will limit the depth of the ocean to 80-100 km. For Ganymede and Callisto, the ice I shell on top of the H₂O-NH₃ ocean will be around 60- to 80-km thick and the oceans may be 200- to 350-km deep. Previous models have suggested that efficient convection in the ice will freeze any existing ocean. The present conclusions are different mainly because they are based on a parameterization of convective heat transport in fluids with strongly temperature dependent viscosity rather than a parameterization derived from constant-viscosity convection models. The present parameterization introduces a conductive stagnant lid at the expense of the thickness of the convecting sublayer, if the latter exists at all. The stagnant lid causes the temperature in the sublayer to be warmer than in a comparable constant-viscosity convecting layer. We have further modified the parameterization to account for the strong increase in homologous temperature, and therefore decrease in viscosity, with depth along an adiabat. This modification causes even thicker stagnant lids and further elevated temperatures in the well-mixed sublayer. It is the stagnant lid and the comparatively large temperature in the sublayer that frustrates ocean freezing.     

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.     

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.     

Onset of Convection,Heat Flow and Thickness of the Europa‘s ice Shell

The observational evidence given by Galileo spacecraft about Europa supports an icy rigid layer of several kilometers over another ductile layer of ice in convection, which floats over an internal ocean of liquid water. Before the onset of convection, heat is transmitted into the crust by conduction. The heat flow analysis in the potentially convective layer gives values higher than those obtained previously by tidal dissipation models, and suggests that the ice may be limited to a thin layer of ∼4 km total thickness. This revised version was published online in July 2006 with corrections to the Cover Date.     

Grooved terrain on Ganymede

The icy crust of Ganymede comprises bright and dark areas. Investigation of Voyager 1 and 2 images has shown that bright terrain is grooved and separates dark polygons of cratered terrain. The grooved terrain contains alternating ridges and grooves in straight and curvilinear sets, which are locally interrupted by smooth patches and swaths. Cratered terrain, where 'it occurs in small wedges and slivers, has a pervasive grain of narrowly spaced furrows, and thus is transitional to grooved terrain. An analysis of the morphology of terrain types, and of superposition and cross-cutting relations, suggests that grooved terrain grew at the expense of cratered terrain, that tracts of cratered terrain were converted into grooved terrain in situ, and that vertical tectonism and shear movements dominated in the restructuring of Ganymede's surface. It is postulated that during a period in the planet's history when the lithosphere was thin, upwelling convection currents caused incipient rifting accompanied by intensive normal faulting; where rifting went to completion, crustal segments separated, locally spread apart, and sheared past one another. In places subduction and compression may have occured, but the evidence is inconclusive. Thus, the grooved terrain on Ganymede may record an early phase of ice-plate tectonics that caused rifting and drifting of the icy lithosphere, but, unlike silicate plate tectonics on Earth, may have resulted in only minor vertical turnover.     

On the resurfacing of Ganymede by liquid-water volcanism

A long-popular model for producing Ganymede's bright terrain involves flooding of low-lying graben with liquid water, slush, or warm, soft ice. The model suffers from major problems, however, including the absence of obvious near-surface heat sources, the negative buoyancy of liquid water, and the lack of a mechanism for confining the flows to graben floors. We present new models for cryovolcanic resurfacing to overcome these difficulties. Tidal heating within an ancient Laplace-like orbital resonance (Showman and Malhotra 1997, Icarus 127, 93; Showman et al., 1997, Icarus 129, 367) provides a plausible heat source and could allow partial melting to occur as shallow as 5-10 km depth. Our favored mechanism for delivering this water to the surface invokes the fact that topography—such as a global set of graben—causes subsurface pressure gradients that can pump water or slush upward onto the floors of topographic lows (graben) despite the negative buoyancy of the liquid. These eruptions can occur only within the topographic lows; furthermore, as the low areas become full, the pressure gradients disappear and the resurfacing ceases. This provides an explanation for the observed straight dark-bright terrain boundaries: water cannot overflow the graben, so resurfacing rarely embays craters or other rough topography. Pure liquid water can be pumped to the surface from only 5-10 km depth, but macroscopic bodies of slush ascending within fractures can reach the surface from much greater depths due to the smaller negative buoyancy of slush. A challenge for these models is the short predicted gravitational relaxation timescale of topographic features at high heat flows; the resurfacing must occur before the graben topography disappears. We also evaluate alternate resurfacing mechanisms, such as pumping of liquid water to the surface by thermal expansion stresses and buoyant rise of water through a silicate-contaminated crust that is denser than liquid water, and conclude that they are unlikely to explain Ganymede's bright terrain.     

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.     

The origin of Ganymede's polar caps

Since their discovery in Voyager images, the origin of the bright polar caps of Ganymede has intrigued investigators. Some models attributed the polar cap formation to thermal migration of water vapor to higher latitudes, while other models implicated plasma bombardment in brightening ice. Only with the arrival of Galileo at Jupiter was it apparent that Ganymede possesses a strong internal magnetic field, which blocks most of the plasma from bombarding the satellite's equatorial region while funneling plasma onto the polar regions. This discovery provides a plausible explanation for the polar caps as related to differences in plasma-induced brightening in the polar and the equatorial regions. In this context, we analyze global color and high resolution images of Ganymede obtained by Galileo, finding a very close correspondence between the observed polar cap boundary and the open/closed field lines boundary obtained from new modeling of the magnetic field environment. This establishes a clear link between plasma bombardment and polar cap brightening. High resolution images show that bright polar terrain is segregated into bright and dark patches, suggesting sputter-induced redistribution and subsequent cold trapping of water molecules. Minor differences between the location of the open/closed field lines boundary and the observed polar cap boundary may be due to interaction of Ganymede with Jupiter's magnetosphere, and our neglect of higher-order terms in modeling Ganymede's internal field. We postulate that leading-trailing brightness differences in Ganymede's low-latitude surface are due to enhanced plasma flux onto the leading hemisphere, rather than darkening of the trailing hemisphere. In contrast to Ganymede, the entire surface of Europa is bombarded by jovian plasma, suggesting that sputter-induced redistribution of water molecules is a viable means of brightening that satellite's surface.     

Resurfacing history of Europa from pole-to-pole geological mapping

We produced geologic maps from two regional mosaics of Galileo images across the leading and trailing hemispheres of Europa in order to investigate the temporal distribution of units in the visible geologic record. Five principal terrain types were identified (plains, bands, ridges, chaos, and crater materials), which are interpreted to result from (1) tectonic fracturing and lineament building, (2) cryovolcanic reworking of surface units, with possible emplacement of sub-surface materials, and (3) impact cratering. The geologic histories of both mapped areas are essentially similar and reflect some common trends: Tectonic resurfacing dominates the early geologic record with the formation of background plains by intricate superposition of lineaments, the opening of wide bands with infilling of inter-plate gaps, and the buildup of ridges and ridge complexes along prominent fractures in the ice. It also appears that lineaments are narrower and more widely spaced with time. The lack of impact craters overprinted by lineaments indicate that the degree of tectonic resurfacing decreased rapidly after ridged plains formation. In contrast, the degree of cryovolcanic resurfacing appears to increase with time, as chaos formation dominates the later parts of the geologic record. These trends, and the transition from tectonic- to cryovolcanic-dominated resurfacing could be attributed to the gradual thickening of Europa's cryosphere during the visible geologic history, that comprises the last 2% or 30-80 Myr of Europa's history: An originally thin, brittle ice shell could be pervasively fractured or melted through by tidal and endogenic processes; the degree of fracturing and plate displacements decreased with time in a thickening shell, and lineaments became narrower and more widely spaced; formation of chaos regions could have occurred where the thickness threshold for solid-state convection was exceeded, and can be aided by preferential tidal heating of more ductile ice. In a long-term context it is not clear at this point whether this inferred thickening trend would reflect a drastic change in the thermal evolution of the satellite, or cyclic or irregular episodes of tectonic and cryovolcanic activity.     

Heat flow, lenticulae spacing, and possibility of convection in the ice shell of europa

Two opposing models to explain the geological features observed on Europa’s surface have been proposed. The thin-shell model states that the ice shell is only a few kilometers thick, transfers heat by conduction only, and can become locally thinner until it exposes an underlying ocean on the satellite’s surface. According to the thick-shell model, the ice shell may be several tens of kilometers thick and have a lower convective layer, above which there is a cold stagnant lid that dissipates heat by conduction. Whichever the case, from magnetic data there is strong support for the presence of a layer of salty liquid water under the ice. The present study was performed to examine whether the possibility of convection is theoretically consistent with surface heat flows of ∼100-200 mW m⁻², deduced from a thin brittle lithosphere, and with the typical spacing of 15-23 km proposed for the features usually known as lenticulae. It was obtained that under Europa’s ice shell conditions convection could occur and also account for high heat flows due to tidal heating of the convective (nearly isothermal) interior, but only if the dominant water ice rheology is superplastic flow (with activation energy of 49 kJ mol⁻¹; this is the rheology thought dominant in the warm interior of the ice shell). In this case the ice shell would be ∼15-50 km thick. Furthermore, in this scenario explaining the origin of the lenticulae related to convective processes requires ice grain size close to 1 mm and ice thickness around 15-20 km.     

The orbital-thermal evolution and global expansion of Ganymede

The tectonically and cryovolcanically resurfaced terrains of Ganymede attest to the satellite's turbulent geologic history. Yet, the ultimate cause of its geologic violence remains unknown. One plausible scenario suggests that the Galilean satellites passed through one or more Laplace-like resonances before evolving into the current Laplace resonance. Passage through such a resonance can excite Ganymede's eccentricity, leading to tidal dissipation within the ice shell. To evaluate the effects of resonance passage on Ganymede's thermal history we model the coupled orbital-thermal evolution of Ganymede both with and without passage through a Laplace-like resonance. In the absence of tidal dissipation, radiogenic heating alone is capable of creating large internal oceans within Ganymede if the ice grain size is 1 mm or greater. For larger grain sizes, oceans will exist into the present epoch. The inclusion of tidal dissipation significantly alters Ganymede's thermal history, and for some parameters (e.g. ice grain size, tidal Q of Jupiter) a thin ice shell (5 to 20 km) can be maintained throughout the period of resonance passage. The pulse of tidal heating that accompanies Laplace-like resonance capture can cause up to 2.5% volumetric expansion of the satellite and contemporaneous formation of near surface partial melt. The presence of a thin ice shell and high satellite orbital eccentricity would generate moderate diurnal tidal stresses in Ganymede's ice shell. Larger stresses result if the ice shell rotates non-synchronously. The combined effects of satellite expansion, its associated tensile stress, rapid formation of near surface partial melt, and tidal stress due to an eccentric orbit may be responsible for creating Ganymede's unique surface features.     

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.     

The evolution of tectonic features on Ganymede

The bands of bright resurfaced terrain on Ganymede are probably broad grabens formed by global expansion and filled with deposits of ice. Grooves within the bands are thought to be extensional features formed during the same episode of expansion. The crust of Ganymede is modeled as a viscoelastic material subjected to extensional strain. With sufficiently high strain rates and stresses, deep normal faulting will occur, creating broad grabens that may then be filled. Continuing deformation at high strain rates and stresses will cause propagation of deep faults up into the flood deposits and normal faulting at the surface, while lower strain rates and stresses will cause formation of open extension fractures or, if the crustal strength is very low, grabens at the surface. The spacing between adjacent fractures may reflect the geothermal gradient at the time of deformation. Surface topography resulting from fracturing and normal faulting will decay with time as a result of viscous relaxation and mass-wasting.     

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Librations of the Galilean satellites: The influence of global internal liquid layers

The four Galilean satellites are thought to harbor one or even two global internal liquid layers beneath their surface layer. The iron core of Io and Ganymede is most likely (partially) liquid and also the core of Europa may be liquid. Furthermore, there are strong indications for the existence of a subsurface ocean in Europa, Ganymede, and Callisto. Here, we investigate whether libration observations can be used to prove the existence of these liquid layers and to constrain the thickness of the overlying solid layers. For Io, the presence of a small liquid core increases the libration of the mantle by a few percent with respect to an entirely solid Io and mantle libration observations could be used to determine the mantle thickness with a precision of several tens of kilometers given that the libration amplitude can be measured with a precision of 1 m. For Europa, Ganymede, and Callisto, the presence of a water ocean close to the surface increases by at least an order of magnitude the ice shell libration amplitude with respect to an entirely solid satellite. The shell libration depends essentially on the shell thickness and to a minor extent on the density difference between the ocean and the ice shell. The possible presence of a liquid core inside Europa and Ganymede has no noticeable influence on their shell libration. For a precision of several meters on the libration measurements, in agreement with the expected accuracy with the NASA/ESA EJSM orbiter mission to Europa and Ganymede, an error on the shell thickness of a few tens kilometers is expected. Therefore, libration measurements can be used to detect liquid layers such as Io’s core or water subsurface oceans in Europa, Ganymede, and Callisto and to constrain the thickness of the overlying solid surface layers.     

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.