Multiple ages of rutile from a single sample of granulite

Durgalakshmi1, Ian S. Williams1, K. Sajeev2

1 Research School of Earth Sciences, Australian National University, Canberra, Australia, 2 Centre for Earth Sciences, Indian Institute of Science, Bengaluru, India

Rutile (TiO2) is a common accessory mineral in hydrothermal and metamorphic rocks that is stable across a wide range of P-T conditions. It can incorporate up to 200 ppm of U, and has a lower closure temperature than zircon, making it a reliable mineral with which to date retrograde metamorphism and low- to medium-grade metamorphic events by the U-Pb technique. Further, Pb and Th are incompatible in rutile, so corrections for its small initial Pb can be made accurately using its 208Pb content.

During prograde metamorphism, rutile forms following the breakdown of biotite or ilmenite as part of a continuous reaction. During retrogression, rutile can be replaced by ilmenite and titanite. Secondary rutile can be formed by hydrothermal alteration, oxidation or exsolution. In high-grade metamorphic rocks, rutile occurs as single crystals in the matrix and/or as inclusions in other minerals such as garnet, pyroxene and amphibole. In low- to medium-grade rocks it usually occurs as needles or polycrystalline aggregates. Using an ion microprobe (SIMS), rutile can be dated in its textural context in thin section, providing age information directly linked to metamorphic reactions. Dating and trace element analysis of separated rutile grains have a wide range of applications from sedimentary provenance studies to dating vein mineralisation and granitic pegmatites that host mineral deposits.

In the Neoarchæan Southern Granulite Terrane, India, rutile preserves a record of the late thermal history that is not provided by other datable minerals such as zircon or monazite. One studied sample of granulite grade felsic gneiss contains at least three distinct generations of rutile that preserve a range of Neoproterozoic to early Palæozoic U-Pb ages. The zircon from the same sample is Neoarchæan, with no evidence of a younger component. The rutile occurs as single ~ 0.5–0.8 mm crystals in the matrix, some of which are rimmed by titanite. Rutile-forming reactions, which can be linked to the metamorphic conditions, have been dated, contributing to unravelling the polymetamorphic and tectonic history of this complex terrane.


Miss Durgalakshmi is a PhD student at Research School of Earth Sciences, ANU. She works on the Archaean rocks of Southern Granulite Terrane, India.

Way out west – does the Arunta Orogen continue westward beneath the Canning Basin?

Kelsey, David E.1,3, Spaggiari, Catherine V.1, Wingate, Michael T.D.1, Lu, Yongjun1, Fielding, Imogen O.H.1, Finch, Emily G.1,2,3

1Geological Survey of Western Australia, 100 Plain St, East Perth, WA 6004, Australia, 2University of South Australia, 101 Currie St, Adelaide, SA 5001, Australia, 3MinEx CRC

Crystalline basement lies beneath the southeastern margin of the Canning Basin and immediately west of the exposed Arunta Orogen, although whether that basement is the continuation of the Arunta Orogen is unknown. The major bounding fault of the Canning Basin coincides with the inferred trace of the Lasseter Shear Zone, which truncates dominant east-trending structures of the orogen and potentially also terminates it. The Top Up Rise prospect is located above a distinct northeast-trending gravity anomaly bound by northeast-trending shear zones coincident with the Lasseter Shear Zone. Five co-funded Exploration Incentive Scheme diamond drillcores from Top Up Rise contain partially melted or melt-injected upper-amphibolite to low-granulite facies basement rocks. These are currently the only drillcores that intersect the Canning Basin basement in this region and provide a means to test its tectonic affinity, and the significance of the Lasseter Shear Zone.

