The formation of Ganymede's grooved terrain: Numerical modeling of extensional necking instabilities

Ganymede's grooved terrain likely formed during an epoch of global expansion, when unstable extension of the lithosphere resulted in the development of periodic necking instabilities. Linear, infinitesimal-strain models of extensional necking support this model of groove formation, finding that the fastest growing modes of an instability have wavelengths and growth rates consistent with Ganymede's grooves. However, several questions remain unanswered, including how nonlinearities affect instability growth at large strains, and what role instabilities play in tectonically resurfacing preexisting terrain. To address these questions we numerically model the extension of an icy lithosphere to examine the growth of periodic necking instabilities over a broad range of strain rates and temperature gradients. We explored thermal gradients up to 45 K km⁻¹ and found that, at infinitesimal strain, maximum growth rates occur at high temperature gradients (45 K km⁻¹) and moderate strain rates (10⁻¹³ s⁻¹). Dominant wavelengths range from 1.8 to 16.4 km (post extension). Our infinitesimal growth rates are qualitatively consistent with, but an order of magnitude lower than, previous linearized calculations. When strain exceeds ∼10% growth rates decrease, limiting the total amount of amplification that can result from unstable extension. This fall-off in growth occurs at lower groove amplitudes for high-temperature-gradient, thin-lithosphere simulations than for low-temperature-gradient, thick-lithosphere simulations. At large strains, this shifts the ideal conditions for producing large amplitude grooves from high temperature gradients to more moderate temperature gradients (15 K km⁻¹). We find that the formation of periodic necking instabilities can modify preexisting terrain, replacing semi-random topography up to 100 m in amplitude with periodic ridges and troughs, assisting the tectonic resurfacing process. Despite this success, the small topographic amplification produced by our model presents a formidable challenge to the necking instability mechanism for groove formation. Success of the necking instability mechanism may require rheological weakening or strain localization by faulting, effects not included in our analysis.     

Constraints on the location, magnitude, and dimensions of Ganymede's mass anomalies

Previously, radio Doppler data, generated with NASA's Galileo spacecraft during its second encounter with Jupiter's moon Ganymede, were used to infer the locations and magnitudes of mass anomalies on Ganymede using point-mass models. However, the point-mass solutions do not provide the vertical and horizontal extent of the anomalous mass concentrations. Here, we provide the results of a new study using spherical cap disks to model Ganymede's mass anomalies. The spherical cap disk models not only provide the locations and magnitudes of the mass anomalies, but also their vertical and horizontal dimensions. The new models show that three disks, a positive mass located at (53.0° N, 127.0° W) and two negative masses located at (22.0° N, 87.0° W) and (49.0° N, 219.0° W), can explain the data. The magnitudes of the mass anomalies are on the order of 10¹⁸ kg. The diameters of the anomalies are a few thousand kilometers. The positive anomaly is about 100 meters thick and both negative anomalies have a thickness of less than a kilometer. We use the additional information provided by the disk models to investigate the viability of mass anomalies at Ganymede's surface by comparing the diameters of the anomalies to the sizes of regiones and sulci and the anomalies' thicknesses to accumulated layers of rock and clean ice on the surface. We find that the dimensions of the mass anomalies could be explained by concentrations of rock in the regio and rock-free ice in the sulci. These results confirm that mass anomalies may reside on or near Ganymede's surface and that positive mass anomalies are contained within areas of dark terrain and negative mass anomalies within bright terrain.     

Resurfacing of Titan by ammonia-water cryomagma

