<div>Building on the national-scale thermochemical tomographic inversions of Haynes & Afonso (2023), we infer regions of subduction-driven metasomatic alteration within the Australia sub-continental lithospheric mantle. Such regions are inferred on the basis of age-corrected magnesium number anomalies for bulk composition of the lithospheric mantle, and the spatial correlation of these features with electrical conductors. This defines a mappable criteria for mineral system conceptual models focused on the transport of melts from re-enriched upper mantle sources. Mapping this feature through stochastic uncertainty propogation of inferred mantle compositions enables us to quantify the level of agreement in the spatial constraints on the feature. Here, we present a voting map that quantifies the relative presence or absence of such features across Australia under any arbitrary model realisation.</div><div><br></div><div>Geoscience Australia's Exploring for the Future program provides precompetitive information to inform decision-making by government, community and industry on the sustainable development of Australia's mineral, energy and groundwater resources. By gathering, analysing and interpreting new and existing precompetitive geoscience data and knowledge, we are building a national picture of Australia's geology and resource potential. This leads to a strong economy, resilient society and sustainable environment for the benefit of all Australians. This includes supporting Australia's transition to a low emissions economy, strong resources and agriculture sectors, and economic opportunities and social benefits for Australia's regional and remote communities. The Exploring for the Future program, which commenced in 2016, is an eight year, $225m investment by the Australian Government.</div>

<div>Raster datasets of inferred magnesium number for the bulk lithospheric mantle across the Australian continent. The magnesium number is an elemental ratio, defined by Mg / (Mg + Fe), which correlates to the relative enrichment or depletion in incompatible elements. Oxide concentrations are inferred in from thermo-chemical inverse modelling of Rayleigh phase velocities, surface heat flow, geoid anomalies, and topography. The magnesium number rasters summarise the results of a Markov-chain Monte Carlo sampling of the posterior model space from an ensemble of plausible candidate models. Model release 'FR23' is developed using primary-mode Rayleigh phase velocity grids adapted from Fishwick & Rawlinson (2012; "3-D structure of the Australian lithosphere from evolving seismic datasets"). Model release 'KY22' is developed using the primary-mode Rayleigh phase velocity grids of Yoshizawa (2014; "Radially anisotropic 3-D shear wave structure of the Australian lithosphere and asthenosphere from multi-mode surface waves"). All models are products of the Exploring for the Future program.</div>

<div>The lithology, geochemistry, and architecture of the continental lithospheric mantle (CLM) underlying the Kimberley Craton of north-western Australia has been constrained using pressure-temperature estimates and mineral compositions for &gt;5,000 newly analyzed and published garnet and chrome (Cr) diopside mantle xenocrysts from 25 kimberlites and lamproites of Mesoproterozoic to Miocene age. Single-grain Cr diopside paleogeotherms define lithospheric thicknesses of 200–250 km and fall along conductive geotherms corresponding to a surface heat flow of 37–40 mW/m 2. Similar geotherms derived from Miocene and Mesoproterozoic intrusions indicate that the lithospheric architecture and thermal state of the CLM has remained stable since at least 1,000 Ma. The chemistry of xenocrysts defines a layered lithosphere with lithological and geochemical domains in the shallow (&lt;100 km) and deep (&gt;150 km) CLM, separated by a diopside-depleted and seismically slow mid-lithosphere discontinuity (100–150 km). The shallow CLM is comprised of Cr diopsides derived from depleted garnet-poor and spinel-bearing lherzolite that has been weakly metasomatized. This layer may represent an early (Meso to Neoarchean?) nucleus of the craton. The deep CLM is comprised of high Cr2O3 garnet lherzolite with lesser harzburgite, and eclogite. The peridotite components are inferred to have formed as residues of polybaric partial mantle melting in the Archean, whereas eclogite likely represents former oceanic crust accreted during Paleoproterozoic subduction. This deep CLM was metasomatized by H2O-rich melts derived from subducted sediments and high-temperature FeO-TiO2 melts from the asthenosphere.</div><div><br></div><div>Geoscience Australia’s Exploring for the Future program provides precompetitive information to inform decision-making by government, community and industry on the sustainable development of Australia's mineral, energy and groundwater resources. By gathering, analysing and interpreting new and existing precompetitive geoscience data and knowledge, we are building a national picture of Australia’s geology and resource potential. This leads to a strong economy, resilient society and sustainable environment for the benefit of all Australians. This includes supporting Australia’s transition to net zero emissions, strong, sustainable resources and agriculture sectors, and economic opportunities and social benefits for Australia’s regional and remote communities. The Exploring for the Future program, which commenced in 2016, is an eight year, $225m investment by the Australian Government.</div><div><br></div><div><strong>Citation:</strong></div><div>Sudholz, Z.J., et al. (2023) Mapping the Structure and Metasomatic Enrichment of the Lithospheric Mantle Beneath the Kimberley Craton, Western Australia,&nbsp;<em><i>Geochemistry, Geophysics, Geosystems</i>,</em>&nbsp;24, e2023GC011040.</div><div>https://doi.org/10.1029/2023GC011040</div>