Greenstone belts in the West African Craton (WAC) are separated by several generations of granitoids intruded at ca. 2.18–1.98 Ga. Simultaneous folding and exhumation play an important role in the formation of greenstone-granitoid belts. However, the overall tectonic regime and origin of granitoids remain controversial. In this study, we present the estimates of the mantle potential temperature (Tp) for the WAC, which yields values of about 1500–1600 °C, pressure estimates of initial and final melting yield values of about 3.7–5.2 GPa and 1–1.3 GPa, respectively. Subsequently, 2D thermo-mechanical models have been constructed to explore the width of volcano-sedimentary basin on spatial-temporal evolution of diapirs that emplaced in the lower-middle crust during compression. The models show that the width of the volcano-sediment layer plays an important role in the formation mechanisms of greenstone-granitoid belts. The lower crust beneath sedimentary sequences is deformed into a buckle fold during the first compressional stage, through which relief uplifts slowly. Subsequently, the buckle fold is further deformed into several individual folds. Diapirs made of lower crust rocks ascend and emplace in the middle-upper crust resulting from instability. Benefitting from the mantle temperature, the pressure estimates and the numerical modelling results, a new geodynamic model was constructed. This model indicates that a series of sheet-like granitoids possibly derived from either subducted mélanges, lower crust and/or mantle melting that are accumulated at depths of the subcontinental mantle would channel along diapirs before feeding the upper crust. When the granitoids arrive at the solidified lids of the diapirs, they would favour migrating horizontally and intrude into the upper crust through weakening zones between the diapirs. Our geodynamic model also suggests an asymmetry of structures between the upper and middle-lower crust, with the dome-like granitoids overlying high-grade sedimentary synforms and high-grade diapirs underlying low-grade greenstone belts.

The Earth’s continental crust represents the outermost envelope of the solid Earth, controlling exchanges within the geosphere and reflecting geodynamics processes. One of the fundamental issues of Earth Science aims to determine crustal thickness in past geodynamic environments in order to discuss the evolution of certain geodynamic processes through time. Despite presenting a continuing challenge, the evolution of crustal thickness during the last 3 billion years can be investigated using indirect clues yielded by the chemical signature of mafic magmas and associated ferromagnesian minerals (pyroxene, amphibole). Here, we present a new statistical assessment of a global database of magmatic and mineral chemical information. Analysis reveals the increasing occurrence of high-temperature pyroxenes and amphiboles growing in Ca-rich, Fe-poor magma since ~1 Ga, which contrasts with lower temperature conditions of minerals crystallization throughout the Meso- and Palaeoproterozoic times. This is interpreted to reflect temporal changes in the control of Earth’s crust on mantle-derived magma composition, related to changes in lithospheric thickness and mantle secular cooling. We propose that thick existing crust is associated with deeper, hotter magmatic reservoirs, potentially elucidating the mineral chemistry and the contrasting iron content between primary and derivative mafic magmas. Based on both the chemical and mineral information of mafic magma, an integrated approach provides qualitative estimates of past crustal thickness and associated magmatic systems. Our findings indicate that the Proterozoic was characterized by thicker crustal sections (>40–50 km) relative to the Phanerozoic and Archean (<35 km). This period of crustal thickening appears at the confluence of major changes on Earth, marked by the onset of mantle cooling and Plate Tectonics and the assembly of Columbia, the first supercontinent.