The Cassini Titan Radar Mapper observed on Titan several large features interpreted as cryovolcanic during the October 26, 2004 pass at high northern latitudes [Lopes, R.M.C., and 43 colleagues, 2007. Icarus 186, 395-412]. To date, models of ammonia-water resurfacing have not been tied to specific events or evolutionary stages of Titan. We propose a model of cryovolcanism that involves cracking at the base of the ice shell and formation of ammonia-water pockets in the ice. As these ammonia-water pockets undergo partial freezing in the cold ice shell, the ammonia concentration in the pockets increases, decreasing the negative buoyancy of the ammonia-water mixture. If the ice shell is contaminated by silicates delivered in impacts, the liquid-solid density difference would be even less. While the liquid cannot easily become buoyant relative to the surrounding ice, these concentrated ammonia-water pockets are sufficiently close to the neutral buoyancy point that large-scale tectonic stress patterns (tides, non-synchronous rotation, satellite volume changes, solid state convection, or subsurface pressure gradients associated with topography) would enable the ammonia to erupt effusively onto the surface. Rather than suggesting steady-state volcanism over the history of the Solar System, we favor a scenario where the cryovolcanic features could have been associated with episodic (potentially late) geological activity.     

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.     

Thermal convection in ice-I shells of Titan and Enceladus

Cassini-Huygens observations have shown that Titan and Enceladus are geologically active icy satellites. Mitri and Showman [Mitri, G., Showman, A.P., 2005. Icarus 177, 447-460] and McKinnon [McKinnon, W.B., 2006. Icarus 183, 435-450] investigated the dynamics of an ice shell overlying a pure liquid-water ocean and showed that transitions from a conductive state to a convective state have major implications for the surface tectonics. We extend this analysis to the case of ice shells overlying ammonia-water oceans. We explore the thermal state of Titan and Enceladus ice-I shells, and also we investigate the consequences of the ice-I shell conductive-convective switch for the geology. We show that thermal convection can occur, under a range of conditions, in the ice-I shells of Titan and Enceladus. Because the Rayleigh number Ra scales with δ³/ηb, where δ is the thickness of the ice shell and ηb is the viscosity at the base of the ice-I shell, and because ammonia in the liquid layer (if any) strongly depresses the melting temperature of the water ice, Ra equals its critical value for two ice-I shell thicknesses: for relatively thin ice shell with warm, low-viscosity base (Onset I) and for thick ice shell with cold, high-viscosity base (Onset II). At Onset I, for a range of heat fluxes, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. Switches between these states can cause large, rapid changes in the ice-shell thickness. For Enceladus, we demonstrate that an Onset I transition can produce tectonic stress of ∼500 bars and fractures of several tens of km depth. At Onset II, in contrast, we demonstrate that zero equilibrium states exist for a range of heat fluxes. For a mean heat flux within this range, the satellite experiences oscillations in surface heat flux and satellite volume with periods of ∼50-800 Myr even when the interior heat production is constant or monotonically declining in time; these oscillations in the thermal state of the ice-I shell would cause repeated episodes of extensional and compressional tectonism.     

Viscous relaxation of craters within the martian south polar layered deposits

The South Polar Layered Deposits (PLD) are of fundamental importance to martian climatology, as they may comprise the largest reservoir of near-surface water on Mars. The South PLD exhibit relatively young crater retention surface ages, which are widely attributed to recent resurfacing. However, we show that both constructional and destructional resurfacing mechanisms (such as dust deposition and water ice sublimation, respectively) are inconsistent with the size, depth, and spatial distributions of South PLD craters. We demonstrate that another process—viscous creep relaxation of dusty water ice—is more compatible with the observed cratering of South PLD surfaces. The results of our finite element relaxation simulations suggest that, despite their apparent youthfulness, the PLD have been stable for at least several hundred million and perhaps even over a billion years. Consequently, our modeling implies that the time scales for the formation and preservation of the layers characteristic of the South (and possibly North) PLD are much longer than generally assumed.     

The topography of Iapetus' leading side

We have used Cassini stereo images to study the topography of Iapetus' leading side. A terrain model derived at resolutions of 4-8 km reveals that Iapetus has substantial topography with heights in the range of −10 km to +13 km, much more than observed on the other middle-sized satellites of Saturn so far. Most of the topography is older than 4 Ga [Neukum, G., Wagner, R., Denk, T., Porco, C.C., 2005. Lunar Planet. Sci. XXXVI. Abstract 2034] which implies that Iapetus must have had a thick lithosphere early in its history to support this topography. Models of lithospheric deflection by topographic loads provide an estimate of the required elastic thickness in the range of 50-100 km. Iapetus' prominent equatorial ridge [Porco, C.C., and 34 colleagues, 2005. Science 307, 1237-1242] reaches widths of 70 km and heights of up to 13 km from their base within the modeled area. The morphology of the ridge suggests an endogenous origin rather than a formation by collisional accretion of a ring remnant [Ip, W.-H., 2006. Geophys. Res. Lett. 33, doi:10.1029/2005GL025386. L16203]. The transition from simple to complex central peak craters on Iapetus occurs at diameters of 11±3 km. The central peaks have pronounced conical shapes with flanking slopes of typically 11° and heights that can rise above the surrounding plains. Crater depths seem to be systematically lower on Iapetus than on similarly sized Rhea, which if true, may be related to more pronounced crater-wall slumping (which widens the craters) on Iapetus than on Rhea. There are seven large impact basins with complex morphologies including central peak massifs and terraced walls, the largest one reaches 800 km in diameter and has rim topography of up to 10 km. Generally, no rings are observed with the basins consistent with a thick lithosphere but still thin enough to allow for viscous relaxation of the basin floors, which is inferred from crater depth-to-diameter measurements. In particular, a 400-km basin shows up-domed floor topography which is suggestive of viscous relaxation. A model of complex crater formation with a viscoplastic (Bingham) rheology [Melosh, H.J., 1989. Impact Cratering. Oxford Univ. Press, New York] of the impact-shocked icy material provides an estimate of the effective cohesion/viscosity at . The local distribution of bright and dark material on the surface of Iapetus is largely controlled by topography and consistent with the dark material being a sublimation lag deposit originating from a bright icy substrate mixed with the dark components, but frost deposits are possible as well.     

Strength of Faults on Mars from MOLA Topography

The stresses which must be maintained on faults bounding the rift topography at Tempe Fossae—the “North Tempe Rift” (NTR)—and Valles Marineris (VM) on Mars are estimated, using a simple elastic model and topographic data from the Mars Orbiter Laser Altimeter (MOLA). The absence of rift-flank uplift at the NTR is consistent with an elastic thickness, T_e, of 20 km or greater at the time of rift formation. The maximum resolved shear stresses on bounding faults due to this topography do not therefore exceed 20 MPa, similar to the inferred strength of terrestrial faults. Elastic thickness estimates at VM are mostly around 50 km or greater. Therefore, for canyon widths of ∼400 km, the bounding faults of VM, if present, must be able to withstand stresses of up to approximately 100 MPa. However, if the fault-controlled sections of the canyons do not exceed 150 km in width, as suggested by geomorphological analysis, the fault strength required is only 20 MPa. Although the maximum resolved shear stresses required to support the topography at VM may need to be greater than the stresses which terrestrial faults can support, at least some faults on Mars are no stronger than similar features on Earth. This observation is consistent with the existence of liquid water in the shallow subsurface of Mars at the time the faults were active. On Venus, plate tectonics is probably prevented by the frictional resistance to motion across strong faults. On Mars, it is more likely that the large thickness of the elastic layer of the lithosphere and the possible positive buoyancy of the crust are responsible for the observed lack of plate tectonics.     

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.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.     

HST BVI photometry of Triton and Proteus

BVI photometry of Triton and Proteus was derived from HST images taken in 1997. The VEGAMAG photometric technique was used. Triton was found to be brighter by a few percent than observations of the 1970's and 1980's, as expected due to the increasingly greater exposure of the bright south polar region. The leading side was also found to be brighter than the trailing side by 0.09 mag in all filters—50% larger than reported by Franz [Franz, O.G., 1981. Icarus 45, 602-606]. Contrary to our previous results [Pascu, D., et al., 1998. Bull. Am. Astron. Soc. 30, 1101], we found no episodic reddening. Our previous conclusions were based on an inaccurate early version of the Charge Transfer Efficiency (CTE) correction. The present result limits the start of the reddening event reported by Hicks and Buratti [Hicks, M.D., Buratti, B.J., 2004. Icarus 171, 210-218]. Our (B-V) result of 0.70±0.01 supports the global blueing described by Buratti et al. [Buratti, B.J., Goguen, J.D., Gibson, J., Mosher, J., 1994. Icarus 110, 303-314]. Our observations of July 1997 agree with the Voyager results and are among the bluest colors seen. We found Proteus somewhat brighter than earlier studies, but in good agreement with the recent value given by Karkoschka [Karkoschka, E., 2003. Icarus 162, 400-407]. A leading/trailing brightness asymmetry was detected for Proteus, with the leading side 0.1 mag brighter. The unique differences in action of the endogenic and exogenic processes on Triton and Proteus provides an opportunity to separate the endogenic and exogenic effects on Triton.     

Preliminary mariner 9 report on the geology of Mars

Mariner 9 pictures indicate that the surface of Mars has been shaped by impact, volcanic, tectonic, erosional and depositional activity. The moonlike cratered terrain, identified as the dominant surface unit from the Mariner 6 and 7 flyby data, has proven to be less typical of Mars than previously believed, although extensive in the mid- and high-latitude regions of the southern hemisphere. Martian craters are highly modified but their size-frequency distribution and morphology suggest that most were formed by impact. Circular basins encompassed by rugged terrain and filled with smooth plains material are recognized. These structures, like the craters, are more modified than corresponding features on the Moon and they exercise a less dominant influence on the regional geology. Smooth plains with few visible craters fill the large basins and the floors of larger craters; they also occupy large parts of the northern hemisphere where the plains lap against higher landforms. The middle northern latitudes of Mars from 90 to 150† longitude contain at least four large shield volcanoes each of which is about twice as massive as the largest on Earth. Steep-sided domes with summit craters and large, fresh-appearing volcanic craters with smooth rims are also present in this region. Multiple flow structures, ridges with lobate flanks, chain craters, and sinuous rilles occur in all regions, suggesting widespread volcanism. Evidence for tectonic activity postdating formation of the cratered terrain and some of the plains units is abundant in the equatorial area from 0 to 120° longitude.Some regions exhibit a complex semiradial array of graben that suggest doming and stretching of the surface. Others contain intensity faulted terrain with broader, deeper graben separated by a complex mosaic of flat-topped blocks. An east-west-trending canyon system about 100–200 km wide and about 2500 km long extends through the Coprates-Eos region. The canyons have gullied walls indicative of extensive headward erosion since their initial formation. Regionally depressed areas called chaotic terrain consist of intricately broken and jumbled blocks and appear to result from breaking up and slumping of older geologic units. Compressional features have not been identified in any of the pictures analyzed to data. Plumose light and dark surface markings can be explained by eolian transport. Mariner 9 has thus revealed that Mars is a complex planet with its own distinctive geologic history and that it is less primitive than the Moon.     

Solid tidal friction above a liquid water reservoir as the origin of the south pole hotspot on Enceladus