Petrography of granitic gneiss and paragneiss indicates distinctly higher metamorphic grade than is observed in exposed rocks of the west Arunta Orogen, consistent with their separation by a significant shear zone. SHRIMP U–Pb zircon geochronology of several granitic gneiss samples has so far revealed a single magmatic protolith crystallization age of c. 1875 Ma, which is distinctly older than known magmatic rocks in the Arunta Orogen. Maximum deposition ages of 1877–1822 Ma for the metasedimentary rocks are similar to or younger than the magmatic protolith ages of the granitic gneiss, suggesting that emplacement of granitic rocks predated deposition of at least some of the sedimentary rocks. Zircon rims in both granitic gneiss and paragneiss samples record high-grade metamorphism at 1622–1604 Ma. Although metamorphism of this age is unknown in the west Arunta Orogen, thermal events of similar age may have occurred in the central Arunta Orogen (Alessio et al., 2020, Lithos,, and granitic magmatism occurred at this time in the Haast Bluff Domain of the Warumpi Province (NTGS Special Publication 5, 2013).

The granitic protolith ages of c. 1875 Ma from Top Up Rise are different to those of granitic suites in the Lamboo Province and Granites–Tanami Orogen; however, felsic volcanic rocks of the lower Halls Creek Group have been dated at c. 1880 Ma (Phillips et al., 2016, GSWA Report 164) and c. 1880 Ma granitic rocks are known in the Arnhem Province (NTGS Record 2017-008). In contrast, detrital zircon age components at c. 1875 Ma are widespread and form significant age components in the Lander Rock Formation of the Arunta Orogen, the Marboo Formation and Tickalara Metamorphics of the Halls Creek Orogen, as well as in Top Up Rise paragneisses. Hence, the newly identified granitic basement in the Top Up Rise drillcores may be representative of a major source component that fed detritus into turbidite fan systems that included the Lander Rock Formation. This basement could represent part of the ‘missing’ Arunta basement, which would support interpretations that the Arunta does continue westwards, or it could be part of a previously unrecognized Proterozoic crustal element that underlies the Canning Basin.


The author team is made up of geologists and geochronologists from the GSWA. Emily Finch is a MinEx CRC Embedded Researcher at the GSWA.

Neodymium and oxygen isotope maps of Western Australia

Lu, Yongjun1, Smithies, RH1, Champion, DC2, Wingate, MTD1, Johnson, SP1, Martin, L3, Jeon, H4, Poujol, M5, Zhao, J6, Maas, R7, Creaser, RA8

1Geological Survey of Western Australia, 100 Plain Street, East Perth, WA 6004, 2Geoscience Australia, GPO Box 378, Canberra ACT 2601, 3Centre for Microscopy, Characterisation, and Analysis, University of Western Australia, Perth, WA 6009, 4Swedish Museum of Natural History, Box 50 007, SE-104 05 Stockholm, Sweden, 5Univ Rennes, CNRS, Géosciences Rennes – UMR 6118, 35000 Rennes, France, 6Radiogenic Isotope Facility, School of Earth Sciences, The University of Queensland, Brisbane, QLD 4072, Australia, 7School of Earth Sciences, University of Melbourne, Parkville, VIC 3010, Australia, 8Dept. Earth & Atmospheric Sciences, University of Alberta, Edmonton, Alberta, Canada

Multi-isotope maps can characterise lithospheric architecture through time, play an increasingly important role in predictive exploration targeting, and are consequently sought-after datasets by industry. We present the first zircon oxygen isotope map and an updated whole-rock Sm–Nd isotope map of Western Australia. These shed new light on crustal evolution and mineral system distributions. The isotope maps were generated from datasets that are subject to ongoing updates as new data are generated and compiled.

Median zircon δ18O values for about 125 igneous rocks have been spatially visualized so far, and coverage currently extends across the Pilbara and Yilgarn Cratons, the Capricorn, Paterson, and Albany-Fraser Orogens and the Eucla basement (Madura and Coompana Provinces). The Pilbara and Yilgarn Cratons are dominated by mantle-like δ18O values (4.7 – 5.9‰), consistent with reworking of igneous material that had not been exposed at the surface. A c. 3.47 Ga diorite, four c. 3.3 Ga hornblende-bearing granitic rocks, and a c. 2.95 Ga hornblende monzogranite in the Pilbara Craton exhibit weakly elevated zircon δ18O values (5.9 – 6.5‰), which together with trace element enrichment were attributed to hydrous sanukitoids or to derivation from a sanukitoid-enriched source. The Capricorn, Paterson, and Albany–Fraser Orogens and the Eucla basement also contain rocks with elevated δ18O values (6.6 – 9.9‰), suggesting significant reworking of upper crustal material during magma genesis. Zircons with sub-mantle δ18O values (<4.7‰) were found for granitic rocks of c. 3.55 Ga in the Sylvania Inlier, of c. 3.44 Ga in the northern Pilbara Craton, and of c. 3.0 and 2.67 Ga in the South West Terrane, suggesting recycling of crustal material subjected to high-temperature hydrothermal alteration, such as observed in post-Archean rift systems or calderas.

Sm–Nd isotopes for about 1120 felsic igneous rocks provide regionally extensive images of crustal architecture. The map of two-stage depleted mantle model ages (TDM2) highlights the distinction between Archean cratons (TDM2 >2.6 Ga) and Proterozoic orogens (TDM2 <2.2 Ga), and isotopic boundaries correlate well with most existing proposed terrane boundaries. However, the isotopic boundary between the South West Terrane and the Youanmi Terrane appears to be about 100 km west of the previously proposed boundary, but correlates well with magnetic and gravity anomaly zones and the distribution of gold mineralization. The crustal residence map highlights predominantly short residence times (<0.5 Ga) for the Pilbara and Yilgarn Cratons, and much longer crustal residence times (>0.8 Ga) in the Paterson, Albany–Fraser, Pinjarra and Capricorn Orogens, suggesting decreased juvenile crust generation in these orogens. 

These maps are directly applicable to metallogeny. For example, most giant gold deposits in WA are located on or near significant isotopic boundaries and tectonic structures. Interestingly, Telfer, Plutonic and giant gold deposits in the Murchison are aligned along a northeast-trending isotopic boundary. Similar boundaries occur between the eastern and western parts of the Pilbara Craton and between the Yilgarn Craton and the Albany–Fraser Orogen. These isotopically defined discontinuities may be important clues to the earliest architectural elements in Western Australia.


Dr. Yongjun Lu is Senior Geochronologist Isotope Geologist at GSWA. Yongjun has over 15 years’ geological experience, extensive collaboration with industry, government and academia. He has made important contributions to mineral systems science, highlighted by being the 52nd recipient of the Waldemar Lindgren Award of Society of Economic Geologists (SEG).

Records of the Earth’s early crust from apatite inclusions in zircon – development and applications of in situ 87Sr/86Sr analysis by SIMS

Gillespie, Jack1, Kinny, Pete1, Martin, Laure2, Kirkland, Christopher1, Nemchin, Alexander1, Cavosie, Aaron1

1The Institute for Geoscience Research (TIGeR), School of Earth and Planetary Sciences, Curtin University, Perth, WA 6845, Australia, 2Center for Microscopy, Characterisation and Analysis (CMCA), University of Western Australia, 35 Stirling Highway, Perth, WA 6009, Australia

Rb-Sr isotopes in geological materials provide a system for tracing crustal differentiation processes and insights into the evolution of planetary bodies. The ingrowth of radiogenic 87Sr from the decay of 87Rb leads to increased 87Sr/86Sr over time, and due to the strong relationship between the Rb/Sr and SiO2 contents of igneous rocks, this provides a time-integrated window into the evolution of geochemical reservoirs. However, the high geological mobility of both Rb and Sr in whole rocks during metamorphism and fluid alteration means that this record becomes progressively less reliable in older rocks that have experienced post-crystallization geological events.

Strontium is easily substituted into the crystal lattice of apatite, occurring as a trace element in concentrations ranging from less than a hundred parts per million to several weight percent. In contrast, apatite nearly entirely excludes Rb (<1 ppm) resulting in negligible radiogenic ingrowth of 87Sr, and consequently the initial 87Sr/86Sr ratio of the melt from which an apatite crystallizes is faithfully recorded by the mineral. Inclusions of apatite within magmatic zircons are particularly valuable as they are armoured by the more robust host mineral, allowing them to survive subsequent events that might otherwise cause isotopic reset or recrystallization. However, the typically very small size of apatite inclusions in zircon and the complex isobaric interferences on the isotopes of Sr during in-situ mass spectrometry have previously limited the information that can be obtained from this archive.

Our recent development of a method to measure the 87Sr/86Sr ratio in apatite by SIMS with a spot size appropriate for accessing typical mineral inclusions in zircon (<15 µm) makes it possible to routinely analyse the commonly occurring inclusions of apatite in zircon. We have applied this method to determine the initial 87Sr/86Sr ratios of various Eo-Meso Archean igneous rocks by analysing the Sr isotope composition of apatite inclusions. High resolution SEM imaging and EPMA analysis illustrate the primary nature of these inclusions. Combining the measured 87Sr/86Sr of apatite inclusions with the U-Pb age and Hf isotopic composition of the co-genetic zircon host allows for the ‘triangulation’ of the Rb/Sr necessary for the ingrowth of radiogenic strontium over the crustal residence interval calculated from the crystallization and Hf model ages. Examples from SW Greenland and the Narryer Gneiss Terrane of Western Australia suggest that these rocks were derived from the melting of ancient crustal material that was on average of intermediate-felsic rather than mafic composition.


Jack Gillespie is a post-doctoral researcher at Curtin University working on developing new methods for understanding the evolution of the early earth

Deriving Quantitative Alteration Mineralogy from TIR Hyperspectral Data in IOCG Systems

1Stromberg, Jessica, 1Schlegel, Tobias, 1Pejcic, Bobby, 1Birchall, Renee, 1Shelton, Tina

1CSIRO Mineral Resources, Kensington, Australia

Identifying alteration mineral zonation around hydrothermal ore systems is critical to the mineral exploration process. Hyperspectral methods are commonly used to map alteration because they are fast, inexpensive, and require little to no sample preparation compared to other mineralogical techniques such as scanning electron microscope (SEM) based mineral mapping, for example. The visible-near and shortwave infrared (VN-SWIR, 350-2500 nm) spectral regions are most used as they are sensitive to hydrated mineral phases including chlorite and white micas. However, in mineral systems which are iron oxide-rich or where key alteration assemblages include abundant anhydrous phases, such as in iron oxide-alkali-calcic alteration systems, using this spectral range can be problematic. In this case, the thermal-infrared (TIR) spectral range (6000-14500 nm) may be more appropriate as it is sensitive to anhydrous silicates such as quartz and feldspars. However, there are inherent challenges in unmixing hyperspectral data for deriving quantitative mineral abundances. In particular, when the quantification of minor or spectrally similar phases are key to the alteration assemblages, such as in the quantification of different feldspars. This can be overcome with the use of a calibration dataset such as quantitative X-ray diffraction or SEM-based quantitative mineralogy in combination with Partial-Least Square (PLS) regression methods. In developing such models, scale is of critical importance. Variability in the sampling area and volume between datasets is one of the greatest challenges in validating hyperspectral data with quantitative mineralogy, and in integrating any geoscience datasets. In this work, we used several hyperspectral and spectroscopic instruments (HyLogger, Agilent 4300 FTIR, Bruker Vertex FTIR, ASD Fieldspec Pro) to evaluate the impact of scale on hyperspectral data validation in IOCG systems. In this process we developed a methodology for creating the first scale-consistent dataset of VNIR-SWIR, TIR, and SEM-based quantitative mineralogy data on drill core samples. This dataset comprises 250 samples from a world-class IOCG deposit in which the key mineral phases, assemblages, and alteration patterns were identified using the SEM-based quantitative mineralogy. Hyperspectral data was processed using The Spectral Geologist Software (TSGTM) software and even with a scale consistent dataset, and a constrained mineral library based on the SEM-based mineralogy, conventional unmixing methods such as the The Spectral Assistant (TSATM) were unable to reproduce key alteration patterns for vectoring towards ore. Using the SEM-based mineralogy data, PLS modelling was applied to derive predictive models for key mineral phases from the TIR hyperspectral data. The resulting models produced quantitative mineralogy with r2 > 0.94 for key phases including quartz, K-feldspar, albite, calcite, and >0.80 for magnetite, biotite, and plagioclase. More importantly, the key mineral assemblages change with distance to ore and the relative abundance feldspar species (albite, K-feldspar, plagioclase) identified by the SEM-based mineral mapping were reproduced using PLS-derived mineral abundances in a validation drill core. This method and the models generated will provide a framework for improving the application of hyperspectral data for mapping alteration in IOGC systems.


Jessica is a Research Scientist in the Mineral Footprints team at CSIRO Mineral Resources where she leads Activity 5 of the NVCL and works on projects across multiple commodities applying the combined use of lab and field-based spectroscopic and geochemical techniques to address industry-specific challenges.

Modelling of petrophysical from hyperspectral drill core data collected from the Osborne Cu-Au deposit, Mount Isa Inlier, Queensland

Laukamp, Dr Carsten1, Francis, Neil1, Hauser, Dr Juerg1, Gopalakrishnan, Dr Suraj2, Mule, Shane1

1CSIRO, Perth, Australia, 2Geologial Survey of Queensland, Brisbane, Australia

The combination of magnetic susceptibility and density allows identification of iron oxide copper-gold (IOCG) mineralisation by estimating proportions of magnetite, sulphide and hematite alteration which can be indicative of IOCGs. In the frame of the National Virtual Core Library project, hyperspectral reflectance spectra acquired from drill core of the Osborne Cu-Au deposit, Mount Isa Inlier, Queensland using a HyLogger3 at GSQ’s Exploration Data Centre were compared with magnetic susceptibility and density measurements. In this study we explore the feasibility of inferring the petrophysical data from the 1) visible-near (VNIR), shortwave (SWIR) and 2) thermal (TIR) infrared wavelength regions. Specifically, we seek to predict magnetic susceptibility and density values in drill core sections where petrophysical data are not available and potentially extrapolate these to other hyperspectral data sets, such as those acquired by field or Earth Observation instruments.

