Numerical Insights into the Formation and Stability of Cratons

<div>Geophysical, geochemical, and geological investigations have attributed the stable behaviour of Earth's continents to the presence of strong and viscous cratons underlying the continental crust. The cratons are underlain by thick and cold mantle keels, which are composed of melt-depleted and low density peridotite residues [1]. Progressive melt extraction increases the magnesium number Mg# in the residual peridotite, thereby making the roots of cratons chemically buoyant [2, 3] and counteracting their negative thermal buoyancy. Recent global models have shown the production of Archean continental crust by two-step mantle differentiation, however this primordial crust gets recycled and no stable continents form [4]. This points to the missing ingredient of cratonic lithosphere in these models, which could act as a stable basement for the crustal material to accumulate on and may also help with the transition of global regime from "vertical tectonics'' to "horizontal tectonics''. Based on the bulk FeO and MgO content of the residual peridotites, it has been proposed that cratonic mantle formed by hot shallow melting with mantle potential temperature, which was higher by 200-300 °C than present-day [5]. We introduce Fe-Mg partitioning between mantle peridotite and melt to track the Mg# variation through melting, and parametrise craton formation using the corresponding P-T formation conditions. Using self-consistent global convection models, we show the dynamic formation of cratons as a result of naturally occurring lateral compression and thickening of the lithosphere, which has been suggested by geochemical and petrological data. To allow for the material to compact and thicken, but prevent it from collapsing under its own weight, a combination of lithospheric strength, plastic yielding, dehydration strengthening, and depletion-induced density reduction of the depleted mantle material is necessary.</div><div> </div><div> [1] Boyd, F. R. High-and low-temperature garnet peridotite xenoliths and their possible relation to the lithosphere- asthenosphere boundary beneath Africa. In Nixon, P. H. (ed.) <em>Mantle Xenolith</em>, 403–412 (John Wiley & Sons Ltd., 1987).</div><div>[2] Jordan, T. H. Mineralogies, densities and seismic velocities of garnet lherzolites and their geophysical implications. In <em>The Mantle Sample: Inclusion in Kimberlites and Other Volcanics</em>, 1–14 (American Geophysical Union, Washington, D. C., 1979).</div><div>[3] Schutt, D. L. & Lesher, C. E. Effects of melt depletion on the density and seismic velocity of garnet and spinel lherzolite. <em>Journal of Geophysical Research </em><strong>111</strong> (2006).</div><div>[4] Jain, C., Rozel, A. B., Tackley, P. J., Sanan, P. & Gerya, T. V. Growing primordial continental crust self-consistently in global mantle convection models. <em>Gondwana Research</em> <strong>73</strong>, 96–122 (2019).</div><div>[5] Lee, C.-T. A. & Chin, E. J. Calculating melting temperatures and pressures of peridotite protoliths: Implications for the origin of cratonic mantle. <em>Earth and Planetary Science Letters</em> <strong>403</strong>, 273–286 (2014)</div>