Flood-formed dunes in Athabasca Valles, Mars: morphology, modeling, and implications

Estimates of discharge for martian outflow channels have spanned orders of magnitude due in part to uncertainties in floodwater height. A methodology of estimating discharge based on bedforms would reduce some of this uncertainty. Such a methodology based on the morphology and granulometry of flood-formed (‘diluvial’) dunes has been developed by Carling (1996b, in: Branson, J., Brown, A.G., Gregory, K.J. (Eds.), Global Continental Changes: The Context of Palaeohydrology. Geological Society Special Publication No. 115, London, UK, 165-179) and applied to Pleistocene flood-formed dunes in Siberia. Transverse periodic dune-like bedforms in Athabasca Valles, Mars, have previously been classified both as flood-formed dunes and as antidunes. Either interpretation is important, as they both imply substantial quantities of water, but each has different hydraulic implications. We undertook photoclinometric measurements of these forms, and compared them with data from flood-formed dunes in Siberia. Our analysis of those data shows their morphology to be more consistent with dunes than antidunes, thus providing the first documentation of flood-formed dunes on Mars. Other reasoning based on context and likely hydraulics also supports the bedforms' classification as dunes. Evidence does not support the dunes being aeolian, although a conclusive determination cannot be made with present data. Given the preponderance of evidence that the features are flood-formed instead of aeolian, we applied Carling's (1996b, in: Branson, J., Brown, A.G., Gregory, K.J. (Eds.), Global Continental Changes: The Context of Palaeohydrology. Geological Society Special Publication No. 115, London, UK, 165-179) dune-flow model to derive the peak discharge of the flood flow that formed them. The resultant estimate is approximately 2×10⁶ m³/s, similar to previous estimates. The size of the Athabascan dunes' in comparison with that of terrestrial dunes suggests that these martian dunes took at least 1-2 days to grow. Their flattened morphology implies that they were formed at high subcritical flow and that the flood flow that formed them receded very quickly.     

Depth profiles of venusian sinuous rilles and valley networks

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

The role of arcuate ridges and gullies in the degradation of craters in the Newton Basin region of Mars

A survey of craters in the vicinity of Newton Basin, using high-resolution images from Mars Global Surveyor and Mars Odyssey, was conducted to find and analyze examples of gullies and arcuate ridges and assess their implications for impact crater degradation processes. In the Phaethontis Quadrangle (MC-24), we identified 225 craters that contain these features. Of these, 188 had gullies on some portion of their walls, 118 had arcuate ridges at the bases of the crater walls, and 104 contained both features, typically on the same crater wall. A major result is that the pole-facing or equator-facing orientation of these features is latitude dependent. At latitudes >44° S, equator-facing orientations for both ridges and gullies are prevalent, but at latitudes <44° S, pole-facing orientations are prevalent. The gullies and arcuate ridges typically occupy craters between ∼2 and 30 km in diameter, at elevations between −1 and 3 km. Mars Orbiter Laser Altimeter (MOLA) elevation profiles indicate that most craters with pole-facing arcuate ridges have floors sloping downward from the pole-facing wall, and some of these craters show asymmetry in crater rim heights, with lower pole-facing rims. These patterns suggest viscous flow of ice-rich materials preferentially away from gullied crater walls. Clear associations exist between gullies and arcuate ridges, including (a) geometric congruence between alcoves and sinuous arcs of arcuate ridges and (b) backfilling of arcuate ridges by debris aprons associated with gully systems. Chronologic studies suggest that gullied walls and patterned crater floor deposits have ages corresponding to the last few high obliquity cycles. Our data appear consistent with the hypothesis that these features are associated with periods of ice deposition and subsequent erosion associated with obliquity excursions within the last few tens of millions of years. Arcuate ridges may form from cycles of activity that also involve gully formation, and the ridges may be in part due to mass-wasted, ice-rich material transported downslope from the alcoves, which then interacts with previously emplaced floor deposits. Most observed gullies may be late-stage features in a degradational cycle that may have occurred many times on a given crater wall.     

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

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

Embayed intermediate volcanoes on Venus: Implications for the evolution of the volcanic plains

Volcanoes on Venus are classified according to size with studies on the stratigraphic position of large volcanoes proposing that most of the large volcanoes postdate the regional volcanic materials. Some studies regarding intermediate volcanoes proposed that some of these volcanic features could be large volcanoes with embayed flow aprons, a situation that would alter the previous stratigraphic considerations about large volcanoes on Venus.In this work I analyze the global population of embayed intermediate-size volcanoes and compare their summits with that of other edifices classified as large volcanoes. Intermediate-size volcanoes are considered embayed when: (1) flows from another source clearly overlap the volcano slopes, and (2) display scarps related to flank-failure processes but with the associated collapse deposits being absent (i.e. interpreted as covered). As result of the survey 88 embayed intermediate-size volcanoes have been catalogued and integrated into a Geographic Information System. These embayed volcanoes have summit sizes and characteristics similar to large volcanoes and, therefore, could be interpreted as possible large volcanoes with their flow aprons embayed. Embayment materials for these volcanoes include all the units present in the history of the volcanic plains and would indicate that this type of central volcanic edifice would occur throughout the geologic history recorded in the venusian plains.     

Millochau crater, Mars: Infilling and erosion of an ancient highland impact crater

The geology and stratigraphy of Millochau crater (21.4° S, 275° W), located in the highlands of Tyrrhena Terra, Mars, are documented through geomorphic analyses and geologic mapping. Crater size-frequency distributions and superposition relationships are used to constrain relative ages of geologic units and determine the timing and duration of the geologic processes that modified Millochau rim materials and emplaced deposits on Millochau's floor. Crater size-frequency distributions show a Middle Noachian age for rim materials and Middle Noachian to Early Hesperian ages for most of the interior deposits. Valley networks and gullies incised within Millochau's rim materials and interior wall, respectively, indicate fluvial activity was an important erosional process. Millochau contains an interior plateau, offset northeast of Millochau's center, which rises up to 400 m above the surrounding crater floor and slopes downward to the south and west. Layers exposed along the northern and eastern scarp boundaries of the plateau are tens to hundreds of meters thick and laterally continuous in MOC images. These layers suggest most materials within Millochau were emplaced by sedimentary processes (e.g., fluvial or eolian), with the potential for lacustrine deposition in shallow transient bodies of water and contributions of volcanic airfall. Mass wasting may have also contributed significant quantities of material to Millochau's interior, especially to the deposits surrounding the plateau. Superposition relationships combined with impact crater statistics indicate that most deposition and erosion of Millochau's interior deposits is ancient, which implies that fluvial activity in this part of Tyrrhena Terra is much older than in the eastern Hellas region. Eolian processes mobilized sediment to form complicated patterns of long- and short-wavelength dunes, whose emplacement is controlled by local topography. These deposits are some of the youngest within Millochau (Amazonian) and eolian modification may be ongoing.     

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.     

The structural control of venusian polygonal impact craters

Pre-impact inhomogeneities of the target material sometimes cause the rim of an impact crater to be composed of several straight segments, instead of being circular. The venusian surface hosts 121 such polygonal impact craters (PICs)>12 km in diameter. Their straight rim segments are often parallel to the orientations of the surrounding tectonic structures, particularly those in tessera terrain and young rift zones, as well as the concentric components of coronae. This match is notably more distinct in distances less than two crater diameters between the PICs and the tectonic structures than further away. Surrounding wrinkle ridges, lineaments or radial components of volcano-tectonic features seem to have very little influence on the orientations of PIC rims. These results imply that the formation of straight segments of venusian PIC rims is controlled by pre-existing tectonic structures of the crust, but not by the apparently most surficial ones. Thus, PICs could be used to provide further constraints on the distribution and orientations of zones of weakness in the venusian crust.     

Beta Regio, Venus: Evidence for uplift, rifting, and volcanism due to a mantle plume

