Depth of faulting and ancient heat flows in the Kuiper region of Mercury from lobate scarp topography

Mercurian lobate scarps are interpreted to be the surface expressions of thrust faults formed by planetary cooling and contraction, which deformed the crust down to the brittle–ductile transition (BDT) depth at the time of faulting. In this work we have used a forward modeling procedure in order to analyze the relation between scarp topography and fault geometries and depths associated with a group of prominent lobate scarps (Santa Maria Rupes and two unnamed scarps) located in the Kuiper region of Mercury for which Earth-based radar altimetry is available. Also a backthrust associated with one of the lobate scarps has been included in this study. We have obtained best fits for depths of faulting between 30 and 39 km; the results are consistent with the previous results for other lobate scarps on Mercury.The so-derived fault depths have been used to calculate surface heat flows for the time of faulting, taking into account crustal heat sources and a heterogeneous surface temperature due to the variable insolation pattern. Deduced surface heat flows are between 19 and 39 mW m^?2 for the Kuiper region, and between 22 and 43 mW m^?2 for Discovery Rupes. Both BDT depths and heat flows are consistent with the predictions of thermal history models for the range of time relevant for scarp formation.     

Accommodation of lithospheric shortening on Mercury from altimetric profiles of ridges and lobate scarps measured during MESSENGER flybys 1 and 2

The Mercury Laser Altimeter on the NASA MESSENGER mission has ranged to several ridges and lobate scarps during two equatorial flybys of the planet Mercury. The tectonic features sampled, like others documented by spacecraft imaging and Earth-based radar, are spatially isolated and have vertical relief in excess of 1 km. The profiles also indicate that the faulting associated with their formation penetrated to tens of kilometers depth into the lithosphere and accommodated substantial shortening. To gain insight into the mechanism(s) of strain accommodation across these structures, we perform analytical and numerical modeling of representative dynamic localization mechanisms. We find that ductile localization due to shear heating is not favored, given our current understanding of thermal gradients and shallow thermal structure of Mercury at the time of ridge and scarp formation, and is likely to be of secondary importance at best. Brittle localization, associated with loss of resistance during fault development or with velocity weakening during sliding on mature faults, is weakly localizing but permits slip to accumulate over geological time scales. The range of shallow thermal gradients that produce isolated faults rather than distributed fault sets under the assumption of modest fault weakening is consistent with previous models for Mercury’s early global thermal history. To be consistent with strain rates predicted from thermal history models and the amount of shortening required to account for the underlying large-offset faults, ridges and scarps on Mercury likely developed over geologically substantial time spans.     

Pleistocene northwards fold propagation of the Jura within the southern Upper Rhine Graben: seismotectonic implications

On the basis of the geophysical (seismic profiles and electric tomography), geomorphic and geological data, we re-evaluate the post-Pliocene structural interpretation of the southern Upper Rhine graben (Basel–Mulhouse area): we demonstrate a Plio–Pleistocene northward propagation of the Jura thrust and fold belt up to Mulhouse proceeding from a succession of four 10 km apart ramps (from north to south Ferrette, Muespach, Magstatt and Rixheim) rooted within the late Triassic evaporitic marls acting as a decollement. This domain was previously considered as having undergone an on-going continuous extension (horst of Mulhouse bounded by the Quaternary Sirentz and Dannemarie grabens).The Quaternary activity of this thin-skinned tectonics induces the growth of a sedimentary wedge whose regional slope, which comprised between 1.4° and 1° to the north, also attests to a low friction basal detachment. More into details, these ramps correspond to 40–50-m high jumps within the forward topographic slope. Pleistocene activity is suggested just above the Muespach ramp by the presence of a 5–10-m north-facing scarp corresponding in depth to a 3-m vertical offset of early Pleistocene alluvial deposits. Farther to the north, a stronger incision of the Rhine Würm terrace can be interpreted as the result of the growth of the Mulhouse–Rixheim frontal ramp.This northward propagation of the Jura thrust and fold belt is strongly controlled by the Oligocene structural inheritage. The development of the frontal ramp in Mulhouse has to be related to the Oligocene significant vertical offset of the Triassic evaporite along the Mulhouse Railway Station fault preventing a propagation of the decollement farther to the north. In the same way, the fold propagation is laterally segmented by the N20°E trending Oligocene fabrics (from East to West, Rhine Valley flexure fault, Allschwil–Istein fault system and Illfurth fault) which acts above the decollement as lateral ramps. To the west, the development of a shallow anticline along the Illfurth fault suggests that the thin-skinned propagation is oblique with respect to the Oligocene fabrics. It results in spacial contrast between a left-lateral-reverse and a right-lateral–normal shallow kinematics along the western and eastern lateral ramps, respectively. In depth to the east, it also induces a vertical contrast between shallow (right-lateral–normal) and deep (left-lateral given by fault plane solutions) kinematics along the Istein–Allschwill–Rhine Valley fault system.Few arguments supporting a nucleation of the Basel-1356 earthquake, the strongest event in NW Europe in the last thousand years, onto the Rhine Valley fault system beneath the decollement have been given. However, we emphasize that the above mentioned coeval thin (aseismic)- and thick (seismic)-skinned tectonics along the Istein–Allschwill–Rhine Valley fault system would make difficult both the identification and the interpretation of the surface rupture of the Basel-1356 earthquake.     