Using The Spectral Geologist (TSGTM) software, partial least squares (PLS) was employed to derive predictive models using 23 unique magnetic susceptibility and density measurements, that were assigned to all nearby spectral measurements (+/- ~10cm). The values of the input magnetic sustainability and density values ranged from 0 to 2.3 K (Si) and 2.7 to 4.9 g/cm3, respectively. The hyperspectral data were not spatially re-sampled to fit with the drill core interval measured for petrophysical data. Instead, the original 1 cm spatial resolution was used to evaluate the variability of hyperspectral data within the petrophysical sample intervals. The correlation between the 23 measured and corresponding modelled magnetic susceptibility (n = 157) for the same 23 depth intervals was high for the VNIR-SWIR (r2 = 0.95) and the TIR (r2 = 0.969), but the PLS-modelled magnetic susceptibility values showed a large variance (± 0.8 and ± 0.5, respectively). Similarly, the correlation between the measured and modelled density was high for the VNIR-SWIR (r2 = 0.958) and the TIR (r2 = 0.989), with the PLS-modelled density values showing a large variance (± 0.4 g/cm3 for both wavelength ranges). However, HyLogger3 high-resolution RGB imagery showed that the predicted value ranges were sufficiently different to discriminate drill core intervals dominated by magnetite-rich rocks, from magnetite-rich breccia and least-altered (non-mineralised) rocks. PLS models based on the VNIR-SWIR wavelength ranges were mainly driven by depth changes of electronic transition absorption features related to iron and copper in the VNIR, which are most intense in the highly altered, magnetite- and/or sulphide rich rocks. PLS models based on the TIR wavelength ranges were highly influenced by the thermal background typically associated with iron oxides and sulphides. Density values modelled from VNIR-SWIR compared to those modelled from TIR showed a good correlation (r2 = 0.729), whereas the correlation between magnetic susceptibility modelled from VNIR-SWIR and TIR was comparably low (r2 = 0.518). While the small amount of data used to infer the models discussed here means that their predictive power needs to be assessed comprehensively, our results nevertheless indicate a high potential for successfully inferring petrophysical from hyperspectral data and cost-effective mapping of IOCG-related alteration. 


