Data in the GEOCHEM database comprises inorganic geochemical analytical data and associated metadata. Geochemical data comprises concentration data (value, error, unit of measure) measured on a range of analytical instruments, for a range of elements of the periodic table. Associated metadata includes information on analytical techniques, analytical methodology, laboratory, analysts, date of analysis, detection limits, accuracy, and precision. The GEOCHEM database also records results for reference standards. Data is specifically for rocks, soils and other unconsolidated geological material and does not include oils, gases or water analyses. Geochemical data may be total rock (i.e., whole rock analysed) or for a variety of fractions of the total rock, e.g., various non-total acid digests, mineral separates, differing size fractions. It also includes quantitative to semi-quantitative data from field measurements, such as portable x-ray fluorescence (XRF). It does not include geochemical data for individual minerals. <b>Value: </b>Geochemical data underpins much geoscientific study, and is used directly to classify, characterise and understand geological material and its formation. It has direct relevance to understanding the formation of the earth, the continents, and the processes that create and shape the surface we live on. For example, this information is used within: both discovering and the understanding of mineral deposits we depend on; the nature, health and sustainability of the soils we live and farm on; as well as providing input into a range of potential geohazards. <b>Scope: </b>The collection includes data from over 60 years of Geoscience Australia (GA) and state/territory partner regional geological projects within Australia, as well as continental-scale and regional geochemical surveys like National Geochemical Survey of Australia (NGSA) and Northern Australia Geochemical Survey (NAGS) (Exploring for the Future- EFTF). It also includes data from other countries that GA has worked with, e.g., Papua New Guinea, Antarctica, Solomon Islands and New Zealand. Explore the <b>Geoscience Australia portal - <a href="https://portal.ga.gov.au/">https://portal.ga.gov.au/</a></b>

