Replication Data for: Overturn of Ilmenite-Bearing Cumulates in a Rheologically Weak Lunar Mantle|月球地质数据集|流变学模拟数据集

The crystallization of the lunar magma ocean (LMO) determines the initial structure of the solid Moon. Near the end of the LMO crystallization, ilmenite-bearing cumulates (IBC) form beneath the plagioclase crust. Being denser than the underlying mantle, IBC are prone to overturn, a hypothesis that explains several aspects of the Moon's evolution. Yet the formation of stagnant lid due to the temperature dependence of viscosity can easily prevent IBC from sinking. To infer the rheological conditions allowing IBC to sink, we calculated the LMO crystallization sequence and performed high-resolution numerical simulations of the overturn dynamics. We assumed a diffusion creep rheology and tested the effects of reference viscosity, activation energy, and compositional viscosity contrast between IBC and mantle. The overturn strongly depends on reference viscosity and activation energy and is facilitated by a low IBC viscosity. For a reference viscosity of 1021 Pa s, characteristic of a dry rheology, IBC overturn cannot take place. For a reference viscosity of 1020 Pa s, the overturn is possible if the activation energy is a factor of 2–3 lower than the values typically assumed for dry olivine. These low activation energies suggest a role for dislocation creep. For lower-reference viscosities associated with the presence of water or trapped melt, more than 95% IBC can sink regardless of the activation energy. Scaling laws for Rayleigh-Taylor instability confirmed these results but also showed the need of numerical simulations to accurately quantify the overturn dynamics. Whenever IBC sink, the overturn occurs via small-scale diapirs.