Carsten Laukamp is a senior research geoscientist at CSIRO Mineral Resources, based in Perth, Australia and is project leader of the National Virtual Core Library. Carsten explores the potential for combined use of reflectance spectroscopy, geochemistry and geophysics for tracing hydrothermal alteration signatures through cover and advancing ore body knowledge.

HIGH RESOLUTION VS. STANDARD RESOLUTION: How an increase in spectral resolution using a field portable spectrometer affects quality of data – a case study on Nickel exploration

Shelton Pieniazek, Lori1

1Spectral Evolution, United States

Standard resolution of a field portable spectrometer at 3nm (UV), 8nm (VIS) and 6nm(NIR) has been used successfully in the geological remote sensing and mining industry for a variety of applications. Alteration mapping has been key for finding different mineral assemblages that indicate an area of interest. Using a higher resolution field spectrometer has proven to be beneficial in understanding and identifying new features in a variety of different mineral groups and subgroups.

A study was conducted using three different field spectrometers, each with different resolutions. As you increase resolution, features that are not visible with a standard resolution spectrometer appear in all three ranges (UV,VIS,NIR) thus yielding a better understanding of alteration changes and geochemical conditions


Mrs. Lori Shelton Pieniazek is the Geologist for Spectral Evolution. She received her degree in Geosciences from Tennessee Technological University. Her main focus is on applying spectral geology applications to the mining and geological remote sensing industries.

