Zircon chronochemistry of high heat-producing granites in Queensland and Europe

High heat-producing granites (HHPGs) are reservoir rocks for enhanced geothermal systems (EGS), yet the origins of their anomalous chemistry remain poorly understood. To gain a better understanding of the characteristic distribution of elemental depletions and enrichments (focussing on U, Th & K) within granite suites of different heritage and tectonic setting, and the processes that lead to these enrichments, we are undertaking a systematic accessory-mineral chronochemical study of two suites of S- and I-type granites in northern Queensland, as well as two archetypal HHPGs in Cornwall, England (S-type) and Soultz-sous- Forêts, France (I-type). Novel zircon LA-ICP-MS chronochemical methods will later be underpinned by a systematic petrographic, scanning electron microscope (SEM), and electron microprobe (EPMA) study of all the REE-Y-Th-U-rich accessory minerals to fully characterise how the composition, textural distributions and associations change with rock chemistry between and among the suites. Preliminary results indicate that zircons with inherited ages do not have anomalously high U (>1000 ppm) & Th (>400 ppm) values (Ahrens, 1965). Instead, enrichment in these HPE is seen in zircons dated to around the time of magmatic emplacement. These results indicate that enrichment arose primarily through fractional crystallisation of the granitic magmas. Our results support the suggestion that a source pre-enriched in the HPEs does not appear to be fundamental for the formation of all HHPGs. Instead fractional crystallisation processes, and the accessory minerals formed in magmas of differing initial compositions, are the key controls on the levels of enrichment observed (e.g. Champion & Chappell, 1992; Chappell & Hine, 2006). One implication is that the most fractionated granites may not be the most enriched in the HPEs and therefore prospective to future EGS development.

Two cycles of voluminous pyroclastic volcanism and sedimentation related to episodic granite emplacement during the late Archean : Eastern Yilgarn Craton, Western Australia

The thick piles of late-Archean volcaniclastic sedimentary successions that overlie the voluminous greenstone units of the eastern Yilgarn Craton, Western Australia, record the important transition from the cessation in mafic-ultramafic volcanism to cratonisation between about 2690 and 2655 Ma. Unfortunately, an inability to clearly subdivide the superficially similar sedimentary successions and correlate them between the various geological terranes and domains of the eastern Yilgarn Craton has led to uncertainty about the timing and nature of the region's palaeogeographic and palaeotectonic evolution. Here, we present the results of some 2025 U–Pb laser-ablation-ICP-MS analyses and 323 Sensitive High-Resolution Ion Microprobe (SHRIMP) analyses of detrital zircons from 14 late-Archean felsic clastic successions of the eastern Yilgarn Craton, which have enabled correlation of clastic successions. The results of our data, together with those compiled from previous studies, show that the post-greenstone sedimentary successions include two major cycles that both commenced with voluminous pyroclastic volcanism and ended with widespread exhumation and erosion associated with granite emplacement. Cycle One commences with an influx of rapidly reworked feldspar-rich pyroclastic debris. These units, here-named the Early Black Flag Group, are dominated by a single population of detrital zircons with an average age of 2690–2680 Ma. Thick (up to 2 km) dolerite bodies, such as the Golden Mile Dolerite, intrude the upper parts of the Early Black Flag Group at about 2680 Ma. Incipient development of large granite domes during Cycle One created extensional basins predominantly near their southeastern and northwestern margins (e.g., St Ives, Wallaby, Kanowna Belle and Agnew), into which the Early Black Flag Group and overlying coarse mafic conglomerate facies of the Late Black Flag Group were deposited. The clast compositions and detrital-zircon ages of the late Black Flag Group detritus match closely the nearby and/or stratigraphically underlying successions, thus suggesting relatively local provenance. Cycle Two involved a similar progression to that observed in Cycle One, but the age and composition of the detritus were notably different. Deposition of rapidly reworked quartz-rich pyroclastic deposits dominated by a single detrital-zircon age population of 2670–2660 Ma heralded the beginning of Cycle Two. These coarse-grained quartz-rich units, are name here the Early Merougil Group. The mean ages of the detrital zircons from the Early Merougil Group match closely the age of the peak in high-Ca (quartz-rich) granite magmatism in the Yilgarn Craton and thus probably represent the surface expression of the same event. Successions of the Late Merougil Group are dominated by coarse felsic conglomerate with abundant volcanic quartz. Although the detrital zircons in these successions have a broad spread of age, the principal sub-populations have ages of about 2665 Ma and thus match closely those of the Early Merougil Group. These successions occur most commonly at the northwestern and southeastern margins of the granite batholiths and thus are interpreted to represent resedimented units dominted by the stratigraphically underlying packages of the Early Merougil Group. The Kurrawang Group is the youngest sedimentary units identified in this study and is dominated by polymictic conglomerate with clasts of banded iron formation (BIF), granite and quartzite near the base and quartz-rich sandstone units containing detrital zircons aged up to 3500 Ma near the top. These units record provenance from deeper and/or more-distal sources. We suggest here that the principal driver for the major episodes of volcanism, sedimentation and deformation associated with basin development was the progressive emplacement of large granite batholiths. This interpretation has important implication for palaeogeographic and palaeotectonic evolution of all late-Archean terranes around the world.

