Articles | Volume 15, issue 4
Data description paper
14 Apr 2023
Data description paper |  | 14 Apr 2023

A strontium isoscape of northern Australia

Patrice de Caritat, Anthony Dosseto, and Florian Dux

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. 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, and digested (hydrofluoric acid + nitric acid followed by aqua regia) to release total Sr. The Sr was then separated by chromatography, and the 87Sr/86Sr ratio was determined by multicollector inductively coupled plasma mass spectrometry. The 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. The 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). Conversely, In proximity to the outcropping Proterozoic metamorphic basement of the Tennant, McArthur, Murphy, and Mount Isa geological regions, radiogenic 87Sr/86Sr values (>0.7655) are observed. A potential correlation between mineralization 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 and dispersion), and, potentially, recognize components of mineral systems. Furthermore, the resulting Sr isoscape and future models derived therefrom can also be utilized in forensic, archaeological, palaeontological, and ecological studies that aim to investigate, for example, 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;

1 Introduction

Strontium isotope ratios (87Sr/86Sr) can be measured in many geological materials, as the trace element strontium (Sr) is relatively abundant and readily substitutes for calcium (Ca) in minerals and organic tissues. The 87Sr/86Sr of a mineral or rock is a function of (1) its initial and unchanging 86Sr content, (2) its initial rubidium (Rb) content (thus Rb/Sr ratio), and (3) time (e.g. McNutt, 2000). Rubidium substitutes readily for potassium (K) in minerals and is, thus, also relatively common. As one of the two naturally occurring Rb isotopes, 87Rb, which accounts for 27.8 % of Rb, decays over time (by emitting a negative beta particle) to stable 87Sr (t1/2=49.6×109 years), and the 87Sr/86Sr ratio of that material slowly increases with time (e.g. Rotenberg et al., 2012; Nebel and Stammeier, 2018). During geological processes, such as mineral dissolution or precipitation, or biological processes, such as bone and tooth growth, the 87Sr/86Sr remains constant, as there is no fractionation (e.g. Gosz et al., 1983; Nebel and Stammeier, 2018). These characteristics make the Sr isotopic system very useful in the geosciences, where it has been used for decades, for instance, in the study of the following:

  • sediment diagenesis and low-grade metamorphism (e.g. Swart et al., 1987; Schultz et al., 1989; Mountjoy et al., 1992; Schaltegger et al., 1994);

  • sediment, dust, and soil provenance (e.g. Douglas et al., 1995; Revel-Rolland et al., 2006; Bathgate et al., 2011; Bataille and Bowen, 2012; De Deckker et al., 2014; Jomori et al., 2017; De Deckker, 2020; de Caritat et al., 2022b);

  • stratigraphy and reservoir/basin analysis (e.g. Rundberg and Smalley, 1989; Smalley et al., 1989, 1992; McArthur et al., 2012);

  • fluid flow (e.g. McNutt et al., 1987; Stueber et al., 1987; Sullivan et al., 1990; Bagheri et al., 2014);

  • weathering and pedogenesis (e.g. Åberg et al., 1989; Faure and Felder, 1981; Wickman and Jacks, 1992; Blum et al., 1994; Bullen et al., 1994, 1997; Quade et al., 1995; Probst et al., 2000; Harrington and Herczeg, 2003; Oliver et al., 2003; Green et al., 2004);

  • terrestrial ecosystems and catchment processes (e.g. Graustein, 1989; Jacks et al., 1989; Lyons et al., 1995; Négrel and Grosbois, 1999; Négrel and Pauwels, 2004; Gosselin et al., 2004; Chadwick et al., 2009; Hagedorn et al., 2011);

  • hydrology and hydrogeology (e.g. Collerson et al., 1988; Andersson et al., 1992, 1994; Douglas et al., 1995; Yang et al., 1996; Grobe et al., 2000; Dogramaci and Herczeg, 2002; Ojiambo et al., 2003; Palmer et al., 2004; Cartwright et al., 2007; Christensen et al., 2018; Shin et al., 2021);

  • environment and environmental change (e.g. Andersson et al., 1990; Åberg, 1995; Åberg et al., 1995; Oishi, 2021);

  • (palaeo)climatology and palaeogeographic reconstructions (e.g. Blum and Erel, 1995; Wei et al., 2018; Flecker et al., 2002);

  • determining the 87Sr/86Sr of seawater through geological time (e.g. Veizer, 1989; Gruszczynski et al., 1992; Denison et al., 1994a, b; Shields, 2007); and

  • ore genesis and mineral exploration (e.g. Le Bas et al., 1997; de Caritat et al., 2005; Daneshvar et al., 2020; Wei et al., 2020; Zhao et al., 2021).

Outside of the geosciences, food tracing and provenancing have also been underpinned by the use of Sr isotopes, although, in this case, generally relying on the bioavailable Sr rather than the total Sr (e.g. Voerkelius et al., 2010; Di Paola-Naranjo et al., 2011; Vinciguerra et al., 2015; Hoogewerff et al., 2019; Moffat et al., 2020). Anthropological studies have relied on 87Sr/86Sr isotope ratios to locate archaeological artefacts or reconstruct ancient human behaviours (e.g. Frei and Frei, 2013; Willmes et al., 2014, 2018; Adams et al., 2019; Pacheco-Forés et al., 2020; Washburn et al., 2021). Animal migration studies have also relied on Sr isotope data (e.g. Price et al., 2017). More recently, large-scale compilations and machine-learning-based predictions of the 87Sr/86Sr variations up to the continental and even global scale have been proposed (e.g. Bataille et al., 2014, 2018, 2020).

Strontium isotope landscape maps (“isoscapes”) provide the fundamental context required for the interpretation of more detailed scientific research about processes or provenance. Despite the plethora of research using Sr isotopes to address various scientific questions, very few Sr isoscapes exist in the Southern Hemisphere, particularly for soils or covering large swathes of the Earth's surface (see Bataille et al., 2020). Two exceptions to this are (1) the work by Adams et al. (2019), which reported 87Sr/86Sr in plant, soil, and biota over ∼300 000 km2 on the Cape York Peninsula in Australia, and (2) the ∼500 000 km2 Sr isoscape of inland southeastern Australia recently published by our team (de Caritat et al., 2022b). The present study affords an opportunity to further redress this deficiency and will reduce the Northern Hemisphere bias in future global 87Sr/86Sr models. It also pertains to a land surface that has not been rejuvenated by recent glaciation, consisting of over 85 % regolith or weathered material (Wilford, 2012), and, as a result, is abundant in minerals such as kaolinite, illite–smectite, goethite, and hematite. The choice of total rather than bioavailable Sr as the focus of this work was driven by an emphasis on geological sources and processes.

2 Setting