HyLoggerTM mineralogy from chips: A RoXplorer pilot study

Moltzen, Jake1

1Geological Survey Of New South Wales, Londonderry, Australia

The RoXplorer® coiled tube drilling system is scheduled to be used for MinEx CRC drilling activities in New South Wales (NSW) from 2022 for National Drilling Initiative (NDI) areas South Cobar, North Cobar, Mundi, Forbes and Dubbo. This recently developed system provides a faster and more cost-effective alternative to conventional diamond and reverse circulation drilling, while also offering safety advantages and reduced environmental impacts (Hillis et al. 2014). Ongoing technology development within MinEx CRC will maximise RoXplorer® production during MinEx CRC drilling campaigns, with 5 km spaced holes planned throughout NDI areas to better understand basement and cover sequence geology and the potential for undiscovered ore deposits buried beneath cover.

In this pilot study, two trial drillholes from the RoXplorer® system were scanned using the HyLogger-3TM instrument to determine the best method for scanning RoXplorer® chips and to assess the quality of visible-near infrared to shortwave infrared (VNIR-SWIR) and thermal infrared (TIR) spectra for routine interpretation. Downhole chip samples produced by the RoXplorer® range in size from ~10 µm to ~4 mm (Tiddy et al. 2019), however a separated coarse fraction with sizes ranging from ~1 mm to ~4 mm was used for both drillholes. Outcomes from this study included:

  • VNIR‒SWIR and TIR spectra were found to be reliable for routine mineral identification, tracking mineral chemistry and identifying lithological variations downhole.
  • Black plastic chip trays produced better spectral results than clear or white trays.
  • Larger 50 mm x 50 mm chip compartments produced better spectral results than standard 25 mm x 50 mm chip compartments.
  • Regular linescan mode produced spectra with higher reflectance and therefore better sensitivity, than chip scanning mode.
  • Rinsing chips with water to remove clinging fines increased TIR reflectance by up to 40% and reduced volume scattering effects.

