The evolution of the martian elastic lithosphere and implications for crustal and mantle rheology

Estimates of the martian elastic lithosphere thickness Te imply that Te increased from around 20 km in the Noachian to about 70 km in the Amazonian period. A phase of rapid lithospheric growth is observed during the Hesperian and we propose that this elastic thickness history is a consequence of the martian crustal rheology and its thermal evolution. A wet crustal rheology is found to generate a mechanically incompetent layer in the lower crust during the early evolution and the rapid growth of Te during the Hesperian results from the disappearance of this layer due to planetary cooling. The incompetent layer and the related rapid lithospheric growth are absent for a dry basaltic crustal rheology, which is therefore incompatible with the observations. Furthermore, we find that the observed elastic thickness evolution is best compatible with a wet mantle rheology, although a dry mantle cannot be ruled out. It therefore seems likely that rheologically significant amounts of water were retained in the Martian crust and mantle after planetary accretion.     

The present-day thermal state of Mars

The present-day thermal state of the martian interior is a very important issue for understanding the internal evolution of the planet. Here, in order to obtain an improved upper limit for the heat flow at the north polar region, we use the lower limit of the effective elastic thickness of the lithosphere loaded by the north polar cap, crustal heat-producing elements (HPE) abundances based on martian geochemistry, and a temperature-dependent thermal conductivity for the upper mantle. We also perform similar calculations for the south polar region, although uncertainties in lithospheric flexure make the results less robust. Our results show that the present-day surface and sublithospheric heat flows cannot be higher than 19 and 12 mW m⁻², respectively, in the north polar region, and similar values might be representative of the south polar region (although with a somewhat higher surface heat flow due to the radioactive contribution from a thicker crust). These values, if representative of martian averages, do not necessarily imply sub-chondritic HPE bulk abundances for Mars (as previously suggested), since (1) chondritic composition models produce a present-day total heat power equivalent to an average surface heat flow of 14-22 mW m⁻² and (2) some convective models obtain similar heat flows for the present time. Regions of low heat flow may even have existed during the last billions of years, in accordance with several surface heat flow estimates of ∼20 mW m⁻² or less for terrains loaded during Hesperian or Amazonian times. On the other hand, there are some evidences suggesting the current existence of regions of enhanced heat flow, and therefore average heat flows could be higher than those obtained for the north (and maybe the south) polar region.     

Crustal recycling, mantle dehydration, and the thermal evolution of Mars

We have reinvestigated the coupled thermal and crustal evolution of Mars taking new laboratory data concerning the flow behavior of iron-rich olivine into account. The low mantle viscosities associated with the relatively higher iron content of the martian mantle as well as the observed high concentrations of heat producing elements in a crust with a reduced thermal conductivity were found to promote phases of crustal recycling in many models. As crustal recycling is incompatible with an early separation of geochemical reservoirs, models were required to show no episodes of crustal recycling. Furthermore, admissible models were required to reproduce the martian crust formation history, to allow for the formation of partial melt under present day mantle conditions and to reproduce the measured concentrations of potassium and thorium on the martian surface. Taking dehydration stiffening of the mantle viscosity by the extraction of water from the mantle into account, we found that admissible models have low initial upper mantle temperatures around 1650 K, preferably a primordial crustal thickness of 30 km, and an initially wet mantle rheology. The crust formation process on Mars would then be driven by the extraction of a primordial crust after core formation, cooling the mantle to temperatures close to the peridotite solidus. According to this scenario, the second stage of global crust formation took place over a more extended period of time, waning at around 3500 Myr b.p., and was driven by heat produced by the decay of radioactive elements. Present-day volcanism would then be driven by mantle plumes originating at the core-mantle boundary under regions of locally thickened, thermally insulating crust. Water extraction from the mantle was found to be relatively efficient and close to 40% of the total inventory was lost from the mantle in most models. Assuming an initial mantle water content of 100 ppm and that 10% of the extracted water is supplied to the surface, this amount is equivalent to a 14 m thick global surface layer, suggesting that volcanic outgassing of H₂O could have significantly influenced the early martian climate and increased the planet’s habitability.     