Structural control of scarps in the Rembrandt region of Mercury

Lobate scarps, thought to be the surface expression of large thrust faults, are the most spectacular contractional tectonic features visible on Mercury. Most lobate scarps follow a general and relatively simple pattern, with a roughly arcuate or linear form in plan view, and an asymmetric cross section characterized by a steeply rising scarp face and a gently declining back scarp. In this work, we study two peculiar and complex scarps in the Rembrandt region of Mercury through MESSENGER imagery. On the one hand, the formation of these scarps resulted in the deformation of features such as impact craters, fractures, extensional faults, and volcanic plains, while on the other hand, the deformed features partly influenced the formation of the scarps. Evidence for structural control on the formation of the scarps includes their orientation, segmentation, bifurcation, change in structural trend and dip orientation, and transition into high-relief ridges or wrinkle ridge morphologies in some cases. Thus, these two lobate scarps provide examples of complex geological relations among other features, expanding the recognized richness of mercurian geology. Also, the southern scarp records a complex history of contraction, suggesting that the development of some mercurian lobate scarps may be more complex than usually thought.     

The equatorial shape and gravity field of Mercury from MESSENGER flybys 1 and 2

On 14 January and 6 October 2008 the MESSENGER spacecraft passed within 200 km of the surface of Mercury. These flybys by MESSENGER provided the first observations of Mercury from a spacecraft since the Mariner 10 flybys in 1974 and 1975. Data from the Mercury Laser Altimeter (MLA) provided new information on the equatorial shape of Mercury, and Doppler tracking of the spacecraft through the flybys provided new data on the planet’s gravity field. The MLA passes were on opposite hemispheres of the planet and span collectively ∼40% of the equatorial circumference. The mean elevation of topography observed during flyby 1, in the longitude range 0-90°E, is greater than that seen during flyby 2 in the longitude range 180-270°E, indicating an offset between centers of mass and figure having a magnitude and phase in general agreement with topography determined by Earth-based radar. Both MLA profiles are characterized by slopes of ∼0.015° downward to the east, which is consistent with a long-wavelength equatorial shape defined by a best-fitting ellipse. The Doppler tracking data show sensitivity to the gravitational structure of Mercury. The equatorial ellipticity of the gravitational field, C_2,2, is well determined and correlates with the equatorial shape. The S_2,2 coefficient is ∼0, as would be expected if Mercury’s coordinate system, defined by its rotational state, is aligned along its principal axes of inertia. The recovered value of the polar flattening of the gravitational potential, J₂, is considerably lower in magnitude than the value obtained from Mariner 10 tracking, a result that is problematic for internal structure models. This parameter is not as well constrained as the equatorial ellipticity because the flyby trajectories were nearly in the planet’s equatorial plane. The residuals from the Doppler tracking data suggest the possibility of mascons on Mercury, but flyby observations are of insufficient resolution for confident recovery. For a range of assumptions on degree of compensation and crustal and mantle densities, the allowable crustal thickness is consistent with the upper limit of about 100 km estimated from the inferred depth of faulting beneath a prominent lobate scarp, an assumed ductile flow law for crustal material, and the condition that temperature at the base of the crust does not exceed the solidus temperature. The MESSENGER value of C_2,2 has allowed an improved estimate of the ratio of the polar moment of inertia of the mantle and crust to the full polar moment (C_m/C), a refinement that strengthens the conclusion that Mercury has at present a fluid outer core.     