A `weighted geometric median’ approach has been used to estimate the median surface reflectance of the barest state (i.e., least vegetation) observed through Landsat-8 OLI observations from 2013 to September 2018 to generate a six-band Landsat-8 Barest Earth pixel composite mosaic over the Australian continent. The bands include BLUE (0.452 - 0.512), GREEN (0.533 - 0.590), RED, (0.636 - 0.673) NIR (0.851 - 0.879), SWIR1 (1.566 - 1.651) and SWIR2 (2.107 - 2.294) wavelength regions. The weighted median approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. The product reduces the influence of vegetation and allows for more direct mapping of soil and rock mineralogy. Reference: Dale Roberts, John Wilford, and Omar Ghattas (2018). Revealing the Australian Continent at its Barest, submitted.

An estimate of the spectra of the barest state (i.e., least vegetation) observed from imagery of the Australian continent collected by the Landsat 5, 7, and 8 satellites over a period of more than 30 years (1983 – 2018). The bands include BLUE (0.452 - 0.512), GREEN (0.533 - 0.590), RED, (0.636 - 0.673) NIR (0.851 - 0.879), SWIR1 (1.566 - 1.651) and SWIR2 (2.107 - 2.294) wavelength regions. The approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. The product reduces the influence of vegetation and allows for more direct mapping of soil and rock mineralogy. This product complements the Landsat-8 Barest Earth which is based on the same algorithm but just uses Landsat8 satellite imagery from 2013-2108. Landsat-8’s OLI sensor provides improved signal-to-noise radiometric (SNR) performance quantised over a 12-bit dynamic range compared to the 8-bit dynamic range of Landsat-5 and Landsat-7 data. However the Landsat 30+ Barest Earth has a greater capacity to find the barest ground due to the greater temporal depth. Reference: Exposed Soil and Mineral Map of the Australian Continent Revealing the Land at its Barest - Dale Roberts, John Wilford and Omar Ghattas Ghattas (2019). Nature Communications, DOI: 10.1038/s41467-019-13276-1. https://www.nature.com/articles/s41467-019-13276-1

A `weighted geometric median' approach has been used to estimate the median surface reflectance of the barest state (i.e., least vegetation) observed through Landsat-8 Operational Land Image (OLI) observations from 2013 to September 2018 to generate a six-band Landsat-8 Barest Earth pixel composite mosaic over the Australian continent. The bands include BLUE (0.452 - 0.512), GREEN (0.533 - 0.590), RED, (0.636 - 0.673) NIR (0.851 - 0.879), SWIR1 (1.566 - 1.651) and SWIR2 (2.107 - 2.294) wavelength regions. The weighted median approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. The product reduces the influence of vegetation and allows for more direct mapping of soil and rock mineralogy. Reference: Dale Roberts, John Wilford, and Omar Ghattas (2018). Revealing the Australian Continent at its Barest, submitted. <b>Value: </b>Has broad application in mapping surface geochemistry and mineralogy of exposed soil and bedrock. Has applications in geological mapping and natural resource management including mapping of soil characteristics. <b>Scope: </b>Two enhanced bare earth products have been generated reflecting different Landsat satellites and acquisition periods. The first only uses Landsat 8 observations from 2013 to 2018. The second incorporates the full 30+ year archive combining Landsat 5, 7, and 8 from 1986 to 2018.

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%. <b>Citation:</b> Wilford, J. and Roberts, D., 2020. Enhanced barest earth Landsat imagery for soil and lithological modelling. In: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A. and Slatter, E. (eds.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, 1–4.

The Sentinel-2 Bare Earth thematic product provides the first national scale mosaic of the Australian continent to support improved mapping of soil and geology. The bare earth algorithm using all available Sentinel-2 A and Sentinel-2 B observations up to September 2020 preferentially weights bare pixels through time to significantly reduce the effect of seasonal vegetation in the imagery. The result are image pixels that are more likely to reflect the mineralogy and/or geochemistry of soil and bedrock. The algorithm uses a high-dimensional weighted geometric median approach that maintains the spectral relationships across all Sentinel-2 bands. A similar bare earth algorithm has been applied to Geoscience Australia’s deeper Landsat time series archive (please search for "Landsat barest Earth". Both bare earth products have spectral bands in the visible near infrared and shortwave infrared region of the electromagnetic spectrum. However, the main visible and near-infrared Sentinel-2 bands have a spatial resolution of 10 meters compared to 30m for the Landsat TM equivalents. The weighted median approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. Not all the sentinel-2 bands have been processed - we have excluded atmospheric bands including 1, 9 and 10. The remaining bands have been re-number 1-10 and these bands correlate to the original bands in brackets below: 1 = blue (2) , 2 = green (3) , 3 = red (4), 4 = vegetation red edge (5), 5 = vegetation red edge (6), 6= vegetation red edge (7), 7 = NIR(8), 8 = Narrow NIR (8a), 9 = SWIR1 (11) and 10 = SWIR2(12). All 10 bands have been resampled to 10 meters to facilitate band integration and use in machine learning.

During the last 10-20 years, Geological Surveys around the world have undertaken a major effort towards delivering fully harmonized and tightly quality controlled low-density multi-element soil geochemical maps and datasets of vast regions including up to whole continents. Concentrations of between 45 and 60 elements commonly have been determined in a variety of different regolith types (e.g., sediment, soil). The multi-element datasets are published as complete geochemical atlases and made available to the general public. Several other geochemical datasets covering smaller areas but generally at a higher spatial density are also available. These datasets may, however, not be found by superficial internet-based searches because the elements are not mentioned individually either in the title or in the keyword lists of the original references. This publication attempts to increase the visibility and discoverability of these fundamental background datasets covering large areas up to whole continents. <b>Citation:</b> P. de Caritat, C. Reimann, D.B. Smith, X. Wang, Chemical elements in the environment: Multi-element geochemical datasets from continental- to national-scale surveys on four continents, <i>Applied Geochemistry</i>, Volume 89, 2018, Pages 150-159, ISSN 0883-2927, https://doi.org/10.1016/j.apgeochem.2017.11.010

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

A predictive model of weathering intensity or the degree of weathering has been generate over the Australian continent. The model has been generated using the Random Forest decision tree machine learning algorithm. The algorithm is used to establish predictive relationships between field estimates of the degree of weathering and a comprehensive suite of covariate or predictive datasets. The covariates used to generate the model include satellite imagery, terrain attributes, airborne radiometric imagery and mapped geology. The weathering intensity model is an estimate of the degree of surface weathering only. The interpretation of the weathering intensity is different for in-situ or residual landscapes compared with transported materials within depositional landscapes. In residual landscapes, weathering process are operating locally whereas in depositional landscapes the model is reflecting the degree of weathering either prior to erosion and subsequent deposition, or weathering of sediments after being deposited. The degree of surface weathering is particularly important in Australia where variations in weathering intensity correspond to the nature and distribution of regolith (weathered bedrock and sediments) which mantles approximately 90% of the Australian continent. The weathering intensity prediction has been generated using the Random Forest decision tree machine learning algorithm. The algorithm is used to establish predictive relationships between field estimates of the degree of weathering and a comprehensive suite of covariate or predictive datasets. The covariates used to generate the model include satellite imagery, terrain attributes, airborne radiometric imagery and mapped geology. Correlations between the training dataset and the covariates were explored through the generation of 300 random tree models. An r-squared correlation of 0.85 is reported using 5 K-fold cross-validation. The mean of the 300 models is used for predicting the weathering intensity and the uncertainty in the weathering intensity is estimated at each location via the standard deviation in the 300 model values. The predictive weathering intensity model is an estimate of the degree of surface weathering only. The interpretation of the weathering intensity is different for in-situ or residual landscapes compared with transported materials within depositional landscapes. In residual landscapes, weathering process are operating locally whereas in depositional landscapes the model is reflecting the degree of weathering either prior to erosion and subsequent deposition, or weathering of sediments after being deposited. The weathering intensity model has broad utility in assisting mineral exploration in variably weathered geochemical landscapes across the Australian continent, mapping chemical and physical attributes of soils in agricultural landscapes and in understanding the nature and distribution of weathering processes occurring within the upper regolith. <b>Value: </b>Weathering intensity is an important characteristic of the earth's surface that has a significant influence on the chemical and physical properties of surface materials. Weathering intensity largely controls the degree to which primary minerals are altered to secondary components including clay minerals and oxides. In this context the weathering intensity model has broad application in understanding geomorphological and weathering processes, mapping soil/regolith and geology. <b>Scope: </b>National dataset which over time can be improved with additional sites for training and thematic datasets for prediction.