Aram Chaos and its constraints on the surface heat flux of Mars

The surface heat flux of a planet is an important parameter to characterize its internal activity and to determine its thermal evolution. Here we report on a new method to constrain the surface heat flux of Mars during the Hesperian. For this, we explore the consequences for the martian surface heat flux from a recently presented new hypothesis for the formation of Aram Chaos (Zegers, T.E., Oosthoek, J.H.P., Rossi, A.P., Blom, J.K., Schumacher, S. [2010]. Earth Planet. Sci. Lett. 297, 496-504. doi:10.1016/j.epsl.2010.06.049.). In this hypothesis the chaotic terrain is thought to have formed by melting of a buried ice sheet. The slow sedimentation and burial of the ice sheet led to an increased thermal insulation of the ice and subsequently to a temperature increase high enough to trigger melting and the formation of the subsurface lake. As these processes highly depend on the thermal properties of the subsurface and especially on the surface heat flux, it is possible to constrain the latter by using numerical simulations. Based on the hypothesis for the formation of Aram Chaos, we conducted an extensive parameter study to determine the parameter settings leading to sufficient melting of the buried ice sheet. We find that the surface heat flux in the Aram Chaos region during the Hesperian was most likely between 20 and 45 mW m⁻² with a possible maximum value of up to 60 mW m⁻².     

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.     

The thermal evolution of Mars as constrained by paleo-heat flows

Lithospheric strength can be used to estimate the heat flow at the time when a given region was deformed, allowing us to constrain the thermal evolution of a planetary body. In this sense, the high (>300 km) effective elastic thickness of the lithosphere deduced from the very limited deflection caused by the north polar cap of Mars indicates a low surface heat flow for this region at the present time, a finding difficult to reconcile with thermal history models. This has started a debate on the current heat flow of Mars and the implications for the thermal evolution of the planet. Here we perform refined estimates of paleo-heat flow for 22 martian regions of different periods and geological context, derived from the effective elastic thickness of the lithosphere or from faulting depth beneath large thrust faults, by considering regional radioactive element abundances and realistic thermal conductivities for the crust and mantle lithosphere. For the calculations based on the effective elastic thickness of the lithosphere we also consider the respective contributions of crust and mantle lithosphere to the total lithospheric strength. The obtained surface heat flows are in general lower than the equivalent radioactive heat production of Mars at the corresponding times, suggesting a limited contribution from secular cooling to the heat flow during the majority of the history of Mars. This is contrary to the predictions from the majority of thermal history models, but is consistent with evidence suggesting a currently fluid core, limited secular contraction for Mars, and recent extensive volcanism. Moreover, the interior of Mars could even have been heating up during part of the thermal history of the planet.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

The influence of mantle melting on the evolution of Mars

We present a parameterized convection model of Mars by incorporating a new heat-flow scaling law for stagnant-lid convection, to better understand how the evolution of Mars may be affected by mantle melting. Melting in the mantle during convection leads to the formation of a compositionally buoyant lithosphere, which may also be intrinsically more viscous by dehydration. The consequences of these melting effects on the evolution of terrestrial planets have not been explored before. The temporal evolution of crust and lithospheric mantle is modeled in a self-consistent manner considering mantle melting, convective instability, and the rewetting of dehydrated lithosphere from below by hydrogen diffusion. Though the effect of compositional buoyancy turns out to be minimal, the introduction of viscosity contrast between wet and dry mantle can considerably slow mantle cooling and sometimes lead to non-monotonic core cooling. Furthermore, with or without dehydration stiffening, our model predicts that the martian mantle must have been degassed more extensively (>80%) than previously suggested (<10%); the loss of such a large amount of water from the mantle to surface has significant implications about the role of water in the early surface and climate evolution of Mars.     

Scenarios for the formation of Chasma Boreale, Mars