The findings of this work will be used to help develop a consistent approach to scanning future RoXplorer® drillholes across the state and territory Geological Survey HyLoggerTM nodes.


Hillis R. R., Giles D., Van Der Wielen S. E., Baensch A., Cleverley J. S., Fabris A., Uvarova Y. (2014). Coiled tubing drilling and real time sensing – enabling prospecting drilling in the 21st Century? Society of Economic Geologists Special Publications, 18, 243–259.

Tiddy C. J., Hill S. M., Giles D., van der Hoek, B. G., Normington V. J., Anand R. R., Baudet E., Custance K., Hill R., Johnson A., McLennan S., Mitchell C., Zivak D., Salama W., Stoate K. Wolff K. (2019). Utilising geochemical data for the identification and characterisation of mineral exploration sample media within cover sequence materials, Australian Journal of Earth Sciences, DOI: 10.1080/08120099.2019.1673484


Jake completed undergraduate studies at Curtin University and honours at the University of Tasmania before spending 4 years working for Mineral Resources Tasmania with the HyLoggerTM team. He is currently employed as a spectral geologist with the Geological Survey of NSW to oversee the delivery of the NVCL project.

Petrographically constrained in situ sulfur isotopes: why the “SEDEX” can’t be used as model for sediment-hosted sulfide deposits in the 1.6 Ga Edmund Basin, Australia

