Next Article in Journal
Editorial for the Special Issue “Platinum-Group Minerals: New Results and Advances in PGE Mineralogy in Various Ni-Cu-Cr-PGE Ore Systems”
Next Article in Special Issue
Copper-Arsenic Nanoparticles in Hematite: Fingerprinting Fluid-Mineral Interaction
Previous Article in Journal
Cryogenic Raman Spectroscopic Studies on Common Ore-forming Fluid Systems
Previous Article in Special Issue
Crystals from the Powellite-Scheelite Series at the Nanoscale: A Case Study from the Zhibula Cu Skarn, Gangdese Belt, Tibet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Zircon at the Nanoscale Records Metasomatic Processes Leading to Large Magmatic–Hydrothermal Ore Systems

1
School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
2
Adelaide Microscopy, The University of Adelaide, Adelaide, SA 5005, Australia
3
BHP Olympic Dam, Adelaide, SA 5000, Australia
4
Key Laboratory of Mineral Resources, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Minerals 2019, 9(6), 364; https://doi.org/10.3390/min9060364
Submission received: 30 April 2019 / Revised: 1 June 2019 / Accepted: 13 June 2019 / Published: 16 June 2019
(This article belongs to the Special Issue Minerals Down to the Nanoscale: A Glimpse at Ore-Forming Processes)

Abstract

:
The petrography and geochemistry of zircon offers an exciting opportunity to better understand the genesis of, as well as identify pathfinders for, large magmatic–hydrothermal ore systems. Electron probe microanalysis, laser ablation inductively coupled plasma mass spectrometry, high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) imaging, and energy-dispersive X-ray spectrometry STEM mapping/spot analysis were combined to characterize Proterozoic granitic zircon in the eastern Gawler Craton, South Australia. Granites from the ~1.85 Ga Donington Suite and ~1.6 Ga Hiltaba Suite were selected from locations that are either mineralized or not, with the same style of iron-oxide copper gold (IOCG) mineralization. Although Donington Suite granites are host to mineralization in several prospects, only Hiltaba Suite granites are considered “fertile” in that their emplacement at ~1.6 Ga is associated with generation of one of the best metal-endowed IOCG provinces on Earth. Crystal oscillatory zoning with respect to non-formula elements, notably Fe and Cl, are textural and chemical features preserved in zircon, with no evidence for U or Pb accumulation relating to amorphization effects. Bands with Fe and Ca show mottling with respect to chloro–hydroxy–zircon nanoprecipitates. Lattice defects occur along fractures crosscutting such nanoprecipitates indicating fluid infiltration post-mottling. Lattice stretching and screw dislocations leading to expansion of the zircon structure are the only nanoscale structures attributable to self-induced irradiation damage. These features increase in abundance in zircons from granites hosting IOCG mineralization, including from the world-class Olympic Dam Cu–U–Au–Ag deposit. The nano- to micron-scale features documented reflect interaction between magmatic zircon and corrosive Fe–Cl-bearing fluids in an initial metasomatic event that follows magmatic crystallization and immediately precedes deposition of IOCG mineralization. Quantification of α-decay damage that could relate zircon alteration to the first percolation point in zircon gives ~100 Ma, a time interval that cannot be reconciled with the 2–4 Ma period between magmatic crystallization and onset of hydrothermal fluid flow. Crystal oscillatory zoning and nanoprecipitate mottling in zircon intensify with proximity to mineralization and represent a potential pathfinder to locate fertile granites associated with Cu–Au mineralization.

1. Introduction

Ore deposits falling within the broad iron-oxide-Cu–Au (IOCG) clan are a fascinating, and in some cases, poorly understood, family of magmatic–hydrothermal ore systems (References [1,2] and references therein). Debate on their genesis includes the relative timing and evolution of fluids during transition from the magmatic to hydrothermal stage. Detailed investigation of magmatic accessory minerals at appropriate scales of observation has enabled fluid–mineral interaction to be fingerprinted (e.g., using apatite [3,4]).
Zircon (ZrSiO4) is a common accessory mineral and has been pivotal for the evolution of U–(Th)–Pb geochronology (References [5] and references therein). Aside from actinides, magmatic zircon will incorporate a range of trace elements (lanthanides, Y, Hf, and Ti) which are valuable tracers for mantle versus crustal origin of melts and for underpinning plate tectonic models during Earth cycles (e.g., [6,7,8]). Contemporary consensus views the incorporation of other, “non-formula” components in zircon, including Ca, Fe, Al, and H2O, as being unrelated to primary magmatic processes (Reference [9], and references therein). These are most often measured within altered zircon (e.g., [10]) and are attributed to metamictization-related effects.
In this contribution, we bridge micron- to nanoscale observations to comprehensively characterize magmatic zircon from Proterozoic granites in the eastern Gawler Craton, South Australia (Figure 1A). Iron-oxide copper gold systems of various sizes are found within granites of different ages in the Olympic Cu–Au Province [11], for example, the Olympic Dam deposit [12] and Wirrda Well prospect, hosted by ~1.6 Ga (Hiltaba Suite) and ~1.85 Ga (Donington Suite) granitoids, respectively. The U–Pb age of hydrothermal hematite from Olympic Dam, a product of Fe-metasomatism, has been constrained by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) to 1589.91 ± 0.94 Ma [13], crystallizing roughly 2–4 Ma after zircon within the granite hosting the deposit (chemical abrasion ID-TIMS; 1593.87 ± 0.21 Ma) [14].
The overarching aims of the present work were two-fold. Firstly, we aimed to document metasomatic, fluid-assisted alteration of zircon that occurs in a magmatic–hydrothermal environment prior to structural damage induced by alpha-decay. Secondly, we tested this hypothesis and its implications using samples from metallogenetically productive and unproductive granites of the same suites across the Gawler Craton. We set out to show that micron- to nanoscale observations of zircon geochemistry and crystal structure can, in combination, be used to trace early mineralizing fluids. Furthermore, we provide atomic-scale resolution images of zircon, offering insights into metasomatic processes and products, such as the formation of zircon nanoprecipitates and structural defects.
Such an approach provides potential qualitative indicators of magma fertility (the ability of a magma to generate hydrothermal fluids that are sufficiently well-endowed with Cu, etc., to form a sizable ore deposit). Such an idea is validated by geochemical analysis of known productive and unproductive granites, with a broad application in IOCG exploration across the Gawler Craton and elsewhere.

2. Background and Rationale

2.1. Zircon Chemistry and Structure

Magmatic zircon has a crystal structure resistant to physical and chemical degradation during secondary processes, properties which have motivated research into its usage as a durable host for nuclear waste storage [15]. However, zircon undergoes self-irradiation damage during U–(Th)–Pb decay, ultimately transitioning from crystalline to an amorphous (metamict) state [10], which can be followed by partial to complete recrystallization. The amorphous domains begin to interconnect due to the radiation damage (α-decay), over a time period determined by U/Th concentrations and annealing rates [16]. The time dependency of structural damage is highly relevant for geochronology, as it can result in discordance of the U–Pb system, most commonly through Pb-loss. Disturbances to zircon U–Pb systematics via metamictization can, however, be selectively eliminated by chemical abrasion of damaged zones, prior to ID-TIMS, permitting high-precision geochronological constraints [17].
The secondary textures commonly observed to crosscut growth zones in magmatic zircon are attributed to re-equilibration of zircon in aqueous fluids and melts (Reference [18] and references therein). Primary oscillatory zoning in zircon can, however, be locally preserved, even in cases when the superimposed reaction temperature exceeds ~600 °C [19]. Structural damage can also create pathways for fluids to enter and precipitate minerals as inclusions within zircon displaying secondary microstructures. For example, Fe–Ti oxides and ZrO2 observed along dislocation and pores in Hadean zircons (Jack Hills, Western Australia) are considered to postdate primary crystallization [20].
Elemental abundance and isotope ratios in zircon, including rare earth element (REE) fractionation patterns and Eu/Eu*, Sr/Y, and V/Sc anomalies, all readily determined by microbeam analysis, can be used to define magma fertility and track its evolution (Reference [21] and references therein). These geochemical pathfinders have been extensively applied to target metallogenically productive (“fertile”) and non-productive (“infertile”) granitic intrusions [22]. Despite the potential for identifying geochemical changes in zircon as a response to metasomatism, the commonly used microanalytical instrumentation, such as laser-ablation inductively-coupled plasma mass-spectrometry (LA-ICP-MS), cannot accurately spatially resolve compositions smaller than a few microns at best, a high threshold relative to typical internal textures and structures commonly found in zircon.
Z-contrast techniques such as high-angle annular dark-field scanning transmission electron microscopy (HAADF STEM; e.g., [23]) allows processes of metamictization and alteration to be addressed down to the nanoscale. For example, fission tracks relating to radiation damage, the nanoscale occurrence of Pb, and evidence for nanoscale U-mobility have been revealed by this technique [24,25,26,27]. Increasing levels of spatial resolution and analytical precision combined with in situ micro-sampling techniques have permitted HAADF STEM study of foils prepared by focused ion beam (FIB)-scanning electron microscopy (SEM) to visualize atomic structure and defects in minerals (e.g., [28,29]).

2.2. The Olympic Cu–Au Province

The Olympic Cu–Au Province, South Australia (Figure 1A) hosts IOCG-style mineralization attributed to hydrothermal activity associated with emplacement of the Gawler Silicic Large Igneous Province at ~1.6 Ga [11], accompanying the onset of Columbia/Nuna supercontinent breakup [30]. Magmatic activity is represented by the Hiltaba Suite and contemporaneous Gawler Range Volcanics [31]. The ~700 km long Olympic Cu–Au Province strikes roughly N–S along the eastern margin of the Gawler Craton from the Mount Woods Inlier in the north, which hosts the Prominent Hill deposit, through Olympic Dam and Carrapateena, to the Moonta and Hillside deposits in the south (Figure 1A). Each major deposit is surrounded by numerous smaller, less explored prospects. Host lithologies range from Gawler Range Volcanics (e.g., Acropolis prospect; Figure 1B), granites of Hiltaba Suite affiliation (e.g., Olympic Dam within Roxby Downs Granite; RDG), or older granitoids (e.g., ~1.85 Ga Donington Suite granites at Wirrda Well and Carrapateena; Figure 1A,B).
Several other styles of mineralization occur on the Eyre Peninsula to the immediate west of the Olympic Province. These include older banded iron formations [32] and younger U-mineralization within cover rocks suprajacent to granites belonging to the Samphire Pluton of Hiltaba Suite affiliation [33]. A further Hiltaba Suite granitoid is the Charleston Granite (Figure 1A). No IOCG-style mineralization has yet been discovered associated with either the Charleston Granite or the Samphire Pluton. Donington Suite granitoids are present throughout the Eyre Peninsula and outcrop as far south as Cape Donington, south of Port Lincoln. Although the potential for an extension of the Olympic Cu–Au Province westwards into the Eyre Peninsula remains largely untested, we note that several of the aforementioned banded iron formation ores feature geochemical signatures suggestive of overprinting by hydrothermal fluids of granitic affiliation [32,34,35].
Formation of IOCG deposits remains highly contested, particularly with respect to sources and timing of fluids and contained metals (References [1,2], and references therein). In the Gawler Craton, lithospheric architecture, mantle metasomatism, and cratonic metallogeny [37,38] have been linked to generic models for IOCG genesis with global applicability [1]. The large volumes of Fe required to form IOCG deposits associated with granitoid-derived fluids or felsic magmatism have recently been proposed to be derived from underplating mafic magma in deeper reservoirs allowing for volatile- and metal-transfer via magma mingling and magnetite flotation [39]. If Fe is sourced from a deep magma reservoir, then zircon, as a ubiquitous magmatic accessory mineral, should record this, either in the magmatic growth stage or as an overprint following interaction with early, low pH hydrothermal fluids released from the locally Fe-enriched magma. Zircon from fertile and barren granitoids should, therefore, differ from one another, both in terms of their geochemical signatures and textures at the micron- to nanoscale.

3. Sample Suite

For this study, we used magmatic zircon from five granitoids of Donington- (~1.85 Ga) and Hiltaba-affiliation (~1.6 Ga) from within and outside the Olympic Cu–Au Province (Figure 1C; Table 1). Except for the example from Cape Donington, all granites considered here displayed macroscopic alteration, albeit to varying degrees, while still retaining granitoid textures and some relict magmatic minerals. The sample suite is summarized in Table 1 with respect to rock type, location, key features, and types of analysis reported in the present contribution.
Zircon from all five case studies were temporally well-constrained by U–Pb methods, either previously published, or given in this study (Figure 1C). We included zircon from RDG samples, collected at locations distal and proximal to the Olympic Dam deposit (Figure 1B), to evaluate whether metasomatic processes analogous to those reported for feldspars [40] and apatite [3] could be tracked. Secondly, we sought to identify whether such processes are recognizable in older zircon from Donington Suite granite hosting the Wirrda Well prospect [41], 25 km SSE of Olympic Dam. Zircon from the Charleston and Cape Donington granites are included to assess the potential of zircon as an indicator of the metallogenic fertility of plutons in the Eyre Peninsula bordering the Olympic Cu–Au Province, where IOCG-style mineralization has not thus far been discovered.

4. Methodology