Geophysical data have led to the interpretation that Beta Regio, a 2000×25000 km wide topographic rise with associated rifting and volcanism, formed due to the rise of a hot mantle diapir interpreted to be caused by a mantle plume. We have tested this hypothesis through detailed geologic mapping of the V-17 quadrangle, which includes a significant part of the Beta Regio rise, and reconnaissance mapping of the remaining parts of this region. Our analysis documents signatures of an early stage of uplift in the formation of the Agrona Linea fracture belts before the emplacement of regional plains and their deformation by wrinkle ridging. We see evidence that the Theia rift-associated volcanism occurred during the first part of post-regional-plains time and cannot exclude that it continued into later time. We also see evidence that Devana Chasma rifting was active during the first and the second parts of post-regional-plains time. These data are consistent with uplift, rifting and volcanism associated with a mantle diapir. Geophysical modeling shows that diapiric upwelling may continue at the present time. Together these data suggest that the duration of mantle diapir activity was as long as several hundred million years. The regional plains north of Beta rise and the area east and west of it were little affected by the Beta-forming plume, but the broader area (at least 4000 km across), whose center-northern part includes Beta Regio, could have experienced earlier uplift as morphologically recorded in formation of tessera transitional terrain.     

Tectonic and stratigraphic implications of the relative ages of venusian plains and wrinkle ridges

A major ongoing controversy concerns the style of crustal evolution on Venus. At one extreme is a directional model that proposes a sequence of depositional and deformational events that occur at specific times in the evolution of the crust and that are global in extent. At the other extreme is a model that argues for different ages of these events in different places on the planet. A test of the directional model is here focused on whether wrinkle ridges formed at a single time in the recorded crustal history of Venus. Where sets of wrinkle ridges intersect it commonly is possible to determine that one set is older than the other. Also, the deformation responsible for wrinkle ridges is, in places, clearly progressive with respect to stratigraphic material units. These observations are not consistent with a specific single time for the formation of wrinkle ridges within the stratigraphic sequence. Within an area including about 1/3 of the surface of Venus 15% of craters that are younger than regional plains are older than wrinkle ridges, 85% are younger than wrinkle ridges. Taking 750 myr as a reasonable mean age for the regional plains, this implies that the mean age of wrinkle ridges is ∼110 myr younger than the mean age of plains. Solomon et al. (1999, Science 286, 87) propose that the emplacement of a large volume of plains lava would lead to a major atmospheric temperature increase. Their model predicts thermal stresses in the lithosphere that, at shallow depth, would reach peak compressive stresses in about 100 myr, a number very similar to the time lag between plains emplacement and wrinkle ridge formation indicated by the crater data. The thermal compressive stresses responsible for wrinkle ridges would be maintained at a level sufficient to deform basalt for at least 100 myr and possibly for as long as 350 myr. These time intervals are not really short compared to the mean age of the plains. Finally, because wrinkle ridges are demonstrably younger than the plains they deform, they cannot be related to the processes that formed the plains and thus should not be used to define a “plains with wrinkle ridges” unit.     

Longevity of fluorine-bearing tremolite on Venus

There is a general belief that hydrous minerals cannot exist on Venus under current surface conditions. This view was challenged when Johnson and Fegley (2000, Icarus 146, 301-306) showed that tremolite (Ca₂Mg₅Si₈O₂₂(OH)₂), a hydrous mineral, is stable against thermal decomposition at current Venus surface temperatures, e.g., 50% decomposition in 4 Ga at 740 K. To further explore hydrous mineral thermal stability on Venus, we experimentally determined the thermal decomposition kinetics of fluorine-bearing tremolite. Fluor-tremolite is thermodynamically more stable than OH-tremolite and should decompose more slowly. However how much slower was unknown. We measured the decomposition rate of fluorine-bearing tremolite and show that its decomposition is several times to greater than ten times slower than that of OH-tremolite. We also show that F-bearing tremolite is depleted in fluorine after decomposition and that fluorine is lost as a volatile species such as HF gas. If tremolite ever formed on Venus, it would probably also contain fluorine. The exceptional stability of F-bearing tremolite strengthens our conclusions that if hydrous minerals ever formed on Venus, they could still be there today.     

Rheological constraints on martian landslides

We use a dynamic finite-difference model to simulate martian landslides in the Valles Marineris canyon system and Olympus Mons aureole using three different modal rheologies: frictional, Bingham, and power law. The frictional and Bingham modes are applied individually. Fluidized rheology is treated as a combination of frictional and power-law modes; general fluidization can include pore pressure contributions, whereas acoustic fluidization does not. We find that general fluidization most often produces slides that best match landslide geometry in the Valles Marineris. This implies that some amount of supporting liquid or gas was present in the material during failure. The profile of the Olympus Mons aureole is not well matched by any landslide model, suggesting an alternative genesis. In contrast, acoustic fluidization produces the best match for a lunar slide, a result anticipated for dry crust with no overlying atmosphere. The presence of pressurized fluid during Valles Marineris landsliding may be due to liquid water beneath a thin cryosphere (<1-2 km) or flash sublimation of CO₂.     