Magnitude of global contraction on Mars from analysis of surface faults: Implications for martian thermal history

Faults provide a record of a planet’s crustal stress state and interior dynamics, including volumetric changes related to long-term cooling. Previous work has suggested that Mars experienced a pulse of large-scale global contraction during Hesperian time. Here we evaluate the evidence for martian global contraction using a recent compilation of thrust faults. Fault-related strains were calculated for wrinkle ridges and lobate scarps to provide lower and upper bounds, respectively, on the magnitude of global contraction from contractional structures observed on the surface of Mars. During the hypothesized pulse of global contraction, contractional strain of −0.007% to −0.13% is indicated by the structures, corresponding to decreases in planetary radius of 112 m to 2.24 km, respectively. By contrast, consideration of all recognized thrust faults regardless of age produces a globally averaged contractional strain of −0.011% to −0.22%, corresponding to a radius decrease of 188 m to 3.77 km since the Early Noachian. The amount of global contraction predicted by thermal models is larger than what is recorded by the faults at the surface, paralleling similar studies for Mercury and the Moon, which suggests that observations of fault populations at the surface may provide tighter bounds on planetary thermal evolution than models alone.     

Despinning plus global contraction and the orientation of lobate scarps on Mercury: Predictions for MESSENGER

Observations by the Mariner 10 spacecraft suggest that the lobate scarps on Mercury, which have been interpreted to record at most 1-2 km of radial contraction of the planet after the end of the Late Heavy Bombardment, possess a global, preferred N-S orientation but lack a strong latitudinal dependence on their surface expression. Here, we reexamine the idea that a decrease in the planetary rotation rate (despinning) coupled with global contraction of at least 3-5.5 km prior to the end of Late Heavy Bombardment resulted in global N-S oriented thrust faults. The surface expression of these faults is assumed to have been erased by the end of the Late Heavy Bombardment, and the faults were subsequently reactivated by later global contraction, producing generally N-S oriented thrust faults from an isotropic stress field. We use the estimate of >3-5.5 km contraction prior to ∼4 Ga as an additional constraint to thermomechanical simulations of the evolution of Mercury, finding that a wide range of models are consistent with this observation. The fact that a wide range of states are consistent with the contraction of Mercury prior to the end of Late Heavy Bombardment but only a restricted set of states are consistent with the at most 1-2 km of subsequent contraction bolsters the idea that there may be hidden strain on Mercury, features unseen by Mariner 10 but likely visible to the MESSENGER spacecraft.     

Ancient heat flow and crustal thickness at Warrego rise, Thaumasia highlands, Mars: Implications for a stratified crust

Heat flow calculations based on geological and/or geophysical indicators can help to constrain the thickness, and potentially the geochemical stratification, of the martian crust. Here we analyze the Warrego rise region, part of the ancient mountain range referred to as the Thaumasia highlands. This region has a crustal thickness much greater than the martian average, as well as estimations of the depth to the brittle-ductile transition beneath two scarps interpreted to be thrust faults. For the local crustal density (2900 kg m⁻³) favored by our analysis of the flexural state of compensation of the local topography, the crustal thickness is at least 70 and 75 km at the scarp locations. However, for one of the scarp locations our nominal model does not obtain heat flow solutions permitting a homogeneous crust as thick as required. Our results, therefore, suggest that the crust beneath the Warrego rise region is chemically stratified with a heat-producing enriched upper layer thinner than the whole crust. Moreover, if the mantle heat flow (at the time of scarp formation) was higher than 0.3 of the surface heat low, as predicted by thermal history models, then a stratified crust rise seems unavoidable for this region, even if local heat-producing element abundances lower than average or hydrostatic pore pressure are considered. This finding is consistent with a complex geological history, which includes magmatic-driven activity.     

Rümker hills: a lunar volcanic dome complex

The Rümker Hills, a volcanic dome-flow complex in the northern Oceanus Procellarum, is characterized by overlapping plains-forming units with lobate scarps, volcanic domes, a 60 km ring, and a scarp which separates the plateau from surrounding mare materials. Plains-forming units are interpreted as fluid volcanic flows, and domes as viscous extrusions. One dome may be a stratovolcano. The ring system is discordant with regional structural trends and probably has a local origin. The Rümker Hills is the closest lunar analog to the large martian shield structures revealed on the Mariner 9 photographs of Mars.     

