Ilya Fomin

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 data sets and thermodynamic modeling. 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 data sets. 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.

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.

The Archean Superior craton was formed by the assemblage of continental and oceanic terranes at ∼2.6 Ga. The craton is surrounded by multiple Proterozoic mobile belts, including the Paleoproterozoic Trans‐Hudson Orogen which brought together the Superior and Rae/Hearne cratons at ∼1.9–1.8 Ga. Despite numerous studies on Precambrian lithospheric formation and evolution, the deep thermochemical structure of the Superior craton and its surroundings remains poorly understood. Here we investigate the upper mantle beneath the region from the surface to 400 km depth by jointly inverting Rayleigh wave phase velocity dispersion data, elevation, geoid height and surface heat flow, using a probabilistic inversion to obtain a (pseudo‐)3D model of composition, density and temperature. The lithospheric structure is dominated by thick cratonic roots (>300 km) beneath the eastern and western arms of the Superior craton, with a chemically depleted signature (Mg# > 92.5), consistent with independent results from mantle xenoliths. Beneath the surrounding Proterozoic and Phanerozoic orogens, the Mid‐continent Rift and Hudson Strait, we observe a relatively thinner lithosphere and more fertile composition, indicating that these regions have undergone lithospheric modification and erosion. Our model supports the hypothesis that the core of the Superior craton is well‐preserved and has evaded lithospheric destruction and refertilization. We propose three factors playing a critical role in the craton's stability: (a) the presence of a mid‐lithospheric discontinuity, (b) the correct isopycnic conditions to sustain a strength contrast between the craton and the surrounding mantle, and (c) the presence of weaker mobile belts around the craton.

The thermochemical structure of the lithosphere controls melting mechanisms in the mantle, as well as the location of volcanism and ore deposits. Obtaining reliable images of the lithosphere structure, and its complex interactions with the asthenosphere, requires the joint inversion of multiple data sets and their associated uncertainties. In particular, the combination of seismic velocity and electrical conductivity, along with proxies for bulk composition and elusive minor phases, represents a crucial step toward fully understanding large‐scale lithospheric structure and melting processes. We apply a novel probabilistic approach for joint inversions of 3D magnetotelluric and seismic data to image the lithosphere beneath southeast Australia. The results show a highly heterogeneous lithosphere 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 (a) boundaries of geological provinces, elucidating the subduction‐accretion process in the region; (b) distribution of leucitite and basaltic magmatism; (c) independent geochemical data, and (d) 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 primary mechanisms for melting in eastern Australia, contrary to the conventional notion of mantle plume activity. Our study presents an integrated lithospheric model for southeastern Australia and provides valuable insight into the mechanisms driving surface geological processes.

The Archean Superior craton was formed by the assemblage of continental and oceanic terranes at ∼2.6 Ga. The craton is surrounded by multiple Proterozoic mobile belts, including the Paleoproterozoic Trans-Hudson Orogen which brought together the Superior and Rae/Hearne cratons at ∼1.9-1.8 Ga. Despite numerous studies on Precambrian lithospheric formation and evolution, the deep thermochemical structure of the Superior craton and its surroundings remains poorly understood. Here we investigate the upper mantle beneath the region from the surface to 400 km depth by jointly inverting Rayleigh wave phase velocity dispersion data, elevation, geoid height and surface heat flow, using a probabilistic inversion to obtain a (pseudo-)3D model of composition, density and temperature. The lithospheric structure is dominated by thick cratonic roots (>300 km) beneath the eastern and western arms of the Superior craton, with a chemically depleted signature (Mg# > 92.5), consistent with independent results from mantle xenoliths. Beneath the surrounding Proterozoic and Phanerozoic orogens, the Mid-continent Rift and Hudson Strait, we observe a relatively thinner lithosphere and more fertile composition, indicating that these regions have undergone lithospheric modification and erosion. Our model supports the hypothesis that the core of the Superior craton is well-preserved and has evaded lithospheric destruction and refertilization. We propose three factors playing a critical role in the craton’s stability: (i) the presence of a mid-lithospheric discontinuity, (ii) the correct isopycnic conditions to sustain a strength contrast between the craton and the surrounding mantle, and (iii) the presence of weaker mobile belts around the craton.

The thermochemical structure of the lithosphere controls melting mechanisms in the mantle, as well as the location of volcanism and ore deposits. Obtaining reliable images of the lithosphere structure, and its complex interactions with the asthenosphere, requires the joint inversion of multiple data sets and their associated uncertainties. In particular, the combination of seismic velocity and electrical conductivity, along with proxies for bulk composition and elusive minor phases, represents a crucial step towards fully understanding large-scale lithospheric structure and melting processes. We apply a novel probabilistic approach for joint inversions of 3D magnetotelluric and seismic data to image the lithosphere beneath southeast Australia. The results show a highly heterogeneous lithosphere 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) boundaries of geological provinces, elucidating 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 primary mechanisms for melting in eastern Australia, contrary to the conventional notion of mantle plume activity. Our study presents an integrated lithospheric model for southeastern Australia and provides valuable insight into the mechanisms driving surface geological processes.

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.

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.

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.