An assessment of evidence for pingos on Mars using HiRISE

Pingos are small hills with cores of ice, formed by injection and freezing of pressurized water. The possibility of pingos on Mars is of particular interest because of the associated implications for liquid water. We have systematically searched for candidate pingos using images from the High Resolution Imaging Science Experiment (HiRISE) camera. Since pingos are expected to develop surface fractures due to extension of the frozen ground over the ice core, we have searched for fractured features and identified a variety of mounds. These features are confined to the martian mid-latitudes, in the bands where gullies are also most common. The observed fractured mounds have a variety of morphologies and are likely of multiple origins. Isolated fractured mounds found on the floors of gullied craters in the southern hemisphere match the general morphologic characteristics of terrestrial pingos and are the best candidates for martian pingos, but there is currently no direct evidence for presence of ice cores and it is difficult to produce the necessary water volumes, so these features should still be interpreted with caution. Other fractured mounds appear more likely to be erosional remnants of an unusual mantling layer or possibly thermokarst structures. Flat-topped mounds in Utopia have some characteristics (fracture pattern and latitudinal distribution) consistent with pingos but differ in other aspects such as shape and setting. While we do not rule out a pingo origin, we prefer an erosional model for these enigmatic features.     

A revised model for cycloid growth mechanics on Europa: Evidence from surface morphologies and geometries

Although a single model currently exists to explain the development of curved Europan cycloids, there have been no systematic studies of the range of morphologies and quantifiable geometric parameters of cycloidal features. We address variations in geometry along individual cycloid segments, characterizing differences in cusp styles and angles, and addressing the morphologic aspects of cycloid segments and cusps. In so doing, we illustrate how geometric and morphologic evidence imply a formation mechanism that differs from the existing model in several aspects. The current model states that cycloids are initiated as tensile fractures that grow in a curved path in response to rotating diurnal tidal stresses on Europa. However, the geometry of a cycloid cusp necessitates that shear stress was resolved onto the existing cycloid segment by the rotating diurnal stresses at the instant of cusp formation. Furthermore, we observe that cycloid cusps have a strikingly similar geometry to tailcracks that developed at the tips of many ridge-like strike-slip faults on Europa in response to shearing at the fault tip. We suggest that this similarity in geometries can be attributed to an identical formation mechanism whereby cycloid cusps form by a tailcracking process. We therefore present a revised, mechanically-based model for cycloid formation that retains the basic premise that crack growth is governed by diurnal stresses, but describes the development of cycloid cusps in response to resolved shear stresses at the tips of existing cycloid segments. The ratio of normal to shear stress at the time of tailcrack formation dictates the cusp angle and, over longer time periods, influences the morphologic evolution of the cycloid segment as it is repeatedly reworked by tidal stresses.     

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

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

Tharsis Tholus: an unusual martian volcano

Tharsis Tholus is unusual martian shield volcano in that the edifice is cut by a series of large normal faults that appear to penetrate the entire volcano. Northeast-trending narrow graben also cut the flank. The large normal faults may be caused by loading of a ductile subsurface layer allowing failure of the edifice; the narrow graben are typical tensional faults. The flank is heavily mantled by aeolian material. Despite the bulbous appearance, the overall morphology of Tharsis Tholus suggests it is a basaltic shield. Crater counts indicate an age of early Hesperian placing Tharsis Tholus in the middle of the period of activity that built the other small Tharsis volcanoes.     

Modeling crater populations on Venus and Titan