The martian polar caps feature large chasmata and smaller trough systems which have no counterpart in terrestrial ice sheets. Chasma Boreale cuts about 500 km into the western part of the north-polar cap, is up to 100 km wide and up to 2 km deep. One possible formation mechanism is by a temporary heat source under the ice due to tectono-thermal or volcanic activity, which melts the ice from below. It is demonstrated by model simulations that this process is feasible, a moderately increased heat flux of 0.5-1 W m⁻², sustained over at least tens of thousands of years, producing a topographic depression which resembles the real chasma. Associated meltwater discharge rates are small (), but can exceed 10 km³ a⁻¹ if a stronger heat flux of 10 W m⁻² is assumed. Local ice-flow velocities during the process of chasma formation can exceed 1 m a⁻¹ at the head and scarps of the chasma. However, if the thermal anomaly shuts down, glacial flow quickly decreases, so that the chasma can stay open for an indefinite amount of time without an ongoing, sustaining process under the climate conditions of the most recent millions of years.     

The deep interior of Venus, Mars, and the Earth: A brief review and the need for planetary surface-based measurements

A comparison of the internal structure of Earth-like planets is unavoidable to understand the formation and evolution of the solar system, and the differences between Earth’s, Mars’, and Venus’ atmospheres, surfaces and tectonic behaviors. Recent studies point at the role of core structure and dynamics in the evolution of the atmosphere, mantle and crust. On Earth, the crust thickness and the radius and physical state of the cores are known for almost one century, since the advent of seismological observations, but the lack of long-term surface-based geodetic, electromagnetic and seismological observations on the other planets, results in very large uncertainties on the crust thickness, on the temperature and composition of their mantle, and on the size and physical state of their cores. According to the currently available geodetic data, Mars’ dimensionless mean-moment-of-inertia ratio is equal to 0.3653±0.0008. When combined with geochemical observations and with the inputs of laboratory experiments on planetary materials at high pressure and high temperature, this result constrains a narrow range of density values for Mars’ mantle and favors a light [6200-6765 kg m⁻³] sulfur-rich core, but it still allows for a 1600-1750 km range for the core radius, i.e. an uncertainty at least ten times larger than the precision obtained in 1913 by Gutenberg for the Earth’s core. Mars’ mantle density distribution may be explained by a large range of temperatures and mineralogical compositions, either olivine- or pyroxene-rich. The unknown mean thickness of Mars’ crust makes necessary a number of working assumptions for the interpretation of gravimetric and magnetic data. The situation is worse for Venus, and the most conservative model of its deep interior is a transposition of the Earth’s structure scaled to Venus’ radius and mass. The temperature conditions at the surface of this planet hardly make possible long-term ground-based measurements, but this is indeed feasible at the surface of Mars. Precise measurements of Mars’ crust thickness, core radius and structure, and the proof of the existence or absence of an inner core, would put tight constraints on mantle dynamics and thermal evolution, and on possible scenarios leading to the extinction of Mars’ magnetic field about 4.0 Ga ago. Long-lasting surface-based geodetic, seismological and magnetic observations would provide this information, as well as the distributions as a function of depth of the density, elastic and anelastic parameters, and electrical conductivity. Current studies on the structure of Earth’s deep interior demonstrate that the latter data set, when constrained by laboratory experiments, may be inverted in terms of temperature, chemical, and mineralogical compositions.     

A top-down origin for martian mantle plumes

The two main volcanic centers on Mars, Tharsis and Elysium, are often interpreted in terms of mantle plume hotspots, even though there are several problems with the plume hypothesis for Mars. We present results of 2D cylindrical shell numerical mantle convection experiments in which we try to ascertain whether flushing of the hot lower mantle could provide a mechanism for the generation of a small number of plume-like features, i.e., localized upwelling of hot material. In this scenario the formation of hot upwellings is driven from the top by cold downwellings rather than from a hot thermal boundary layer at the CMB. First we construct a range of Mars interior structure models consistent with observations in order to demonstrate that the presence of a thin lower mantle in the martian interior is a viable scenario. Then we use a series of numerical convection experiments to investigate the effects of solid-state phase transitions, different stratified and temperature-dependent viscosity models, and the presence of a thick southern hemisphere crust on the operation of such a mechanism. Our results show that it is possible to generate hot strong localized upwellings from top-down dynamics if the lithosphere is thin or actively involved in the convective pattern. The presence of a thick, immobile, insulating southern hemisphere crust reduces the number of upwellings, and the perovskite phase transition causes a focusing of the upwellings. Further experiments demonstrate that an initial 500 Myr phase of mobile lid is sufficient to start this process create an upwelling which is stable for billions of years.     

Influence of ice rheology and dust content on the dynamics of the north-polar cap of Mars

The evolution and dynamics of the north-polar cap (residual-ice-cap/layered-deposits complex) of Mars is simulated with a thermomechanical ice-sheet model. We consider a scenario with ice-free initial conditions at 5 Ma before present due to the large obliquities which prevailed prior to this time. The north-polar cap is then built up to its present shape, driven by a parameterized climate forcing (surface temperature, surface mass balance) based on the obliquity and eccentricity history. The effects of different ice rheologies and different dust contents are investigated. It is found that the build-up scenarios require an accumulation rate of approximately 0.15-0.2 mm a⁻¹ at present. The topography evolution is essentially independent of the ice dynamics due to the slow ice flow. Owing to the uncertainties associated with the ice rheology and the dust content, flow velocities can only be predicted within a range of two orders of magnitude. Likely present values are of the order of 0.1-1 mm a⁻¹, and a strong variation over the climatic cycles is found. For all cases, computed basal temperatures are far below pressure melting.     

Non-uniform seasonal defrosting of subpolar dune field on Mars

Images of the sub-polar regions of Mars show existence of various seasonal-depending features, including mosaics of dark patches called Dalmatian spots. We interpret formation of dark patches on dune fields by means of the non-uniform sublimation of winter CO₂ ice deposit.     

Exploration of hydrothermal targets on Mars

Based on various lines of geologic, geomorphic, topographic, geophysical, spectral, and elemental evidence, we conclude that hydrothermal environments have certainly existed on Mars and are likely to still exist. Here, we present candidate targets of endogenic- and exogenic-driven hydrothermal environments on Mars based on a set of selection criteria and suggest strategies for the detection of such targets. This includes a re-evaluation of potential targets using both existing and yet-to-be-released remote information provided by the instruments onboard the Mars orbiters and rovers. We also provide terrestrial analogs for possible martian hydrothermal environments to highlight the implications of these targets for potential martian life. This compilation and synthesis of data from martian localities indicating hydrothermal activity is timely and a first step towards prioritizing candidate targets for further investigation, which will likely add more targets to this list. Future in situ exploration will have to focus on the most promising of the hydrothermal targets and investigate them utilizing a novel integrated multi-tier, multi-agent reconnaissance mission architecture.     

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.     

Constraints on the composition of the martian south polar cap from gravity and topography

The polar caps of Mars have long been acknowledged to be composed of unknown proportions of water ice, solid CO₂ (dry ice), and dust. Gravity and topography data are here analyzed over the southern cap to place constraints on its density, and hence composition. Using a localized spectral analysis combined with a lithospheric flexure model of ice cap loading, the best fit density of the volatile-rich south polar layered deposits is found to be 1271 kg m⁻³ with 1-σ limits of 1166 and 1391 kg m⁻³. The best fit elastic thickness of this geologically young deposit is 140 km, though any value greater than 102 km can fit the observations. The best fit density implies that about 55% dry ice by volume could be sequestered in these deposits if they were completely dust free. Alternatively, if these deposits were completely free of solid CO₂, the dust content would be constrained to lie between about 14 and 28% by volume. The bulk thermal conductivity of the polar cap is not significantly affected by these maximum allowable concentrations of dust. However, even if a moderate quantity of solid CO₂ were present as horizontal layers, the bulk thermal conductivity of the polar cap would be significantly reduced. Reasonable estimates of the present day heat flow of Mars predict that dry ice beneath the thicker portions of the south polar cap would have melted. Depending on the quantity of solid CO₂ in these deposits today, it is even possible that water ice could melt where the cap is thickest. If independent estimates for either the dust or CO₂ content of the south polar cap could be obtained, and if radar sounding data could determine whether this polar cap is presently experiencing basal melting or not, it would be possible to use these observations to place tight constraints on the present day heat flow of Mars.     

