Imaging the Earth’s thermochemical structure is crucial for understanding its dynamics and evolution. Moreover, the increased demand for critical minerals and geothermal energy driven by the energy transition has intensified the need for reliable subsurface models. Multi-Observable Thermochemical Tomography (MTT) is a simulation-based, probabilistic inversion platform designed to harness the combined sensitivities of multiple geophysical datasets and thermodynamic modelling. It produces internally-consistent estimates of the Earth’s interior as probability distributions, offering a powerful means for uncertainty quantification. Here, we present an updated MTT formalism and assess its benefits and limitations to image the thermochemical structure of the lithosphere-asthenosphere system. Individual and combined sensitivities of different observables to parameters of interest (e.g. temperature, composition, crustal architecture) are explored using challenging synthetic models. Our findings demonstrate that a judicious combination of observables can retrieve complex thermochemical structures relevant to greenfields exploration. We then apply MTT to study two cratonic regions of geological and economic significance. In the Superior Craton, we jointly invert receiver functions, gravity anomalies, gravity gradients, geoid anomalies, Rayleigh-wave dispersion curves, absolute elevation and surface heat flow. In the North Australian Craton, we incorporate new data from the Ausarray and add teleseismic P- and S-phase travel times to the datasets. The imaged lithospheric architectures provide new insights into the tectonic evolution of these two regions and the physical meaning of geophysical signatures. Additionally, these models offer unique proxies to guide exploration efforts for clean energy and critical minerals and serve as reference models for future high-resolution studies.

A document by María Constanza Manassero. Click on the document to view its contents.

The thermochemical structure of the lithosphere exerts control on melting mechanisms in the mantle as well as the location of volcanism and ore deposits. Imaging the complex interactions between the lithosphere and asthenospheric mantle requires the joint inversion of multiple data sets and their uncertainties. In particular, the combination of seismic velocity and electrical conductivity with data proxies for bulk composition and elusive minor phases is a crucial step towards fully understanding large-scale lithospheric structure and melting. We apply a novel probabilistic approach for joint inversions of 3D magnetotelluric and seismic data to image the lithosphere beneath southeast Australia. Results show a highly heterogeneous lithospheric structure with deep conductivity anomalies that correlate with the location of Cenozoic volcanism. In regions where the conductivities have been at odds with sub-lithospheric temperatures and seismic velocities, we observe that the joint inversion provides conductivity values consistent with other observations. The results reveal a strong relationship between metasomatized regions in the mantle and i) the limits of geological provinces in the crust, which elucidates the subduction-accretion process in the region; ii) distribution of leucitite and basaltic magmatism; iii) independent geochemical data, and iv) a series of lithospheric steps which constitute areas prone to generating small-scale instabilities in the asthenosphere. This scenario suggests that shear-driven upwelling and edge-driven convection are the dominant melting mechanisms in eastern Australia rather than mantle plume activity, as conventionally conceived. Our study offers an integrated lithospheric model for southeastern Australia and provides insights into the feedback mechanism driving surface processes.