Stereo topographic models of Mercury after three MESSENGER flybys

From photogrammetric analysis of stereo images of Mercury obtained during three MESSENGER flybys, we have produced three digital terrain models (DTMs) that have a grid spacing of 1 km and together cover 30% of the planet's surface. The terrain models provide a rich source of information on the morphology of Mercury's surface, including details of tectonic scarp systems as well as impact craters and basins. More than 400 craters larger than 15 km in diameter are included in the models. Additionally, the models provide important test cases for the analysis of stereo image data to be collected during MESSENGER's orbital mission phase. Small lateral offsets and differences in trends between stereo DTMs and laser altimeter profiles may be due to remaining errors in spacecraft position, instrument pointing, or Mercury coordinate knowledge. Such errors should be resolved during the orbital mission phase, when more joint analyses of data and detailed orbit modeling will be possible.     

Exposure of spectrally distinct material by impact craters on Mercury: Implications for global stratigraphy

MESSENGER’s Mercury Dual Imaging System (MDIS) obtained multispectral images for more than 80% of the surface of Mercury during its first two flybys. Those images have confirmed that the surface of Mercury exhibits subtle color variations, some of which can be attributed to compositional differences. In many areas, impact craters are associated with material that is spectrally distinct from the surrounding surface. These deposits can be located on the crater floor, rim, wall, or central peak or in the ejecta deposit, and represent material that originally resided at depth and was subsequently excavated during the cratering process. The resulting craters make it possible to investigate the stratigraphy of Mercury’s upper crust. Studies of laboratory, terrestrial, and lunar craters provide a means to bound the depth of origin of spectrally distinct ejecta and central peak structures. Excavated red material (RM), with comparatively steep (red) spectral slope, and low-reflectance material (LRM) stand out prominently from the surrounding terrain in enhanced-color images because they are spectral end-members in Mercury’s compositional continuum. Newly imaged examples of RM were found to be spectrally similar to the relatively red, high-reflectance plains (HRP), suggesting that they may represent deposits of HRP-like material that were subsequently covered by a thin layer (∼1 km thick) of intermediate plains. In one area, craters with diameters ranging from 30 km to 130 km have excavated and incorporated RM into their rims, suggesting that the underlying RM layer may be several kilometers thick. LRM deposits are useful as stratigraphic markers, due to their unique spectral properties. Some RM and LRM were excavated by pre-Tolstojan basins, indicating a relatively old age (>4.0 Ga) for the original emplacement of these deposits. Detailed examination of several small areas on Mercury reveals the complex nature of the local stratigraphy, including the possible presence of buried volcanic plains, and supports sequential buildup of most of the upper ∼5 km of crust by volcanic flows with compositions spanning the range of material now visible on the surface, distributed heterogeneously across the planet. This emerging picture strongly suggests that the crust of Mercury is characterized by a much more substantial component of early volcanism than represented by the phase of mare emplacement on Earth’s Moon.     

The planet Mercury: Synthesis of resolved images of unknown part in the longitude range 250-290°W