The study area in northern Australia stretches across Western Australia, the Northern Territory, and Queensland (hereafter abbreviated to WA, NT, and Qld, respectively), north of 21.5 S (Fig. 1). It is surrounded by the Indian Ocean to the west, the Timor and Arafura seas and the Gulf of Carpentaria to the north, and the Coral Sea (Pacific Ocean) to the east. The main climate zones in the area are described as “Hot humid summer” in the north and along the coast and “Warm humid summer” further south and inland (Bureau of Meteorology, 2022a). The vegetation zones are dominated by “Tropical savanna” and “Hot grassland with winter drought” (Bureau of Meteorology, 2022b). The 10-year (1996–2005) average minimum and maximum temperatures mostly range from 15 to 24 C and from 27 to 36 C, respectively (Bureau of Meteorology, 2022c). Average annual rainfall over the 4-year period to November 2009 (when the bulk of the sampling was completed) mostly ranges from 400 to 1500 mm yr−1 and is strongly seasonal (summer rain) (Bureau of Meteorology, 2022d). Physiographically, the study area includes, from west to east, the Kimberley, North Australian Plateaus, Barkly–Tanami Plains, Carpentaria Fall, Carpentaria Lowlands, Peninsular Uplands, and Burdekin Uplands provinces (Pain et al., 2011). Topographic altitude ranges from sea level to 1622 m a.s.l. (above sea level) at Mount Bartle Frere (Qld's highest point, located just south of Cairns), and the mean altitude is ∼260 m a.s.l. (Geoscience Australia, 2008).

The soil types most commonly encountered in the study area are, according to the Australian Soil Classification scheme (Isbell and National Committee on Soil and Terrain, 2021; Australian Soil Resource Information System, 2022a), Vertosol, Kandosol, or Tenosol (∼24 % of the sample sites each), followed by Rudosol, Hydrosol, or Sodosol (7 %–9 %) and the rarer Calcarosol, Chromosol, Dermosol, or Podosol (<2 %). The major river basins that divide the area belong, from west to east, to the Pilbara, Sandy Desert–Mackay, Fitzroy, North Kimberley, Ord, Victoria, Tennant, Daly, Darwin, Melville, Alligator, Roper, Arnhem, McArthur, Channel Country, Burketown, Leichhardt, Flinders, Gilbert, Weipa, Mitchell, Princess Charlotte Bay, Burdekin, Barron, Cooper Creek, and Whitsunday water regions (Geoscience Australia, 1997). Land use over the area is overwhelmingly “Grazing natural vegetation”, followed by “Other protected areas (inc. indigenous uses)” and “Nature conservation” (Australian Soil Resource Information System, 2022b).

Figure 1The northern Australia Sr isotope study area (see the inset for location: WA – Western Australia, NT – Northern Territory, Qld – Queensland, NSW – New South Wales, ACT – Australian Capital Territory, Vic – Victoria, Tas – Tasmania, and SA – South Australia) and the National Geochemical Survey of Australia (NGSA) and Northern Australia Geochemical Survey (NAGS) 87Sr/86Sr sample locations (small black squares) shown with towns (large orange squares) and places (medium orange squares), main roads (solid green lines in panel a), secondary roads (dashed grey lines in panel a), and the NGSA catchment boundaries (medium blue lines in panel b). Map projection: Albers equal area.

The main geological region groups (Blake and Kilgour, 1998) over the study area are, from west to east, the Pilbara, Canning, Kimberley, Bonaparte, Ord, Victoria, Pine Creek, Victoria, Wiso, Arunta–Tanami, Tennant, Money, McArthur, Arafura, Georgina, Isa, Carpentaria, Eromanga, Georgetown–Coen, Quinkan, and Paleozoic-Qld groups (see Fig S1 in the Supplement). The higher-level geological region domains (Blake and Kilgour, 1998) over the study area are, from west to east, the Pilbara, Paleozoic, North Australian Craton, Proterozoic, and Meso-Cenozoic domains (see Fig S2 in the Supplement). Thus, the study area displays a complex and protracted geological history spanning over 3.6 Gyr (Ollier, 1988; Braun et al., 1998; Betts et al., 2002; Blewett, 2012; Withnall et al., 2013; Ahmad and Scrimgeour, 2013), with two cratonic nuclei in the west and centre (the Archean Pilbara Craton and Palaeoproterozoic–Mesoproterozoic North Australian Craton) flanked by younger Proterozoic and Palaeozoic orogenic belts and basins, including the Phanerozoic Tasmanides to the east. Mesozoic and Cenozoic sedimentary sequences of the Eromanga and Carpentaria basins conceal much of the basement terrain over the eastern third of the study area. The whole region has experienced extensive weathering, resulting in a ubiquitous and locally thick regolith mantle.

Figure 2a shows that the decidedly dominant bedrock unit intersected by the sample sites is “Sedimentary and low-grade metamorphics” (89 % of sites), followed by the much less frequent “Felsic igneous intrusive” (4 %) and “Mafic to ultramafic igneous volcanic” (3 %) lithologies. Representing 2 % or less each are the “Medium-grade metamorphics”, “Felsic to mafic igneous volcanic”, “High-grade metamorphics”, “Mafic igneous intrusive”, and “Felsic to intermediate igneous volcanic” lithologies. Further investigation of the dominant bedrock unit “Sedimentary and low-grade metamorphics” reveals that it is comprised of “Regolith” (79 %), “Sedimentary siliciclastic” (4 %), “Feldspar- or lithic-rich arenite to rudite” (2 %), “Sedimentary carbonate” (2 %), and <1 % each of “Argillaceous detrital sediment”, “Metasedimentary siliciclastic”, and “Quartz-rich arenite to rudite”. Further, the “Regolith” class is reported to mainly consist of “Alluvium” (51 %), “Sand plain” (9 %), “Estuarine and delta deposit” (4 %), “Colluvium” (4 %), “Calcrete” (3 %), “Undifferentiated sediment” (3 %), “Lake deposit” (2 %), and <2 % each of “Black soil plain”, “Dune”, and “Ferruginous duricrust”.

The bedrock ages (periods) intersected at the sample sites are overwhelmingly Tertiary–Quaternary (58 %), followed by much less frequent Mesoproterozoic (8 %), Cenozoic (6 %), Cambrian (5 %), Palaeoproterozoic (5 %), Palaeoproterozoic–Mesoproterozoic (3 %), Carboniferous–Permian (3 %), Permian (3 %), Devonian–Carboniferous (2 %), Archaean (2 %), and Cambrian–Ordovician (2 %) ages. Representing <2 % are the Devonian, Jurassic–Cretaceous, Ordovician, Silurian–Devonian, Jurassic, Neoproterozoic–Cambrian, and Triassic periods.

Numerous mineral occurrences are found in northern Australia, and Fig. 2b shows the most important ones classified as “Mineral Deposits” (i.e. those with an inferred resource), “Operating Mines” (currently producing), or “Developing Projects” (approved but not yet producing). As mineralization is the end point of geologically “unusual” processes taking place regionally in “mineral systems” (Wyborn et al., 1994), it is useful to investigate if mineral system processes leave a recognizable 87Sr/86Sr fingerprint in the country rock and sediment derived therefrom. The study area is particularly rich in “Base Metals” deposits (both Pb–Zn- and Cu-dominated deposits, e.g. the Mount Isa, McArthur River, and George Fisher deposits) and is host to the “North Australian Zinc Belt”, the world's largest Zn–Pb province (Huston et al., 2022). In addition, notable occurrences of “Battery or Alloy Metals” (Ni, Co, Mn, V, Mo, and Mg, e.g. the Spinifex Ridge, Ripon Hills, Julia Creek, and Richmond deposits) and “Precious Metals” (Au and Ag, e.g. the Mount Isa and George Fisher deposits) are also found here. Further, “Other Metals” (Sn, Sb, W, Ta, and Nb, e.g. the Mount Carbine deposit), “Rare Earth Elements” (e.g. Thunderbird and Nolans Bore deposits), “Platinum Group Elements” (Pt, Pd, and Rh, e.g. the Munni Munni deposit), and “Heavy Mineral Sands” (e.g. Thunderbird deposit) occurrences are also catalogued. Thus, the area is significant for underpinning Australia's critical mineral resources supply now and into the future (e.g. Department of Industry, Science, Energy and Resources, 2022), without which a global transition to a lower-carbon economy will be challenging.

Figure 2The northern Australia Sr isotope study area and National Geochemical Survey of Australia (NGSA) and Northern Australia Geochemical Survey (NAGS) 87Sr/86Sr sample locations (small black squares) shown with (a) surface geology from Raymond et al. (2012) and (b) NGSA catchment boundaries (medium blue) and mineral occurrences (Geoscience Australia, 2022a). Map projection: Albers equal area.

3 Material and methods

3.1 Material

This study makes use of archive “catchment outlet sediment” samples collected during the National Geochemical Survey of Australia (NGSA), which covered ∼80 % of Australia (de Caritat and Cooper, 2011a, 2016; de Caritat, 2022). The sampling philosophy of the NGSA was to collect naturally mixed and fine-grained fluvial/alluvial sediments from large catchments, thereby obtaining a representative average of the main soil and rock types contributing sediment through weathering. This allowed an ultralow sampling density (approximately one sample per 5200 km2) that was still representative of large-scale natural variations (de Caritat and Cooper, 2011b). Catchment outlet sediments are similar to floodplain sediments in the sense that they are deposited during receding floodwaters outside the riverbanks; however there is an added complexity in that, in Australia, many areas can also experience the addition (or loss) of aeolian material. The sampled floodplain geomorphological entities are typically vegetated and biologically active (e.g. plants, worms, and ants), thereby making the collected materials true soils, albeit soils developed on transported alluvium parent material.

The sampling medium and density were both strategically chosen in the NGSA project to prioritize coverage over resolution. This was justified by the fact that the NGSA was Australia's first and, to date, only fully integrated, internally consistent geochemical survey with a truly national scope. In terms of the present study area, it is clear that these choices have implications with respect to the granularity of the patterns revealed by the Sr isoscape; as the collection of Sr isotope data in Australia using NGSA samples grows in the future (e.g. de Caritat et al., 2022b, and this contribution), it is hoped that the value of coverage will prevail over a relative low resolution of detailed features.

The NGSA collected samples at two depths: a “top outlet sediment” (TOS) from a shallow (0.1 m) soil pit approximately 0.8 m × 0.8 m in area and a “bottom outlet sediment” (BOS) from a minimum of three auger holes that were generally drilled within ∼10 m of the TOS pit. The auger holes were drilled as deep as possible (to refusal or to a maximum depth of 1 m), and the BOS sample was collected from an average depth of 0.6–0.8 m from all drilled holes. A field manual was compiled to record all sample collection method details, including site selection (Lech et al., 2007). Sampling for the NGSA took place between July 2007 and November 2009, and the field data were recorded in Cooper et al. (2010). In the laboratory, the samples were air-dried at 40 C for a minimum of 48 h (or to constant mass) before being further prepared (see de Caritat et al., 2009); for the comprehensive geochemical analysis programme of the NGSA, the reader is referred to de Caritat et al. (2010). An aliquot of a minimum of ∼1 g of sample was milled to a fine powder using a carbon steel ring mill or, for a few samples only, an agate micromill. The main sample type selected for the present Sr isotope study was NGSA BOS <2 mm in order to be as representative as possible of the geogenic background unaffected by modern land use practices and inputs (e.g. fertilizers). A few NGSA TOS <2 mm samples, prepared in an identical fashion, were also analysed.

Several additional samples from the Northern Australia Geochemical Survey (NAGS; Bastrakov and Main, 2020) were included in this project, as they provide a higher-density coverage over part of the study area. The NAGS project used the same sampling philosophy and sample collection, preparation, and analysis methods as the NGSA, with a higher sampling density of one sample per ∼500 km2. The NAGS project collected only TOS (0–0.1 m depth) samples, and the NAGS TOS <2 mm aliquots were prepared and analysed as per the NGSA samples. Sampling for the NAGS took place in May and June 2017.

Overall, 326 NGSA BOS <2 mm, 18 NGSA TOS <2 mm (including 15 with BOS also analysed), and 28 NAGS TOS samples, resulting in a total of 372 analyses from 357 samples, were analysed for Sr isotopes as detailed in Sect. 3.2 below. Given that there are ∼10 % field duplicates in the NGSA and the smaller NAGS catchments are nested within NGSA catchments, all of those samples originate from within 307 NGSA catchments, which together cover 1.536 ×106 km2 of northern Australia (see Fig. 1).

3.2 Methods

Samples were prepared and analysed for Sr isotopes (87Sr/86Sr) at the Wollongong Isotope Geochronology Laboratory (WIGL). Approximately 50 mg of sample was weighed and digested in a 2:1 mixture of hydrofluoric and nitric acids. All reagents used were SEASTAR™ BASELINE® grade, with Sr concentrations typically <10 ng kg−1. Following digestion, samples were redissolved in aqua regia (twice if needed) in order to eliminate any fluorides, followed by nitric acid twice. Finally, samples were redissolved in 2 M nitric acid prior to ion exchange chromatography. Strontium was isolated from the sample matrix using an automated, low-pressure chromatographic system (Elemental Scientific prepFAST MC™) and a 1 mL Sr–Ca column (Eichrom™) (Romaniello et al., 2015). The Sr elutions were redissolved in 0.3 M nitric acid. Strontium isotope analysis was performed on a Neptune Plus™ (ThermoFisher Scientific) multicollector inductively coupled plasma mass spectrometer (MC-ICP-MS) at WIGL. The sample introduction system consists of Apex-ST PFA MicroFlow nebulizer (Elemental Scientific, Inc.) with an uptake rate of ∼0.1 mL min−1, a quartz SSI dual cyclonic spray chamber, jet sample, and X-skimmer cones. Measurements were performed in low-resolution mode. The instrument was tuned at the start of each session with a 20 µg kg−1 Sr solution, and sensitivity for 88Sr was typically around 4 V. Masses 88, 87, 86, 85, 84, and 83 were collected on Faraday cups. Instrumental mass bias was internally corrected using measured 87Sr/86Sr. Masses 85 and 83 were used to correct for the isobaric interference of 87Rb and 86Kr, respectively. Maps were prepared using the open software QGIS® (version 3.16.14 Hannover) and applying an Albers equal area projection. Symbology for displaying 87Sr/86Sr data here either classified point data in eight equal quantile classes (12.5 % of the data each, using a gradient from green denoting low to red denoting high) at the sampling site or attributed the same value and colour to the whole catchment from which the outlet sediment comes, reflecting the sampling medium, catchment outlet sediment, being a representative sample of the average materials in the catchment (see Sect. 3.1).

3.3 Quality assessment

National Institute of Standards and Technology (NIST) strontium carbonate isotope Standard Reference Material SRM987 was used as a secondary standard and measured after every five samples to assess accuracy during analysis. The accuracy of the whole procedure was assessed by processing United States Geological Survey (USGS) reference material basalt from the Columbia River standard BCR-2 (Plumlee, 1998). The mean ±2 standard errors 87Sr/86Sr value for BCR-2 in this study is 0.704961±35 (n=13), which is within the error of the value in Jweda et al. (2016) (0.704500±11). Total procedure blanks ranged between 0.025 and 0.245 ng Sr (n=12). A total of 20 field duplicate sample pairs (collected at a median distance of ∼80 m from one another on the same landscape unit; see Lech et al., 2007) were analysed for 87Sr/86Sr in the BOS <2 mm sample, and they returned a median relative standard deviation of 0.17 %. The relative standard deviation from field duplicates includes natural variability (mineralogical/chemical heterogeneity of the alluvial deposit) as well as sample collection, preparation, and analysis uncertainties.

Overall, we feel that the quality of the 87Sr/86Sr data presented herein is adequate for the purpose of regional mapping and that reporting 87Sr/86Sr data to the third decimal place with an indicative fourth decimal place is appropriate for this work. This relatively low precision obtained for field duplicates is attributed to the heterogeneity of the alluvial deposits, as precision relating to sample preparation and analysis for Sr isotopes is at the fifth decimal place (see the results for BCR-2 above).

3.4 Data analysis

Data management was performed using Microsoft Excel®, graphing and visualization was carried out using IMDEX ioGas®, and spatial analysis and mapping was undertaken with QGIS® (an open-source geographical information system). For the purposes of generating the Sr isoscape from combined BOS and TOS samples, calculated “BOS-equivalent” values of TOS were derived from regression analysis (see Sect. 5.2). Catchment-based Sr isoscapes were constructed by assigning the 87Sr/86Sr value of a catchment's outlet sediment sample to each corresponding NGSA catchment or, if more than one sample was analysed per catchment (e.g. field duplicates or higher-resolution NAGS samples), by assigning the mean of those multiple samples. All maps are shown in Albers equal area projection.

4 Results

The soil 87Sr/86Sr values reported herein range from 0.7048 to 1.0330 (range = 0.3282). The median is 0.7405 and the mean is 0.7532 (standard deviation = 0.0480; kurtosis = 8.1602; skewness = 2.4640). Figure 3 illustrates the univariate structure of the new data. Spatially, the 87Sr/86Sr values define large-scale, coherent patterns with multi-point low and high regions (Fig. 4). The main high-value (radiogenic) regions are found in the central and northern parts of the NT, most of Cape York, and along some of the coast (including the WA part of the study). Prominent low-87Sr/86Sr-value (unradiogenic) regions include a large, central, northwest–southeast-trending elongated area in the NT, a smaller and similarly trending area in central Qld, and the northernmost and easternmost parts of the Qld study area.

Figure 3Univariate distribution of the 87Sr/86Sr data (n=357) from the northern Australia Sr isotope study area: (a) histogram (20 bins, 0.016 wide); (b) cumulative frequency plot; and (c) Tukey box plot (the solid white dot shows the mean, circles denote outliers, and triangles represent extreme outliers) (Tukey, 1977). The sample medium is the <2 mm fraction of NGSA bottom outlet sediment (BOS) or equivalent (see the text for further information).


Figure 4Strontium isoscape for the northern Australia study area with (a) data points classed by quantiles at the sampling sites and (b) as whole catchments coloured by same colour ramp. The sample medium is the <2 mm fraction of NGSA bottom outlet sediment (BOS) or equivalent (see the text for further information). Map projection: Albers equal area.

5 Discussion

5.1 Comparison with other datasets

A comparison of the present results with selected soil 87Sr/86Sr datasets from around the world is offered in Table 1. All discoverable data from Australia and the Southern Hemisphere comprising more than a handful of sites were included, whereas only selected Northern Hemisphere datasets focusing on recent and large-scale datasets were included. Table 1 combines both bioavailable/exchangeable and bulk/total 87Sr/86Sr data, representing variable soil grain-size fractions (where reported), parent materials, and land uses. Despite this variability, one can make several observations. Firstly, the northern Australia Sr isotope dataset has the highest maximum (1.0330) and largest range (0.3282) of 87Sr/86Sr values amongst those compiled. Secondly, there is no observed standard protocol for collecting or preparing soils for either bioavailable or total 87Sr/86Sr determination. Thirdly, total 87Sr/86Sr datasets tend to have a wider range and greater variance than bioavailable ones, at least partly accounting for the higher values encountered here. These issues complicate data compilation and integration across projects/countries but do not preclude them. Indeed, through contemporaneous data analytics, including artificial intelligence and machine learning, it is likely that the relationships between bioavailable and total Sr isotope values and a host of other environmental variables (including from climatic, topographic, biotic, and geoscientific categories) can be teased out and that high-spatial-resolution models/predictions of bioavailable or total Sr isotope distribution can be derived (e.g. Bataille et al., 2020). This is indeed a future research direction that we propose/support for isoscape science in general as well as in an Australian Sr isoscape in particular.

Table 1Comparison of 87Sr/86Sr data from this study with selected datasets from around the world, with an emphasis on Australia and the Southern Hemisphere (count, minimum, median, mean, standard deviation, maximum, and range). The regions are abbreviated as follows: Aus – Australia, Qld – Queensland, SA – South Australia, SE – southeastern, Vic – Victoria, and WA – Western Australia. The digestions are abbreviated as follows: AcAc – acetic acid (CH3COOH), AmAc – ammonium acetate (CH3COONH4), AmNt – ammonium nitrate (NH4NO3), and AR – aqua regia.

Download Print Version | Download XLSX

5.2 Top–bottom relationship

Based on the 15 sites where both TOS and BOS samples were analysed, a strong correlation (Fig. 5) between the aforementioned two depths is found:

(1) ( 87 Sr / 86 Sr ) BOS = 0.9913 × ( 87 Sr / 86 Sr ) TOS + 0.0069 ( r 2 = 0.97 ; p < 0.001 ; n = 15 ) .

Equation (1) was used to infer BOS 87Sr/86Sr from those sites where only TOS samples could be obtained (i.e. 3 NGSA and 28 NAGS sites), effectively deriving a BOS-equivalent 87Sr/86Sr value for the purposes of performing internally consistent statistical and spatial analysis.

As a result of the robust correlation between TOS and BOS, the BOS (or subsurface) Sr isoscapes of Fig. 4 can be recalculated to TOS (or surface) Sr isoscapes by replacing the legend values for the eight quantile class boundaries with the values 0.7050–0.7156–0.7215–0.7293–0.7401–0.7522–0.7674–0.7910–1.0276. A TOS map and values would be useful when studying surface processes, such as plant and animal uptake or the effect of land use on sediment composition and dynamics.

Figure 5Scatterplot of the BOS (deep) versus TOS (surface) 87Sr/86Sr data (n=15) from the northern Australia Sr isotope study area with total uncertainty derived from field duplicates (0.17 %; grey plus symbols) and linear square regression line (dashed line).


5.3 Relationship to bedrock

When grouped by age of their respective geological region, the new 87Sr/86Sr data show a general trend of increasing 87Sr/86Sr with increasing bedrock age, as illustrated by the box and scatter plots in Fig. 6. In Fig. 6, average numerical ages (mid-points between beginning and end ages) are attributed to each geological period group. The observed trend is consistent with that of the known controls on mineral (and rock) 87Sr/86Sr values, namely their age.

Figure 6Distribution of the 87Sr/86Sr data from the northern Australia Sr isotope study area by geological region age: (a) notched Tukey box plots (Tukey, 1977) of, left to right, the whole dataset first (grey; n=357) and then sorted by average age (early Tertiary, n=6; Middle Jurassic to Holocene, n=80; Jurassic to Tertiary, n=5; Cretaceous, n=4; Palaeozoic to Mesozoic, n=28; Silurian to Triassic, n=6; Palaeozoic, n=61; Neoproterozoic to early Permian, n=7; Neoproterozoic to Devonian, n=49; Mesoproterozoic to Permian, n=2; Mesoproterozoic to Neoproterozoic, n=2; Mesoproterozoic, n=4; Palaeoproterozoic to Permian, n=7; Palaeoproterozoic to Neoproterozoic, n=13; Palaeoproterozoic to Mesoproterozoic, n=59; Palaeoproterozoic, n=7; Archaean to Palaeoproterozoic, n=12; and Archaean, n=5); (b) a scatterplot against the average age of the geological region with the group averages (crosses) and linear square regression line (dashed line).


Another known control on Sr isotopic values is the Rb content, which, over time, contributes 87Sr by radiogenic decay. Closer examination of the geology in the catchments with the 10 highest 87Sr/86Sr values indicates that they contain lithologies likely to be enriched in Rb relative to Sr (e.g. felsic, felsic-to-intermediate igneous rocks, and medium- to high-grade metamorphic rocks), as shown using detailed maps in the Supplement. Of the 5 top 10 87Sr/86Sr values in the northern NT, 3 are from the Pine Creek geological region. The two northernmost of these (sample nos. 2007190099 and 2007191390) with 87Sr/86Sr values of 0.9413 and 0.9900, respectively, are in catchments downstream of 1780–1790 Ma I-type granite intrusions that have been reported to have 87Sr/86Sr values from 0.71 to 8.13, 87Rb/86Sr values from ∼1 to 290, Rb concentrations of up to 956 mg kg−1, and Sr concentrations of up to 459 mg kg−1 (Riley, 1980, their Table 1).

The elevated samples from the WA coast are also generally proximal to, or downstream of, felsic igneous lithologies (including the granitoid and gneiss underlying the Mallina Basin; Van Kranendonk et al., 2002) and are consistent with 87Sr/86Sr values reported for clay fraction sediments from the region by De Deckker (2019), as shown in the Supplement. The northernmost samples along the coast of WA overlap with the King Leopold Orogen, abutting the western flank of the Kimberley geological region. Few 87Sr/86Sr data are available for this region, but those that exist within the Hooper Complex (Griffin et al., 2000) indicate radiogenic 87Sr/86Sr values: 0.7708–0.9620 for the Whitewater Volcanics, 0.7336–0.8524 for the Kongorow Granite, and 0.7353–0.7602 for the Lennard Granite (Bennett and Gellatly, 1970, their Table 2). Thus, observed elevated 87Sr/86Sr values of 0.8080, 0.7709, and 0.7923 in three sediment samples (sample nos. 200719103, 2007190317, and 2007190966, respectively) near Derby, WA, are entirely compatible with these bedrock source rocks.

Page et al. (1980, their Table 5) published Rb and Sr data from Alligator Rivers uranium field rocks, including from the Nimbuwah Complex granitoid close to Cooper Creek, <20 km north of the Nabarlek uranium deposit (their sample no. 7212.4063). For this sample, they reported 87Rb/86Sr values of 1.25–1.47, 87Sr/86Sr values of 0.7387–0.7442, Rb concentrations of 146–164 mg kg−1, and Sr concentrations of 320–336 mg kg−1. Our sample (sample no. 2007190710) from this same catchment has a 87Sr/86Sr value of 0.8019 and is located on a felsic granite polygon.

Black et al. (1983, their Tables 3 and 5) reported calculated initial rock 87Sr/86Sr ratios in the range from 0.76 to 0.92 from granites, including the Wangala and the Haverson granites, in the northern Arunta geological region (southern NT). These intrusions intersect three sampled NGSA catchments for which we obtained elevated 87Sr/86Sr values, namely 0.8140, 0.7763, and 0.8651 (sample nos. 2007190727, 2007190562, and 2007190123, respectively). These catchments are also just upstream of some of our top 10 87Sr/86Sr values (between 0.8935 and 1.0330).

In the Georgetown geological region of northwestern Qld, broad consistency between whole-rock 87Sr/86Sr and catchment sediment 87Sr/86Sr is also observed. For instance, our sample no. 2007190858, from a catchment that drains mostly the Esmeralda Granite (average whole-rock initial 87Sr/86Sr of 0.7314±0.0116; Black, 1973, their Table 2) on the western edge of the Georgetown Inlier, has a relatively elevated 87Sr/86Sr value of 0.7467. Conversely, the three catchments that mostly directly drain the less radiogenic Newcastle Range Volcanics in the centre of the Georgetown Inlier, which have a lower average whole-rock initial 87Sr/86Sr of 0.7152±0.0011 (Black, 1973, Table 2), have correspondingly lower sediment 87Sr/86Sr values of 0.7318 (sample no. 2007190572), 0.7268 (sample no. 2007190348), and a field duplicate pair of 0.7258 (sample no. 2007190548) and 0.7282 (sample no. 2007190548).

The floodplain sediment Sr isotopic values recorded in areas dominated by sedimentary carbonate and mafic/ultramafic igneous rocks are usually within an intermediate range between the more radiogenic and unradiogenic end-members discussed above. Indeed, samples sited within 0.1 (∼10 km) of lithologies recorded as sedimentary carbonate and mafic/ultramafic igneous rocks have median values of 0.7387 (n=96) and 0.7422 (n=42), respectively.

Figure 7Distribution of the 87Sr/86Sr data from the northern Australia Sr isotope study area by mineral occurrence hosted in same catchment. The categories shown, from left to right, are as follows: all (n=97); Base Metals – Cu (Zn, Pb, Ag, Au) (n=5); Base Metals – Zn, Pb (Cu, Ag) (n=19); Battery/Alloy Metals – Ni, Co, Mn, V, Mo, Mg (n=17); Bauxite (n=3); Coal (n=2); Diamond (n=3); Fertilizer Elements – P, K (n=12); Graphite (n=2); Heavy Mineral Sands (n=1); Iron Ore (n=5); Other Metals – Sn, Sb, W, Ta, Nb (n=6); Platinum Group Elements (n=2); Precious Metals – Au, Ag (n=13); Rare Earth Elements (n=6); and Uranium (n=1).


5.4 Relationship to mineralization

Figure 7 illustrates the range of 87Sr/86Sr values of catchments that contain known mineralization of various types. A total of 44 NGSA catchments host 97 mineral occurrences as catalogued by Geoscience Australia (2022a). These resources include 19 “Base Metals – Zn, Pb (Cu, Ag)”; 17 “Battery/Alloy Metals – Ni, Co, Mn, V, Mo, Mg”; 13 “Precious Metals – Au, Ag”; 12 “Fertilizer Elements – P, K”; 6 “Rare Earth Elements”; 6 “Other Metals – Sn, Sb, W, Ta, Nb”; 5 “Iron Ore”; and 5 “Base Metals – Cu (Zn, Pb, Ag, Au)” occurrences. Whilst the vast majority of mineral occurrences are found in catchments with a sediment outlet 87Sr/86Sr signature ≪0.84, six outlier occurrences are associated with higher catchment outlet sediment 87Sr/86Sr values: 0.8434 (the Coronation Hill Pt–Pd deposit in the catchment containing sample no. 2007191552); 0.9373 (Nolans Bore hard-rock rare earth element (REE) and Mount Peake vanadium deposits in the catchment containing sample no. 2007191112); 0.9748 (Batman Au deposit in the catchment containing sample no. 2007191329); and, highest of all, 0.9900 (the Browns Pb–Zn and the Browns Co deposits in the catchment containing sample no. 2007191387). Some mineral occurrence types are not found in catchments with 87Sr/86Sr signature less radiogenic than 0.75 in this study: the minima reported for the Platinum Group Elements, Rare Earth Elements/Heavy Mineral Sands, and Diamond commodity groups are 0.7529, 0.7606, 0.7606, and 0.7567, respectively. Whilst it is unlikely that mineral deposits themselves impart a significant control on the 87Sr/86Sr of sediment collected down-catchment given their usually limited size, the hosting country rock may, however, potentially be more widely affected by mineral system processes, and, upon further investigation, appears to record some Sr isotopic effect that may be useful for mineral vectoring. One of the most radiogenic 87Sr/86Sr signatures (0.9373) is found for the catchment hosting both the Mount Peake vanadium deposit (Simandl and Paradis, 2022) and the Nolans Bore hydrothermal REE deposit (Schoneveld et al., 2015). Hydrothermal or magmatic REE deposits are not usually characterized by highly radiogenic Sr signatures (e.g. 0.7029–0.7262 at Bayan Obo – Le Bas et al., 1997; 0.701–0.708 in global review of carbonatite deposits – Bolonin, 2019; and 0.7037–0.7062 at Caotan – Wei et al., 2020) and nor is Nolans Bore specifically (0.7054–0.7079; Huston et al., 2016). Thus, the elevated 87Sr/86Sr reported here from this catchment is attributed to the radiogenic bedrock occurring in the study area and perhaps in upstream catchments immediately to the south of it.

6 Data availability

The new spatial Sr isotope dataset for the northern Australia region is publicly available from (de Caritat et al., 2022a) and via the Geoscience Australia portal (last access: 27 February 2023; Geoscience Australia, 2023). Metadata are also available via the Geoscience Australia portal (, last access: 27 February 2023) (Geoscience Australia, 2022b).

7 Conclusions

A total of 372 new strontium (Sr) isotopic compositions (87Sr/86Sr) are reported from 357 catchment outlet sediment samples from northern Australia (north of 21.5 S). The analysed material originates from the sample archives of the National Geochemical Survey of Australia (NGSA; n=344) and Northern Australia Geochemical Survey (NAGS; n=28) projects, both of which targeted overbank or floodplain landforms near the outlet of large catchments. The sampled catchments together cover 1.536 ×106 km2 of northern Australia. For the most part, bottom outlet sediment (BOS) samples, retrieved mostly by drilling (with an auger) to, on average, 0.6 to 0.8 m depth, were analysed. However, a few top outlet sediment (TOS) samples, collected from the top 0.1 m of soil, were also included, notably all of the NAGS samples and 18 NGSA samples. Total digestion of milled <2 mm grain-size fractions from these sediments yielded a wide range of 87Sr/86Sr values from a minimum of 0.7048 to a maximum of 1.0330.

To date, the present study represents the largest Sr isoscape in the Southern Hemisphere, a region that is critically under-represented in terms of isotopic data worldwide. We found a very strong correlation between the 87Sr/86Sr values in the BOS and TOS samples, allowing us to confidently infer BOS-equivalent values where only TOS samples were analysed. A map of the 87Sr/86Sr distribution (isoscape) across northern Australia reveals spatial patterns reflecting the ages and lithologies of the source material for the sediment, which is principally carried down catchment by fluvial processes. There is an overall increase in 87Sr/86Sr observed with increasing age of the geological region intersected at the sampling points. In areas of outcropping or sub-cropping felsic igneous rocks, relatively radiogenic Sr signatures are observed, whereas the inverse is noted for areas of mafic igneous rocks or marine sediments. The new Sr isotopic data are also examined in terms of the mineral occurrences found in their catchment. Whilst most mineral occurrences in the region are found in a catchment with an 87Sr/86Sr signature ≪0.84, six outlier occurrences are associated with higher values. Some mineral occurrence types (Platinum Group Elements, Rare Earth Elements, Heavy Mineral Sands, and Diamond) are not found in catchments with 87Sr/86Sr signature less radiogenic than 0.75 in this study. Whether mineral system or regional petrogenetic processes are responsible for the association of relatively elevated 87Sr/86Sr with mineral deposits is a subject for further investigation.

Although we have focused the discussion of the new 87Sr/86Sr data on sediment sources in terms of rock ages and types, potential applications of the present isoscape and modelling derived therefrom could be extended to studies of mineralization, hydrology, food tracing, dust provenancing or sourcing, and historic migrations of people and animals.


The supplement related to this article is available online at:

Author contributions

PdC provided the concept, samples, and funding, curated and analyzed the data, created the figures, and wrote and edited the manuscript. AD provided technical guidance, resources, and supervision, curated the data, and edited the manuscript. FD carried out the analyses, provided technical support, and curated the data.

Competing interests

The contact author has declared that none of the authors has any competing interests.


This paper is published with the permission of the Chief Executive Officer, Geoscience Australia.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The National Geochemical Survey of Australia (NGSA) project would not have been possible without Commonwealth funding through the “Onshore Energy Security Program” (, last access: 27 February 2023), and Geoscience Australia appropriation. The Northern Australia Geochemical Survey (NAGS) project was funded under the “Exploring for the Future” programme (, last access: 27 February 2023) and Geoscience Australia appropriation. Collaboration with the geoscience agencies of all states and the Northern Territory is gratefully recognized. We acknowledge all land owners and custodians, whether private, corporate, and/or traditional, for granting access to the field sites for the purposes of sampling. We are also grateful to Geoscience Australia laboratory staff for assistance with preparing and analysing the samples. We thank Geoscience Australia reviewers, Kathryn Waltenberg and David Huston; journal reviewers, Ian Moffat, Malte Willmes, and Jodie Pritchard; and the journal editor, Attila Demény, for their detailed and constructive critique of our work.

Financial support

This research has been supported by the Australian Government (Exploring for the Future,, last access: 27 February 2023). Geoscience Australia's Exploring for the Future programme provides pre-competitive 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 pre-competitive 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 programme, which commenced in 2016, is an 8-year, AUD 225 million investment by the Australian Government.

Review statement

This paper was edited by Attila Demény and reviewed by Jodie Pritchard, Ian Moffat, and Malte Willmes.


Åberg, G.: The use of natural strontium isotopes as tracers in environmental studies, Wat. Air Soil Poll., 79, 309–322,, 1995. 

Åberg, G., Jacks, G., and Hamilton, P. J.: Weathering rates and 87Sr/86Sr ratios: an isotopic approach, J. Hydrol., 109, 65–78,, 1989. 

Åberg, G., Wickman, T., and Mutvei, H.: Strontium isotope ratios in mussel shells as indicators of acidification, Ambio, 24, 265–268, 1995. 

Adams, S., Grün, R., McGahan, D., Zhao, J.-X., Feng, Y., Nguyen, A., Willmes, M., Quaresimin, M., Lobsey, B., Collard, M., and Westaway, M. C.: A strontium isoscape of north-east Australia for human provenance and repatriation, Geoarchaeol., 34, 231–251,, 2019. 

Ahmad, M. and Scrimgeour, I. R.: Geological framework, chap. 2, in: Geology and Mineral Resources of the Northern Territory, edited by: Ahmad, M. and Munson, T. J., Northern Territory Geol. Surv., Special Publication, 5, 2:1–2:16, (last access: 27 February 2023), 2013. 

Andersson, P., Lofvendahl, R., and Åberg, G.: Major element chemistry, δ2H, δ18O and 87Sr/86Sr in a snow profile across central Sweden, Atmos. Environ., 24A, 2601–2608,, 1990. 

Andersson, P. S., Wasserburg, G. J., and Ingri, J.: The sources and transport of Sr and Nd isotopes in the Baltic Sea, Earth Planet. Sc. Lett., 113, 459–472,, 1992. 

Andersson, P. S., Wasserburg, G. J., Ingri, J., and Stordal, M. C. Strontium, dissolved and particulate loads in fresh and brackish waters: the Baltic Sea and Mississippi Delta, Earth Planet. Sc. Lett., 124, 195–210,, 1994. 

Australian Soil Resource Information System (ASRIS): Digital Atlas of Australian Soils, ASRIS [data set], (last access: 27 February 2023), 2022a. 

Australian Soil Resource Information System (ASRIS): ASRIS Map Viewer Tool, Landcover Layer, Land Use 2001–2002, (last access: 27 February 2023), 2022b. 

Bagheri, R., Nadri, A., Raeisi, E., Eggenkamp, H. G. M., Kazemi, G. A., and Montaseri, A.: Hydrochemical and isotopic (δ18O, δ2H, 87Sr/86Sr, δ37Cl and δ81Br) evidence for the origin of saline formation water in a gas reservoir, Chem. Geol., 384, 62–75,, 2014. 

Bastrakov, E. N. and Main, P. T.: Northern Australia Geochemical Survey: a review of regional soil geochemical patterns, in: Exploring for the Future: Extended Abstracts, edited by: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A., and Slatter, E., Geoscience Australia, Canberra, 1–4,, 2020. 

Bataille, C. P. and Bowen, G. J.: Mapping 87Sr/86Sr variations in bedrock and water for large scale provenance studies, Chem. Geol., 304, 39–52,, 2012. 

Bataille, C. P., Brennan, S. R., Hartmann, J., Moosdorf, N., Wooller, M. J., and Bowen, G. J.: A geostatistical framework for predicting variability in strontium concentrations and isotope ratios in Alaskan rivers, Chem. Geol., 389, 1–15,, 2014. 

Bataille, C. P., von Holstein, I. C. C., Laffoon, J. E., Willmes, M., Liu, X.-M., and Davies, G. R.: A bioavailable strontium isoscape for Western Europe: a machine learning approach, PLoS ONE, 13, e0197386,, 2018. 

Bataille, C. P., Crowley, B. E., Wooller, M. J., and Bowen, G. J.: Advances in global bioavailable strontium isoscapes, Palaeogeogr. Palaeocl. Palaeoecol., 555, 109849,, 2020. 

Bathgate, H., Lovett, A., Clarke, J., Ueckermann, H., and Hoogewerff, J.: The use of 87Sr/86Sr ratios, trace elements and biological profiling for forensic provenancing of soils, Eur. Geoph. Union Gen. Assembly 2011, 3–8 April 2011, Vienna, Austria, Geophys. Res. Abstr., 13, EGU2011-6592, 2011. 

Bennett, R. and Gellatly, D. C.: Rb-Sr age determinations of some rocks from the West Kimberley Region, Western Australia, Bur. Min. Res. Geol. Geoph. Record, 1970/20, 14 pp., (last access: 27 February 2023), 1970. 

Betts, P. G., Giles, D., Lister, G. S., and Frick, L. R.: Evolution of the Australian lithosphere, Austral. J. Earth Sci., 49, 661–695,, 2002. 

Black, L. P.: Tables of Isotopic Ages from the Georgetown Inlier, North Queensland, Bur. Min. Res. Geol. Geoph. Record, 1973/50, 2 pp., (last access: 27 February 2023), 1973. 

Black, L. P., Shaw, R. D., and Stewart, A. J.: Rb-Sr geochronology of Proterozoic events in the Arunta inlier, central Australia, Bur. Min. Res. J. Austral. Geol. Geoph., 8, 129–137, 1983. 

Blake, D. H. and Kilgour, B.: Geological Regions of Australia 1:5 000 000 scale, Geosci. Austral., Canberra [data set], (last access: 27 February 2023), 1998. 

Blewett, R. (Ed.): Shaping a Nation – A Geology of Australia, Geosci. Austral. and ANU E Press, Canberra,, 2012. 

Blum, J. D. and Erel, Y.: A silicate weathering mechanism linking increases in marine 87Sr/86Sr with global glaciation, Nature, 373, 415–418,, 1995. 

Blum, J. D., Erel, Y., and Brown, K.: 87Sr/86Sr ratios of Sierra Nevada stream waters: implications for relative mineral weathering rates, Geochim. Cosmochim. Ac., 57, 5019–5025,, 1994. 

Bolonin, A.: Oxygen and carbon isotope composition in primary carbonatites of the world: data summary and linear trends. Open J. Geol., 9, 424–439,, 2019. 

Braun, J., Dooley, J., Goleby, B., van der Hilst, R., and Klootwijk, C. (Eds.): Structure and Evolution of the Australian Plate, Am. Geoph. Union, Geodyn. Series, 26, ISBN 9780875905280,, 1998. 

Bullen, T., White, A., Blum, A., Harden, J., and Schulz, M.: Chemical weathering of a soil chronosequence on granitoid alluvium: II. Mineralogic and isotopic constraints on the behavior of strontium, Geochim. Cosmochim. Ac., 61, 291–306,, 1997. 

Bullen, T. D., Blum, A. E., White, A. F., and Schulz, M. S.: A model for the temporal evolution of 87Sr/86Sr of the cation exchange pool in a simple granitoid weathering system: approach and implications, Eos, 75, 280–281, 1994. 

Bureau of Meteorology (BOM): Climate classification maps, (last access: 27 February 2023), 2022a. 

Bureau of Meteorology (BOM): Map of Climate Zones of Australia, (last access: 27 February 2023), 2022b. 

Bureau of Meteorology (BOM): Decadal and multi-decadal temperature, (last access: 27 February 2023), 2022c. 

Bureau of Meteorology (BOM): Recent and historical rainfall maps, (last access: 27 February 2023), 2022d. 

Cartwright, I., Weaver, T., and Petrides, B.: Controls on 87Sr/86Sr ratios of groundwater in silicate-dominated aquifers: SE Murray Basin, Australia, Chem. Geol., 246, 107–123,, 2007. 

Chadwick, O. A., Derry, L. A., Bern, C. R., and Vitousek, P. M.: Changing sources of strontium to soils and ecosystems across the Hawaiian Islands, Chem. Geol., 267, 64–76,, 2009. 

Christensen, J. N., Dafflon, B., Shiel, A. E., Tokunaga, T. K., Wan, J., Faybishenko, B., Dong, W., Williams, K. H., Hobson, C., Brown, S. T., and Hubbard, S. S.: Using strontium isotopes to evaluate the spatial variation of groundwater recharge, Sci. Total Environ., 637–638, 672–685,, 2018. 

Collerson, K. D., Ullman, W. J., and Torgersen, T. Ground waters with unradiogenic 87Sr/86Sr ratios in the Great Artesian Basin, Australia, Geology, 16, 59–63,<0059:GWWUSS>2.3.CO;2, 1988. 

Cooper, M., de Caritat, P., Burton, G., Fidler, R., Green, G., House, E., Strickland, C., Tang, J., and Wygralak, A.: National Geochemical Survey of Australia: Field Data, Record, 2010/18, Geosci. Austral., Canberra,, 2010. 

Daneshvar, N., Azizi, H., Asahara, Y., Tsuboi, M., and Hosseini, M.: Rare earth elements and Sr isotope ratios of large apatite crystals in Ghareh Bagh mica mine, NW Iran: tracing for petrogenesis and mineralization, Minerals, 10, 833,, 2020. 

Dart, R. C., Barovich, K. M., Chittleborough, D. J., and Hill, S. M.: Calcium in regolith carbonates of central and southern Australia: its source and implications for the global carbon cycle, Palaeogeogr. Palaeocl. Palaeoecol., 249, 322–334,, 2007. 

de Almeida, B. S.: 87Sr/86Sr Isotopic Characterization as a Tool for the Designation of Origin and Geographical Indication: Application to Volcanic Rocks, Soils, Grapes and Wines from Brazil and Italy, Doctoral Thesis, University Frederic II of Naples, 106 pp.,, 2021. 

de Caritat, P.: The National Geochemical Survey of Australia: review and impact. Geochemistry: Exploration, Environment, Analysis, geochem2022-032,, 2022. 

de Caritat, P. and Cooper, M.: National Geochemical Survey of Australia: The Geochemical Atlas of Australia, Record, 2011/20, Geosci. Austral., Canberra,, 2011a. 

de Caritat, P. and Cooper, M.: National Geochemical Survey of Australia: Data Quality Assessment, Record, 2011/21, Geosci. Austral., Canberra, (last access: 27 February 2023), 2011b. 

de Caritat, P. and Cooper, M.: A continental-scale geochemical atlas for resource exploration and environmental management: the National Geochemical Survey of Australia, Geochem. Explo. Env. Anal., 16, 3–13,, 2016. 

de Caritat, P., Kirste, D., Carr, G., and McCulloch, M.: Groundwater in the Broken Hill region, Australia: recognising interaction with bedrock and mineralisation using S, Sr and Pb isotopes, Appl. Geochem., 20, 767–787,, 2005. 

de Caritat, P., Cooper, M., Lech, M., McPherson, A., and Thun, C.: National Geochemical Survey of Australia: Sample Preparation Manual, Record, 2009/08, Geosci. Austral., Canberra, (last access: 27 February 2023), 2009. 

de Caritat, P., Cooper, M., Pappas, W., Thun, C., and Webber, E.: National Geochemical Survey of Australia: Analytical Methods Manual, Record, 2010/15, Geosci. Austral., Canberra, (last access: 27 February 2023), 2010. 

de Caritat, P., Dosseto, A., and Dux, F.: A strontium isoscape of northern Australia, Geosci. Austral., Canberra [data set],, 2022a. 

de Caritat, P., Dosseto, A., and Dux, F.: A strontium isoscape of inland southeastern Australia, Earth Syst. Sci. Data, 14, 4271–4286,, 2022b. 

De Deckker, P.: OZ-3 northern Western Australia sampling trip and analysis: strontium, rubidium and neodymium isotopes, Pangaea [data set],, 2019. 

De Deckker, P.: Airborne dust traffic from Australia in modern and Late Quaternary times, Global Planet. Change, 184, 103056,, 2020. 

De Deckker, P., Munday, C. I., Brocks, J., O'loingsigh, T., Allison, G. E., Hope, J., Norman, M., Stuut, J.-B. W., Tapper, N. J., and van der Kaars, S.: Characterisation of the major dust storm that traversed over eastern Australia in September 2009; a multidisciplinary approach, Aeol. Res., 15, 133–149,, 2014. 

Denison, R. E., Koepnick, R. B., Burke, W. H., Hetherington, E. A., and Fletcher, A.: Construction of the Mississippian, Pennsylvanian and Permian seawater 87Sr/86Sr curve, Chem. Geol., 112, 145–167,, 1994a. 

Denison, R. E., Koepnick, R. B., Fletcher, A., Howell, M. W., and Callaway, W. S.: Criteria for the retention of original seawater 87Sr/86Sr in ancient shelf limestones, Chem. Geol., 112, 131–143,, 1994b. 

Department of Industry, Science, Energy and Resources (DISER): 2022 Critical Minerals Strategy, Aust. Gov., Canberra,, last access: 27 February 2023, 2022. 

Di Paola-Naranjo, R. D., Baroni, M. V., Podio, N. S., Rubinstein, H. R., Fabani, M. P., Badini, R. G., Inga, M., Ostera, H. A., Cagnoni, M., Gallegos, E., Gautier, E., Peral-Garcia, P., Hoogewerff, J., and Wunderlin, D. A.: Fingerprints for main varieties of Argentinean wines: terroir differentiation by inorganic, organic, and stable isotopic analyses coupled to chemometrics, J. Ag. Food Chem., 59, 7854–7865,, 2011. 

Dogramaci, S. S. and Herczeg, A. L.: Strontium and carbon isotope constraints on carbonate-solution interactions and inter-aquifer mixing in groundwaters of the semi-arid Murray Basin, Australia, J. Hydrol., 262, 50–67,, 2002. 

Douglas, G. B., Gray, C. M., Hart, B. T., and Beckett, R.: A strontium isotopic investigation of the origin of suspended particulate matter (SPM) in the Murray-Darling River system, Australia, Geochim. Cosmochim. Ac., 59, 3799–3815,, 1995. 

Evans, J. A., Montgomery, J., Wildman, G., and Boulton, N.: Spatial variations in biosphere 87Sr/86Sr in Britain, J. Geol. Soc., 167, 1–4,, 2010. 

Faure, G. and Felder, R. P.: Isotopic composition of strontium and sulfur in secondary gypsum crystals, Brown Hills, Transantarctic Mountains, J. Geoch. Explo., 14, 265–270,, 1981. 

Flecker, R., de Villiers, S., and Ellam, R. M.: Modelling the effect of evaporation on the salinity-87Sr/86Sr relationship in modern and ancient marginal-marine systems: the Mediterranean Messinian Salinity Crisis, Earth Planet. Sc. Lett., 203, 221–233,, 2002. 

Frei, R. and Frei, K. M.: The geographic distribution of Sr isotopes from surface waters and soil extracts over the island of Bornholm (Denmark) – A base for provenance studies in archeology and agriculture, Appl. Geochem., 38, 147–160,, 2013. 

Geoscience Australia (GA): Australia's River Basins 1997 – Product User Guide, Geosci. Austral., Canberra, (last access: 27 February 2023), 1997. 

Geoscience Australia (GA): 9 Second Digital Elevation Model of Australia Version 3, Geosci. Austral., Canberra [data set], (last access: 27 February 2023), 2008. 

Geoscience Australia (GA): Australian Operating Mines Map 2021 Data, Geosci. Austral., Canberra [data set], (last access: 27 February 2023), 2022a. 

Geoscience Australia: Geoscience Australia Portal: Geochronology and Isotopes – Isotopes – Rb-Sr Isotope – Points, Australian Government [data set], (last access: 27 February 2023), 2022b. 

Geoscience Australia: Geochronology and Isotopes Data Portal,, Geoscience Australia [data set], last access: 27 February 2023. 

Gingele, F. X. and De Deckker, P.: Clay mineral, geochemical and Sr-Nd-isotopic fingerprinting of sediments in the Murray-Darling fluvial system, southeast Australia, Aust. J. Earth Sci., 52, 965–974,, 2005. 

Gosselin, D. C., Harvey, F. E., Frost, C., Stotler, R., and Macfarlane, P. A.: Strontium isotope geochemistry of groundwater in the central part of the Dakota (Great Plains) aquifer, USA, Appl. Geochem., 19, 359–377,, 2004. 

Gosz, J. R., Brookins, D. G., and Moore, D. I.: Using strontium isotope ratios to estimate inputs to ecosystems, Bioscience, 33, 23–30,, 1983. 

Graustein, W. C.: 87Sr/86Sr ratios measure the sources and flow of strontium in terrestrial ecosystems, chap. 28, in: Stable Isotopes in Ecological Research, edited by: Rundel, P. W., Ehleringer, J. R., and Nagy, K. A., Springer-Verlag, New York, 491–512,, 1989. 

Green, G. P., Bestland, E. A., and Walker, G. S.: Distinguishing sources of base cations in irrigated and natural soils: evidence from strontium isotopes, Biogeochem., 68, 199–225,, 2004. 

Griffin, T. J., Page, R., Sheppard, S., and Tyler, I.: Tectonic implications of Palaeoproterozoic post-collisional, high-K felsic igneous rocks from the Kimberley region of northwestern Australia, Precambrian Res., 101, 1–23,, 2000. 

Grobe, M., Machel, H. G., and Heuser, H.: Origin and evolution of saline groundwater in the Münsterland Cretaceous Basin, Germany: oxygen, hydrogen, and strontium isotope evidence, J. Geoch. Explo., 69–70, 5–9,, 2000. 

Gruszczynski, M., Hoffman, A., Krzysztof, M., and Veizer, J.: Seawater strontium isotopic perturbation at the Permian-Triassic boundary, west Spitsbergen, and its implications for the interpretation of strontium isotopic data, Geology, 20, 779–782,<0779:SSIPAT>2.3.CO;2, 1992. 

Hagedorn, B., Cartwright, I., Raveggi, M., and Maas, R.: Rare earth element and strontium geochemistry of the Australian Victorian Alps drainage system: Evaluating the dominance of carbonate vs. aluminosilicate weathering under varying runoff, Chem. Geol., 284, 105–126,, 2011. 

Harrington, G. A. and Herczeg, A. L.: The importance of silicate weathering of a sedimentary aquifer in arid Central Australia indicated by very high 87Sr/86Sr ratios, Chem. Geol., 199, 281–292,, 2003. 

Hoogewerff, J. A., Reimann, C., Ueckermann, H., Frei, R., Frei, K. M., van Aswegen, T., Stirling, C., Reid, M., Clayton, A., Ladenberger, A., and The GEMAS Project Team: Bioavailable 87Sr/86Sr in European soils: a baseline for provenancing studies, Sci. Total Environ., 672, 1033–1044,, 2019. 

Huston, D. L., Maas, R., Cross, A., Hussey, K. J., Mernagh, T. P., Fraser, G., and Champion, D. C.: The Nolans Bore rare-earth element-phosphorus-uranium mineral system: geology, origin and post-depositional modifications, Miner. Deposita, 51, 797–822,, 2016. 

Huston, D. L., Champion, D. C., Czarnota, K., Duan, J., Hutchens, M., Paradis, S., Hoggard, M., Ware, B., Gibson, G. M., Doublier, M. P., Kelley, K., McCafferty, A., Hayward, N., Richards, F., Tessalina, S., and Carr, G.: Zinc on the edge–isotopic and geophysical evidence that cratonic edges control world-class shale-hosted zinc-lead deposits, Miner. Deposita, 58, 707–729,, 2022. 

Isbell, R. F. and National Committee on Soil and Terrain: The Australian Soil Classification, third edn., CSIRO Publishing, Melbourne, Victoria, 181 pp., (last access: 27 February 2023), 2021. 

Jacks, G., Åberg, G., and Hamilton, P. J.: Calcium budgets for catchments as interpreted by strontium isotopes, Hydrol. Res., 20, 85–96,, 1989. 

Jomori, Y., Minami, M., Sakurai-Goto, A., and Ohta, A.: Comparing the 87Sr/86Sr of the bulk and exchangeable fractions in stream sediments: implications for 87Sr/86Sr mapping in provenance studies, Appl. Geochem., 86, 70–83,, 2017. 

Jweda, J., Bolge, L., Class, C., and Goldstein, S. L.: High precision Sr-Nd-Hf-Pb isotopic compositions of USGS reference material BCR-2, Geostand. Geoanal. Res., 40, 101–115,, 2016. 

Kennedy, M. J., Hedin, L. O., and Derry, L. A.: Decoupling of unpolluted temperate forests from rock nutrient sources revealed by natural 87Sr/86Sr and 84Sr tracer addition, P. Nat. Acad. Sci. USA, 99, 9639–9644,, 2002. 

Knudson, K. J., Webb, E., White, C., and Longstaffe, F. J.: Baseline data for Andean paleomobility research: a radiogenic strontium isotope study of modern Peruvian agricultural soils, Archaeol. Anthropol. Sci., 6, 205–219,, 2014. 

Le Bas, M., Spiro, B., and Yang, X.: Oxygen, carbon and strontium isotope study of the carbonatitic dolomite host of the Bayan Obo Fe-Nb-REE deposit, Inner Mongolia, N China, Mineral. Mag., 61, 531–541,, 1997. 

Lech, M. E., de Caritat, P., and Mcpherson, A. A.: National Geochemical Survey of Australia: Field Manual, Record, 2007/08, Geosci. Austral., Canberra, (last access: 27 February 2023), 2007. 

Lugli, F., Cipriani, A., Bruno, L., Ronchetti, F., Cavazzuti, C., and Benazzi, S.: A strontium isoscape of Italy for provenance studies, Chem. Geol., 587, 120624,, 2022. 

Lyons, W. B., Tyler, S. W., Gaudette, H. E., and Long, D. T.: The use of strontium isotopes in determining groundwater mixing and brine fingering in a playa spring zone, Lake Tyrrell, Australia, J. Hydrol., 167, 225–239,, 1995. 

McArthur, J. M., Howarth, R. J., and Shields, G. A.: Strontium Isotope Stratigraphy, Chap. 7 in: The Geologic Time Scale, edited by: F. M., Ogg, J. G., Schmitz, M., and Ogg, G., Elsevier BV, 127–144,, 2012. 

McNutt, R. H.: Strontium Isotopes, in: Environmental Tracers in Subsurface Hydrology, edited by: Cook P. G. and Herczeg A. L., Springer, Boston, MA,, 2000. 

McNutt, R. H., Frape, S. K., and Dollar, P.: A strontium, oxygen and hydrogen isotopic composition of brines, Michigan and Appalachian Basins, Ontario and Michigan, Appl. Geochem., 2, 495–505,, 1987. 

Moffat, I., Rudd, R., Willmes, M., Mortimer, G., Kinsley, L., McMorrow, L., Armstrong, R., Aubert, M., and Grün, R.: Bioavailable soil and rock strontium isotope data from Israel, Earth Syst. Sci. Data, 12, 3641–3652,, 2020. 

Mountjoy, E. W., Qing, H., and McNutt, R. H.: Strontium isotopic composition of Devonian dolomites, Western Canada Sedimentary Basin: significance of sources of dolomitizing fluids, Appl. Geochem., 7, 59–75, 1992. 

Nebel, O. and Stammeier, J. A.: Strontium Isotopes, in: Encycl. Geochem., Encycl. Earth Sci. Series, edited by: White W. M., Springer, Cham,, 2018. 

Négrel, P. and Grosbois, C.: Changes in chemical and 87Sr/86Sr signature distribution patterns of suspended matter and bed sediments in the upper Loire river basin (France), Chem. Geol., 156, 231–249,, 1999. 

Négrel, P. and Pauwels, H.: Interaction between different groundwaters in Brittany catchments (France): characterizing multiple sources through strontium- and sulphur isotope tracing, Wat. Air Soil Poll., 151, 261–285,, 2004. 

Oishi, Y.: Is the Sr isotope ratio of mosses a good indicator for Asian dust (Kosa)?, Landscape Ecol. Eng., 18, 11–17,, 2021. 

Ojiambo, S. B., Lyons, W. B., Welch, K. A., Poreda, R. J., and Johannesson, K. H.: Strontium isotopes and rare earth elements as tracers of groundwater-lake water interactions, Lake Naivasha, Kenya, Appl. Geochem., 18, 1789–1805,, 2003. 

Oliver, L., Harris, N., Bickle, M., Chapman, H., Dise, N., and Horstwood, M.: Silicate weathering rates decoupled from the 87Sr/86Sr ratio of the dissolved load during Himalayan erosion, Chem. Geol., 201, 119–139,, 2003. 

Ollier, C. D.: The regolith in Australia, Earth-Sci. Rev., 25, 355–361,, 1988. 

Pacheco-Forés, S. I., Gordon, G. W., and Knudson, K. J.: Expanding radiogenic strontium isotope baseline data for central Mexican paleomobility studies, PLoS ONE, 15, e0229687,, 2020. 

Page, R. W., Needham, R. S., and Compston, W.: Geochronology and evolution of the Late-Archaean basement and Proterozoic rocks in the Alligator Rivers uranium field, Northern Territory, Australia, in: Proc. Int. Uranium Symp., Pine Creek Geosyncline (Sydney, New South Wales, 4–8 June 1979), edited by: Ferguson, J. and Goleby, A. B., Int. Atom. En. Agency, Vienna, Austria, 39–68,, last access: 27 February 2023, 1980. 

Pain, C., Gregory, L., Wilson, P., and McKenzie, N.: The Physiographic Regions of Australia – Explanatory Notes, Australian Collaborative Land Evaluation Program (ACLEP) and National Committee on Soil and Terrain (NCST), Canberra, 30 pp., (last access: 27 February 2023), 2011. 

Palmer, M. R., Helvaci, C., and Fallick, A. E.: Sulphur, sulphate oxygen and strontium isotope composition of Cenozoic Turkish evaporates, Chem. Geol., 209, 341–356,, 2004. 

Plumlee, G.: Basalt, Columbia River, BCR-2, Prelim. U.S. Geol. Surv. Cert. Anal., (last access: 27 February 2023), 1998. 

Price, G. J., Ferguson, K. J., Webb, G. E., Feng, Y. X., Higgins, P., Nguyen, A. D., Zhao, J. X., Joannes-Boyau, R., and Louys, J.: Seasonal migration of marsupial megafauna in Pleistocene Sahul (Australia-New Guinea), P. Roy. Soc. B-Biol. Sci., 284, 20170785,, 2017. 

Probst, A., El Gh'Mari, A., Aubert, D., Fritz, B., and McNutt, R.: Strontium as a tracer of weathering processes in a silicate catchment polluted by acid atmospheric inputs, Strengbach, France, Chem. Geol., 170, 203–219,, 2000. 

Quade, J., Chivas, A. R., and McCulloch, M. T.: Strontium and carbon isotope tracers and the origins of soil carbonate in South Australia and Victoria, Palaeogeogr. Palaeocl. Palaeoecol., 113, 103–117,, 1995. 

Raymond, O. L., Gallagher, R., Shaw, R., Yeates, A. N., Doutch, H. F., Palfreyman, W. D., Blake, D. H., and Highet, L.: Surface Geology of Australia 1 : 2.5 million scale dataset 2012 edition, edited by: Raymond, O. L. and Gallagher, R., Geosci. Austral., Canberra [data set],, 2012. 

Revel-Rolland, M., De Deckker, P., Delmonte, B., Hesse, P. P., Magee, J. W., Basile-Doelsch, I., Grousset, F., and Bosch, D.: Eastern Australia: a possible source of dust in East Antarctica interglacial ice, Earth Planet. Sc. Lett., 249, 1–13,, 2006. 

Riley, G. H.: Granite ages in the Pine Creek Geosyncline, in: Proc. Int. Uranium Symp., Pine Creek Geosyncline, 4–8 June 1979, Sydney, New South Wales, edited by: Ferguson, J. and Goleby, A. B., Int. Atom. En. Agency, Vienna, Austria, 69–72,, last access: 27 February 2023, 1980. 

Romaniello, S. J., Field, M. P., Smith, H. B., Gordon, G. W., Kim, M. H., and Anbar, A. D.: Fully automated chromatographic purification of Sr and Ca for isotopic analysis, J. Anal. Atomic Spectro., 30, 1906–1912,, 2015. 

Rotenberg, E., Davis, D. W., Amelin, Y., Ghosh, S., and Bergquist, B. A.: Determination of the decay-constant of 87Rb by laboratory accumulation of 87Sr, Geochim. Cosmochim. Ac., 85, 41–57,, 2012. 

Rundberg, Y. and Smalley, P. C.: High-resolution dating of Cenozoic sediments from northern North Sea using 87Sr/86Sr stratigraphy, Am. Assoc. Petrol. Geol. Bull., 73, 298–308,, 1989. 

Schaltegger, U., Stille, P., Rais, N., Pique, A., and Clauer, N.: Neodymium and strontium isotopic dating of diagenesis and low-grade metamorphism of argillaceous sediments, Geochim. Cosmochim. Ac., 58, 1471–1481,, 1994. 

Schoneveld, L., Spandler, C., and Hussey, K.: Genesis of the central zone of the Nolans Bore rare earth element deposit, Northern Territory, Australia, Contrib. Mineral. Petrol., 170, 11,, 2015. 

Schultz, J. L., Boles, J. R., and Tilton, G. R.: Tracking calcium in the San Joaquin basin, California: a strontium isotopic study of carbonate cements at North Coles Levee, Geochim. Cosmochim. Ac., 53, 1991–1999,, 1989. 

Shewan, L. G., Armstrong, R. A., and O'Reilly, D.: Baseline bioavailable strontium isotope values for the investigation of residential mobility and resource-acquisition strategies in prehistoric Cambodia, Archaeom., 62, 810–826,, 2020. 

Shields, G. A.: A normalised seawater strontium isotope curve and the Neoproterozoic-Cambrian chemical weathering event, eEarth Discuss., 2, 69–84, 2007. 

Shin, W.-J., Ryu, J.-S., Kim, R.-H., and Min, J.-S.: First strontium isotope map of groundwater in South Korea: applications for identifying the geographical origin, Geosci. J., 25, 173–181,, 2021. 

Simandl, G. J. and Paradis, S.: Vanadium as a critical material: economic geology with emphasis on market and the main deposit types, Appl. Earth Sci., 131, 218–236,, 2022. 

Smalley, P. C., Råheim, A., Rundberg, Y., and Johansen, H.: Strontium-isotope stratigraphy: applications in basin modelling and reservoir correlation, chap. 3, in: Correlation in Hydrocarbon Exploration, Norwegian Petrol. Soc., Graham and Trotman, 23–31,, 1989. 

Smalley, P. C., Lønøy, A., and Råheim, A.: Spatial 87Sr/86Sr variations in formation water and calcite from the Ekofisk chalk oil field: implications for reservoir connectivity and fluid composition, Appl. Geochem., 7, 341–350,, 1992. 

Stueber, A. M., Pushkar, P., and Hetherington, E. A.: A strontium isotopic study of formation waters from the Illinois basin, USA, Appl. Geochem., 2, 477–494,, 1987. 

Sullivan, M. D., Haszeldine, R. S., and Fallick, A. E.: Linear coupling of carbon and strontium isotopes in Rotliegend Sandstone, North Sea: evidence for cross-formational fluid flow, Geology, 18, 1215–1218,<1215:LCOCAS>2.3.CO;2, 1990. 

Swart, P. K., Ruiz, J., and Holmes, C. W.: Use of strontium isotopes to constrain the timing and mode of dolomitization of upper Cenozoic sediments in a core from San Salvador, Bahamas, Geology, 15, 262–265,<262:UOSITC>2.0.CO;2, 1987. 

Tukey, J. W.: Exploratory Data Analysis, Addison-Wesley Publishing Company, Reading, MA, 506 pp., 1977. 

Van Kranendonk, M. J., Hickman, A. H., Smithies, R. H., Nelson, D. R., and Pike, G.: Geology and tectonic evolution of the Archean North Pilbara Terrain, Pilbara Craton, Western Australia, Econ. Geol., 97, 695–732,, 2002. 

Veizer, J.: Strontium isotopes in seawater trough time, Ann. Rev. Earth Planet. Sci., 17, 141–167,, 1989. 

Vinciguerra, V., Stevenson, R., Pedneault, K., Poirer, A., Hélie, J.-F., and Widory, D.: Strontium isotope characterization of Wines from the Quebec (Canada) Terroir, Proc. Earth Planet. Sci., 13, 252–255,, 2015. 

Voerkelius, S., Lorenz, G. D., Rummel, S., Quétel, C. R., Heiss, G., Baxter, M., Brach-Papa, C., Deters-Itzelsberger, P., Hoelzl, S., Hoogewerff, J., Ponzevera, E., Bocxstaele, M., and Van Ueckermann, H.: Strontium isotopic signatures of natural mineral waters, the reference to a simple geological map and its potential for authentication of food, Food Chem., 118, 933–940,, 2010. 

Vorster, C., Greeff, L., and Coetzee, P. P.: The determination of 11B/10B and 87Sr/86Sr isotope ratios by quadrupole-based ICP-MS for the fingerprinting of South African wine, South Afric. J. Chem., 63, 207–214, 2010. 

Washburn, E., Nesbitt, J., Ibarra, B., Fehren-Schmitz, L., and Oelze, V. M.: A strontium isoscape for the Conchucos region of highland Peru and its application to Andean archaeology, PLOS ONE 16, e0248209,, 2021. 

Wei, C.-W., Xu, C., Chakhmouradian, A. R., Brenna, M., Kynicky, J., and Song, W.-L.: Carbon-strontium isotope decoupling in carbonatites from Caotan (Qinling, China): implications for the origin of calcite carbonatite in orogenic settings, J. Petrol., 61, egaa024,, 2020. 

Wei, X., Wang, S., Ji, H., and Shi, Z.: Strontium isotopes reveal weathering processes in lateritic covers in southern China with implications for paleogeographic reconstructions, PLOS ONE, 13, e0191780,, 2018. 

Wickman, T. and Jacks, G.: Strontium isotopes in weathering budgeting. In: Water-Rock Interaction, Proc. 7th Int. Symp. Water-Rock Interaction, 13–18 July 1992, Park City, Utah, edited by: Kharaka, Y. K. and Maest, A. S., Balkema, Rotterdam, 1, 611–614, ISBN 9054100753, 1992. 

Wilford, J.: A weathering intensity index for the Australian continent using airborne gamma-ray spectrometry and digital terrain analysis, Geoderma, 183–184, 124–142,, 2012. 

Willmes, M., McMorrow, L., Kinsley, L., Armstrong, R., Aubert, M., Eggins, S., Falguères, C., Maureille, B., Moffat, I., and Grün, R.: The IRHUM (Isotopic Reconstruction of Human Migration) database – bioavailable strontium isotope ratios for geochemical fingerprinting in France, Earth Syst. Sci. Data, 6, 117–122,, 2014. 

Willmes, M., Bataille, C. P., James, H. F., Moffat, I., McMorrow, L., Kinsley, L., Armstrong, R. A., Eggins, S., and Grün, R.: Mapping of bioavailable strontium isotope ratios in France for archaeological provenance studies, Appl. Geochem., 90, 75–86,, 2018. 

Withnall, I. W., Hutton, L. J., Armit, R. J., Betts, P. G., Blewett, R. S., Champion, D. C., and Jell, P. A.: North Australian Craton. In: Geology of Queensland, edited by: Jell, P. A., Geol. Surv. Queensland, Brisbane, 23–112, ISBN 9781921489761, (last access: 27 February 2023), 2013.  

Wyborn, L. A. I., Heinrich, C. A., and Jaques, A. L.: Australian Proterozoic mineral systems: essential ingredients and mappable criteria, Austral. Inst. Mining Metall. (AusIMM) Ann. Conf. Proc., 5–9 August 1994, Darwin, Northern Territory, 109–115, (last access: 27 February 2023), 1994. 

Yang, C., Telmer, K., and Veizer, J.: Chemical dynamics of the “St. Lawrence” riverine system: δDH2O, δ18OH2O, δ13CDIC, δ34Ssulfate, and dissolved 87Sr/86Sr, Geochim. Cosmochim. Ac., 60, 851–866,, 1996. 

Zhao, Z., Leach, D. L., Wei, J., Liang, S., and Pfaff, K.: Origin of the Xitieshan Pb-Zn deposit, Qinghai, China: evidence from petrography and S-C-O-Sr isotope geochemistry, Ore Geol. Rev., 139A, 104429,, 2021. 

Short summary
This new, extensive (~1.5×106 km2) dataset from northern Australia contributes considerable new information on Australia's strontium (Sr) isotope coverage. The data are discussed in terms of lithology and age of the source areas. This dataset will reduce Northern Hemisphere bias in future global Sr isotope models. Other potential applications of the new data include mineral exploration, hydrology, food tracing, dust provenancing, and examining historic migrations of people and animals.
Final-revised paper