In search of a hidden long-term isolated sub-chondritic 142Nd/144Nd reservoir in the deep mantle : implications for the Nd isotope systematics of the Earth

Here we search for evidence of the existence of a sub-chondritic 142Nd/144Nd reservoir that balances the Nd isotope chemistry of the Earth relative to chondrites. If present, it may reside in the source region of deeply sourced mantle plume material. We suggest that lavas from Hawai’i with coupled elevations in 186Os/188Os and 187Os/188Os, from Iceland that represent mixing of upper mantle and lower mantle components, and from Gough with sub-chondritic 143Nd/144Nd and high 207Pb/206Pb, are favorable samples that could reflect mantle sources that have interacted with an Early-Enriched Reservoir (EER) with sub-chondritic 142Nd/144Nd. High-precision Nd isotope analyses of basalts from Hawai’i, Iceland and Gough demonstrate no discernable 142Nd/144Nd deviation from terrestrial standards. These data are consistent with previous high-precision Nd isotope analysis of recent mantle-derived samples and demonstrate that no mantle-derived material to date provides evidence for the existence of an EER in the mantle. We then evaluate mass balance in the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd. The Nd isotope systematics of EERs are modeled for different sizes and timing of formation relative to ε143Nd estimates of the reservoirs in the μ142Nd = 0 Earth, where μ142Nd is ((measured 142Nd/144Nd/terrestrial standard 142Nd/144Nd)−1 * 10−6) and the μ142Nd = 0 Earth is the proportion of the silicate Earth with 142Nd/144Nd indistinguishable from the terrestrial standard. The models indicate that it is not possible to balance the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd unless the μ142Nd = 0 Earth has a ε143Nd within error of the present-day Depleted Mid-ocean ridge basalt Mantle source (DMM). The 4567 Myr age 142Nd–143Nd isochron for the Earth intersects μ142Nd = 0 at ε143Nd of +8 ± 2 providing a minimum ε143Nd for the μ142Nd = 0 Earth. The high ε143Nd of the μ142Nd = 0 Earth is confirmed by the Nd isotope systematics of Archean mantle-derived rocks that consistently have positive ε143Nd. If the EER formed early after solar system formation (0–70 Ma) continental crust and DMM can be complementary reservoirs with respect to Nd isotopes, with no requirement for significant additional reservoirs. If the EER formed after 70 Ma then the μ142Nd = 0 Earth must have a bulk ε143Nd more radiogenic than DMM and additional high ε143Nd material is required to balance the Nd isotope systematics of the Earth.

Were intercalated komatiites and dacites at the Black Swan nickel sulphide mine, Yilgarn Craton, Western Australia, emplaced as extrusive lavas or intrusive bodies? The significance of breccia textures and contact relationships