East-west trending magnetic anomalies in the Southern Hemisphere of Mars: Modeling analysis and interpretation

Maps of the vector components of the martian crustal magnetic field over the strongly magnetized Terra Cimmeria/Sirenum region are constructed using Mars Global Surveyor magnetometer data. Although pronounced east-west trending anomalies are present on the radial and north field component maps at the mapping altitude (∼360-380 km), these trends are much less prominent at the lower aerobraking altitude (∼90-150 km). Comparisons with similar maps produced using artificial data at the aerobraking altitude indicate that elongated sources in this region may have maximum lengths along the martian surface of ∼500 km and maximum aspect ratios of ∼2. Iterative forward modeling of several relatively isolated anomalies in the mapped region yields paleomagnetic pole positions consistent with those estimated in previous studies of other anomalies using mapping phase and science phasing orbit data. On this basis, it is inferred that sources in the studied region are most probably magnetized primarily in northward or southward directions. Using this additional constraint, iterative forward modeling is then applied to determine a magnetization distribution that is consistent with data at both the aerobraking altitude and the mapping altitude. The model magnetization distribution, which includes 41 discrete sources, again indicates no highly elongated sources. An examination of surface geology in the region as well as a consideration of the global distribution of anomalies suggests that magmatic intrusions (e.g., subsurface dike swarms), cooling in the presence of water, are the most likely sources of the magnetic anomalies.     

Antarctic Dry Valleys and indigenous weathering in Mars meteorites: Implications for water and life on Mars

The Dry Valleys of Antarctica are an excellent analog of the environment at the surface of Mars. Soil formation histories involving slow processes of sublimation and migration of water-soluble ions in polar desert environments are characteristic of both Mars and the Dry Valleys. At the present time, the environment in the Dry Valleys is probably the most similar to that in the mid-latitudes on Mars although similar conditions may be found in areas of the polar regions during their respective Mars summers. It is thought that Mars is currently in an interglacial period, and that subsurface water ice is sublimating poleward. Because the Mars sublimation zones seem to be the most similar to the Antarctic Dry Valleys, the Dry Valleys-type Mars climate is migrating towards the poles. Mars has likely undergone drastic obliquity changes, which means that the Dry Valleys analog to Mars may be valid for large parts of Mars, including the polar regions, at different times in geologic history. Dry Valleys soils contain traces of silicate alteration products and secondary salts much like those found in Mars meteorites. A martian origin for some of the meteorite secondary phases has been verified previously; it can be based on the presence of shock effects and other features which could not have formed after the rocks were ejected from Mars, or demonstrable modification of a feature by the passage of the meteorite through Earth's atmosphere (proving the feature to be pre-terrestrial). The martian weathering products provide critical information for deciphering the near-surface history of Mars. Definite martian secondary phases include Ca-carbonate, Ca-sulfate, and Mg-sulfate. These salts are also found in soils from the Dry Valleys of Antarctica. Results of earlier Wright Valley work are consistent with what is now known about Mars based on meteorite and orbital data. Results from recent and current Mars missions support this inference. Aqueous processes are active even in permanently frozen Antarctic Dry Valleys soils, and similar processes are probably also occurring on Mars today, especially at the mid-latitudes. Both weathering products and life in Dry Valleys soils are distributed heterogeneously. Such variations should be taken into account in future studies of martian soils and also in the search for possible life on Mars.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).