Lampinen, Dr Heta1, LaFlamme, Dr Crystal2, Occhipinti, Dr Sandra1, Fiorentini, Dr Marco3, Spinks, Dr Sam1

1CSIRO, Kensigton, Australia, 2Université Laval, , Canada, 3Centre for Exploration Targeting, School of Earth Sciences, University of Western Australia, Crawley, Australia

A common foundation for sediment-hosted massive sulfide (SHMS) deposit sulfur isotope data interpretation is the assumption of sedimentary exhalative “SEDEX” model. The model presumes synsedimentary sulfide precipitation and the sulfur mainly sourced from the contemporaneous ocean via bacterial sulfate reduction, which can be further interpreted to reflect the evolution of ancient hydrosphere. However, synsedimentary SEDEX model has been challenged or disproven for many SHMS deposits, including ones in the McArthur Basin, Australia. Many SHMS deposits also contain multiple coexisting sulfide generations and/or express geospatial associations between the isotope signature and distance from the hydrothermal vent. Due to the internal complexity of SHMS systems, unravelling their sulfur isotope architecture requires both a robust paragenetic framework and a well-known geological context for the data. In situ secondary ion mass spectrometry (SIMS) sulfur isotope analysis has this capability.

Petrographically constrained in situ sulfur isotope SIMS analysis was applied to pyrite and chalcopyrite (n=135) to investigate the spatial and temporal sulfur isotope architecture of replacement and synsedimentary-style SHMS deposits at four sites (including the Abra deposit) in the ca. 1680-1455 Ma Edmund Basin, Western Australia. From this data, the sulfur isotope fractionation associated with the hydrothermal mineral systems, and representativeness for the secular evolution interpretations of the seawater sulfate through the Proterozoic Eon was evaluated.