Intercalated Archean komatiites and dacites sit above a thick footwall dacite unit in the host rock succession at the Black Swan Nickel Mine, north of Kalgoorlie in the Yilgarn Craton, Western Australia. Both lithofacies occur in units that vary in scale from laterally extensive at the scale of the mine lease to localized, thin, irregular bodies, from > 100 m thick to only centimetres thick. Some dacites are only slightly altered and deformed, and are interpreted to post-date major deformation and alteration (late porphyries). However, the majority of the dacites display evidence of deformation, especially at contacts, and metamorphism, varying from silicification and chlorite alteration at contacts to pervasive low grade regional metamorphic alteration represented by common assemblages of chlorite, sericite and albite. Texturally, the dacites vary from entirely massive and coherent to partially brecciated to totally brecciated. Strangely, some dacites are coherent at the margins and brecciated internally. Breccia textures vary from cryptically defined, to blocky, closely packed, in situ jig-saw fit textures with secondary minerals in fractures between clasts, to more apparent matrix rich textures with round clast forms, giving apparent conglomerate textures. Some clast zones have multi-coloured clasts, giving the impression of varied provenance. Strangely however, all these textural variants have gradational relationships with each other, and no bedding or depositional structures are present. This indicates that all textures have an in situ origin. The komatiites are generally altered and pervasively carbonate veined. Preservation of original textures is patchy and local, but includes coarse adcumulate, mesocumulate, orthocumulate, crescumulate-harrisite and occasionally spinifex textures. Where original contacts between komatiites and dacites are preserved intact (i.e. not sheared or overprinted by alteration), the komatiites have chilled margins, whereas the dacites do not. The margins of the dacites are commonly silicified, and inclusions of dacite occur in komatiite, even at the top contacts of komatiite units, but komatiite clasts do not occur in the dacites. The komatiites therefore were emplaced as sills into the dacites, and the intercalated relationships are interpreted as intrusive. The brecciation and alteration in the dacites are interpreted as being largely due to hydraulic fracturing and alteration induced by contact metamorphic effects and hydrothermal alteration deriving from the intrusion of komatiites into the felsic pile. The absence of autobreccia and hyaloclastite textures in the dacites suggest that they were emplaced as an earlier intrusive (sill?) complex at a high level in the crust.

A new database compilation of whole-rock chemical and geochronological data of igneous rocks in Queensland : a new resource for HDR geothermal resource exploration

The Geothermal industry in Australia and Queensland is in its infancy and for hot dry rock (HDR) geothermal energy, it is very much in the target identification and resource definition stages. As a key effort to assist the geothermal industry and exploration for HDR in Queensland, we are developing a comprehensive and new integrated geochemical and geochronological database on igneous rocks. To date, around 18,000 igneous rocks have been analysed across Queensland for chemical and/or age information. However, these data currently reside in a number of disparate datasets (e.g., Ozchron, Champion et al., 2007, Geological Survey of Queensland, journal publications, and unpublished university theses). The goal of this project is to collate and integrate these data on Queensland igneous rocks to improve our understanding of high heat producing granites in Queensland, in terms of their distribution (particularly in the subsurface), dimensions, ages, and controlling factors in their genesis.

Evaluating the role of deep granitic rocks in generating anomalous temperatures in SW-Queensland

Across central Australia and south-west Queensland, a large (~800,000km2) subsurface temperature anomaly occurs (Figure 1). Temperatures are interpreted to be greater than 235°C at 5km depth, ca. 85°C higher than the average geothermal gradient for the upper continental crust (Chopra & Holgate, 2005; Holgate & Gerner, 2011). This anomaly has driven the development of Engineered Geothermal Systems (EGS) at Innamincka, where high temperatures have been related to the radiogenic heat production of High Heat Producing Granites (HHPG) at depth, below thermally insulative sedimentary cover (Chopra & Holgate, 2005; Draper & D’Arcy, 2006; Meixner & Holgate, 2009). To evaluate the role of granitic rocks at depth in generating the broader temperature anomaly in SW-Queensland, we sampled 25 granitic rocks from basement intervals of petroleum drill cores below thermal insulative cover along two transects (WNW–ESE and NNE–SSW — Figure 1) and performed a multidisciplinary study involving petrography, whole-rock chemistry, zircon dating and thermal conductivity measurements.

QGECE research on delineation of Australian geothermal resources

