A combined geodynamical-geochemical modelling approach to investigating the Lu-Hf isotopic evolution of the terrestrial mantle and crust

The chemical and isotopic evolution of the Earth's mantle is largely influenced by the formation of oceanic and continental crust at spreading ridges and through arc volcanism, and the subsequent recycling of this crust back into the mantle via subduction. In this study we use a combined geodynamical-geochemical modelling approach to investigate the Lu-Hf isotopic evolution of the terrestrial mantle and crust. We utilise the geodynamic mantle convection model developed by Brandenburg et al., 2008. This model satisfies the geophysical constraints of oceanic heat flow and average plate velocities, as well as geochemical observations such as 40Ar in the atmosphere. It has also been shown to reproduce the observed geochemical distributions in multiple isotope systems (U-Th-Pb, Rb-Sr, Sm-Nd, and Re-Os) that define the DMM, HIMU and EM1 mantle endmembers. We go on to extend this application to investigate the Lu-Hf isotope system, specifically in combination with Sm-Nd. The model has been updated to include a self-consistent reorganisation of the plates with regions of up-/down-wellings. The model is initiated at 4.55 Ga, assumes continental crust is produced from 4 Ga and that a transition from `dry' to `wet' subduction occurs at 3 Ga. The results of the geodynamic model suggest that the ƐHf composition and evolution of the upper mantle can be generated through the extraction and recycling of oceanic crust, which creates an enriched and radiogenic reservoir at the core-mantle boundary. The formation of continental crust, which is extracted at each time-step from the oceanic crust to imitate subduction zone processes, and the recycling of this continental crust as sediments, plays a lesser role. Depending on the selected partition coefficients DMM, FOZO and HIMU mantle endmember compositions are also produced via the simple extraction and recycling of oceanic crust. The formation of continental crust produces spread in the ƐNd vs. ƐHf array and extends the model values into the HIMU region of the terrestrial array. We go on to use this geodynamic-geochemical model to investigate different models of continental growth, by observing the effects on the coupled crustal-mantle reservoirs. Brandenburg, J.P., Hauri, E.H., van Keken, P.E., Ballentine, C.J., 2008. Earth and Planetary Science Letters 276, 1-13.