Plans to send orbiter missions to Mercury (e.g., NASA's Messenger and ESA's BepiColombo) have prompted renewed efforts to investigate the surface of Mercury using ground-based remote sensing. While the highest resolution instrumentation optical telescopes (e.g. HST) cannot be used at small angular distances (<45°) from the Sun (Mercury's elongation never exceeds 28° seen from Earth), advanced ground-based astronomical techniques and modern processing software can be used to construct resolved images of the poorly known part of Mercury. Our observations of the planet presented here were carried out mainly in April and May, 2002, at evening elongation of the planet, at the Skinakas astrophysical observatory of Heraklion University (Crete, Greece). A synthesis of the acquired images of the hemisphere of Mercury, which was not observed by the Mariner 10 mission (1974-1975), is presented. A double rim basin with an internal diameter of about 1000 km and an external rim about 2000 km is suggested by the data. We present the observational method, the data analysis approach, and the resulting images.     

The transition from complex crater to peak-ring basin on the Moon: New observations from the Lunar Orbiter Laser Altimeter (LOLA) instrument

Impact craters on planetary bodies transition with increasing size from simple, to complex, to peak-ring basins and finally to multi-ring basins. Important to understanding the relationship between complex craters with central peaks and multi-ring basins is the analysis of protobasins (exhibiting a rim crest and interior ring plus a central peak) and peak-ring basins (exhibiting a rim crest and an interior ring). New data have permitted improved portrayal and classification of these transitional features on the Moon. We used new 128 pixel/degree gridded topographic data from the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter, combined with image mosaics, to conduct a survey of craters >50 km in diameter on the Moon and to update the existing catalogs of lunar peak-ring basins and protobasins. Our updated catalog includes 17 peak-ring basins (rim-crest diameters range from 207 km to 582 km, geometric mean = 343 km) and 3 protobasins (137-170 km, geometric mean = 157 km). Several basins inferred to be multi-ring basins in prior studies (Apollo, Moscoviense, Grimaldi, Freundlich-Sharonov, Coulomb-Sarton, and Korolev) are now classified as peak-ring basins due to their similarities with lunar peak-ring basin morphologies and absence of definitive topographic ring structures greater than two in number. We also include in our catalog 23 craters exhibiting small ring-like clusters of peaks (50-205 km, geometric mean = 81 km); one (Humboldt) exhibits a rim-crest diameter and an interior morphology that may be uniquely transitional to the process of forming peak rings. A power-law fit to ring diameters (D_ring) and rim-crest diameters (D_r) of peak-ring basins on the Moon [D_ring = 0.14 ± 0.10(D_r)^1.21±0.13] reveals a trend that is very similar to a power-law fit to peak-ring basin diameters on Mercury [D_ring = 0.25 ± 0.14(D_rim)^1.13±0.10] [Baker, D.M.H. et al. [2011]. Planet. Space Sci., in press]. Plots of ring/rim-crest ratios versus rim-crest diameters for peak-ring basins and protobasins on the Moon also reveal a continuous, nonlinear trend that is similar to trends observed for Mercury and Venus and suggest that protobasins and peak-ring basins are parts of a continuum of basin morphologies. The surface density of peak-ring basins on the Moon (4.5 × 10⁻⁷ per km²) is a factor of two less than Mercury (9.9 × 10⁻⁷ per km²), which may be a function of their widely different mean impact velocities (19.4 km/s and 42.5 km/s, respectively) and differences in peak-ring basin onset diameters. New calculations of the onset diameter for peak-ring basins on the Moon and the terrestrial planets re-affirm previous analyses that the Moon has the largest onset diameter for peak-ring basins in the inner Solar System. Comparisons of the predictions of models for the formation of peak-ring basins with the characteristics of the new basin catalog for the Moon suggest that formation and modification of an interior melt cavity and nonlinear scaling of impact melt volume with crater diameter provide important controls on the development of peak rings. In particular, a power-law model of growth of an interior melt cavity with increasing crater diameter is consistent with power-law fits to the peak-ring basin data for the Moon and Mercury. We suggest that the relationship between the depth of melting and depth of the transient cavity offers a plausible control on the onset diameter and subsequent development of peak-ring basins and also multi-ring basins, which is consistent with both planetary gravitational acceleration and mean impact velocity being important in determining the onset of basin morphological forms on the terrestrial planets.     

Importance of aeolian processes in the origin of the north polar chasmata, Mars

A detailed examination of the location and orientation of sand dunes and other aeolian features within the north polar chasmata indicates that steep scarps strongly influence the direction and intensity of prevailing winds. These steep scarps are present at the heads and along the margins of the north polar chasmata. Topographic profiles of the arcuate head scarps and equator-facing wall of Chasma Boreale reveal unusually steep polar slopes ranging from ∼6°-30°. The relatively steep-sloped (∼8°), sinuous scarp at the head of two smaller chasmata, located west of Chasma Boreale, exhibits an obvious resistant cap-forming unit. Scarp retreat is occurring in places where the cap unit is actively being undercut by descending slope winds. Low-albedo surfaces lacking sand dunes or dust mantles are present at the base of the polar scarps. A ∼100-300 m deep moat, located at the base of the scarps, corresponds with these surfaces and indicates an area of active aeolian scour from descending katabatic winds. Small local dust storms observed along the equator-facing wall of Chasma Boreale imply that slope wind velocities in Chasma Boreale are sufficient to mobilize dust and sand-sized particles in the Polar Layered Deposits (PLD). Two amphitheater forms, located above the cap-forming unit of the sinuous scarp and west of Chasma Boreale, may represent an early stage of polar scarp and chasma formation. These two forms are developing within a younger section of polar layered materials. The unusually steep scarps associated with the polar chasmata have developed where resistant layers are present in the PLD, offering resistance during the headward erosion and poleward retreat of the scarps. Steep slopes that formed under these circumstances enhance the flow of down-scarp katabatic winds. On the basis of these observations, we reject the fluvial outflood hypothesis for the origin of the north polar chasmata and embrace a wind erosion model for their long-term development. In the aeolian model, off-pole katabatic winds progressively remove materials from the steep slopes below chasmata scarps, undermining resistant layers at the tops of scarps and causing retreat by headward erosion. Assuming a minimum age for the onset of formation of Chasma Boreale (10⁵ yr) results in a maximum volumetric erosion rate of . Removal of this volume of material from the equator-facing wall and head scarps of chasma would require a rate for scarp retreat of .     

East-west faults due to planetary contraction

Contraction, expansion and despinning have been common in the past evolution of Solar System bodies. These processes deform the lithosphere until it breaks along faults. Their characteristic tectonic patterns have thus been sought for on all planets and large satellites with an ancient surface. While the search for despinning tectonics has not been conclusive, there is good observational evidence on several bodies for the global faulting pattern associated with contraction or expansion, though the pattern is seldom isotropic as predicted. The cause of the non-random orientation of the faults has been attributed either to regional stresses or to the combined action of contraction/expansion with another deformation (despinning, tidal deformation, reorientation). Another cause of the mismatch may be the neglect of the lithospheric thinning at the equator or at the poles due either to latitudinal variation in solar insolation or to localized tidal dissipation. Using thin elastic shells with variable thickness, I show that the equatorial thinning of the lithosphere transforms the homogeneous and isotropic fault pattern caused by contraction/expansion into a pattern of faults striking east-west, preferably formed in the equatorial region. By contrast, lithospheric thickness variations only weakly affect the despinning faulting pattern consisting of equatorial strike-slip faults and polar normal faults. If contraction is added to despinning, the despinning pattern first shifts to thrust faults striking north-south and then to thrust faults striking east-west. If the lithosphere is thinner at the poles, the tectonic pattern caused by contraction/expansion consists of faults striking north/south. I start by predicting the main characteristics of the stress pattern with symmetry arguments. I further prove that the solutions for contraction and despinning are dual if the inverse elastic thickness is limited to harmonic degree two, making it easy to determine fault orientation for combined contraction and despinning. I give two methods for solving the equations of elasticity, one numerical and the other semi-analytical. The latter method yields explicit formulas for stresses as expansions in Legendre polynomials about the solution for constant shell thickness. Though I only discuss the cases of a lithosphere thinner at the equator or at the poles, the method is applicable for any latitudinal variation of the lithospheric thickness. On Iapetus, contraction or expansion on a lithosphere thinner at the equator explains the location and orientation of the equatorial ridge. On Mercury, the combination of contraction and despinning makes possible the existence of zonal provinces of thrust faults differing in orientation (north-south or east-west), which may be relevant to the orientation of lobate scarps.     

Sodium winds on Mercury

Solar radiation acceleration imparts anti-sunward velocities to sodium atoms in the Mercury exosphere. The Earthward-directed vectors of the Sun-accelerated atom velocities can be observed from Earth as small Doppler shifts, either added to, or subtracted from the Earth-Mercury Doppler shifts. We measured these small Doppler shifts using high resolution spectrographs capable of detecting sodium velocity differences as small as 0.1 km/s. We report here four sets of observations performed at different Mercury true anomaly angles. For these measurements, the spectrograph slit was oriented first east-west, and then north-south on the planet so as to get east-west and north-south transects of the velocities. The velocity patterns in east-west transects could be explained in terms of sodium flows outwards from the subsolar point, except for unexpectedly large Earthward velocities observed above the dawn terminator, which we interpreted to be the result of evaporation of sodium as the cold surface is heated by the rising Sun. North-south transects also showed a general pattern consistent with sodium flows outwards from the subsolar point. However, in all cases, the velocities were higher in one hemisphere relative to the other. For two cases, excess sodium emission was observed in the same hemisphere as the velocity excess. We interpreted these results to mean that there existed sources of sodium at high latitudes, which could appear in either hemisphere.     

Mechanical modeling of thrust faults in the Thaumasia region, Mars, and implications for the Noachian heat flux

Insight into the state of the early martian lithosphere is gained by modeling the topography above surface breaking thrust faults in the southern Thaumasia region. Crater counts of key surface units associated with the faulting indicate a scarp emplacement in the late Noachian-early Hesperian periods between 4.0 and 3.7 Gyr. The seismogenic layer thickness at the time of faulting is constrained to 27-35 km and 21-28 km for the two scarps investigated, implying paleo geothermal gradients of 12-18 and 15-23 K km⁻¹, corresponding to heat flows of 24-36 and 30-46 mW m⁻². The heat flow values obtained in this study are considerably lower than those derived from rift flank uplift at the close-by Coracis Fossae for a similar time period, indicating that surface heat flow is a strong function of regional setting. If viewed as representative for magmatically active and inactive regions, the thermal gradients at rifts and scarps span the range of admissible global mean values. This implies , with the true value probably being closer to the lower bound.     

Ice-lubricated gravity spreading of the Olympus Mons aureole deposits

Gravity sliding and spreading at low strain rates can account for the general morphology and structure of the aureoles and basal scarp of Olympus Mons. Detachment sliding could have occurred around the volcano if either pore-fluid pressures were exceptionally high (greater than 90%) or the rocks had very low resistance to shear (about 1 × 10⁵ Pa or 1 bar). Because of the vast areal extent and probable shallow depth of the detachment zone, development of ubiquitous, high pore-fluid pressures beneath aureole-forming material was unlikely. However, a zone of sufficiently weak material consisting of about 10% interstitial or interbedded ice could have been present. If so, a simple rheologic model for the aureole deposits can be applied that consists of a thin ductile layer overlain by a thicker brittle layer. According to this model, extensional deformation would have occurred near the shield and compressional deformation in its distal parts. Proximal grabens and distal corrugations on aureole surfaces support this model. A submarine slide at Kitimat Arm, British Columbia, is a valid qualitative analogy for the observed features and inferred emplacement style of the aureole deposits. Ground-ice processes have been considered the cause of many geologic features on Mars; a 3% average concentration of ground ice in the regolith is predicted by theoretical models for the ice budget and cryosphere. Ice may have been deposited in higher concentrations below the aureole-forming material; the source of the ice could have been juvenile water circulated hydrothermally by Olympus Mons volcanism. The basal scarp of Olympus Mons apparently demarcates the transition between the upper, stable part of the shield and its lower part that decoupled and formed the aureole deposits. This transition may reflect a change in the bulk shear strength of the shield, caused either by a radial dependence in the abundance of ice or fluid in the shield materials or by the concentration of intrusive dikes within the volcano. Other Martian volcanoes exhibit virtually no evidence of similar large-scale gravity spreading and basal scarps. Perhaps such evidence, if it existed, has been buried by lava flows, or perhaps the smaller size of other volcanoes did not permit the development of these features.     

Observations of the sodium tail of Mercury

Cross-sections of the sodium emission tail of Mercury were measured at various distances down the tail when Mercury was moving away from the Sun (true anomaly angles <180°), and again when Mercury was moving towards the Sun (true anomaly angles >180°). As predicted in early modeling studies, significant differences were expected between these two cases, as the result of Doppler shifts to higher solar intensity in the former case, and to lower solar intensity for the latter case. For observations with Mercury moving away from the Sun, the sodium tail was observed out to about 40,000 kilometers (16 Mercury radii, RM) downstream, expanding, on average, at a rate of 1.9±0.3 km/s. The source rates for sodium generation from Mercury into the tail were found to be in the range 2-5×¹⁰²³ atoms/s, corresponding to between 1 and 10% of the estimated total sodium production rate on the planet. The limiting value of radiation acceleration required to produce an observable sodium tail was estimated to be 112±24 cm/s². For observations where Mercury was moving towards the Sun, the emission intensity in the sodium tail decreased very rapidly with distance downstream, disappearing entirely beyond 12,000 (6 RM) kilometers for radiation accelerations of 128.7 and 135.4 cm/s². For smaller radiation accelerations, the sodium tail was not detectable at all, yielding a limiting value for tail generation of about 122±2 cm/s². Interpretation of the limiting radiation acceleration values suggests that the process that generates the sodium tail yields atoms with energies greater than 3 eV. Particle sputtering is the most reasonable source process.