Earth, Jupiter's moon Io and Saturn's tiny moon Enceladus are the only solid objects in the Solar System to be sufficiently geologically active for their internal heat to be detected by remote sensing. Interestingly, the endogenic activity on Enceladus is only located on a specific region at the south pole, from which jets of water vapor and ice particles have been observed [Spencer, J.R., and 9 colleagues, 2006. Science 311, 1401-1405; Porco, C.C., and 24 colleagues, 2006. Science 311, 1393-1401]. The current polar location of the thermal anomaly can possibly be explained by diapir-induced reorientation of the satellite [Nimmo, F., Pappalardo, R.T., 2006. Nature 441, 614-616], but the thermal anomaly triggering and the heat power required to sustain it over geological timescales remain problematic. Using a three-dimensional viscoelastic numerical model simulating the response of Enceladus to tidal forcing, we explore the effect of a low-viscosity anomaly in the ice shell, localized to the south polar region, on the tidal dissipation patterns. We demonstrate that only interior models with a liquid water layer at depth can explain the observed magnitude of dissipation rate and its particular location at the south pole. Moreover, we show that tidally-induced heat must be generated over a relatively broad region in the southern hemisphere, and it is then transferred toward the south pole where it is episodically released during relatively short resurfacing events. As large tidal dissipation and internal melting cannot be induced in the south polar region in the absence of a pre-existing liquid decoupling layer, we propose that liquid water must have been present in the interior for a very long period of time, and possibly since the satellite formation. Owing to the orbital equilibrium requirement [Meyer, J., Wisdom, J., 2007. Icarus 188, 535-539], sustaining some liquid water at depth is impossible if heat is continuously emitted at a rate of 4-8 GW at the south pole. Based on that requirement, we propose that the current thermal emission is not in equilibrium with the heat production, and that the thermal emission rate is abnormally high at present time. Alternatively, continuous dissipation at a rate of 1-2 GW in the ice shell at the south pole should be sufficient to induce internal melting and it could sustain a layer of liquid water at depth over geologic timescales.     

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.     

Topography and geomorphology of the Huygens landing site on Titan

The Descent Imager/Spectral Radiometer (DISR) aboard the Huygens Probe took several hundred visible-light images with its three cameras on approach to the surface of Titan. Several sets of stereo image pairs were collected during the descent. The digital terrain models constructed from those images show rugged topography, in places approaching the angle of repose, adjacent to flatter darker plains. Brighter regions north of the landing site display two styles of drainage patterns: (1) bright highlands with rough topography and deeply incised branching dendritic drainage networks (up to fourth order) with dark-floored valleys that are suggestive of erosion by methane rainfall and (2) short, stubby low-order drainages that follow linear fault patterns forming canyon-like features suggestive of methane spring-sapping. The topographic data show that the bright highland terrains are extremely rugged; slopes of order of 30° appear common. These systems drain into adjacent relatively flat, dark lowland terrains. A stereo model for part of the dark plains region to the east of the landing site suggests surface scour across this plain flowing from west to east leaving ∼100-m-high bright ridges. Tectonic patterns are evident in (1) controlling the rectilinear, low-order, stubby drainages and (2) the “coastline” at the highland–lowland boundary with numerous straight and angular margins. In addition to flow from the highlands drainages, the lowland area shows evidence for more prolific flow parallel to the highland–lowland boundary leaving bright outliers resembling terrestrial sandbars. This implies major west to east floods across the plains where the probe landed with flow parallel to the highland–lowland boundary; the primary source of these flows is evidently not the dendritic channels in the bright highlands to the north.     

Models of high velocity impacts into dust-covered ice: Application to Martian northern lowlands

The detection of fresh impact craters with bright floors and ejecta (arising from fresh clean water ice) in the northern lowlands of Mars (Byrne et al., 2009b, Science 325, 1674), together with observations of polygonal structures and evidence from the Phoenix probe, suggests that there are substantial water ice deposits just below the surface over large areas. Specifically in cases of the larger craters observed, the impacts themselves may have raised the temperature and the pressure of the water ice deposits locally to values which allow phase changes. In this paper, we use smoothed particle hydrodynamics to model hyper-velocity impacts. We estimate peak shock pressures in a solid water ice target covered by a layer of loose material, modeled by pre-damaged dunite. In addition, we account for the possibility of a thin layer of sub-surface water ice by using a three-layer model where the ice is surrounded by dunite. We find that the peak shock pressures reached in the simulated events are high enough to produce several 100-1000 kg of liquid water depending upon the impact parameters and the exact shock pressure needed for the phase change. A difficulty remains however in determining whether liquid is generated or whether a type of fluidized ice is produced (or indeed some combination of the two). We also note that the process can become rather complex as the number of layers increases because of reflections of the shock at sub-surface boundaries—a process which should lead to increased fluidization.     