We describe a model for crater populations on planets and satellites with dense atmospheres, like those of Venus and Titan. The model takes into account ablation (or mass shedding), pancaking, and fragmentation. Fragmentation is assumed to occur due to the hydrodynamic instabilities promoted by the impactors’ deceleration in the atmosphere. Fragments that survive to hit the ground make craters or groups thereof. Crater sizes are estimated using standard laws in the gravity regime, modified to take into account impactor disruption. We use Monte Carlo methods to pick parameters from appropriate distributions of impactor mass, zenith angle, and velocity. Good fits to the Venus crater populations (including multiple crater fields) can be found with reasonable values of model parameters. An important aspect of the model is that it reproduces the dearth of small craters on Venus: this is due to a cutoff on crater formation we impose, when the expected crater would be smaller than the (dispersed) object that would make it. Hydrodynamic effects alone (ablation, pancaking, fragmentation) due to the passage of impactors through the atmosphere are insufficient to explain the lack of small craters. In our favored model, the observed number of craters (940) is produced by ∼5500 impactors with masses , yielding an age of (1-σ uncertainty) for the venusian surface. This figure does not take into account any uncertainties in crater scaling and impactor population characteristics, which probably increase the uncertainty to a factor of two in age.We apply the model with the same parameter values to Titan to predict crater populations under differing assumptions of impactor populations that reflect present conditions. We assume that the impactors (comets) are made of 50% porous ice. Predicted crater production rates are ≈190 craters . The smallest craters on Titan are predicted to be in diameter, and ≈5 crater fields are expected. If the impactors are composed of solid ice (density ), crater production rates increase by ≈70% and the smallest crater is predicted to be in diameter. We give cratering rates for denser comets and atmospheres 0.1 and 10 times as thick as Titan's current atmosphere. We also explicitly address leading-trailing hemisphere asymmetries that might be seen if Titan's rotation rate were strictly synchronous over astronomical timescales: if that is the case, the ratio of crater production on the leading hemisphere to that on the trailing hemisphere is ≈4:1.     

Boulders and ponds on the Asteroid 433 Eros

There are ∼300 features on the Asteroid 433 Eros that morphologically resemble ponds (flat-floored and sharply embaying the bounding depression in which they sit). Because boulders on Eros are apparently eroding in place and because ponds with associated boulders tend to be larger than ponds without blocks, we propose that ponds form from thermally disaggregated and seismically flattened boulder material, under the assumption that repeated day/night cycling causes material fatigue that leads to erosion of the boulders. Results from a simple boulder emplacement/thermal erosion model with boulders emplaced in a few discrete events (i.e., large impacts) match well the observed size distribution. Under this scenario, the subtle color differences of ponds (somewhat bluer than the rest of the surface) might be due to some combination of less space-weathered material and density stratification of silicate-rich chondrules and more metal-rich matrix from a disaggregated boulder. Volume estimates of ponds derived from NEAR Laser Rangefinder profiles are consistent with what can be supplied by boulders. Ponds are also observed to be concentrated in regions of low slope and high elevation, which suggests the presence of a less mobile regolith and thus a contrast in the resistance to seismic shaking between the pond material and the material that makes up the bounding depression. Future tests include shake-table experiments and temperature cycling (fatigue) of ordinary chondrites to test the thermal erosion mechanism.     

Layering stratigraphy of eastern Coprates and northern Capri Chasmata, Mars

Distinct competent layers are observed in the slopes of eastern Coprates Chasma, part of the Valles Marineris system on Mars. Our observations indicate that the stratigraphy of Coprates Chasma consists of alternating thin strong layers and thicker sequences of relatively weak layers. The strong, competent layers maintain steeper slopes and play a major role in controlling the overall shape and geomorphology of the chasmata slopes. The topmost competent layer in this area is well preserved and easy to identify in outcrops on the northern rim of Coprates Chasma less than 100 m below the southern Ophir Planum surface. The volume of the topmost emplaced layer is at least 70 km³ and may be greater than 2100 km³ if the unit underlies most of Ophir Planum. The broad extent of this layer allows us to measure elevation offsets within the north rim of the chasma and in a freestanding massif within Coprates Chasma where the layer is also observed. Rim outcrop morphology and elevation differences between Ophir and Aurorae Plana may be indicative of the easternmost extent of the topmost competent layer. These observations allow an insight into the depositional processes that formed the stratigraphic stack into which this portion of the Valles Marineris is carved, and they present a picture of some of the last volcanic activity in this area. Furthermore, the elevation offsets within the layer are evidence of significant subsidence of the massif and surrounding material.