The geology/reservoir program of the Queensland Geothermal Energy Centre of Excellence (QGECE) has the mission to improve the existing knowledge and develop new innovative scientific approaches for the identification of geothermal resources in Australia, with a particular focus on Queensland. Specifically, the QGECE geology/reservoir program is currently (1) producing a comprehensive geochemical dataset for high heat producing rocks, (2) conducting detailed mineralogical and geochronological studies of granites and hydrothermal alteration minerals, and ; (3) investigating the Cooper Basin representing a superb natural laboratory for understanding of radiogenic heat enrichment process and possible involvement of mantle heat flow. Seven research projects have been established, which are being conducted largely as PhD studies. In the preliminary studies, high quality and valuable results were obtained to address the research topics of understanding the causes and timing of heat producing element enrichment.

Heat-producing crust regulation of subsurface temperatures : a stochastic model re-evaluation of the geothermal potential in western Queensland, Australia

A large subsurface, elevated temperature anomaly is well documented in Central Australia. High Heat Producing Granites (HHPGs) intersected by drilling at Innamincka are often assumed to be the dominant cause of the elevated subsurface temperatures, although their presence in other parts of the temperature anomaly has not been confirmed. Geological controls on the temperature anomaly remain poorly understood. Additionally, methods previously used to predict temperature at 5 km depth in this area are simplistic and possibly do not give an accurate representation of the true distribution and magnitude of the temperature anomaly. Here we re-evaluate the geological controls on geothermal potential in the Queensland part of the temperature anomaly using a stochastic thermal model. The results illustrate that the temperature distribution is most sensitive to the thermal conductivity structure of the top 5 km. Furthermore, the results indicate the presence of silicic crust enriched in heat producing elements between and 40 km.

Timing and source constraints on the relationship between mafic and felsic intrusions in the Emeiashan large igneous province

Several I- and A-type granite, syenite plutons and spatially associated, giant Fe–Ti–V deposit-bearing mafic ultramafic layered intrusions occur in the Pan–Xi(Panzhihua–Xichang) area within the inner zone of the Emeishan large igneous province (ELIP). These complexes are interpreted to be related to the Emeishan mantle plume. We present LA-ICP-MS and SIMS zircon U–Pb ages and Hf–Nd isotopic compositions for the gabbros, syenites and granites from these complexes. The dating shows that the age of the felsic intrusive magmatism (256.2 ± 3.0–259.8 ± 1.6 Ma) is indistinguishable from that of the mafic intrusive magmatism (255.4 ± 3.1–259.5 ± 2.7 Ma) and represents the final phase of a continuous magmatic episode that lasted no more than 10 Myr. The upper gabbros in the mafic–ultramafic intrusions are generally more isotopically enriched (lower eNd and eHf) than the middle and lower gabbros, suggesting that the upper gabbros have experienced a higher level of crustal contamination than the lower gabbros. The significantly positive eHf(t) values of the A-type granites and syenites (+4.9 to +10.8) are higher than those of the upper gabbros of the associated mafic intrusion, which shows that they cannot be derived by fractional crystallization of these bodies. They are however identical to those of the mafic enclaves (+7.0 to +11.4) and middle and lower gabbros, implying that they are cogenetic. We suggest that they were generated by fractionation of large-volume, plume-related basaltic magmas that ponded deep in the crust. The deep-seated magma chamber erupted in two stages: the first near a density minimum in the basaltic fractionation trend and the second during the final stage of fractionation when the magma was a low density Fe-poor, Si-rich felsic magma. The basaltic magmas emplaced in the shallowlevel magma chambers differentiated to form mafic–ultramafic layered intrusions accompanied by a small amount of crustal assimilation through roof melting. Evolved A-type granites (synenites and syenodiorites) were produced dominantly by crystallization in the deep crustal magma chamber. In contrast, the I-type granites have negative eNd(t) [-6.3 to -7.5] and eHf(t) [-1.3 to -6.7] values, with the Nd model ages (T Nd DM2) of 1.63-1.67 Ga and Hf model ages (T Hf DM2) of 1.56-1.58 Ga, suggesting that they were mainly derived from partial melting of Mesoproterozoic crust. In combination with previous studies, this study also shows that plume activity not only gave rise to reworking of ancient crust, but also significant growth of juvenile crust in the center of the ELIP.

