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.