<div>Strontium isotopes (87Sr/86Sr) are useful to trace processes in the Earth sciences as well as in forensic, archaeological, palaeontological, and ecological sciences. As very few large-scale Sr isoscapes exist in Australia, we have identified an opportunity to determine 87Sr/86Sr ratios on archived fluvial sediment samples from the low-density National Geochemical Survey of Australia (www.ga.gov.au/ngsa; last access: 15 December 2022). The present study targeted the northern parts of Western Australia, the Northern Territory and Queensland, north of 21.5 °S. The samples were taken mostly from a depth of ~60-80 cm in floodplain deposits at or near the outlet of large catchments (drainage basins). A coarse (< 2 mm) grain-size fraction was air-dried, sieved, milled then digested (hydrofluoric acid + nitric acid followed by aqua regia) to release <em>total</em> Sr. The Sr was then separated by chromatography and the 87Sr/86Sr ratio determined by multicollector-inductively coupled plasma mass spectrometry. Results demonstrate a wide range of Sr isotopic values (0.7048 to 1.0330) over the survey area, reflecting a large diversity of source rock lithologies, geological processes and bedrock ages. Spatial distribution of 87Sr/86Sr shows coherent (multi-point anomalies and smooth gradients), large-scale (> 100 km) patterns that appear to be broadly consistent with surface geology, regolith/soil type, and/or nearby outcropping bedrock. For instance, the extensive black clay soils of the Barkly Tableland define a > 500 km-long northwest-southeast-trending unradiogenic anomaly (87Sr/86Sr < 0.7182). Where sedimentary carbonate or mafic/ultramafic igneous rocks dominate, low to moderate 87Sr/86Sr values are generally recorded (medians of 0.7387 and 0.7422, respectively). In proximity to the outcropping Proterozoic metamorphic basement of the Tennant, McArthur, Murphy and Mount Isa geological regions, conversely, radiogenic 87Sr/86Sr values (> 0.7655) are observed. A potential correlation between mineralisation and elevated 87Sr/86Sr values in these regions needs to be investigated in greater detail. Our results to-date indicate that incorporating soil/regolith Sr isotopes in regional, exploratory geoscience investigations can help identify basement rock types under (shallow) cover, constrain surface processes (e.g. weathering, dispersion), and, potentially, recognise components of mineral systems. Furthermore, the resulting Sr isoscape and future models derived therefrom can also be utilised in forensic, archaeological, paleontological and ecological studies that aim to investigate, e.g., past and modern animal (including humans) dietary habits and migrations. The new spatial Sr isotope dataset for the northern Australia region is publicly available (de Caritat et al., 2022a; https://dx.doi.org/10.26186/147473; last access: 15 December 2022).</div> <b>Citation:</b> de Caritat, P., Dosseto, A., and Dux, F.: A strontium isoscape of northern Australia, <i>Earth Syst. Sci. Data</i>, 15, 1655–1673, https://doi.org/10.5194/essd-15-1655-2023, <b>2023</b>.

<div>Strontium isotopes (87Sr/86Sr) are useful in the earth sciences (e.g. recognising geological provinces, studying geological processes) as well in archaeological (e.g. informing on past human migrations), palaeontological/ecological (e.g. investigating extinct and extant taxa’s dietary range and migrations) and forensic (e.g. validating the origin of drinks and foodstuffs) sciences. Recently, Geoscience Australia and the University of Wollongong have teamed up to determine 87Sr/86Sr ratios in fluvial sediments selected mostly from the low-density National Geochemical Survey of Australia (www.ga.gov.au/ngsa), with a few additional Northern Australia Geochemical Survey infill samples. The present study targeted the northern parts of Western Australia, the Northern Territory and Queensland in Australia, north of 21.5 °S. The samples were taken mostly from a depth of ~60-80 cm depth in floodplain deposits at or near the outlet of large catchments (drainage basins). A coarse grain-size fraction (&lt;2 mm) was air-dried, sieved, milled then digested (hydrofluoric acid + nitric acid followed by aqua regia) to release total strontium. Preliminary results demonstrate a wide range of strontium isotopic values (0.7048 &lt; 87Sr/86Sr &lt; 1.0330) over the survey area, reflecting a large diversity of source rock lithologies, geological processes and bedrock ages. Spatial distribution of 87Sr/86Sr shows coherent (multi-point anomalies and smooth gradients), large-scale (&gt;100 km) patterns that appears to be consistent, in many places, with surface geology, regolith/soil type and/or nearby outcropping bedrock. For instance, the extensive black clay soils of the Barkly Tableland define a &gt;500 km-long northwest-southeast-trending low anomaly (87Sr/86Sr &lt; 0.7182). Where carbonate or mafic igneous rocks dominate, a low to moderate strontium isotope signature is observed. In proximity to the outcropping Proterozoic metamorphic provinces of the Tennant, McArthur, Murphy and Mount Isa geological regions, conversely, high 87Sr/86Sr values (&gt; 0.7655) are observed. A potential link between mineralisation and elevated 87Sr/86Sr values in these regions needs to be investigated in greater detail. Our results to-date indicate that incorporating soil/regolith strontium isotopes in regional, exploratory geoscience investigations can help identify basement rock types under (shallow) cover, constrain surface processes (e.g. weathering, dispersion), and, potentially, recognise components of mineral systems. Furthermore, the resulting strontium isoscape and model derived therefrom can also be utilised in archaeological, paleontological and ecological studies that aim to investigate past and modern animal (including humans) dietary habits and migrations. &nbsp;The new spatial dataset is publicly available through the Geoscience Australia portal https://portal.ga.gov.au/.</div>

<div>Soil is a complex and spatially variable material that has a demonstrated potential to be a useful evidence class in forensic casework and intelligence operations. Here, the capability to spatially constrain searches and prioritise resources by triaging areas as low and high interest is advantageous. Conducted between 2017 and 2021, a forensically relevant topsoil survey (0-5 cm depth; 1 sample per 1 km2) has been carried out over Canberra, Australia, with the aims of documenting the distribution of chemical elements in an urban/suburban environment, and of acting as a testbed for investigating various aspects of forensic soil provenancing. Geochemical data from X-Ray Fluorescence (XRF; for total major oxides) and Inductively Coupled Plasma-Mass Spectrometry (ICP-MS; for trace elements) following a total digestion (HF + HNO3) were obtained from the survey’s 685 topsoil samples (plus 138 additional quality control samples and six “Blind” simulated evidentiary samples). Using those “Blind” samples, we document a likelihood ratio approach where for each grid cell the analytical similarity between the grid cell and evidentiary sample is attributed from a measure of overlap between both Cauchy distributions, including appropriate uncertainties. Unlike existing methods that base inclusion/exclusion on an arbitrary threshold (e.g., ± three standard deviations), our approach is free from strict binary or Boolean thresholds, providing an unconstrained gradual transition dictated by the analytical similarity. Using this provenancing model, we present and evaluate a new method for upscaling from a fine (25 m x 25 m) interpolated grid to a more appropriate coarser (500 m x 500 m) grid, in addition to an objective method using Random Match Probabilities for ranking individual variables to be used for provenancing prior to receiving evidentiary material. Our results show this collective procedure generates more consistent and robust provenance maps between two different interpolation algorithms (e.g., inverse distance weighting, and natural neighbour), grid placements (e.g., grid shifts to the north or east) and theoretical users (e.g., different computer systems, or forensic geoscientists).</div> <b>Citation:</b> Michael G. Aberle, Patrice de Caritat, James Robertson, Jurian A. Hoogewerff, A robust interpolation-based method for forensic soil provenancing: A Bayesian likelihood ratio approach,<i> Forensic Science International</i>, Volume 353, 2023, 111883, ISSN 0379-0738. https://doi.org/10.1016/j.forsciint.2023.111883.

Major oxides provide valuable information about the composition, origin, and properties of rocks and regolith. Analysing major oxides contributes significantly to understanding the nature of geological materials and processes (i.e. physical and chemical weathering) – with potential applications in resource exploration, engineering, environmental assessments, agriculture, and other fields. Traditionally most measurements of oxide concentrations are obtained by laboratory assay, often using X-ray fluorescence, on rock or regolith samples. To expand beyond the point measurements of the geochemical data, we have used a machine learning approach to produce seamless national scale grids for each of the major oxides. This approach builds predictive models by learning relationships between the site measurements of an oxide concentration (sourced from Geoscience Australia’s OZCHEM database and selected sites from state survey databases) and a comprehensive library of covariates (features). These covariates include: terrain derivatives; climate surfaces; geological maps; gamma-ray radiometric, magnetic, and gravity grids; and satellite imagery. This approach is used to derive national predictions for 10 major oxide concentrations at the resolution of the covariates (nominally 80 m). The models include the oxides of silicon (SiO2), aluminium (Al2O3), iron (Fe2O3tot), calcium (CaO), magnesium (MgO), manganese (MnO), potassium (K2O), sodium (Na2O), titanium (TiO2), and phosphorus (P2O5). The grids of oxide concentrations provided include the median of multiple models run as the prediction, and lower and upper (5th and 95th) percentiles as measures of the prediction’s uncertainty. Higher uncertainties correlate with greater spreads of model values. Differences in the features used in the model compared with the full feature space covering the entire continent are captured in the ‘covariate shift’ map. High values in the shift model can indicate higher potential uncertainty or unreliability of the model prediction. Users therefore need to be mindful, when interpreting this dataset, of the uncertainties shown by the 5th-95th percentiles, and high values in the covariate shift map. Details of the modelling approach, model uncertainties and datasets are describe in an attached word document “Model approach uncertainties”. This work is part of Geoscience Australia’s Exploring for the Future program that 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. These data are published with the permission of the CEO, Geoscience Australia.

<div>Strontium isotopes (87Sr/86Sr) are useful in the earth sciences (e.g. recognising geological provinces, studying geological processes) as well in archaeological (e.g. informing on past human migrations), palaeontological/ecological (e.g. investigating extinct and extant taxa’s dietary range and migrations) and forensic (e.g. validating the origin of drinks and foodstuffs) sciences. Recently, Geoscience Australia and the University of Wollongong have teamed up to determine 87Sr/86Sr ratios in fluvial sediments selected mostly from the low-density National Geochemical Survey of Australia (NGSA; www.ga.gov.au/ngsa). The present study targeted the Yilgarn geological region in southwestern Australia. The samples were mostly taken from a depth of ~60-80 cm (Bottom Outlet Sediments, BOS) in floodplain deposits at or near the outlet of large catchments (drainage basins). A small number of surface (0-10 cm) samples (Top Outlet Sediments, TOS) were also included in the study. For all, a coarse grain-size fraction (<2 mm) was air-dried, sieved, milled then digested (hydrofluoric acid + nitric acid followed by aqua regia) to release total strontium. Overall, 107 NGSA BOS < 2 mm and 13 NGSA TOS < 2 mm were analysed for Sr isotopes. Given that there are ~10 % field duplicates in the NGSA, all those samples originate from within 97 NGSA catchments, which together cover 533 000 km2 of southwestern Australia. Preliminary results for the BOS samples demonstrate a wide range of strontium isotopic values (0.7152 < 87Sr/86Sr < 1.0909) over the survey area, reflecting a large diversity of source rock lithologies, geological processes and bedrock ages. Spatial distribution of 87Sr/86Sr shows coherent (multi-point anomalies and smooth gradients), large-scale (>100 km) patterns that appear to be consistent, in many places, with surface geology, regolith/soil type and/or nearby outcropping bedrock. For instance, catchments in the western and central Yilgarn dominated by felsic intrusive basement geology have radiogenic 87Sr/86Sr signatures in the floodplain sediments consistent with published whole-rock data. Similarly, unradiogenic signatures in sediments in the eastern Yilgarn are in agreement with published whole-rock data. Our results to-date indicate that incorporating soil/regolith strontium isotopes in regional, exploratory geoscience investigations can help identify basement rock types under (shallow) cover, constrain surface processes (e.g. weathering, dispersion), and, potentially, recognise components of mineral systems. Furthermore, the resulting strontium isoscape and model derived therefrom can also be utilised in archaeological, paleontological and ecological studies that aim to investigate past and modern animal (including humans) dietary habits and migrations.&nbsp; The new spatial dataset is publicly available through the Geoscience Australia portal https://portal.ga.gov.au/.</div>

<div>With a higher demand for lithium (Li), a better understanding of its concentration and spatial distribution is important to delineate potential anomalous areas. This study uses a digital soil mapping framework to combine data from recent geochemical surveys and environmental covariates to predict and map Li content across the 7.6 million km2 area of Australia. Soil samples were collected by the National Geochemical Survey of Australia at a total of 1315 sites, with both top (0–10 cm depth) and bottom (on average 60–80 cm depth) catchment outlet sediments sampled. We developed 50 bootstrap models using a Cubist regression tree algorithm for both depths. The spatial prediction models were validated on an independent Northern Australia Geochemical Survey dataset, showing a good prediction with an RMSE of 3.82 mg kg-1 for the top depth. The model for the bottom depth has yet to be validated. The variables of importance for the models indicated that the first three Landsat bands and gamma radiometric dose have a strong impact on Li prediction. The bootstrapped models were then used to generate digital soil Li prediction maps for both depths, which could select and delineate areas with anomalously high Li concentrations in the regolith. The map shows high Li concentration around existing mines and other potentially anomalous Li areas. The same mapping principles can potentially be applied to other elements.&nbsp;</div> <b>Citation:</b> Ng, W., Minasny, B., McBratney, A., de Caritat, P., and Wilford, J.: Digital soil mapping of lithium in Australia, <i>Earth Syst. Sci. Data</i>, 15, 2465–2482, https://doi.org/10.5194/essd-15-2465-2023, <b>2023</b>.

<div>Alluvial sediments have long been used in geochemical surveys as their compositions are assumed to be representative of areas upstream. Overbank and floodplain sediments, in particular, are increasingly used for regional to continental-scale geochemical mapping. However, during downstream transport, sediments from heterogeneous source regions are carried away from their source regions and mixed. Consequently, using alluvial sedimentary geochemical data to generate continuous geochemical maps remains challenging. In this study we demonstrate a technique that numerically unmixes alluvial sediments to make a geochemical map of their upstream catchments. The unmixing approach uses a model that predicts the concentration of elements in downstream sediments, given a map of the drainage network and element concentrations in the source region. To unmix sedimentary chemistry, we seek the upstream geochemical map that, when mixed downstream, best fits geochemical observations downstream. To prevent overfitting we penalise the roughness of the geochemical model. To demonstrate our approach we apply it to alluvial samples gathered as part of the Northern Australia Geochemical Survey. This survey gathered samples collected over a ∼ 500,000 km2 area in northern Australia. We first validate our approach for this sample distribution with synthetic tests, which indicate that we can resolve geochemical variability at scales greater than 0.5 – 1◦ in size. We proceed to invert real geochemical data from the total digestion of fine-grained fraction of alluvial sediments. The resulting geochemical maps for two elements of potential economic interest, Cu and Nd, are evaluated in detail. We find that in both cases, our predicted downstream concentrations match well against a held-out, unseen subset of the data, as well as against data from an independent geochemical survey. By performing principal component analysis on maps generated for all 46 available elements we produce a synthesis map showing the significant geochemical domains of this part of northern Australia. This map shows strong spatial similarities to the underlying lithological map of the area. Finally, we compare the results from our approach to a geochemical map produced by kriging. We find that, unlike the method presented here, kriging generates geochemical maps that are both dampened relative to expected magnitude, as well as being spatially distorted. We argue that the unmixing approach is the most appropriate method for generating geochemical maps from regional-scale alluvial surveys.&nbsp;</div> <b>Citation:</b> Alex G. Lipp, Patrice de Caritat, Gareth G. Roberts, Geochemical mapping by unmixing alluvial sediments: An example from northern Australia, <i>Journal of Geochemical Exploration,</i> Volume 248, <b>2023</b>, 107174, ISSN 0375-6742, https://doi.org/10.1016/j.gexplo.2023.107174. (https://www.sciencedirect.com/science/article/pii/S0375674223000213)