A new period of volcanogenic massive sulfide formation in the Yilgarn: A volcanological study of the ca 2.76 Ga Hollandaire VMS deposit, Yilgarn Craton, Western Australia

The Archean Hollandaire volcanogenic massive sulfide deposit is a felsic–siliciclastic VMS deposit located in the Murchison Domain of the Youanmi Terrane, Yilgarn Craton, Western Australia. It is hosted in a succession of turbidites, mudstones and coherent rhyodacite sills and has been metamorphosed to upper greenschist/lower amphibolite facies and includes a pervasive S1 deformational fabric. The coherent rhyodacitic sills are interpreted as syndepositional based on geochemical similarities with well-known VMS-associated felsic rocks and similar foliations to the metasediments. We offer several explanations for the absence of textural evidence (e.g. breccias) for syn-depositional origins: 1) the subaqueous sediments were dehydrated by long-lived magmatism such that no pore-water remained to drive quench fragmentation; 2) pore-space occlusion by burial and/or, 3) alteration overprinting and obscuring of primary breccias at contact margins. Mineralisation occurs by sub-seafloor replacement of original host rocks in two ore bodies, Hollandaire Main (~125 x >500 m and ~8 m thick) and Hollandaire West (~100 x 470 m and ~5 m thick), and occurs in three main textural styles, massive sulfides, which are exclusively hosted in turbidites and mudstones, and stringer and disseminated sulfides, which are also hosted in coherent rhyodacite. Most sulfides have textures consistent with remobilisation and recrystallisation. Hydrothermal metamorphism has altered the hangingwall and footwall to similar degrees, with significant gains in Mg, Mn and K and losses in Na, Ca and Sr. Garnet and staurolite porphyryoblasts also exhibit a footprint around mineralisation, extending up to 30 m both above and below the ore zone. High precision thermal ionisation mass spectrometry of zircons extracted from the coherent rhyodacite yield an age of 2759.5 ± 0.9 Ma, which along with geochemical comparisons, places the succession within the 2760–2735 Ma Greensleeves Formation of the Polelle Group of the Murchison Supergroup. Geochemical and geochronological evidence link the coherent rhyodacite sills to the Peter Well Granodiorite pluton ~2 km to the W, which acted as the heat engine driving hydrothermal circulation during VMS mineralisation. This study highlights the importance of both: detailed physical volcanological studies from which an accurate assessment of timing relationships, particularly the possibility of intrusions dismembering ore horizons, can be made; and identifying synvolcanic plutons and other similar suites, for VMS exploration targets in the Youanmi Terrane and worldwide.

The largest Au deposits in the St Ives Goldfield (Yilgarn Craton, Western Australia) may be located in a major Neoarchean volcano-sedimentary depo-centre

The largest Neoarchean gold deposits in the world-class St Ives Goldfield, Western Australia, occur in an area known as the Argo-Junction region (e.g. Junction, Argo and Athena). Why this region is so well endowed with large deposits compared with other parts of the St Ives Goldfield is currently unclear, because gold deposits at St Ives are hosted by a variety of lithologic units and were formed during at least three different deformational events. This paper presents an investigation into the stratigraphic architecture and evolution of the Argo-Junction region to assess its implications for gold metallogenesis. The results show that the region's stratigraphy may be subdivided into five regionally correlatable packages: mafic lavas of the Paringa Basalt; contemporaneously resedimented feldspar-rich pyroclastic debris of the Early Black Flag Group; coarse polymictic volcanic debris of the Late Black Flag Group; thick piles of mafic lavas and sub-volcanic sills of the Athena Basalt and Condenser Dolerite; and the voluminous quartz-rich sedimentary successions of the Early Merougil Group. In the Argo-Junction region, these units have an interpreted maximum thickness of at least 7,130 m, and thus represent an unusually thick accumulation of the Neoarchean volcano-sedimentary successions. It is postulated that major basin-forming structures that were active during deposition and emplacement of the voluminous successions later acted as important conduits during mineralisation. Therefore, a correlation exists between the location of the largest gold deposits in the St Ives Goldfield and the thickest parts of the stratigraphy. Recognition of this association has important implications for camp-scale exploration.