Can a fractionally crystallized magma ocean explain the thermo-chemical evolution of Mars?

Abstract

The impact heat accumulated during the late stage of planetary accretion can melt a significant part or even the entire mantle of a terrestrial body, giving rise to a global magma ocean. The subsequent cooling of the interior causes the magma ocean to freeze from the core-mantle boundary (CMB) to the surface due to the steeper slope of the mantle adiabat compared to the slope of the solidus.

Assuming fractional crystallization of the magma ocean, dense cumulates are produced close to the surface, largely due to iron enrichment in the evolving magma ocean liquid. A gravitationally unstable mantle thus forms, which is prone to overturn. We investigate the cumulate overturn and its influence on the thermal evolution of Mars using mantle convection simulations in 2D cylindrical geometry. We present a suite of simulations using different initial conditions and a strongly temperature-dependent viscosity. We assume that all radiogenic heat sources have been enriched during the freezing-phase of the magma ocean in the uppermost 50 km and that the initial steam-atmosphere created by the degassing of the freezing magma ocean was rapidly lost, implying that the surface temperature is set to present-day values. In this case, a stagnant lid quickly forms on top of the convective interior preventing the uppermost dense cumulates to sink, even when allowing for a plastic yielding mechanism. Below this dense stagnant lid, the mantle chemical gradient settles to a stable configuration. The convection pattern is dominated by small-scale structures, which are difficult to reconcile with the large-scale volcanic features observed over Mars' surface and partial melting ceases in less than 900 Ma. Assuming that the stagnant lid can break because of additional mechanisms and allowing the uppermost dense layer to overturn, a stable density gradient is obtained, with the densest material and the entire amount of heat sources lying above the CMB. This stratification leads to a strong overheating of the lowermost mantle, whose temperature increases to values that exceed the liquidus. The iron-rich melt would most likely remain trapped in the lower part of the mantle. The upper mantle in that scenario cools rapidly and only shows partial melting during the first billion year of evolution. Therefore a fractionated global and deep magma ocean is difficult to reconcile with observations. Different scenarios assuming, for instance, a hemispherical or shallow magma ocean, or a crystallization sequence resulting in a lower density gradient than that implied by pure fractional crystallization will have to be considered.