<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)

<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.

Plutonium (Pu) interactions in the environment are highly complex. Site-specific variables play an integral role in determining the chemical and physical form of Pu, and its migration, bioavailability, and immobility. This paper aims to identify the key variables that can be used to highlight regions of radioecological sensitivity and guide remediation strategies in Australia. Plutonium is present in the Australian environment as a result of global fallout and the British nuclear testing program of 1952 – 1958 in central and west Australia (Maralinga and Monte Bello islands). We report the first systematic measurements of 239+240Pu and 238Pu activity concentrations in distal (≥1,000 km from test sites) catchment outlet sediments from Queensland, Australia. The average 239+240Pu activity concentration was 0.29 mBq.g -1 (n = 73 samples) with a maximum of 4.88 mBq.g -1. 238Pu/239+240Pu isotope ratios identified a large range (0.02 – 0.29 (RSD: 74%)) which is congruent with the heterogeneous nuclear material used for the British nuclear testing programme at Maralinga and Montebello Islands. The use of a modified PCA relying on non-linear distance correlation (dCorr) provided broader insight into the impact of environmental variables on the transport and migration of Pu in this soil system. Primary key environmental indicators of Pu presence were determined to be actinide/lanthanide/heavier transition metals, elevation, electrical conductivity (EC), CaO, SiO2, SO3, landform, geomorphology, land use, and climate explaining 81.7% of the variance of the system. Overall this highlighted that trace level Pu accumulations are associated with the coarse, refractive components of Australian soils, and are more likely regulated by the climate of the region and overall soil type. <b>Citation:</b> Megan Cook, Patrice de Caritat, Ross Kleinschmidt, Joёl Brugger, Vanessa NL. Wong, Future migration: Key environmental indicators of Pu accumulation in terrestrial sediments of Queensland, Australia,<i> Journal of Environmental Radioactivity</i>, Volumes 223–224, 2020, 106398,ISSN 0265-931X, https://doi.org/10.1016/j.jenvrad.2020.106398

<div>The soil gas database table contains publicly available results from Geoscience Australia's organic geochemistry (ORGCHEM) schema and supporting oracle databases for gas analyses undertaken by Geoscience Australia's laboratory on soil samples taken from shallow (down to 1 m below the surface) percussion holes. Data includes the percussion hole field site location, sample depth, analytical methods and other relevant metadata, as well as the molecular and isotopic compositions of the soil gas with air included in the reported results. Acquisition of the molecular compounds are by gas chromatography (GC) and the isotopic ratios by gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS). The concentrations of argon (Ar), carbon dioxide (CO₂), nitrogen (N₂) and oxygen (O₂) are given in mole percent (mol%). The concentrations of carbon monoxide (CO), helium (He), hydrogen (H₂) and methane (C₁, CH₄) are given in parts per million (ppm). Compound concentrations that are below detection limit (BDL) are reported as the value -99999. The stable carbon (¹³C/¹²C) and nitrogen (¹⁵N/¹⁴N) isotopic ratios are presented in parts per mil (‰) and in delta notation as δ¹³C and δ¹⁵N, respectively.</div><div><br></div><div>Determining the individual sources and migration pathways of the components of natural gases found in the near surface are useful in basin analysis with derived information being used to support exploration for energy resources (petroleum and hydrogen) and helium in Australian provinces. These data are collated from Geoscience Australia records with the results being delivered in the Soil Gas web services on the Geoscience Australia Data Discovery portal at https://portal.ga.gov.au which will be periodically updated.</div>

<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>.

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 (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>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>The Exploring for the Future program, led by Geoscience Australia, was a $225 million Australian Government investment over 8 years, focused on revealing Australia’s mineral, energy, and groundwater potential by characterising geology.&nbsp;&nbsp;This report provides an overview of activities, results, achievements and impacts from the Exploring for the Future program, with a particular focus on the last four years (2020-2024). &nbsp;</div>

<b>Please Note:</b> The data related to this Abstract can be obtained by contacting <a href = "mailto: clientservices@ga.gov.au">Manager Client Services</a> and quoting Catalogue number 144231. The data are arranged by regions, so please download the Data Description document found in the Downloads tab to determine your area of interest. Remotely sensed datasets provide fundamental information for understanding the chemical, physical and temporal dynamics of the atmosphere, lithosphere, biosphere and hydrosphere. Satellite remote sensing has been used extensively in mapping the nature and characteristics of the terrestrial land surface, including vegetation, rock, soil and landforms, across global to local-district scales. With the exception of hyper-arid regions, mapping rock and soil from space has been problematic because of vegetation that either masks the underlying substrate or confuses the spectral signatures of geological materials (i.e. diagnostic mineral spectral features), making them difficult to resolve. As part of the Exploring for the Future program, a new barest earth Landsat mosaic of the Australian continent using time-series analysis significantly reduces the influence of vegetation and enhances mapping of soil and exposed rock from space. Here, we provide a brief background on geological remote sensing and describe a suite of enhanced images using the barest earth Landsat mosaic for mapping surface mineralogy and geochemistry. These geological enhanced images provide improved inputs for predictive modelling of soil and rock properties over the Australian continent. In one case study, use of these products instead of existing Landsat TM band data to model chromium and sodium distribution using a random forest machine learning algorithm improved model performance by 28–46%.