Regolith transport in craters on Eros

Images of Eros from the NEAR Shoemaker spacecraft reveal bright and dark albedo features on steep crater walls unlike markings previously observed on asteroids. These features have been attributed to the downslope movement of space-weathered regolith, exposing less weathered material (Science 292 (2001) 484; Meteor. Planet. Sci. 36 (2001) 1617; Icarus 155 (2002) 145). Here we present observations of the interiors of large craters (>1 km in diameter) to test this hypothesis and constrain the origin of the features. We find that bright regions in these craters correspond to steep slopes, consistent with previous work. The geographic distribution of craters with albedo variations shows no pattern and does not resemble the distribution of ponds, another phenomenon on Eros attributed to regolith movement. Shadows and other indications of topography are not observed at feature boundaries, implying that the transported layer is ?1 m thick. The presence of multiple bright and dark units on long slopes with sharp boundaries between them suggests that mobilized regolith may be halted by frictional or other effects before reaching the foot of the slope. Features on crater walls should darken at the same rate as bright ejecta deposits from crater formation; the lack of observed, morphologically fresh craters with bright interiors or ejecta suggests that the albedo patterns are younger than the most recently formed craters greater than about 100 m in diameter. Smaller or micrometeorite impacts, which would not necessarily leave evident deposits of bright ejecta, remain possible causes of albedo patterns. Although their effectiveness is difficult to assess, electrostatic processes and thermal creep are also candidates.     

The production of Ganymede's magnetic field

One of the great discoveries of NASA's Galileo mission was the presence of an intrinsically produced magnetic field at Ganymede. Generation of the relatively strong (750 nT) field likely requires dynamo action in Ganymede's metallic core, but how such a dynamo has been maintained into the present epoch remains uncertain. Using a one-dimensional, three layer thermal model of Ganymede, we find that magnetic field generation can only occur if the sulfur mass fraction in Ganymede's core is very low (?3%) or very high (?21%), and the silicate mantle can cool rapidly (i.e. it has a viscosity like wet olivine). However, these requirements are not necessarily compatible with cosmochemical and physical models of the satellite. We therefore investigate an alternative scenario for producing Ganymede's magnetic field in which passage through an eccentricity pumping Laplace-like resonance in Ganymede's past enables present day dynamo action in the metallic core. If sufficient tidal dissipation occurs in Ganymede's silicate mantle during resonance passage, silicate temperatures can undergo a runaway which prevents the core from cooling until the resonance passage ends. The rapid silicate and core cooling that follows resonance escape triggers dynamo action via thermal and/or compositional convection. To test the feasibility of this mechanism we couple our thermal model with an orbital evolution model to examine the effects of resonance passage on Ganymede's silicate mantle and metallic core. We find that, contrary to expectations, there are no physically plausible scenarios in which tidal heating in the silicates is sufficient to cause the thermal runaway necessary to prevent core cooling. These findings are robust to variations in the silicate rheology, tidal dissipation factor of Jupiter (QJ), structure of the ice shell, and the inclusion of partial melting in the silicate mantle. Resonance passage therefore appears unlikely to explain Ganymede's magnetic field and we must appeal to the special conditions described above to explain the presence of the field.     

The Circum-Hellas Volcanic Province,Mars: Overview