The epigenetic replacement-style SHMS systems in the Edmund Basin yield δ34S from +24 to +54‰ from pyrite and chalcopyrite. The relatively 34S depleted pyrite were associated with ore fluid composition in main hydrothermal channels. The bulk isotopic composition of the ore fluid can be used as proxy for sulfate in the underlying sediments. The extremely 34S enriched were found in pyrite in distal parts of the deposit hydrothermal footprint. This 34S enrichment was possibly caused by deficiency of iron relative to sulfur in low permeability rocks, which decelerates the formation of pyrite allows the mass-dependent Rayleigh distillation of sulfur isotopes to reach extreme residual fraction. The systems with syn-sedimentary sulfide precipitation yield δ34S from +1 to +22‰, which can be associated with seawater sulfate and bacterial activity in the basin.

In situ sulfur isotope analysis offered the capacity to link isotopic data to a comprehensive spatially and temporally constrained framework representative of the stratigraphic and geodynamic context. The results of this study also highlight the importance of using tailored geological constraints and a mineral system model as a framework for isotope chemistry – not a generic SEDEX. Tailored geological constraints and deposit model are particularly important for the data are intended for evaluation of hydrosphere over time.

1 CSIRO Mineral Resources, 26 Dick Perry Avenue, Kensington, WA 6151, Australia

2 Département de géologie et de génie géologique, Université Laval, Pavillon Adrien-Pouliot 1065, av. de la Médecine, Québec, QC G1V 0A6, Canada.

3 Centre for Exploration Targeting, ARC Centre of Excellence for Core to Crust Fluid Systems (CCFS), School of Earth Sciences, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia.

Corresponding author:



Heta hails from Kuortane, Finland and has a MSc from University of Turku and a PhD from University of Western Australia. Her research focuses on delineation of multi-scale hydrothermal mineral footprints of undercover ore deposits using integrated geological, hyperspectral, geochemical and geophysical data.

Regional scale folding in the Arthur Metamorphic Complex: structural constraints for the Keith River – Lyons River area, NW Tasmania

Cumming, G1, Jackman C.1, Everard, J. L.1 and Gray, D.2

 1 Geological Survey Branch – Mineral Resources Tasmania – Geological Survey Branch, Rosny Park, Australia, 2 Consultant Structural Geologist for Mineral Resources Tasmania – Geological Survey Branch, Rosny Park, Australia

New geological mapping in the Keith River-Lyons River area in NW Tasmania has provided insight into the structural framework of the northern Arthur Metamorphic Complex (AMC) in NW Tasmania. The AMC is flanked to the east by the Oonah Formation and to the west by the Rocky Cape Group. The main lithological units of the AMC share transitional metamorphic, interpreted low angle fault, and both conformable and unconformable contacts. At a regional scale a significant north-plunging synform, or a large north-tilted block, is contained within the high-strain core of the AMC. Five deformation episodes can be observed throughout the area at outcrop-scale. Three early deformational episodes are likely related to the Middle Cambrian Tyennan orogeny, manifested as early, high strain events which caused isoclinal folding and development of schistose axial planar fabric. A rotational shear component, apparent as shear bands, suggests north over south or sinistral transport. Subsequent D3 deformation within the AMC occurred during the later stage of the Tyennan Orogeny. This event folded and tightened the various stacked lithostratigraphic units to form non-cylindrical asymmetric folds. These were later subject to generally northwest-directed, potentially Devonian compression (D4) and tilting. At a regional scale, late-stage north-plunging folds are inferred along the highest strain zone of the northern AMC. This area is a locus for Mesozoic or early Cainozoic faults, and a half-graben also extends along this zone, which is also in the core of the AMC. A late stage (D5) folding event may be partly related to Devonian compression, although the timing and nature of this folding event is largely unclear.


Grace currently works as a geologist for the Geological Survey Branch at Mineral Resources Tasmania and has spent the last 8 years undertaking mapping work to compile numerous 1:25,000 geological maps of North West Tasmania. 


About the GSA

The Geological Society of Australia was established as a non-profit organisation in 1952 to promote, advance and support Earth sciences in Australia.

As a broadly based professional society that aims to represent all Earth Science disciplines, the GSA attracts a wide diversity of members working in a similarly broad range of industries.