Four granite samples (3 one-inch polished blocks and one thin-section) and one polished block containing separated zircon (sample PL8) were studied in backscattered electron (BSE) and secondary electron (SE) modes using an FEI Quanta 450 SEM. Cathodoluminescence imaging was undertaken using an FEI Quanta 600 SEM equipped with a tungsten filament electron source. The purpose was to identify zircon grains with representative characteristics, particularly primary growth with oscillatory zoning textures and overprinting thereof. Such grains were analyzed by electron probe microanalysis (EPMA; Cameca SX-Five) and LA-ICP-MS (RESOlution-LR excimer laser microprobe coupled to an Agilent 7900x Quadrupole ICP-MS) for minor to trace element and U–Pb compositions using both spot analysis and mapping. All instruments are housed at Adelaide Microscopy, University of Adelaide. Extended details of instrumentation and analytical parameters are provided in Supplementary Materials A. All geochemical data are tabulated in Supplementary Materials B.
For subsequent S/TEM study, seven foils (Table 1; Figure 2A–E, Figure 3B,C, and Supplementary Materials C, Figure S1) were extracted in-situ, thinned to <100 nm by ion beam (Ga+) milling and attached to Cu grids on an FEI Helios Nanolab 600 FIB-SEM platform. Of the seven foils, five were cut across oscillatory zonation observed on grain surfaces in zircon to target their development at depth, allowing correlation between micron- and nanoscale characteristics. In the case of the Wirrda Well Donington sample (LCD17), foil preparation targeted areas that yielded concordant (Foil #5) and discordant (Foil #4) U/Pb data points. In the Olympic Dam (OD) distal sample (LCD47), two zircon grains were cut, one with primary oscillatory zoning (Foil #1) and a second grain (Foil #2) displaying secondary, highly altered domains lacking oscillatory zoning.
Nanoscale study, HAADF STEM imaging, and EDS spot analysis/mapping were performed using an ultra-high-resolution, probe-corrected, FEI Titan Themis S/TEM operated at 200 kV. This instrument is equipped with the X-FEG Schottky source and Super-X EDS geometry. The Super-X EDS detector provides geometrically symmetric EDS detection with an effective solid angle of 0.8 Sr. Probe correction delivered sub-Ångstrom spatial resolution and an inner collection angle greater than 50 mrad was used for HAADF experiments using the Fischione HAADF detector.
Diffraction indexing was performed using Winwulff© 1.5.2 software and publicly available data from the American Mineralogist Crystal Structure Database. Crystal structure models were obtained using CrystalMaker® version 10.1.1 and STEM for xHREMTM version 4.1 software.

5. Results

5.1. Zircon Petrography, U–Pb Dating, and Selection of Grains for Nanoscale Study

Screening of zircon populations in each sample shows a range of grain sizes from 30 to ~300 μm, displaying both primary and secondary overprinting textures (Figure 2A–D). Zircons in two samples of unknown age but with assumed affiliation were dated by LA-ICP-MS to confirm their Hiltaba and Donington Suite ages.
Uranium–Pb dating of Wirrda Well zircon (LCD17) produced an upper intercept age of 1848 ± 10 Ma with concordant and highly discordant data points (Figure 3A), corresponding to published Donington Suite ages [43,44]. Dating of the OD proximal zircon (LCD13) gives an upper intercept age of 1572 ± 37 Ma, with analyses displaying both a high degree of concordance and discordancy. Full U–Pb data for analyzed zircons and reference material is tabulated in Supplementary Materials B, Table S3.
Selection of grains for nanoscale study was based upon several criteria: (i) well-constrained U–Pb ages in which case the FIB cut was placed next to a LA-ICP-MS crater (Figure 3); (ii) euhedral morphology with elongation along the c-axis; and (iii) large, least-fractured, inclusion-free grains. Subhedral zircon with marginal dissolution was also included to assess the most intense alteration observed within the sample suite (e.g., OD distal zircon; Figure 2B,D). In all five samples, zircon displays crystal oscillatory zoning expressed as bright and dark bands on BSE images (Figure 2 and Figure 3 and Supplementary Materials C, Figure S1A). The intensity contrast of the banding ranged from subtle in zircon in granites from the Eyre Peninsula (PL8 and Charleston; Figure 2E and Supplementary Materials C, Figure S1A) to strong in granites from the Olympic Cu–Au Province. Furthermore, the dark bands could display mottled textures with respect to what appears to be inclusions, varying from sub-micron (in the banding) to ~1–2 µm (in the domains obliterating the zoning) (Figure 2C,D).
Spot EDS-SEM analysis does not show any compositional difference between the “inclusions” and host band/domain suggesting these could instead be pores. However, secondary electron imaging only rarely shows the presence of pores, generally in the scalloped areas with coarser mottled textures. Otherwise, the nature of the sub-micron “inclusions” cannot be resolved with SEM imaging alone.
Except for the Cape Donington sample (PL8; Figure 2E), all grains were fractured to varying degrees and were associated with marginal corrosion and occurrence of darker still domains with irregular and/or scalloped morphologies superimposed on crystal zoning (Figure 2D). Radial fractures, typically interpreted as accompanying metamictization, were present in such cases, notably in the OD distal zircon (Figure 2B,D). Additionally, veinlets and overgrowths of xenotime (Figure 3B), as well as tiny U-bearing phases, were also present in the RDG samples.

5.2. Trace and Minor Element Concentrations

Electron microprobe datasets for analyzed zircon (Table 2 and Supplementary Materials B, Table S4) showed steady, consistent concentrations of HfO2 (1.5–2.0 wt.%), Y2O3, (HREE)2O3 and P2O5, minor UO2, PbO, and ThO2 in most but not all analyses, and sporadic concentrations of Na2O, K2O, CaO, MgO, MnO, and Al2O3 above the minimum detection limit (mdl) values across the sample suite. TiO2 was measured consistently (up to 0.06 wt.%) in the Charleston Granite sample only. Sc2O3 was below mdl in most analyses.
“Bright” and “dark” bands in zircon from Cape Donington zircon (PL8) were indistinguishable in terms of composition but the latter had lower analytical totals before recalculation. “Dark” bands in zircon from the Wirrda Well Donington Suite and Charleston Granite samples (LCD17 and MB158) also had totals ~1 wt.% lower than the “bright” bands, generally higher Fe2O3 and minor but measurable Cl. Analogous trends were seen in samples the two OD samples (LCD47 and LCD13), in which the “dark” bands contained 1.0–1.9 wt.% Fe2O3 and ~0.1 wt.% Cl. Analytical totals were 1–2 wt.% lower than for “bright” bands. Additionally, fluorine was measurable (up to 0.15 wt.%) in a minority of analytical spots from the “dark” bands. Variability of concentrations within each sample, and between “bright” and “dark” bands in the same zircon grain, was a feature of the dataset; the degree of intra-grain variability increased with the total of non-formula elements. Importantly, measured concentrations of Y2O3, (HREE)2O3, P2O5, UO2, and ThO2 were statistically identical in “bright” and “dark” bands across the dataset.
Laser ablation inductively coupled plasma mass spectrometry datasets (Supplementary Materials B, Table S5) corroborated and confirmed the accuracy and inherent variability of the EPMA data. Zircon from the Charleston granite was not considered due to the small grain sizes. Other elements (e.g., Sn, Nb, Ta) were not present at concentrations exceeding (at most) a few tens of ppm. Even using a small diameter spot (19 μm), the LA-ICP-MS method did not permit selective analysis of “bright” and “dark” bands.
Cape Donington zircon showed oscillatory patterns correlating with U, Pb, REE, and Y in LA-ICP-MS maps (Figure 2F; Supplementary Materials C, Figure S2). Chondrite-normalized fractionation trends constructed from LA-ICP-MS data (Figure 4A) showed reproducible patterns for zircon in each granite, although with systematic differences in the LREE portion between Donington and Hiltaba Suite samples.
Principal component analysis (PCA) applied to the centered logratio-transformed LA-ICP-MS dataset provided additional insights. A projection of PC1 versus PC2 (Figure 4B,C) showed two distinct groupings: (1) Y, REE (except Eu and Pr) and P; and (2) Hf, U, Pb, Th, Ti, Nb, Eu, and Pr. On both plots, these groupings are distinct from a third group comprising Fe, Mn, Al, and Ca. These patterns underscore interpretation of the EPMA data in which the geochemical signature associated with the metasomatic overprint was defined by Fe, Mn, Al, and Ca, and is independent of any magmatic signature or pattern of distribution at the grain-scale involving element groups (1) and (2). A ternary Fe–Al–Hf diagram (Figure 4D) allows a visual expression of chemical change from least-altered zircon (sample PL8), effectively without Fe or Al, through moderately altered zircon in sample LCD17, to the sub-populations of zircon in samples LCD13 and LCD17 from Olympic Dam, which contained the highest amounts of substituting non-formula (NF) elements.

5.3. Assessment of Zircon Stoichiometry and Substitution Model

Recalculation of EPMA data was undertaken following a template used for hydrogarnet and underlying assumptions/substitutions as provided by Supplementary Materials A. A key objective was an estimation of H2O content based on allocation of all minor elements (Fe, Al, Mn, and Ca) to structural sites, bringing analytical totals close to 100 wt.% (from ~99.2). The data (Supplementary Materials B, Table S4) showed a correlation between atom per formula unit (apfu) values for (REY + Sc)3+ and P5+ on the plot of these two variables (Figure 4E), albeit with a deficiency in the latter—a feature widely described in several publications (see Supplementary Materials A) and taken to indicate more complex substitution mechanisms. Electron microprobe datasets implicitly carry a margin of analytical error, which in the case of the zircon analyzed here were compounded by concentrations of many elements around or below mdl, and by assumptions in the formula calculation that were in turn based on equivocal evidence given in the literature (see Supplementary Materials A). As such, derived stoichiometries and calculated OH contents represent only a best-possible estimation based on the data available. Nonetheless, using mean compositions and a simplified formulation in which Zr* + REE* + Fe + NF = 1, where Zr * = Zr + Hf + Th + U; REE* = Sc + Y + (La to Lu); NF = Mg + Ca + Mn + Nb + Pb + Al, B = Si + Al + P + Ti = 1, and I = Na + K (interstitial), mean compositions (Table 3) can be defined for zircon from the bright and dark bands within each sample.
The lack of Fe or other non-formula elements and Cl is of note in the Cape Donington zircon, implying these elements were not of magmatic origin. The relatively high OH content in these otherwise less altered zircons may be an artefact of the calculation rather than real and lies within analytical precision of the EPMA method. Critical for the narrative of this contribution, zircon from all three altered granites from the Olympic Cu–Au Province (two of Hiltaba age, one of Donington age) display concomitant enrichment in Fe and other NF elements and in (OH, Cl, F) in the darker bands and over those in unaltered zircon (PL8), in full agreement with the nanoscale observations. Concentrations of other minor (magmatic) elements remain unchanged (within analytical error).
A plot of the main “non-formula” components, and Fe versus OH + Cl + F (Figure 4F) showed a strong correlation for all altered zircon, except Charleston. Mean compositions and empirical formulae for each sub-population are given in Supplementary Materials B, Table S4.

5.4. Distribution of Minor Elements: Micron- to Nanoscale Patterns

The presence of Fe and Cl in zircon was also assessed by EPMA mapping of zircon in two samples: Wirrda Well (LCD17) and OD-proximal (LCD 13). The maps show excellent positive correlation between the two elements in the crystal zoning of Wirrda Well zircon (Figure 5A,B) but antipathetic patterns for OD proximal zircon (Figure 5C,D). The OD proximal zircon displays a combination of primary and overprinting textures in which the oscillatory zoning is only preserved on the rims, whereas the middle part of the grain is affected by fractures and re-crystallization obliterating primary zoning (Figure 3B and Figure 5C,D).
Imaging in HAADF STEM mode depicts oscillatory zoning as bright and dark bands relating to non-formula elements in all foils shown in this study (Table 1), except the Cape Donington zircon (Foil #6). Superimposed mottled textures in the darker bands were only observed in zircon from the Olympic Cu–Au Province. Mapping using STEM EDS showed Fe banding in all zircons except Cape Donington. In the latter, the oscillatory zoning was related to trace elements typical of magmatic zircon (REY, actinides), the concentrations of which were below the detection limits of the EDS STEM method but readily mappable by LA-ICP-MS (Figure 2F; Supplementary Materials C, Figure S2). In contrast, the Wirrda Well–Donington zircon displayed an excellent correlation between micron and nanoscale oscillatory zoning patterns with respect to Fe, as shown by both EPMA and STEM EDS mapping (Figure 5A,B and Figure 6). The Fe distribution, however, did not reveal the mottled texture in the darker Fe-rich bands or veinlets crosscutting the brighter bands, as observed in the HAADF STEM images (Figure 6C).
The most varied chemistry mapped at the nanoscale was observed in zircon from the Charleston Granite (Figure 7A). Banding with respect to Fe, Ca, Ti, Al, Y, and Th was recorded outside of a magmatic domain, itself defined by faint oscillatory banding expressed by variation in major element concentration as shown by O, Si, and Zr maps. No mottled textures were present throughout such chemically complex bands. The EDS STEM maps also showed that Fe displayed a sharp boundary against the magmatic domain, contrasting with the diffuse boundaries shown by the other trace elements. All these elements, including U (but not Pb), were enriched within a fracture around the magmatic domain, and within a fine particle inclusion hosted at the diffuse boundary (Figure 7A). The fracture, leading towards this U-bearing fine particle, was represented by a rectangular network of dark nanodomains that hosted pores (Figure 7B–D). These domains were nonetheless crystalline, albeit disordered and with slightly different orientation relative to one another and to host zircon.
The best example of mottled textures throughout oscillatory banding with respect to Fe (Ca, Cl) was shown by the least altered grain in OD distal zircon (Foil #1; Figure 8A). The HAADF STEM image of the foil mimics the BSE image (Figure 2A,C) with respect to the mottled textures throughout the banding. Relative to the generally homogenous distribution of Fe and Ca throughout a given band, the Cl map shows distinct spots of higher concentration. This was exemplified by the signals for Ca and Cl along a profile (Figure 8B). Likewise, variable correlations were seen between the high Cl areas and concentration of major elements, O and Zr, as shown on the EDS STEM maps obtained from mottled areas (Figure 8C). The highly altered zircon grain (Foil #2) showed two distinct domains: (i) weakly-banded and mottled and (ii) cleaner domains that were strongly enriched in Ca, Fe, and Al. The latter surrounds veinlets of xenotime displaying enrichment in HREE and U (Supplementary Materials C, Figure S3).
Olympic Dam proximal zircon (Foil #3) displayed a combination of primary and overprinting textures correlating with minor trace element patterns (Figure 9). In depth, the foil revealed a subsurface 5 μm-diameter inclusion of hematite (Figure 9B–D), as well as oscillatory banding with mottled textures of variable density. As in the OD distal zircon, Fe was evenly distributed throughout the bands, whereas Cl appeared more heterogeneous (Figure 9E,F). In detail, the Fe-poor bands were crosscut by Fe-rich veinlets (Figure 9F). The largest areas of Cl-enrichment were observed in parts of the Fe-rich bands but typically did not correspond to the crosscutting Fe-rich trails. Such patterns of Fe and Cl explain the inverse trends mapped at the micron-scale in the same grain (Figure 5C,D).

5.5. Mottled Areas: Nanoparticles to Fine Particles of Chloro–Hydroxy–Zircon

The nature of the mottled areas (identified in foils as listed in Table 1) relative to Cl-enrichment and substitution mechanisms for zircon was addressed by EDS spot analysis and mapping at various resolutions (Figure 10 and Figure 11). Higher resolution maps of one Cl-rich area selected from Figure 8C reproduced the same variability between major elements and Cl (and remarkably little variation in terms of minor elements) in the particles observed as dark nanodomains on HAADF STEM image (Figure 10A). An inverse correlation between O and Zr corresponded to either diffuse or stronger Cl concentration, whereas the strongest Cl concentrations were not reflected by major element variation.
The highest Cl concentration was recorded from the smallest particles (10–30 nm-diameter) of grey color, (Foil #1; Figure 10B). A profile across one of the darker particles shows that Cl was concentrated at the edge of the particle (Figure 10C). The HAADF signal intensity and O decreased across the boundary of the particle, whereas major and minor elements displayed a slight increase. These results can also be associated with variable particle thickness, their position relative to the plane of imaging, or the presence of voids/pores, particularly for the coarser and darkest domains. Comparative imaging of the same area in both HAADF and bright-field STEM modes illustrates these features (Supplementary Materials C, Figure S1). Such voids could, however, also have been produced via the plucking out of coarser particles during ion milling.
The morphology and size of the particles throughout samples with mottled textures (Foils #1–5) varies from several nm up to 100 nm (nanoparticles) to hundreds of nm (fine particles) (Figure 11A). The largest particles, with well-defined geometrical shapes, can include pores, whereas cross-cutting darker veins will show numerous defects when imaged at high resolution (see below). Energy-dispersive X-ray maps of one medium-sized particle showed strong (lower and higher) variation in O and Zr concentrations relative to host zircon; no variation in Si was recorded (Foil #1; Figure 11B). Weak enrichment in Cl and U was present whereas Ca-rich nanoparticles occurred just outside, along the edge of the mapped particle. Overall, these results suggest that the main substitution involved replacement of O by OH and/or Cl leading to chloro–hydroxy–zircon formation. Although substitution of (SiO4)4− tetrahedra by 4(OH,Cl) can also be invoked, such a relationship is not obvious from the EDS mapping.

5.6. Atomic-Scale HAADF STEM Imaging—Zircon Crystallinity and Lattice-Scale Defects

Four main zone axes in zircon were imaged throughout the seven foils (Table 1; Figure 12A–D). In each case, imaging showed the same overall orientation across zircon in any given foil, indicating that the chemically heterogenous domains were an integral part of the same single crystal. Despite observing evidence of lattice disorder, zircon in all samples was found to be crystalline with no amorphous areas identified. This absence of amorphous domains was observed in both the brighter bands of magmatic zircon from Cape Donington, throughout the darker or mottled bands in all other samples, and in the proximity of fractures containing newly formed phases such as xenotime.
Atomic arrangements on HAADF STEM images for the four zone axes shown in Figure 12 were interpreted using crystal structure models obtained after indexing of the Fast Fourier Transform (FFT) obtained for each image and also STEM simulations (Figure 13). These show that the brightest spots represent an overlap between Zr and Si atoms on (001) and ( 1 1 ¯ 0 ) zone axes (Figure 13A,B), whereas both Zr (brightest, with dumbbell arrangement) and Si (darker) atoms can be resolved on the (100) and ( 1 1 ¯ 1 ) zone axes (Figure 13C,D). There was a relatively good fit between the HAADF STEM images and simulations in each case. Such image assessment of individual zone axes by HAADF STEM was necessary to identify the presence of lattice-scale defects.
Lattice-scale defects on ( 1 1 ¯ 1 ) zone axis in zircon were observed along a veinlet crosscutting a fine particle in Foil #1 (darker color on HAADF STEM images; Figure 14A–D). One of the defects had a hexagonal shape (Figure 14D) and showed clear modifications in both size and intensity of individual atoms relative to “normal” zircon on this zone axis (Figure 12D and Figure 13D). An intensity profile showed changes across this defect which were interpretable as substitutions of individual Zr and Si atoms along the (011) direction (Figure 14E). Signals for Zr decrease within the defect suggesting substitution by lighter cations (e.g., Ca, Al) in variable proportions. One of the Si sites in the defect showed a higher signal relatively to host zircon, whereas the other was barely visible (correlating with darkest areas on the image). Substitution within the (SiO4)4− tetrahedra can be invoked in which Me5+ (e.g., P5+) replaces Si4+, necessitating charge balance via substitution of (OH,Cl) for O2−.
Lattice-scale defects were imaged on all zone axes in zircon except ( 1 1 ¯ 0 ) (Figure 15). In magmatic zircon from Cape Donington (Foil #7), lattice disorder was found in the brighter U-richer bands (Supplementary Materials B, Table S4; Figure 2F). Square-shaped defects on the (100)zircon zone axis appear darker and showed arrays of single Zr (and Si—although these were less well resolved on the image) columns along the c-axis, which were doubled in the b direction leading to a smaller, squarish arrangement (outlined in Figure 15A).
Linear defects were observed along the edge of xenotime veinlets (e.g., Wirrda Well–Donington zircon; Figure 15B). Atom displacement across the two sides of the defect can be interpreted as resulting from a simple screw dislocation considering the “atom loop” observed on the defect plane. The “loop” is illustrated by the atoms placed between the two planes with different heights relative to one another across the defect (smaller size atoms highlighted in Figure 15B). Linear defects along kink faults were imaged on (001) zone axis in Charleston zircon (Figure 15C). In this case, the fault steps/widths changed irregularly, as shown by the presence or absence of “atom loops” on the defect plane (overlays on Figure 15C). Wider planar defects, tens of nm in width, also occurred in the Charleston zircon (Figure 15D). This image can also be interpreted as a simple screw dislocation (yellow lines in Figure 15D). Interstitial atoms were inferred within the defect area only. Their size and brightness changed relative to the Zr and Si atoms on the (001) zone axis in zircon outside the defect area. Such changes could reflect differences in atom heights between the two planes on either side of the defect, or more speculatively, could be attributed to the presence of foreign atoms occurring interstitially in the zircon lattice (mixed site Ti, Ca, Th, Y, etc.). The latter supposition can be correlated with the presence of mixed sites of Ti, Ca, Th, Y, etc., as is known from zirconolite (CaZrTi₂O₇).
The most complex defects were imaged in Fe-rich bands from OD distal zircon in the [ 1 1 ¯ 1 ] zone axis (Foil #1; Figure 15E). In this case, changes in the atomic arrangement along the (011) can be interpreted as displacements with anti-phase modulation along fault planes (yellow lines in Figure 15E. In contrast, the defect within the fine particle described above (Figure 14D and Figure 15F) can be interpreted as the result of local fluid–mineral interaction leading to an increase in OH/halogen content in the nanoparticles/fine particles.

6. Discussion

6.1. Zircon Metasomatism Down to the Nanoscale: Are Iron and Chlorine of Magmatic or Hydrothermal Origin?

It is difficult to unequivocally distinguish magmatic (primary) from secondary patterns in zircon displaying crystal oscillatory zoning with respect to Fe and Cl due to the lack of crosscutting textures and because such elements could also have been components in either granitic melts or hydrothermal fluids. For example, a magmatic origin for divalent Fe was suggested by Reference [45] based on observed substitutions. In the present study, we were able to rule out a magmatic origin by analysis of zircon from a granite belonging to the same ~1.85 Ga Donington Suite which showed the presence of Fe and Cl only in parts of the craton where they are associated with or host younger (~1.6 Ga) IOCG-type mineralization (Wirrda Well). Moreover, the concentrations of Fe and Cl in the Hiltaba Suite zircons increased with proximity to more intense mineralization. Oscillatory zoning patterns can be produced via metasomatism (ion exchange between a mineral and a percolating fluid moving though pores in a rock) in the absence of open fractures. Therefore, we stipulate that Fe–Cl zoning represents a first stage of zircon alteration via metasomatism that is quite distinct from late-stage fracture-infill or textures superimposed onto the zoning (xenotime veinlets, U-bearing nanoparticles; Figure 6, Figure 7 and Figure 8).
Metasomatism allows partial or complete pseudomorphic replacement of one mineral by another. In some cases, the new mineral may be the same species but differs from the parent in terms of minor/trace elements, as was shown for hematite from Olympic Dam [46]. Such processes are “metasomatic” in a broader sense and occur when replacement reactions couple dissolution with (re)precipitation rates (CDRR). This interpretation concurs with empirical evidence from other studies of Hiltaba Suite granites [47], showing local alteration of magmatic feldspars or accessories via CDRR. Variation among the geochemical patterns recorded by zircon (e.g., Ca and Cl, as well as Ti, Al, Th, and Y in Charleston zircon) can be attributed to local mineralogical changes in the granite during early alkali–calcic alteration, e.g., release of cations during replacement of magmatic plagioclase by sericite or biotite by chlorite [42,47].
Metasomatism was accompanied by development of nanoscale structures, such as formation of chloro–hydroxy–zircon nanoprecipitates in all samples from the Olympic Dam district, but significantly, not in those from the Eyre Peninsula (Table 1).
Defects were observed throughout the Fe-rich bands, as well as in magmatic zircon from Cape Donington, but none show amorphization. Linear defects associated with fracture opening can be attributed to structure swelling along fission tracks (Figure 15B). This allows trace element remobilization, in turn leading to formation of U-nanoparticles or xenotime like those shown in zircon in Archean granites from Jack Hills [24,25]. Likewise, lattice expansion via stretching or screw-dislocation defects (Figure 16) can be considered part of the “swelling” induced by radiation damage.
The substitution model for metasomatized zircon shows cation exchange in which loss of 0.4–2 wt.% Si can be compensated by incorporation of OH– and halogens without significant removal of primary magmatic trace elements (Supplementary Materials A). Zirconium can, however, be replaced by Fe and other non-formulae cations. Chemical changes, affecting both Zr and Si, observed in defects along fractures crosscutting fine particles (Figure 14) were concordant with chloro–hydroxy–zircon of variable chemistry, as assessed by EPMA data. However, heterogeneity was observed at the nanoscale, including localized concentrations of other elements, e.g., Ca or Cl-bearing nanoparticles (Figure 10 and Figure 11), thus explaining much of the scatter in the geochemical data obtained at the EPMA-scale.

6.2. Timing of Zircon Alteration: Is Metasomatism Pre- or Post-Metamictization?

The incorporation of Fe and associated non-formula elements relating to metasomatism as introduced above, has been extensively discussed in terms of zircon alteration linked to self-irradiation. The overarching idea is that zircon should not be able to incorporate such elements at wt.% concentration unless it accumulates radiation damage, a process resulting in swelling and amorphization, and thus, creating pathways for fluid percolation (Reference [18] and references therein).
Quantification of the radiation damage necessary to facilitate major element substitution (the “First Percolation Point”) was defined as a dose of ~2 × 1018 α-decays/g (hereafter, Dc), and representing ~30–40% amorphization as a network of interconnected aperiodic domains within the crystal structure [48,49]. The link between the first percolation point and a sudden increase in non-formula elements, particularly Ca, was shown for altered zircon exhibiting U–Pb discordancy [18,19]. A high-porosity-to-spongy appearance is commonly reported, leading to deficient EPMA totals and darkening on BSE images [50]. The latter authors also contend that H2O incorporation is promoted by pores. Nanoscale Fe- or Pb-bearing inclusions were shown in Archean zircon [24], while pores with Ca- and U-rich margins have been documented from a pegmatite zircon [26].
Using molecular dynamic simulations of alpha-recoil cascades produced via actinide decay, models of self-irradiation damage in zircon show domains of disordered material as small as a few nm [18]. However, imaging of zircon in this study (Figure 12, Figure 14, and Figure 15) did not unequivocally support this model, even for samples within the Dc range required for the crystalline-to-amorphous transition. This may relate to thermal annealing following metasomatism, e.g., associated with emplacement of the Gairdner dyke swarm at ~820 Ma, the most recent major tectono-thermal event recognized throughout the Gawler Craton [51].
Nevertheless, the Dc accumulation since ~820 Ma should have produced recognizable disorder at the nanoscale. The range of defects discussed here (Figure 14, Figure 15 and Figure 16) are the closest expression of lattice disorder, although not truly amorphization. The positive correlation reported in other studies between high non-formula elements and Dc (e.g., [19]) is not as well-defined in this dataset. Sigmoidal trendlines shown in previous studies are comparable to that shown here, although there is a large spread in the data here, with points displaying high Dc and low Fe, and vice versa (Figure 17A). Therefore, the correlation between increasing non-formula elements and Dc may not necessarily be the result of high-U zircon creating their own pathways for increased fluid interaction. Furthermore, high-U zircon from this study indicate the Dc required to reach the first percolation point would not have been met until ~100 Ma after zircon crystallization. This postdates the measured timing of hematite crystallization, itself a result of Fe-metasomatism [13]. The presence of hematite inclusions with an HFSE signature within Olympic Dam zircon, is further evidence of Fe–Cl-metasomatism impacting on zircon (Figure 9B) occurring in the early, magmatic-to-hydrothermal transition stage [52], rather than post-metamictization. As highlighted by trace element analysis (Figure 4; Supplementary Materials Table S4), Donington Suite data plots at lower Dc and Fe2O3 values than Hiltaba Suite zircon (Figure 17A) reflect the lower U-contents and lesser alteration of zircon. Even in the case of zircon in Donington Suite granite from Wirrda Well, superimposed Fe–Cl-metasomatism is recognized due to the presence of adjacent Hiltaba Suite granites linked to the same batholith that generated the RDG, host to mineralization at Olympic Dam.
The Dc can affect U–Pb systematics in zircon and be manifest as discordance of data points via Pb loss (Figure 3). Accurate (CA)-ID-TIMS ages for magmatic crystallization can, nevertheless, be obtained [14], despite the presence of micron- to nanoscale domains, which in the present dataset do not appear to concentrate remobilized U or Pb. Such phenomena have been shown in previous studies [24,27], where they relate to alteration associated with highly metamict zircons.

6.3. Zircon Alteration Model and Magma “Fertility”

Concentrations of Cl within zircon containing Fe and other non-formula elements was measured and imaged for the first time in natural zircon. Electron microprobe maps show consistent distributions of Cl correlating with Fe in micron-scale banding (Figure 5). At the nanoscale, however, Cl was present as nanoparticles within Fe–Ca-metasomatized mottled bands (Figure 8, Figure 9, Figure 10 and Figure 11). Such an association between Fe and Cl would be expected in early IOCG fluids exsolved after granite crystallization [52]. Cl-rich nanoparticles became volumetrically more significant in samples from the vicinity of the Olympic Dam deposit, while smaller amounts of Cl were also measured in the nearby Wirrda Well Donington zircon as well as Charleston zircon. In contrast, the Cape Donington zircon displayed a complete absence of Cl. Therefore, we concluded that Cl was inherent to hydrothermal fluids derived from the Hiltaba Suite granites. Moreover, such Fe–Cl-bearing fluids also affected nearby older granites (Donington Suite), host to ~1.6 Ga IOCG mineralization at Wirrda Well.
Zircon metasomatism can be attributed to a low-fluid, mineral-buffered diffusion-reaction process in the beginning (only crystal zoning), followed by formation of chloro–hydroxy–zircon nanoprecipitates when fluids become highly focused (Figure 17B). Iron-rich nano-veinlets begin to interconnect across crystal zones previously metasomatized during increase in the flux of Fe-rich fluids (formation of hematite inclusions; Figure 9). Comparable phenomena in which a mineral can host nanoscale inclusions of the same mineral but with subtle yet distinct chemical differences, are known from Fe oxides, notably Si–Fe nanoprecipitates within silician magnetite. These are documented both from banded iron formation deposits [53] and the Olympic Dam deposit [54].
The preservation of pre-existing crystal orientation in zircon nanoprecipitates and host zircon is characteristic of reactions via a sharp interface, typical of CDRR replacement. Moreover, transient porosity developed during CDRR progression provides sites for deposition of nanoprecipitates, as suggested here and elsewhere (e.g., [55,56,57]). Hydrothermal fluids interacting with zircon can be considered acidic based on the hydrated-substitution model for zircon (Supplementary Materials A and B, Table S4). Fluid transport of Fe, from the granite-derived fluid and Ca + Al from breakdown of minerals in the granite via alkali–calcic alteration is likely to occur as chloride complexes and can lead to Si and/or Zr loss, according to the reactions:
Z r S i O 4 + 4 x H + Z r ( S i O 4 ) 1 x ( O H ) 4 x + x S i 4 +
Z r S i O 4 + x ( F e n + + n C l ) F e x Z r 1 x S i O 4 x n C l x n + x Z r O 2 + x ( n 2 ) O 2 ,
where n = 2 or 3, and x = 0 to 1.
The Fe–Cl-rich signature of post-magmatic fluids is recorded within the mottled textures of zircon and represents the onset of IOCG mineralization during the magmatic-to-hydrothermal transition. From ID-TIMS dating of hematite, iron metasomatism is known to have occurred at Olympic Dam shortly after granite crystallization (~2–4 Ma later) [13]. The present data support early, as well as prolonged, zircon alteration, illustrated here from areas that are overprinting primary growth (Figure 2B,D). The Fe–Cl-rich metasomatized zones within zircon were also imaged directly adjacent to magmatic relict domains (Figure 8A). Recrystallization of zircon surrounding xenotime veinlets was observed down to the nanoscale (Supplementary Materials C, Figure S3). Whereas zircon metasomatism overprinting primary growth zoning can be related to the earliest hydrothermal fluids exsolved at the depth of granite emplacement, later cycles of zircon alteration could happen during and/or after granite uplift as the IOCG mineralizing system evolved.
Overall, the increase in Cl associated with abundant chloro–hydroxy–zircon nanoprecipitates correlates with proximity to the Olympic Dam orebody. Such features can, thus, be considered a direct indicator of magma “fertility” (the ability of a magma to generate hydrothermal fluids that are sufficiently well-endowed with Cu, U, etc., to form a sizable ore deposit) in Hiltaba Suite granites. Furthermore, the presence of Cl in zircon from other Hiltaba Suite granites such as Charleston, indicates potential undiscovered IOCG mineralization using the zircon pathfinder.

7. Conclusions and Implications

1. From the case studies presented and comparison between them, zircon may be substituted with non-formula elements, notably Fe and Cl, released ~2–4 Ma after magmatic crystallization, during a widespread metasomatic event, prior to fluid pathway development during metamictization.
2. Chloro–hydroxy–zircon nanoprecipitates record the Fe–Cl-signature of fluids released during the magmatic-to-hydrothermal transition in IOCG systems, which volumetrically increase in granites of different age with proximity to Olympic Dam. Such nanoprecipitates represent a potential pathfinder to mineralization associated with fertile granites. Although we recommend that the data should initially be tested in other IOCG terranes, there exists potential for similar features to be identified within zircon from other deposit types involving penetrative metasomatic alteration.
3. Relationships between micron- and nanoscale observations in zircon are a largely untapped research field. With ever increasing instrumental resolution and analytical precision, nanoscale studies can provide new insights into processes carrying implications for geochronology, ore deposit genesis and mineral exploration.
4. In summary, this study has shown that early metasomatic alteration can be recorded in magmatic zircon and that these patterns can be preserved over geological time, even when superimposed by the effects of metamictization. This implies that nanoscale studies of magmatic zircon, when coupled with careful micron-scale characterization of textures and compositions, could represent a valuable new tool in mineral exploration.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-163X/9/6/364/s1, Supplementary Materials A: Extended details of methodology (Tables S1 and S2); Supplementary Materials B: Expanded data tables (Tables S3–S5); Supplementary Materials C: Additional figures (Figures S1, S2, and S3A–S3E).

Author Contributions

L.C.-D. and C.L.C. conceived this contribution. Analytical work was conducted L.C-.D., C.L.C., M.R.V.-I., M.D., W.K., J.X. and U.I.D., with A.S. operating the Titan Themis instrument and B.P.W. assisting with EPMA set-up. The manuscript was written by L.C.-D., C.L.C., M.R.V.-I. and N.J.C., with contributions from K.E., A.K.-C. and B.P.W.

Funding

This research was funded from the project “Trace elements in iron-oxides: deportment, distribution and application in ore genesis, geochronology, exploration and mineral processing”, supported by BHP Olympic Dam and the South Australian Government Mining and Petroleum Services Centre of Excellence. N.J.C. and K.E. acknowledge funding from the ARC Research Hub for Australian Copper–Uranium (Grant IH130200033).

Acknowledgments

We acknowledge the AMMRF for instrument access. Sarah Gilbert assisted with LA-ICP-MS analysis and Animesh Basak assisted with FIB-SEM operation. Chris Kalnins is thanked for preparing two of the TEM foils presented. We kindly thank two anonymous reviewers for comments that helped improve this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Groves, D.I.; Bierlein, F.P.; Meinert, L.D.; Hitzman, M.W. Iron oxide copper-gold (IOCG) deposits through Earth history: Implications for origin, lithospheric setting, and distinction from other epigenetic iron oxide deposits. Econ. Geol. 2010, 105, 641–654. [Google Scholar] [CrossRef]
  2. Barton, M.D. Iron oxide (-Cu-Au-REE-P-Ag-U-Co) systems. In Treatise on Geochemistry, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2014; Volume 13, pp. 515–541. [Google Scholar]
  3. Krneta, S.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Kontonikas-Charos, A. Apatite at Olympic Dam, South Australia: A petrogenetic tool. Lithos 2016, 262, 470–485. [Google Scholar] [CrossRef]
  4. Krneta, S.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.J. Numerical modeling of REE fractionation patterns in fluorapatite from the Olympic Dam deposit (South Australia). Minerals 2018, 8, 342. [Google Scholar] [CrossRef]
  5. Condon, D.J.; Schmitz, M.D. One hundred years of geochronology. Elements 2013, 9, 3–80. [Google Scholar] [CrossRef]
  6. Harley, S.L.; Kelly, N.M. Zircon Tiny but Timely. Elements 2007, 3, 13–18. [Google Scholar] [CrossRef]
  7. Rubatto, D.; Hermann, J. Zircon behaviour in deeply subducted rocks. Elements 2007, 3, 31–35. [Google Scholar] [CrossRef]
  8. Scherer, E.E.; Whitehouse, M.J.; Munker, C. Zircon as a Monitor of Crustal Growth. Elements 2007, 3, 19–24. [Google Scholar] [CrossRef]
  9. Hoskin, P.W.O.; Schaltegger, U. The composition of zircon and igneous and metamorphic petrogenesis. Rev. Mineral. Geochem. 2003, 53, 27–62. [Google Scholar] [CrossRef]
  10. Ewing, R.C.; Meldrum, A.; Wang, L.M.; Weber, W.J.; Corrales, L.R. Radiation damage in zircon. Rev. Mineral. Geochem. 2003, 53, 387–425. [Google Scholar] [CrossRef]
  11. Skirrow, R.G.; Bastrakov, E.N.; Barovich, K.; Fraser, G.L.; Creaser, R.A.; Fanning, C.M. Timing of iron oxide Cu-Au-(U) hydrothermal activity and Nd isotope constraints on metal sources in the Gawler craton, South Australia. Econ Geol. 2007, 102, 1441–1470. [Google Scholar] [CrossRef]
  12. Ehrig, K.; McPhie, J.; Kamenetsky, V. Geology and Mineralogical Zonation of the Olympic Dam Iron Oxide Cu-U-Au-Ag Deposit, South Australia. In Geology and Genesis of Major Copper Deposits and Districts of the World: A Tribute to Richard H. Sillitoe; Hedenquist, J.W., Harris, M., Camus, F., Eds.; Society of Economic Geologists Special Publication 16; Society of Economic Geologists: Littleton, CO, USA, 2012; pp. 237–267. [Google Scholar]
  13. Courtney-Davies, L.; Tapster, S.; Ciobanu, C.; Cook, N.J.; Verdugo-Ihl, M.R.; Ehrig, K.J.; Kennedy, A.K.; Gilbert, S.E.; Condon, D.J.; Wade, B.P. A multi-technique evaluation of hydrothermal hematite U-Pb isotope systematics: Implications for ore deposit geochronology. Chem. Geol. 2019, 513, 54–72. [Google Scholar] [CrossRef]
  14. Cherry, A.R.; Ehrig, K.; Kamenetsky, V.S.; McPhie, J.; Crowley, J.L.; Kamenetsky, M.B. Precise geochronological constraints on the origin, setting and incorporation of ca. 1.59 Ga surficial facies into the Olympic Dam Breccia Complex, South Australia. Precambr. Res. 2018, 315, 162–178. [Google Scholar] [CrossRef]
  15. Ewing, R.C.; Lutze, W.; Weber, W.J. Zircon—A Host-Phase for the Disposal of Weapons Plutonium. J. Mater. Res. 1995, 10, 243–246. [Google Scholar] [CrossRef]
  16. Holland, H.D.; Gottfried, D. The effect of nuclear radiation on the structure of zircon. Acta Cryst. 1955, 8, 291–300. [Google Scholar] [CrossRef]
  17. Mattinson, J. Zircon U–Pb chemical abrasion (“CA-TIMS”) method: Combined annealing and multi-step partial dissolution analysis for improved precision and accuracy of zircon ages. Chem. Geol. 2005, 220, 47–66. [Google Scholar] [CrossRef]
  18. Geisler, T.; Schaltegger, U.; Tomaschek, F. Re-equilibration of zircon in aqueous fluids and melts. Elements 2007, 3, 43–50. [Google Scholar] [CrossRef]
  19. Geisler, T.; Rashwan, A.A.; Rahn, M.K.W.; Poller, U.; Zwingmann, H.; Pidgeon, R.T.; Schleicher, H.; Tomaschek, F. Low temperature hydrothermal alteration of natural metamict zircons from the Eastern Desert, Egypt. Mineral. Mag. 2003, 67, 485–508. [Google Scholar] [CrossRef]
  20. Tang, F.; Taylor, J.M.; Einsle, J.F.; Borlina, C.S.; Fu, R.R.; Weiss, B.P.; Williams, H.M.; Williams, W.; Nagy, L.; Midgley, P.; et al. Secondary magnetite in ancient zircon precludes analysis of a Hadean geodynamo. Proc. Natl. Acad. Sci. USA 2019, 116, 407–412. [Google Scholar] [CrossRef]
  21. Lu, Y.J.; Loucks, R.R.; Fiorentini, M.L.; McCuaig, T.C.; Evans, N.J.; Yang, Z.M.; Hou, Z.Q.; Kirkland, C.L.; Parra-Avila, L.A.; Kobussen, A. Zircon compositions as a pathfinder for porphyry Cu ± Mo ± Au deposits. Econ. Geol. 2016, 19, 329–347. [Google Scholar]
  22. Belousova, E.A.; Griffin, W.L.; O’Reilly, S.Y.; Fisher, N.I. Igneous zircon: Trace element composition as an indicator of source rock type. Contrib. Mineral. Petrol. 2002, 143, 602–622. [Google Scholar] [CrossRef]
  23. Utsunomiya, S.; Ewing, R.C. Application of high-angle annular dark field scanning transmission electron microscopy—Energy dispersive X-ray spectrometry, and energy-filtered transmission electron microscopy to the characterization of nanoparticles in the environment. Environ. Sci. Technol. 2003, 37, 786–791. [Google Scholar] [CrossRef]
  24. Utsunomiya, S.; Palenik, C.S.; Valley, J.W.; Cavosie, A.J.; Wilde, S.A.; Ewing, R.C. Nanoscale occurrence of Pb in an Archean zircon. Geochim. Cosmochim. Acta 2004, 68, 4679–4686. [Google Scholar] [CrossRef]
  25. Utsunomiya, S.; Valley, J.W.; Cavosie, A.J.; Wilde, S.A.; Ewing, R.C. Radiation damage and alteration of zircon from a 3.3 Ga porphyritic granite from the Jack Hills, Western Australia. Chem. Geol. 2007, 236, 92–111. [Google Scholar] [CrossRef]
  26. Seydoux-Guillaume, A.M.; Bingen, B.; Paquette, J.L.; Bosse, V. Nanoscale evidence for uranium mobility in zircon and the discordance of U-Pb chronometers. Earth Plan. Sci. Lett. 2015, 409, 43–48. [Google Scholar] [CrossRef]
  27. Kusiak, M.A.; Dunkley, D.J.; Wirth, R.; Whitehouse, M.J.; Wilde, S.A.; Marquardtd, K. Metallic lead nanospheres discovered in ancient zircons. Proc. Natl. Acad. Sci. USA 2015, 112, 4958–4963. [Google Scholar] [CrossRef] [Green Version]
  28. Ciobanu, C.L.; Cook, N.J.; Maunders, C.; Wade, B.P.; Ehrig, K. Focused Ion Beam and Advanced Electron Microscopy for Minerals: Insights and Outlook from Bismuth Sulphosalts. Minerals 2016, 6, 112. [Google Scholar] [CrossRef]
  29. Cook, N.J.; Ciobanu, C.L.; Ehrig, K.; Slattery, A.; Verdugo-Ihl, M.R.; Courtney-Davies, L.; Gao, W. Advances and opportunities in ore mineralogy. Minerals 2017, 7, 233. [Google Scholar] [CrossRef]
  30. Goldfarb, R.J.; Bradley, D.; Leach, D.L. Secular variation in economic geology. Econ. Geol. 2010, 105, 459–465. [Google Scholar] [CrossRef]
  31. Allen, S.R.; McPhie, J.; Ferris, G.; Simpson, C. Evolution and architecture of a large felsic igneous province in western Laurentia: The 1.6 Ga Gawler Range Volcanics, South Australia. J. Volcanol. Geotherm. Res. 2008, 172, 132–147. [Google Scholar] [CrossRef]
  32. Keyser, W.M.; Ciobanu, C.L.; Cook, N.J.; Johnson, G.; Feltus, F.; Johnson, S.; Dmitrijeva, M.; Ehrig, K.; Nguyen, P. Petrography and trace element signatures of iron-oxides in deposits from the Middleback Ranges, South Australia: From banded iron formation to ore. Ore Geol. Rev. 2018, 93, 337–360. [Google Scholar] [CrossRef]
  33. Domnick, U.; Cook, N.J.; Bluck, R.; Brown, C.; Ciobanu, C.L. Petrography of granitoids from the Samphire Pluton, South Australia: Implications for uranium mineralization in overlying sediments. Lithos 2018, 300–301, 1–19. [Google Scholar] [CrossRef]
  34. Keyser, W.; Ciobanu, C.L.; Cook, N.J.; Dmitrijeva, M.; Courtney-Davies, L.; Feltus, H.; Gilbert, S.; Johnson, G.; Ehrig, K. Iron-oxides constrain BIF evolution in terranes with protracted geological histories: The Iron Count prospect, Middleback Ranges, South Australia. Lithos 2019, 324, 20–38. [Google Scholar] [CrossRef]
  35. Dmitrijeva, M.; Metcalfe, A.V.; Ciobanu, C.L.; Cook, N.J.; Frenzel, M.; Keyser, W.M.; Johnson, G.; Ehrig, K. Discrimination and variance structure of trace element signatures in hematite: A case study of BIF-mineralization from the Middleback Ranges, South Australia. Math. Geosci. 2018, 50, 381–415. [Google Scholar] [CrossRef]
  36. Creaser, R.A.; Fanning, C.M. A U-Pb zircon study of the Mesoproterozoic Charleston Granite, Gawler Craton, South Australia. Aust. J. Earth Sci. 1993, 40, 519–526. [Google Scholar] [CrossRef]
  37. Skirrow, R.G.; van der Wielen, S.E.; Champion, D.C.; Czarnota, K.; Thiel, S. Lithospheric architecture and mantle metasomatism linked to iron oxide Cu-Au ore formation: Multidisciplinary evidence from the Olympic Dam Region, South Australia. Geochem. Geophys. Geosyst. 2018, 19, 2673–2705. [Google Scholar] [CrossRef]
  38. Wade, C.E.; Payne, J.L.; Barovich, K.M.; Reid, A.J. Heterogeneity of the sub-continental lithospheric mantle and ‘non-juvenile’ mantle additions to a Proterozoic silicic large igneous province. Lithos 2019, 340–341, 87–107. [Google Scholar] [CrossRef]
  39. Ovalle, J.T.; La Cruz, N.L.; Reich, M.; Barra, F.; Simon, A.C.; Konecke, B.A.; Rodriguez-Mustafa, M.A.; Deditius, A.P.; Childress, T.M.; Morata, D. Formation of massive iron deposits linked to explosive volcanic eruptions. Sci. Rep. 2018, 8, 14855. [Google Scholar] [CrossRef]
  40. Kontonikas-Charos, A.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Krneta, S.; Kamenetsky, V.S. Feldspar evolution in the Roxby Downs Granite, host to Fe-oxide Cu-Au-(U) mineralisation at Olympic Dam, South Australia. Ore Geol. Rev. 2017, 80, 838–859. [Google Scholar] [CrossRef]
  41. Krneta, S.; Cook, N.J.; Ciobanu, C.L.; Ehrig, K.; Kontonikas-Charos, A. The Wirrda Well and Acropolis prospects Gawler Craton, South Australia: Insights into evolving fluid conditions through apatite chemistry. J. Geochem. Explor. 2017, 181, 276–291. [Google Scholar] [CrossRef]
  42. Kontonikas-Charos, A.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Krneta, S.; Kamenetsky, V.S. Rare earth element geochemistry of feldspars: Examples from Fe-oxide Cu-Au systems in the Olympic Cu-Au Province, South Australia. Mineral. Petrol. 2018, 112, 145–172. [Google Scholar] [CrossRef]
  43. Jagodzinski, E.A. Compilation of SHRIMP U-Pb geochronological data Olympic Domain, Gawler Craton, South Australia, 2001–2003. Geosci. Aust. Rec. 2005, 20, 2005. [Google Scholar]
  44. Keyser, W.; Ciobanu, C.L.; Cook, N.J.; Courtney-Davies, L.; Kennedy, A.; Wade, B.P.; Ehrig, K.; Dmitrijeva, M.; Kontonikas-Charos, A.; Feltus, H.; et al. Petrographic and geochronological constraints on the granitic basement to the Middleback Ranges, South Australia. Precambr. Res. 2019, 324, 170–193. [Google Scholar] [CrossRef]
  45. Pérez-Soba, C.; Villaseca, C.; Gonzáles del Tánago, J.; Nasdala, L. The composition of zircon in the peraluminous Hercynian granites of the Spanish Central System batholith. Can. Mineral. 2007, 45, 509–527. [Google Scholar] [CrossRef]
  46. Verdugo-Ihl, M.R.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Courtney-Davies, L.; Gilbert, S. Textures and U-W-Sn-Mo signatures in hematite from the Cu-U-Au-Ag orebody at Olympic Dam, South Australia: Defining the archetype for IOCG deposits. Ore Geol. Rev. 2017, 91, 173–195. [Google Scholar] [CrossRef]
  47. Kontonikas-Charos, A.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Ismail, R.; Krneta, S.; Basak, A. Feldspar mineralogy and rare earth element (re)mobilization in iron-oxide copper gold systems from South Australia: A nanoscale study. Mineral. Mag. 2018, 82, S173–S197. [Google Scholar] [CrossRef]
  48. Salje, E.K.H.; Chrosch, J.; Ewing, R.C. Is “metamictization” of zircon a phase transition? Am. Mineral. 1999, 84, 1107–1116. [Google Scholar] [CrossRef]
  49. Ríos, S.; Malcherek, T.; Salje, E.K.H.; Domeneghetti, C. Localized defects in radiation-damaged zircon. Acta Cryst. 2000, 56, 947–952. [Google Scholar] [CrossRef]
  50. Nasdala, L.; Kronz, A.; Wirth, R.; Váczi, T.; Pérez-Soba, C.; Willner, A.; Kennedy, A.K. Alteration of radiation-damaged zircon and the related phenomenon of deficient electron microprobe totals. Geochim. Cosmochim. Acta 2009, 73, 1637–1650. [Google Scholar] [CrossRef]
  51. Wingate, M.T.D.; Campbell, I.H.; Compston, W.; Gibson, G.M. Ion microprobe U–Pb ages for Neoproterozoic basaltic magmatism in south-central Australia and implications for the breakup of Rodinia. Precambr. Res. 1998, 87, 135–159. [Google Scholar] [CrossRef]
  52. Verdugo-Ihl, M.R.; Ciobanu, C.L.; Cook, N.J.; Ehrig, K.; Courtney-Davies, L. Defining early stages of IOCG systems: Evidence from iron-oxides in the outer shell of the Olympic Dam deposit, South Australia. Mineral. Depos. 2019. [Google Scholar] [CrossRef]
  53. Xu, H.; Shen, Z.; Konishi, H. Si-magnetite nano-precipitates in silician magnetite from banded iron formation: Z-contrast imaging and ab initio study. Am. Mineral. 2014, 99, 2196–2202. [Google Scholar] [CrossRef]
  54. Ciobanu, C.L.; Verdugo-Ihl, M.R.; Slattery, A.; Cook, N.J.; Ehrig, K.; Courtney-Davies, L.; Wade, B.P. Silician magnetite: Si-Fe-nanoprecipitates and other mineral inclusions in magnetite from the Olympic Dam deposit, South Australia. Minerals 2019, 9, 311. [Google Scholar] [CrossRef]
  55. Tooth, B.; Ciobanu, C.L.; Green, L.; O’Neill, B.; Brugger, J. Bi-melt formation and gold scavenging from hydrothermal fluids: An experimental study. Geochim. Cosmochim. Acta 2011, 75, 5423–5443. [Google Scholar] [CrossRef]
  56. Verdugo-Ihl, M.R.; Ciobanu, C.L.; Cook, N.J.; Slattery, A.; Ehrig, K. Copper nanoparticles along fluid inclusion trails in hematite. In Proceedings of the Goldschmidt 2018, Boston, MA, USA, 12–17 August 2018. [Google Scholar]
  57. Owen, N.D.; Ciobanu, C.L.; Cook, N.J.; Slattery, A.; Basak, A. Nanoscale study of clausthalite-bearing symplectites in Cu-Au-(U) ores: Implications for ore genesis. Minerals 2018, 8, 67. [Google Scholar] [CrossRef]
  58. Meldrum, A.; Boatner, L.A.; Weber, W.J.; Ewing, R.C. Radiation damage in zircon and monazite. Geochim. Cosmochim. Acta 1988, 62, 2509–2520. [Google Scholar] [CrossRef]
Figure 1. (A) Geological map of the eastern Gawler Craton and (B) the Olympic Dam district; data were taken from the Geological Survey of South Australia (https://map.sarig.sa.gov.au/). Locations of the five granite samples addressed in this study and the outline of the Olympic Cu–Au district are marked. (C) Table of ages and corresponding uncertainties (colored bars) for each granite sample studied here, with geochronology conducted either in this study (LA-ICP-MS) or in previous published work as indicated. The orange and red uncertainty bars represent the Donington Suite and Hiltaba Intrusive Suite, respectively, and their location within either the Eyre Peninsula or Olympic Cu–Au Province. Corresponding hand specimen photographs of each sample are shown below the age diagram, indicating degree of granite alteration. References: Keyser et al. (2019) [34]; Creaser and Fanning (1993) [36]; Cherry et al. (2018) [14].
Figure 1. (A) Geological map of the eastern Gawler Craton and (B) the Olympic Dam district; data were taken from the Geological Survey of South Australia (https://map.sarig.sa.gov.au/). Locations of the five granite samples addressed in this study and the outline of the Olympic Cu–Au district are marked. (C) Table of ages and corresponding uncertainties (colored bars) for each granite sample studied here, with geochronology conducted either in this study (LA-ICP-MS) or in previous published work as indicated. The orange and red uncertainty bars represent the Donington Suite and Hiltaba Intrusive Suite, respectively, and their location within either the Eyre Peninsula or Olympic Cu–Au Province. Corresponding hand specimen photographs of each sample are shown below the age diagram, indicating degree of granite alteration. References: Keyser et al. (2019) [34]; Creaser and Fanning (1993) [36]; Cherry et al. (2018) [14].
Minerals 09 00364 g001
Figure 2. BSE images of grains #1 (A and C) and #2 (B and D) from OD distal zircon (sample LCD47). Note the mottled dark areas were either bands within crystal oscillatory zoning (C) or scalloped domains superimposed onto crystal zoning (D). (E) BSE (top) and HAADF STEM image (bottom) showing crystal oscillatory zoning for Cape Donington zircon (PL8). (F) Cathodoluminescence image (left) and LA-ICP-MS element maps (right) showing oscillatory zoned patterns for Y and U in coarser zircon grains from PL8. Scale in 10n counts per second (extended maps are given in Supplementary Materials C, Figure S2).
Figure 2. BSE images of grains #1 (A and C) and #2 (B and D) from OD distal zircon (sample LCD47). Note the mottled dark areas were either bands within crystal oscillatory zoning (C) or scalloped domains superimposed onto crystal zoning (D). (E) BSE (top) and HAADF STEM image (bottom) showing crystal oscillatory zoning for Cape Donington zircon (PL8). (F) Cathodoluminescence image (left) and LA-ICP-MS element maps (right) showing oscillatory zoned patterns for Y and U in coarser zircon grains from PL8. Scale in 10n counts per second (extended maps are given in Supplementary Materials C, Figure S2).
Minerals 09 00364 g002
Figure 3. Concordia diagram (A) for OD proximal zircon (LCD13, red ellipses) and Wirrda Well zircon (LCD17, orange ellipses; U–Pb data are presented in Supplementary Materials B, Table S3). Circled letters correspond to LA-ICP-MS craters in BSE images of the dated grains (B and C).
Figure 3. Concordia diagram (A) for OD proximal zircon (LCD13, red ellipses) and Wirrda Well zircon (LCD17, orange ellipses; U–Pb data are presented in Supplementary Materials B, Table S3). Circled letters correspond to LA-ICP-MS craters in BSE images of the dated grains (B and C).
Minerals 09 00364 g003
Figure 4. (A) LA-ICP-MS chondrite-normalized REY (REE + Y) fractionation trends for Donington Suite (light and dark blue) and Hiltaba Intrusive Suite (red and orange) samples. (B and C) PCA of LA-ICP-MS data, used to distinguish associated element signatures, and (D) a ternary Fe–Al–Hf diagram demonstrating the degree of zircon metasomatism (LA-ICP-MS data; Supplementary Materials B, Table S5). (E) Plot of (REY+Sc)3+ versus P5+ (EPMA data). (F) Plot of “non-formula” components Σ(Fe2O3 + Al2O3 + MnO + CaO) versus apfu Σ(OH + Cl + F), respectively. The EPMA diagrams illustrate the relationship between substitution of both non-formulae elements and OH in zircon. Open and closed data points represent bright and dark bands in zircon, respectively; color codes correspond to samples in AD (EPMA data; Supplementary Materials B, Table S4). See main text and Supplementary Materials A for additional explanation.
Figure 4. (A) LA-ICP-MS chondrite-normalized REY (REE + Y) fractionation trends for Donington Suite (light and dark blue) and Hiltaba Intrusive Suite (red and orange) samples. (B and C) PCA of LA-ICP-MS data, used to distinguish associated element signatures, and (D) a ternary Fe–Al–Hf diagram demonstrating the degree of zircon metasomatism (LA-ICP-MS data; Supplementary Materials B, Table S5). (E) Plot of (REY+Sc)3+ versus P5+ (EPMA data). (F) Plot of “non-formula” components Σ(Fe2O3 + Al2O3 + MnO + CaO) versus apfu Σ(OH + Cl + F), respectively. The EPMA diagrams illustrate the relationship between substitution of both non-formulae elements and OH in zircon. Open and closed data points represent bright and dark bands in zircon, respectively; color codes correspond to samples in AD (EPMA data; Supplementary Materials B, Table S4). See main text and Supplementary Materials A for additional explanation.
Minerals 09 00364 g004
Figure 5. (A) BSE image of oscillatory-zoned Wirrda Well–Donington Suite zircon (sample LCD17) that was dated (LA-ICP-MS craters) and used for nanoscale study (foil location as marked). (B). EPMA maps for Cl and Fe2O3 (in wt.%), showing that both elements appear homogeneously enriched and positively correlated with one another within oscillatory bands. (C,D) BSE image and EPMA maps for Cl and Fe2O3 (wt.%) of OD proximal zircon (sample LCD 13; foil marked at the top). Note the high concentrations of both elements within the middle part which was affected by fractures and re-crystallization. In contrast, crystal zoning on the rim shows inverse correlation between the two elements (e.g., banding intersecting the FIB cut).
Figure 5. (A) BSE image of oscillatory-zoned Wirrda Well–Donington Suite zircon (sample LCD17) that was dated (LA-ICP-MS craters) and used for nanoscale study (foil location as marked). (B). EPMA maps for Cl and Fe2O3 (in wt.%), showing that both elements appear homogeneously enriched and positively correlated with one another within oscillatory bands. (C,D) BSE image and EPMA maps for Cl and Fe2O3 (wt.%) of OD proximal zircon (sample LCD 13; foil marked at the top). Note the high concentrations of both elements within the middle part which was affected by fractures and re-crystallization. In contrast, crystal zoning on the rim shows inverse correlation between the two elements (e.g., banding intersecting the FIB cut).
Minerals 09 00364 g005
Figure 6. (A) HAADF STEM images (A and C) and EDS STEM maps for Fe (B and D) for Wirrda Well-Donington zircon showing crystal oscillatory zoning (A and B). The dark bands (even numbers) were enriched in Fe and show mottled texture in the image, but not on the Fe map (B). Note that the oscillatory banding was crosscut by a xenotime-bearing veinlet (extended map shown in Supplementary Materials C, Figure S3). (C,D) Higher resolution maps showing the interface between light and dark bands confirmed that the mottled texture (C) was not reflected in the Fe distribution (D).
Figure 6. (A) HAADF STEM images (A and C) and EDS STEM maps for Fe (B and D) for Wirrda Well-Donington zircon showing crystal oscillatory zoning (A and B). The dark bands (even numbers) were enriched in Fe and show mottled texture in the image, but not on the Fe map (B). Note that the oscillatory banding was crosscut by a xenotime-bearing veinlet (extended map shown in Supplementary Materials C, Figure S3). (C,D) Higher resolution maps showing the interface between light and dark bands confirmed that the mottled texture (C) was not reflected in the Fe distribution (D).
Minerals 09 00364 g006
Figure 7. (A) HAADF STEM image and EDS STEM element maps of Charleston zircon (Foil #6) showing enrichment in trace elements outside of a magmatic domain as marked. Note the sharp boundary (dashed lines) for Fe but not the other elements, which instead display a diffuse boundary. Trace element (re)mobilization (all elements except Fe) along this boundary was observed along the adjacent fracture and the U-bearing fine particle. (BD) HAADF STEM images showing zircon details along the fracture as marked. The U-bearing fine particle (D) displayed a dark rim and was beam-sensitive, indicating that it was most likely amorphous (it could not be imaged at high-resolution due to the thickness of the foil at this location). The image in (B) is the inset box indicated on (A), and image (C) shows the inset box indicated on (B). Image (D) shows detail of the inclusion from image (A).
Figure 7. (A) HAADF STEM image and EDS STEM element maps of Charleston zircon (Foil #6) showing enrichment in trace elements outside of a magmatic domain as marked. Note the sharp boundary (dashed lines) for Fe but not the other elements, which instead display a diffuse boundary. Trace element (re)mobilization (all elements except Fe) along this boundary was observed along the adjacent fracture and the U-bearing fine particle. (BD) HAADF STEM images showing zircon details along the fracture as marked. The U-bearing fine particle (D) displayed a dark rim and was beam-sensitive, indicating that it was most likely amorphous (it could not be imaged at high-resolution due to the thickness of the foil at this location). The image in (B) is the inset box indicated on (A), and image (C) shows the inset box indicated on (B). Image (D) shows detail of the inclusion from image (A).
Minerals 09 00364 g007
Figure 8. Trace element distributions in OD distal zircon from Foil #1. (A) HAADF STEM image and EDS STEM maps for Fe, Ca, and Cl showing oscillatory zoning and mottled textures. Enrichment in Ca and Cl occurred throughout the mottled bands. (B) Profile across the oscillatory zoning displayed positive correlation between Ca and Cl across the most intensively mottled bands. Iron distributions broadly correlated with Ca and Cl across the same bands but showed small-amplitude variation reflecting the finest-scale oscillatory banding. The HAADF intensity signal showed a negative correlation with the Ca- and Cl-enriched bands. (C) Image and EDS STEM maps (extended maps in Supplementary Materials C, Figure S3) of a Cl-rich area as marked in (A). Note the inconsistency in correlation between Cl with O and Zr (compare the areas outlined on the figure), e.g., high Zr and low-O where Cl-enrichment is moderate, but no changes for O and Zr for nanodomains with the highest Cl content.
Figure 8. Trace element distributions in OD distal zircon from Foil #1. (A) HAADF STEM image and EDS STEM maps for Fe, Ca, and Cl showing oscillatory zoning and mottled textures. Enrichment in Ca and Cl occurred throughout the mottled bands. (B) Profile across the oscillatory zoning displayed positive correlation between Ca and Cl across the most intensively mottled bands. Iron distributions broadly correlated with Ca and Cl across the same bands but showed small-amplitude variation reflecting the finest-scale oscillatory banding. The HAADF intensity signal showed a negative correlation with the Ca- and Cl-enriched bands. (C) Image and EDS STEM maps (extended maps in Supplementary Materials C, Figure S3) of a Cl-rich area as marked in (A). Note the inconsistency in correlation between Cl with O and Zr (compare the areas outlined on the figure), e.g., high Zr and low-O where Cl-enrichment is moderate, but no changes for O and Zr for nanodomains with the highest Cl content.
Minerals 09 00364 g008
Figure 9. Trace element distributions in OD proximal zircon (Foil #3). (A) HAADF STEM image showing oscillatory banding, mottled textures, and an inclusion of hematite in zircon. High field strength element (HFSE) enrichment in hematite shown on the EDS STEM map in (B). Identification of hematite from the atomic-scale resolution image in (C) and Fast Fourier Transform (FFT) in (D) as marked. (E,F) HAADF STEM images and EDS STEM maps obtained from areas as marked on (A) (rotated 90 degrees clockwise) showing Fe and Cl enrichment throughout the mottled bands. Note the irregular distribution of Cl relative to Fe in both (E) and (F). Importantly, the Fe map in (F) shows Fe-rich veinlets crosscutting brighter bands.
Figure 9. Trace element distributions in OD proximal zircon (Foil #3). (A) HAADF STEM image showing oscillatory banding, mottled textures, and an inclusion of hematite in zircon. High field strength element (HFSE) enrichment in hematite shown on the EDS STEM map in (B). Identification of hematite from the atomic-scale resolution image in (C) and Fast Fourier Transform (FFT) in (D) as marked. (E,F) HAADF STEM images and EDS STEM maps obtained from areas as marked on (A) (rotated 90 degrees clockwise) showing Fe and Cl enrichment throughout the mottled bands. Note the irregular distribution of Cl relative to Fe in both (E) and (F). Importantly, the Fe map in (F) shows Fe-rich veinlets crosscutting brighter bands.
Minerals 09 00364 g009
Figure 10. Aspects of nanoscale inclusions + pores throughout the mottled bands in the OD distal zircon (Foil #1). (A) Image and EDS STEM maps for a Cl-rich area selected from Figure 8C showing that the highest Cl was recorded in the smallest particles with no differences on the major or trace elements maps. Inverse correlation between O and Zr was recorded for the coarser, darkest inclusions with moderate enrichment in Cl. Note the weak variation in Fe and Ca, but not in Pb. (B) EDS spectrum for Cl-richest nanoparticle circled in (A). (C) EDS STEM maps (left) and profiles (right) across a particle with medium Cl-enrichment in (A). Profiles show high concentrations of Cl at the particle margin with increases in both major and minor elements but a strong HAADF signal decrease.
Figure 10. Aspects of nanoscale inclusions + pores throughout the mottled bands in the OD distal zircon (Foil #1). (A) Image and EDS STEM maps for a Cl-rich area selected from Figure 8C showing that the highest Cl was recorded in the smallest particles with no differences on the major or trace elements maps. Inverse correlation between O and Zr was recorded for the coarser, darkest inclusions with moderate enrichment in Cl. Note the weak variation in Fe and Ca, but not in Pb. (B) EDS spectrum for Cl-richest nanoparticle circled in (A). (C) EDS STEM maps (left) and profiles (right) across a particle with medium Cl-enrichment in (A). Profiles show high concentrations of Cl at the particle margin with increases in both major and minor elements but a strong HAADF signal decrease.
Minerals 09 00364 g010
Figure 11. (A) HAADF STEM images of inclusions representative of all samples with mottled textures showing their variation in size and morphology as marked. Note the presence of rounded pores in the coarsest particles. (B) Image and EDS STEM maps of a single particle from the OD distal zircon (Foil #1). Note the strong variation in O and Zr concentration but not Si relative to host zircon and the weak relative enrichment in Cl, U, and Pb. Ca-rich nanoparticles occur just outside, at the lower edge of the mapped particle.
Figure 11. (A) HAADF STEM images of inclusions representative of all samples with mottled textures showing their variation in size and morphology as marked. Note the presence of rounded pores in the coarsest particles. (B) Image and EDS STEM maps of a single particle from the OD distal zircon (Foil #1). Note the strong variation in O and Zr concentration but not Si relative to host zircon and the weak relative enrichment in Cl, U, and Pb. Ca-rich nanoparticles occur just outside, at the lower edge of the mapped particle.
Minerals 09 00364 g011
Figure 12. (AD) High-resolution HAADF STEM images of zircon on four main zone axes as marked.
Figure 12. (AD) High-resolution HAADF STEM images of zircon on four main zone axes as marked.
Minerals 09 00364 g012
Figure 13. From top to bottom, FFT, HAADF image, STEM simulation, and crystal structure model for each of the main zone axes (A–D, as marked) imaged in zircon (Figure 12). The images and simulations were interpreted in agreement with crystal structure models on zone axes as obtained from indexing FFTs obtained from images in Figure 11. The two zone axes with dumbbell atomic arrangements (C,D) cannot be resolved using the ideal space group for zircon, I41/amd, due to the presence of forbidden reflections (00l; l ≠ 4n and hh0; l ≠ 4n; marked in yellow). Indexing, crystal models, and STEM simulations were obtained using the I212121 space group. The models and simulations did not change when using the I41/amd space group.
Figure 13. From top to bottom, FFT, HAADF image, STEM simulation, and crystal structure model for each of the main zone axes (A–D, as marked) imaged in zircon (Figure 12). The images and simulations were interpreted in agreement with crystal structure models on zone axes as obtained from indexing FFTs obtained from images in Figure 11. The two zone axes with dumbbell atomic arrangements (C,D) cannot be resolved using the ideal space group for zircon, I41/amd, due to the presence of forbidden reflections (00l; l ≠ 4n and hh0; l ≠ 4n; marked in yellow). Indexing, crystal models, and STEM simulations were obtained using the I212121 space group. The models and simulations did not change when using the I41/amd space group.
Minerals 09 00364 g013
Figure 14. (AD) Images showing defects along a fracture crosscutting one of the fine particles (Foil #1). The FFT pattern obtained from the image (inset D) shows a square sub-pattern with weaker reflections superimposed onto the main reflections. Yellow, dashed lines mark the defect shape. (E) Intensity profile across the defect in (D) showing modifications in the Zr and Si signals as marked.
Figure 14. (AD) Images showing defects along a fracture crosscutting one of the fine particles (Foil #1). The FFT pattern obtained from the image (inset D) shows a square sub-pattern with weaker reflections superimposed onto the main reflections. Yellow, dashed lines mark the defect shape. (E) Intensity profile across the defect in (D) showing modifications in the Zr and Si signals as marked.
Minerals 09 00364 g014
Figure 15. Atomic-scale resolution HAADF STEM images of defects in zircon. (A) Square-shaped defect (yellow outline) in U-richer bands from Cape Donington (Foil #7), tilted on the (100) zone axis. The doubling of Zr atoms observed along the c-axis suggests atom “stretching) along b induced by strain, attributable to α-recoil damage during U-decay. (B) Simple linear defect in Wirrda Well–Donington zircon (Foil #4) observed at the edge of a xenotime veinlet (blue overlay); sample tilted on the (001) zone axis in zircon. Screw dislocations (dashed lines) are marked by atom loops in the defect plane (yellow dots) relative to Zr atoms (white dots) outside the defect. (CD) Screw dislocations (dashed lines) along a fault or at the edge of a planar defect on the (001) zone axis in Charleston zircon (Foil #6). Note irregular atom loops along the fault plane (blue overlays). (EF) Defects in the Fe-rich, mottled bands from OD distal zircon (Foil #1) with sample tilted on the [ 1 1 ¯ 1 ] zone axis in zircon. Antiphase atom displacement (dashed line) suggesting lattice stretching along (011) directions in (E). (F) Details of the defect (dotted line) in Figure 14D showing the size and intensity modifications in Zr and Si atoms along the (011) direction (colors as in Figure 14E). The image shows the inferred missing Si atoms as dark, regular spots along this direction.
Figure 15. Atomic-scale resolution HAADF STEM images of defects in zircon. (A) Square-shaped defect (yellow outline) in U-richer bands from Cape Donington (Foil #7), tilted on the (100) zone axis. The doubling of Zr atoms observed along the c-axis suggests atom “stretching) along b induced by strain, attributable to α-recoil damage during U-decay. (B) Simple linear defect in Wirrda Well–Donington zircon (Foil #4) observed at the edge of a xenotime veinlet (blue overlay); sample tilted on the (001) zone axis in zircon. Screw dislocations (dashed lines) are marked by atom loops in the defect plane (yellow dots) relative to Zr atoms (white dots) outside the defect. (CD) Screw dislocations (dashed lines) along a fault or at the edge of a planar defect on the (001) zone axis in Charleston zircon (Foil #6). Note irregular atom loops along the fault plane (blue overlays). (EF) Defects in the Fe-rich, mottled bands from OD distal zircon (Foil #1) with sample tilted on the [ 1 1 ¯ 1 ] zone axis in zircon. Antiphase atom displacement (dashed line) suggesting lattice stretching along (011) directions in (E). (F) Details of the defect (dotted line) in Figure 14D showing the size and intensity modifications in Zr and Si atoms along the (011) direction (colors as in Figure 14E). The image shows the inferred missing Si atoms as dark, regular spots along this direction.
Minerals 09 00364 g015
Figure 16. Schematic representation of atom arrangements in defects indicating lattice stretching and screw dislocations leading to expansion of the zircon structure. Such aspects are attributable to self-induced irradiation damage effects.
Figure 16. Schematic representation of atom arrangements in defects indicating lattice stretching and screw dislocations leading to expansion of the zircon structure. Such aspects are attributable to self-induced irradiation damage effects.
Minerals 09 00364 g016
Figure 17. (A) Plot of Fe2O3 (wt.%) versus Dc (×1018 α-decays/g) for all zircon samples as marked and using EPMA data from Supplementary Materials B, Table S4 (four analyses with Dc > 25 were excluded). A sigmoid through the data was fit by non-linear least squares regression using the software “R”. Although the first percolation point (FPP) for zircon was marked at Dc = 2 (as defined by Salje et al. [48]), a sharp increase in wt.% Fe2O3 was seen only at Dc = ~3. The FPP in Olympic Dam zircon was only reached at least 100 Ma after crystallization, using the dose concentration calculation of Meldrum et al. [58]. Also, note the scatter in the dataset concordant with the weak correlation between the highest measured non-formula element (Fe) and Dc. (B) Schematic diagram illustrating evolution of magmatic zircon from granites in the Olympic Cu–Au Province. Colored areas in (A) correspond to the main evolutionary stages in (B). Note that Cape Donington zircon, with the lowest non-formula elements, plots in a magmatic domain that straddles the FPP line. In contrast, analyses with high U concentration corresponded to intensively altered stages (green), featuring superimposed textures with all the characteristic of metamict zircon, possibly resulting from interaction with late-stage, U-rich hydrothermal fluids. Metasomatic textures and association with abundant, Cl-rich nanoprecipitates are typical features of fertile granites. The association between Fe and Cl recorded in metasomatized zircon during magmatic-to-hydrothermal transition is a diagnostic signature of fertile IOCG systems.
Figure 17. (A) Plot of Fe2O3 (wt.%) versus Dc (×1018 α-decays/g) for all zircon samples as marked and using EPMA data from Supplementary Materials B, Table S4 (four analyses with Dc > 25 were excluded). A sigmoid through the data was fit by non-linear least squares regression using the software “R”. Although the first percolation point (FPP) for zircon was marked at Dc = 2 (as defined by Salje et al. [48]), a sharp increase in wt.% Fe2O3 was seen only at Dc = ~3. The FPP in Olympic Dam zircon was only reached at least 100 Ma after crystallization, using the dose concentration calculation of Meldrum et al. [58]. Also, note the scatter in the dataset concordant with the weak correlation between the highest measured non-formula element (Fe) and Dc. (B) Schematic diagram illustrating evolution of magmatic zircon from granites in the Olympic Cu–Au Province. Colored areas in (A) correspond to the main evolutionary stages in (B). Note that Cape Donington zircon, with the lowest non-formula elements, plots in a magmatic domain that straddles the FPP line. In contrast, analyses with high U concentration corresponded to intensively altered stages (green), featuring superimposed textures with all the characteristic of metamict zircon, possibly resulting from interaction with late-stage, U-rich hydrothermal fluids. Metasomatic textures and association with abundant, Cl-rich nanoprecipitates are typical features of fertile granites. The association between Fe and Cl recorded in metasomatized zircon during magmatic-to-hydrothermal transition is a diagnostic signature of fertile IOCG systems.
Minerals 09 00364 g017
Table 1. Summary of the samples and foils investigated in this study.
Table 1. Summary of the samples and foils investigated in this study.
RegionSample IDGranite Petrography and Age(s), ReferencesLocality
(Drill Hole)
Micron-scale AnalysisNanoscale Analysis: HAADF STEM Imaging and EDS Spot/MappingFigures Showing Zircon Images and Maps
Analytical MethodsFoil No. (#)Key Zircon FeaturesMain Zone Axes Imaged (See Figures 12 and 13)
LA-ICP-MSEPMA MapNF Element Zoning
(STEM Maps)
Mottling/
CHZ-Nano-precipitates
Others
U–Pb DatingSpot AnalysisMap
Olympic Cu–Au ProvinceLCD47HS-Roxby Downs granite
[14,40,42], this study
Distal from OD
(RD2488)
x #1Fe, Ca, ClXscrew dislocation   ( 1 1 ¯ 1 ] )   Figures 2, 8, 10, 11, 14, 15, Supplementary C, Figures S1 and S3
#2 *Fe, Ca, Al, <Cl, YxXtm veinlet, pores ( 001 ] )   Figure 2, Supplementary Figure S3
LCD13Proximal to OD (RD2786A)xx Fe, Cl#3Fe, Cl, <CaxHm inclusions ( 001 ] )   Figures 3, 5, 9, Supplementary Figure S3
LCD17DS granite
[42,43], this study
Wirrda Well (WW46)xx Fe, Cl#4Fe, CaxXtm veinlet, screw dislocation ( 001 ] ) Figures 3, 5, 6, 15, Supplementary Figure S3
#5 **Not mappedxFe detected from spot analysis [ 100 ] ; [ 1 1 ¯ 0 ]   Figures 3 and 5
Eyre PeninsulaMB158HS-Charleston granite
[36,44]
South of Middleback Ranges #6Fe, Ca, Ti, Al, Y, Th 2-D screw dislocation, fracture
U-NP
[ 001 ] Figures 7, 15, Supplementary Figure S1
PL8DS granite
[44]
Cape Donington xx #7None—primary magmatic stretching defects [ 100 ] ; [ 1 1 ¯ 0 ]   Figures 2, 15, Supplementary Figure S2
All samples were characterized by EPMA spot analysis. Sample MB158 for Charleston granite was not analyzed by LA-ICP-MS due to the small size of zircon grains, and particularly the grain from which Foil #6 was extracted (see Figure S1). U–Pb dating was carried out for two samples (including the grains from which the foils were extracted) from locations that were not previously dated. * Foil #2 differs from Foil #1 in containing secondary rather than primary textures. ** Foils #4 and #5 were cut perpendicular and parallel to grain elongation (c-axis), respectively. Upper case X represents the most intense mottling/nanoprecipitates. Abbreviations: CHZ—chloro–hydroxy–zircon; NF—non-formula elements; Hm—hematite; Xtm—xenotime; OD—Olympic Dam; HS—Hiltaba Suite; DS—Donington Suite; U-NP—uranium-bearing nanoparticles.
Table 2. EPMA zircon spot analyses. Gr1: grain 1; Gr2; grain2.
Table 2. EPMA zircon spot analyses. Gr1: grain 1; Gr2; grain2.
LCD13 (OD Proximal)Na2OMgOAl2O3SiO2P2O5K2OCaOTiO2MnOFe2O3Y2O3ZrO2HfO2PbOThO2UO2ClSumH2O *Sum *
Bright 0.0460.016<mdl32.570.0390.0110.0760.010.0350.8220.02363.8342.3550.036<mdl0.1630.116100.600.11100.69
Bright 0.030.011<mdl32.1650.064<mdl0.0420.010.0370.460.14564.1041.4540.0890.080.3320.04899.6400.2399.86
Meann = 20.0380.014<mdl32.3680.052<mdl0.0590.010.0360.6410.08463.9691.9050.063<mdl0.2480.082100.120.17100.275
Dark 0.0830.0470.01132.2910.0630.0110.1530.0150.0290.9720.08562.5071.5320.0320.0780.1780.11998.5990.00298.574
Dark 0.1230.0280.4131.6370.5650.0310.1380.0360.1171.1521.05160.7142.8180.0830.090.8880.075101.2540.197101.434
Dark 0.1970.0190.18631.0560.390.0170.2110.0240.111.3860.86859.6132.4850.0970.0730.8270.06698.8420.46599.292
Dark 0.0580.0220.00832.460.1450.0070.1230.0120.0490.90.48462.7491.1480.0480.1960.250.15699.5470.04999.561
Dark <mdl0.0210.02132.2180.0820.0110.1850.0080.0710.9560.1263.2051.3850.0480.1220.210.37999.7050.14999.769
Dark 0.0670.0130.00731.6010.0340.0160.0850.0080.0680.9120.02962.2042.8550.048<mdl0.2030.16398.7800.35399.096
Dark 0.0640.0210.05631.2390.0660.010.3810.0070.1141.1610.16660.5572.0490.070.0920.7420.08697.6750.43898.052
Meann = 70.0990.0240.131.7860.1920.0150.1820.0160.081.0630.461.652.0390.0610.1090.4710.14999.200.23699.397
LCD47 (OD Distal)Na2OMgOAl2O3SiO2P2O5K2OCaOTiO2MnOFe2O3Y2O3ZrO2HfO2PbOThO2UO2ClSumH2O *Sum *
Bright (Gr2) <mdl0.0170.13231.6450.1130.0060.0320.0430.0590.2270.31767.0911.3790.030.1030.1450.008101.9631.194103.156
Bright (Gr2) 0.0740.0230.12231.2940.0720.0080.5350.0240.1531.8070.21161.7551.3990.0420.1170.2190.07698.3990.87799.259
Meann = 2<mdl0.020.12731.470.0930.0070.2840.0340.1061.0170.26464.4231.3890.0360.110.1820.042100.1811.036101.208
Dark (Gr2) 0.0360.0270.1830.9650.1150.0060.4270.0090.1431.440.16761.611.4020.030.0510.1110.14797.2330.65197.851
Dark (Gr2) 0.0460.030.19631.2730.0540.0060.4350.0380.1741.3330.04361.6921.7360.0370.0230.120.08997.5650.46798.012
Dark (Gr2) 0.0590.0330.83830.8270.1360.0120.7920.0180.2621.9480.55358.4561.225<mdl0.0890.1250.11196.5200.08696.532
Meann = 30.0470.030.40531.0220.1020.0080.5510.0220.1931.5740.25460.5861.4540.0340.0540.1190.11697.1060.40197.465
Bright (Gr1) <mdl0.0110.00933.0170.11<mdl0.010.011<mdl0.910.11665.9511.2630.0360.0360.0480.008101.6410.143101.783
Bright (Gr1) 0.0530.0230.0930.8460.0730.0060.47<mdl0.1381.1770.11361.4421.6180.040.0580.1880.09996.7400.64597.335
Bright (Gr1) 0.1110.0320.14331.5610.2660.0110.3290.0070.2021.3480.48660.9031.6020.0420.0420.2860.06698.0950.14998.193
Bright (Gr1) <mdl0.0160.03431.7910.050.0060.0150.0370.0120.4960.01765.211.2470.030.0170.051<mdl99.2140.46599.679
Bright (Gr1) <mdl0.0120.03532.3760.086<mdl0.0240.0070.0190.3810.05765.0941.4160.0250.0240.0370.00999.7800.07299.850
Bright (Gr1) <mdl0.40.44732.7140.1180.0060.030.0350.030.690.16565.9021.4770.0410.0270.0610.012102.4670.375102.839
Bright (Gr1) <mdl0.020.132.1990.108<mdl0.0330.0170.060.40.1565.1791.3940.0290.0410.0520.012100.0990.253100.349
Meann = 5<mdl0.0730.12332.0720.1160.0070.130.0190.0770.7720.15864.241.4310.0350.0350.1030.03499.7190.3100.004
Dark (Gr1) 0.0340.0290.46730.8970.1290.0140.4290.0130.2141.450.76460.6991.3580.0370.0740.1470.0997.6730.64898.301
Dark (Gr1) 0.0380.0190.02731.1190.0710.0070.2440.0210.1321.2390.0662.2061.4140.0230.0350.0820.09897.2330.61397.823
Dark (Gr1) 0.0230.0270.21431.3590.0720.0060.4550.0480.1631.1750.15961.7161.4480.0290.0530.1160.12297.6600.38698.018
Dark (Gr1) 0.0280.020.22831.430.0670.0080.4340.010.1451.9090.13861.7361.2320.0190.0550.0940.12598.0300.5898.582
Dark (Gr1) 0.0320.030.21631.4540.0640.0080.3740.0770.1861.6660.08161.8941.4050.0170.0350.1050.08598.0140.46298.457
Dark (Gr1) 0.0670.0210.15931.7810.069<mdl0.4420.0060.1711.4130.08762.0241.2580.0250.0420.0850.06598.0420.23298.233
Dark (Gr1) 0.0980.0280.26231.6350.2880.0140.5120.0290.1891.740.93760.7341.3490.0320.1010.2760.06799.3570.46899.748
Dark (Gr1) 0.060.0390.33530.7470.0610.0090.4780.0830.1911.7740.09463.4461.4760.0380.030.1340.06999.3111.397100.692
Dark (Gr1) 0.0940.0230.07330.4810.12<mdl0.280.0080.1430.9680.25262.0811.270.0350.0590.0820.12496.4520.84397.267
Meann = 90.0530.0260.2231.2110.1050.0090.4050.0330.171.4820.28661.8371.3570.0280.0540.1250.09497.9750.62598.569
MB158 (Charleston)Na2OMgOAl2O3SiO2P2O5K2OCaOTiO2MnOFe2O3Y2O3ZrO2HfO2PbOThO2UO2ClSumH2OSum *
Bright <mdl<mdl0.03633.0010.0990.0080.1750.040.1291.2920.22465.1851.3530.0510.1480.130.016102.3810.032102.41
Bright 0.0320.016<mdl31.7080.092<mdl0.6740.0520.1882.0820.14763.2831.320.0490.0810.120.046100.1530.005100.147
Bright 0.0580.0070.0231.930.0870.0070.660.060.0871.7040.10863.2291.3930.0460.0420.1070.03499.980099.973
Bright 0.0480.0530.16432.1520.0960.0050.9330.0310.1441.2240.27362.3781.5360.0380.0380.1540.06399.8420.02299.85
Bright 0.0330.0120.01332.3990.0830.0070.2530.0270.1751.020.11464.1561.5240.0360.0380.1480.016100.3660.005100.368
Meann = 50.0430.0220.05832.2380.0910.0070.5390.0420.1451.4640.17363.6461.4250.0440.0690.1320.035100.5440.013100.55
Dark 0.0440.0190.03531.6460.0950.0071.1150.0290.1331.7840.20662.4981.210.0570.1240.1430.0699.6570.01699.659
Dark 0.040.0160.07332.0770.124<mdl1.3990.0550.1532.2030.35561.5281.160.0450.2320.1470.064100.1750.029100.189
Dark <mdl0.0160.01331.920.082<mdl0.6770.0640.1762.2110.1962.891.3750.0320.1210.160.047100.2930.024100.306
Dark 0.0620.0120.0131.3720.07<mdl0.9810.0140.2191.3060.10261.8791.3970.0360.0460.150.05298.0480098.037
Dark 0.0330.0170.08331.4030.0810.0071.0330.0140.1631.2810.11361.8411.523<mdl0.0840.1540.04998.3170.01498.32
Meann = 50.0450.0160.04331.6840.09<mdl1.0410.0350.1691.7570.19362.1271.3330.0430.1210.1510.05499.2980.01799.302
PL8 (Cape. Don)Na2OMgOAl2O3SiO2P2O5K2OCaOTiO2MnOFe2O3Y2O3ZrO2HfO2PbOThO2UO2ClSumH2O *Sum *
Bright <mdl0.008<mdl32.5180.107<mdl<mdl0.005<mdl<mdl0.1565.7161.2330.0230.0270.035<mdl100.060.012100.071
Bright <mdl0.010.00832.8010.092<mdl<mdl0.0070.0120.0260.03965.941.3920.0190.0150.021<mdl100.5890.003100.592
Bright <mdl<mdl<mdl32.6460.097<mdl0.01<mdl<mdl<mdl0.07565.8911.2750.0380.030.064<mdl100.3960.009100.405
Bright <mdl<mdl<mdl32.7590.074<mdl<mdl0.012<mdl<mdl0.03465.8631.3340.0260.020.045<mdl100.2490100.249
Meann = 4<mdl<mdl<mdl32.6810.093<mdl<mdl0.008<mdl<mdl0.07565.8531.3090.0270.0230.041<mdl100.3240.006100.329
Dark <mdl0.009<mdl31.9970.077<mdl<mdl0.01<mdl0.0320.02865.141.3260.03<mdl0.012<mdl98.9070.14299.049
Dark <mdl0.008<mdl32.4230.09<mdl<mdl0.006<mdl<mdl0.03865.8061.41<mdl<mdl0.016<mdl100.1140.092100.206
Dark <mdl<mdl<mdl32.0050.065<mdl<mdl0.008<mdl0.0410.02565.5021.3260.023<mdl0.015<mdl99.2330.23999.472
Meann = 3<mdl0.009<mdl32.1420.077<mdl<mdl0.008<mdl0.0370.0365.4831.3540.027<mdl0.014<mdl99.4180.15899.576
LCD17 (Wirrda Well)Na2OMgOAl2O3SiO2P2O5K2OCaOTiO2MnOFe2O3Y2O3ZrO2HfO2PbOThO2UO2ClSumH2O *Sum *
Bright <mdl0.011<mdl32.7580.104<mdl0.0180.006<mdl0.1970.09364.5291.2880.0310.0260.051<mdl99.5480.02199.568
Bright <mdl0.013<mdl32.770.172<mdl<mdl0.009<mdl0.1610.10864.5981.2650.0280.0270.048<mdl99.5150.00699.521
Bright <mdl0.012<mdl32.8740.125<mdl<mdl0.008<mdl0.1280.09265.0571.3430.0340.0270.0450.00999.9700.00799.975
Bright <mdl<mdl<mdl32.7890.121<mdl<mdl0.01<mdl0.2360.08365.1591.3890.0340.0250.047<mdl100.1790.009100.188
Bright <mdl0.015<mdl32.7350.1290.0070.0120.006<mdl0.1050.08965.1591.310.0190.0310.059<mdl100.0280.013100.041
Meann = 5<mdl0.013<mdl32.7850.13<mdl<mdl0.008<mdl0.1650.09364.91.3190.0290.0270.05<mdl99.8480.01199.859
Dark 0.0270.0310.2231.3170.1880.0080.260.0140.0851.3380.36862.1091.2820.0380.0660.1080.08198.1930.49898.673
Dark 0.0360.017<mdl32.7240.1380.0060.0460.0060.0380.4710.20263.71.383<mdl0.050.0770.05299.3830.01699.388
Dark 0.0350.013<mdl32.4210.1390.0070.0290.0080.0420.7430.11264.2421.3340.030.0330.0810.02799.7750.06299.831
Dark <mdl0.010<mdl33.0890.0840.0070.010.009<mdl0.2720.05365.8461.2450.026<mdl0.037<mdl100.8930.006100.9
Dark <mdl0.008<mdl32.2590.152<mdl0.0110.0060.0120.1620.11664.6841.2040.0340.0160.054<mdl99.0620.01199.072
Meann = 50.0330.016<mdl32.3620.1400.0070.0710.0090.0440.5970.1764.1161.290.0320.0410.0710.05399.4610.11999.573
average mdl0.0230.0070.0060.0210.0100.0050.0090.0040.0110.0240.0080.0490.1010.0140.0100.0110.007
Sum *: Sum of element totals including calculated H2O (H2O*) and corrected for F and Cl. Sc, Nb, and F were measured but are generally below the mdl. Full values for these elements and REY elements are included in Supplementary Materials B, Table S4. Formulae calculation values and α-decays/mg values are in Supplementary Materials B, Table S4.
Table 3. Mean zircon compositions in bright and dark bands in each studied sample.
Table 3. Mean zircon compositions in bright and dark bands in each studied sample.
Lithology and LocationSample IDZircon Composition
Cape Donington, Donington Suite PL8Bright: (Zr*0.994REE*0.001)∑=0.995(Si1.002P0.002O4);
Dark: (Zr*0.996Fe0.001REE*0.003)∑=1(Si0.986P0.002O3.95)(OH)0.049.
Wirrda Well, Donington Suite LCD17Bright: Zr*0.979Fe0.004REE*0.006NF0.001)∑=0.989(Si1.006P0.005O3.966)(OH)0.001;
Dark: i0.002(Zr*0.943Fe0.031REE*0.014NF0.013)∑=1(Si0.962Al0.008P0.005O3.862)(OH,Cl)0.106.
OD-distal Hiltaba Suite LCD47Bright: (Zr*0.986Fe0.009REE*0.003NF0.002)∑=1(Si0.993Al0.001P0.002O3.978)(OH,Cl)0.015;
Dark: i0.004(Zr*0.934Fe0.04REE*0.004NF0.022)∑=1(Si0.915Al0.012P0.002Ti0.002O3.672)(OH,Cl)0.281.
OD-proximal Hiltaba Suite LCD13Bright: i0.003(Zr*0.97Fe0.019REE*0.007NF0.004)∑=1(Si0.993P0.001O3.959)(OH,Cl)0.028;
Dark: i0.004(Zr*0.942Fe0.027REE*0.013NF0.018)∑=1(Si0.972Al0.002P0.002O3.864)(OH,F,Cl)0.105.
Charleston Granite, Hiltaba SuiteMB158Bright:(Zr*0.962Fe0.019REE*0.009NF0.009)∑=0.999(Si0.985Al0.001P0.003Ti0.001Fe0.01O3.969)(OH,Cl)0.007;
Dark: i0.003(Zr*0.941Fe0.009REE*0.009NF0.041)∑=1(Si0.964Al0.001P0.003Fe0.031O3.933)(OH,Cl)0.006.
Note: All granite samples, with the exception of PL8, show various degrees of alteration (see above). The presented compositions are for grains that were micro-sampled for S/TEM analysis. * To obtain a good stoichiometry between the A and B sites, we used a slightly modified procedure for calculating the formulae of the Charleston zircon, in which Fe was allocated in both A and B sites.

Share and Cite

MDPI and ACS Style

Courtney-Davies, L.; Ciobanu, C.L.; Verdugo-Ihl, M.R.; Slattery, A.; Cook, N.J.; Dmitrijeva, M.; Keyser, W.; Wade, B.P.; Domnick, U.I.; Ehrig, K.; et al. Zircon at the Nanoscale Records Metasomatic Processes Leading to Large Magmatic–Hydrothermal Ore Systems. Minerals 2019, 9, 364. https://doi.org/10.3390/min9060364

AMA Style

Courtney-Davies L, Ciobanu CL, Verdugo-Ihl MR, Slattery A, Cook NJ, Dmitrijeva M, Keyser W, Wade BP, Domnick UI, Ehrig K, et al. Zircon at the Nanoscale Records Metasomatic Processes Leading to Large Magmatic–Hydrothermal Ore Systems. Minerals. 2019; 9(6):364. https://doi.org/10.3390/min9060364

Chicago/Turabian Style

Courtney-Davies, Liam, Cristiana L. Ciobanu, Max R. Verdugo-Ihl, Ashley Slattery, Nigel J. Cook, Marija Dmitrijeva, William Keyser, Benjamin P. Wade, Urs I. Domnick, Kathy Ehrig, and et al. 2019. "Zircon at the Nanoscale Records Metasomatic Processes Leading to Large Magmatic–Hydrothermal Ore Systems" Minerals 9, no. 6: 364. https://doi.org/10.3390/min9060364

APA Style

Courtney-Davies, L., Ciobanu, C. L., Verdugo-Ihl, M. R., Slattery, A., Cook, N. J., Dmitrijeva, M., Keyser, W., Wade, B. P., Domnick, U. I., Ehrig, K., Xu, J., & Kontonikas-Charos, A. (2019). Zircon at the Nanoscale Records Metasomatic Processes Leading to Large Magmatic–Hydrothermal Ore Systems. Minerals, 9(6), 364. https://doi.org/10.3390/min9060364

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop