A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy

A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51-71 km thick with a southern crust thicker by 18-28 km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat-producing elements, the lid of a one-plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well-mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.

The History of Water in Martian Magmas From Thorium Maps

Water inventories in Martian magmas are poorly constrained. Meteorite-based estimates range widely, from 10² to >10⁴ ppm H₂O, and are likely variably influenced by degassing. Orbital measurements of H primarily reflect water cycled and stored in the regolith. Like water, Th behaves incompatibly during mantle melting, but unlike water Th is not prone to degassing and is relatively immobile during aqueous alteration at low temperature. We employ Th as a proxy for original, mantle-derived H₂O in Martian magmas. We use regional maps of Th from Mars Odyssey to assess variations in magmatic water across major volcanic provinces and through time. We infer that Hesperian and Amazonian magmas had ∼100-3,000 ppm H₂O, in the lower range of previous estimates. The implied cumulative outgassing since the Hesperian, equivalent to a global H₂O layer ∼1-40 m deep, agrees with Mars' present-day surface and near-surface water inventory and estimates of sequestration and loss rates.

Thickness and structure of the martian crust from InSight seismic data

A planet's crust bears witness to the history of planetary formation and evolution, but for Mars, no absolute measurement of crustal thickness has been available. Here, we determine the structure of the crust beneath the InSight landing site on Mars using both marsquake recordings and the ambient wavefield. By analyzing seismic phases that are reflected and converted at subsurface interfaces, we find that the observations are consistent with models with at least two and possibly three interfaces. If the second interface is the boundary of the crust, the thickness is 20 ± 5 kilometers, whereas if the third interface is the boundary, the thickness is 39 ± 8 kilometers. Global maps of gravity and topography allow extrapolation of this point measurement to the whole planet, showing that the average thickness of the martian crust lies between 24 and 72 kilometers. Independent bulk composition and geodynamic constraints show that the thicker model is consistent with the abundances of crustal heat-producing elements observed for the shallow surface, whereas the thinner model requires greater concentration at depth.

Groundwater production from geothermal heating on early Mars and implication for early martian habitability

In explaining extensive evidence for past liquid water, the debate on whether Mars was primarily warm and wet or cold and arid 4 billion years (Ga) ago has continued for decades. The Sun's luminosity was ~30% lower 4 Ga ago; thus, most martian climate models struggle to elevate the mean surface temperature past the melting point of water. Basal melting of ice sheets may help resolve that paradox. We modeled the thermophysical evolution of ice and estimate the geothermal heat flux required to produce meltwater on a cold, arid Mars. We then analyzed geophysical and geochemical data, showing that basal melting would have been feasible on Mars 4 Ga ago. If Mars were warm and wet 4 Ga ago, then the geothermal flux would have even sustained hydrothermal activity. Regardless of the actual nature of the ancient martian climate, the subsurface would have been the most habitable region on Mars.

Radiolytic H2 Production in Martian Environments

Hydrogen, produced by water radiolysis, has been suggested to support microbial communities on Mars. We quantitatively assess the potential magnitude of radiolytic H2 production in wet martian environments (the ancient surface and the present subsurface) based on the radionuclide compositions of (1) eight proposed Mars 2020 landing sites, and (2) three sites that individually yield the highest or lowest calculated radiolytic H2 production rates on Mars. For the proposed landing sites, calculated H2 production rates vary by a factor of ∼1.6, while the three comparison sites differ by a factor of ∼6. Rates in wet martian sediment and microfractured rock are comparable with rates in terrestrial environments that harbor low concentrations of microbial life (e.g., subseafloor basalt). Calculated H2 production rates for low-porosity (<35%), fine-grained martian sediment (0.12-1.2 nM/year) are mostly higher than rates for South Pacific subseafloor basalt (∼0.02-0.6 nM/year). Production rates in martian high-porosity sediment (>35%) and microfractured (1 μm) hard rock (0.03 to <0.71 nM/year) are generally similar to rates in South Pacific basalt, while yields for larger martian fractures (1 and 10 cm) are one to two orders of magnitude lower (<0.01 nM/year). If minerals or brine that amplify radiolytic H2 production rates are present, H2 yields exceed the calculated rates.

Present-day heat flow model of Mars

Until the acquisition of in-situ measurements, the study of the present-day heat flow of Mars must rely on indirect methods, mainly based on the relation between the thermal state of the lithosphere and its mechanical strength, or on theoretical models of internal evolution. Here, we present a first-order global model for the present-day surface heat flow for Mars, based on the radiogenic heat production of the crust and mantle, on scaling of heat flow variations arising from crustal thickness and topography variations, and on the heat flow derived from the effective elastic thickness of the lithosphere beneath the North Polar Region. Our preferred model finds heat flows varying between 14 and 25 mW m^-2, with an average value of 19 mW m^-2. Similar results (although about ten percent higher) are obtained if we use heat flow based on the lithospheric strength of the South Polar Region. Moreover, expressing our results in terms of the Urey ratio (the ratio between total internal heat production and total heat loss through the surface), we estimate values close to 0.7-0.75, which indicates a moderate contribution of secular cooling to the heat flow of Mars (consistent with the low heat flow values deduced from lithosphere strength), unless heat-producing elements abundances for Mars are subchondritic.