Numerical Insights into the Formation and Stability of Cratons

Geophysical, geochemical, and geological investigations have attributed the stable behavior 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 (Boyd, 1987). Progressive melt extraction increases the magnesium number Mg# in the residual peridotite, thereby making the roots of cratons chemically buoyant (Jordan, 1979; Schutt and Lesher, 2006) and counteracting their negative thermal buoyancy. Recent global models have shown the self-consistent production of Archean continental crust by two-step mantle differentiation (Jain et al., 2019). These models exhibit intense recycling of crust with delamination and eclogitic dripping in the first 500 million years and this behavior is similar to the "plutonic-squishy lid" that has been suggested for the early Earth. However, no stable continents form and no major regime transition from "vertical tectonics" towards "horizontal tectonics" is observed. 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 help initiate plate 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 (Lee and Chin, 2014). We will introduce Fe-Mg partitioning between mantle peridotite and melt to track the Mg# variation through melting, and parametrize craton formation using the corresponding P-T formation conditions. Grain-size evolution, which has been shown to influence mantle rheology (Hall and Parmentier, 2003) is another mechanism that may contribute towards cratonic strength and will be explored using self-consistent global geodynamic models.