Building on previous studies of volcanoes around the Hellas basin with new studies of imaging (High-Resolution Stereo Camera (HRSC), Thermal Emission Imaging System (THEMIS), Mars Orbiter Camera (MOC), High-Resolution Imaging Science Experiment (HiRISE), Context Imager (CTX)), multispectral (HRSC, Observatoire pour la Minéralogie, l’Eau, les Glaces et l’Activité (OMEGA)), topographic (Mars Orbiter Laser Altimeter (MOLA)) and gravity data, we define a new Martian volcanic province as the Circum-Hellas Volcanic Province (CHVP). With an area of >2.1 million km², it contains the six oldest central vent volcanoes on Mars, which formed after the Hellas impact basin, between 4.0 and 3.6 Ga. These volcanoes mark a transition from the flood volcanism that formed Malea Planum ～3.8 Ga, to localized edifice-building eruptions. The CHVP volcanoes have two general morphologies: (1) shield-like edifices (Tyrrhena, Hadriaca, and Amphitrites Paterae), and (2) caldera-like depressions surrounded by ridged plains (Peneus, Malea, and Pityusa Paterae). Positive gravity anomalies are found at Tyrrhena, Hadriaca, and Amphitrites, perhaps indicative of dense magma bodies below the surface. The lack of positive-relief edifices and weak gravity anomalies at Peneus, Malea, and Pityusa suggest a fundamental difference in their formation, styles of eruption, and/or compositions. The northernmost volcanoes, the ～3.7–3.9 Ga Tyrrhena and Hadriaca Paterae, have low slopes, well-channeled flanks, and smooth caldera floors (at tens of meters/pixel scale), indicative of volcanoes formed from poorly consolidated pyroclastic deposits that have been modified by fluvial and aeolian erosion and deposition. The ～3.6 Ga Amphitrites Patera also has a well-channeled flank, but it and the ～3.8 Ga Peneus Patera are dominated by scalloped and pitted terrain, pedestal and ejecta flow craters, and a general ‘softened’ appearance. This morphology is indicative not only of surface materials subjected to periglacial processes involving water ice, but also of a surface composed of easily eroded materials such as ash and dust. The southernmost volcanoes, the ～3.8 Ga Malea and Pityusa Paterae, have no channeled flanks, no scalloped and pitted terrain, and lack the ‘softened’ appearance of their surfaces, but they do contain pedestal and ejecta flow craters and large, smooth, bright plateaus in their central depressions. This morphology is indicative of a surface with not only a high water ice content, but also a more consolidated material that is less susceptible to degradation (relative to the other four volcanoes). We suggest that Malea and Pityusa (and possibly Peneus) Paterae are Martian equivalents to Earth's giant calderas (e.g., Yellowstone, Long Valley) that erupted large volumes of volcanic materials, and that Malea and Pityusa are probably composed of either lava flows or ignimbrites. HRSC and OMEGA spectral data indicate that dark gray to slightly red materials (often represented as blue or black pixels in HRSC color images), found in the patera floors and topographic lows throughout the CHVP, have a basaltic composition. A key issue is whether this dark material represents concentrations of underlying basaltic material eroded by various processes and exposed by aeolian winnowing, or if the material was transported from elsewhere on Mars by regional winds. Understanding the provenance of these dark materials may be the key to understanding the volcanic diversity of the Circum-Hellas Volcanic Province.     

Mars Exploration Rover candidate landing sites as viewed by THEMIS

The analysis of six landing sites that were candidates for the two NASA Mars Exploration Rovers (MER) benefited from recently available image data from the Thermal Emission Imaging Spectrometer (THEMIS) onboard the 2001 Mars Odyssey spacecraft. The combination of daytime and nighttime thermal infrared images from THEMIS supplemented by additional data sets has lead to new or expanded insights into the nature of each landing site. In Meridiani Planum, a layer of lighter-toned, higher thermal inertia material is observable just below the hematite-bearing layer. Gusev Crater displays a more complex stratigraphy than previously observed, including an upper layer with lobate margins. The highest inertia unit of southern Isidis Planitia is confined to topographic lows in the rim/basin margin and does not appear to be due to highland material transported onto the basin floor. The enigmatic, ovoid, blocky terrain on the floor of Melas Chasma displays higher thermal inertia than its surroundings, an indication that it contains coarser or more indurated material than the adjacent aeolian bedforms. The myriad channels of Athabasca Valles display distinctive thermal signatures despite the presence of a bright layer of dust covering the region. The presence of alluvial fans produced from spur-and-gulley erosion of the walls of Eos Chasma demonstrates that mass movements have occurred following the canyon scouring floods.