Deep geothermal: The ‘Moon Landing’ mission in the unconventional energy and minerals space

Deep geothermal from the hot crystalline basement has remained an unsolved frontier for the geothermal industry for the past 30 years. This poses the challenge for developing a new unconventional geomechanics approach to stimulate such reservoirs. While a number of new unconventional brittle techniques are still available to improve stimulation on short time scales, the astonishing richness of failure modes of longer time scales in hot rocks has so far been overlooked. These failure modes represent a series of microscopic processes: brittle microfracturing prevails at low temperatures and fairly high deviatoric stresses, while upon increasing temperature and decreasing applied stress or longer time scales, the failure modes switch to transgranular and intergranular creep fractures. Accordingly, fluids play an active role and create their own pathways through facilitating shear localization by a process of time-dependent dissolution and precipitation creep, rather than being a passive constituent by simply following brittle fractures that are generated inside a shear zone caused by other localization mechanisms. We lay out a new theoretical approach for the design of new strategies to utilize, enhance and maintain the natural permeability in the deeper and hotter domain of geothermal reservoirs. The advantage of the approach is that, rather than engineering an entirely new EGS reservoir, we acknowledge a suite of creep-assisted geological processes that are driven by the current tectonic stress field. Such processes are particularly supported by higher temperatures potentially allowing in the future to target commercially viable combinations of temperatures and flow rates.

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.

Regulation of carotenoid biosynthesis in banana fruit

This thesis investigated the basis for massive differences in provitamin A carotenoid content in banana fruits. Rather than gene expression levels, carotenoid storage capacity and product degradation explained much of the differences. Such information should provide important insights for future developments in the biofortification of banana. A high carotenoid-containing cultivar, 'Asupina' and a popular commercial but low carotenoid-containing cultivar, 'Cavendish' were used in the investigations.

Copula-based spatial modelling of geometallurgical variables

The most important aspect of modelling a geological variable, such as metal grade, is the spatial correlation. Spatial correlation describes the relationship between realisations of a geological variable sampled at different locations. Any method for spatially modelling such a variable should be capable of accurately estimating the true spatial correlation. Conventional kriged models are the most commonly used in mining for estimating grade or other variables at unsampled locations, and these models use the variogram or covariance function to model the spatial correlations in the process of estimation. However, this usage assumes the relationships of the observations of the variable of interest at nearby locations are only influenced by the vector distance between the locations. This means that these models assume linear spatial correlation of grade. In reality, the relationship with an observation of grade at a nearby location may be influenced by both distance between the locations and the value of the observations (ie non-linear spatial correlation, such as may exist for variables of interest in geometallurgy). Hence this may lead to inaccurate estimation of the ore reserve if a kriged model is used for estimating grade of unsampled locations when nonlinear spatial correlation is present. Copula-based methods, which are widely used in financial and actuarial modelling to quantify the non-linear dependence structures, may offer a solution. This method was introduced by Bárdossy and Li (2008) to geostatistical modelling to quantify the non-linear spatial dependence structure in a groundwater quality measurement network. Their copula-based spatial modelling is applied in this research paper to estimate the grade of 3D blocks. Furthermore, real-world mining data is used to validate this model. These copula-based grade estimates are compared with the results of conventional ordinary and lognormal kriging to present the reliability of this method.

Hydrostratigraphic and hydrochemical characterisation of aquifers, aquitards and coal seams in the Galilee and Eromanga basins, Central Queensland, Australia

Inter-aquifer mixing studies are usually made carrying out hydrochemical and isotopic techniques only. In this thesis these techniques have been integrated with three-dimensional geological modelling proving to be a better approach for inter—aquifer mixing assessment in regional areas, and also highlighting the influence of faulting in the understanding of groundwater and gas migration, which could not be possible using the two fist techniques alone. The results are of particular interest for coal seam gas basins and can even be used as exploration tools as areas of higher permeability and gas migration were identified.

Tectonometamorphic evolution of the Åreskutan Nappe-Caledonian history revealed by SIMS U–Pb zircon geochronology

Secondary ionization mass spectrometry (SIMS) U–Pb dating of zircons from the Åreskutan Nappe in the central part of the Seve Nappe Complex of western central Jämtland provides new constraints on the timing of granulite–amphibolite-facies metamorphism and tectonic stacking of the nappe during the Caledonian orogeny. Peak-temperature metamorphism in garnet migmatites is constrained to c. 442 ± 4 Ma, very similar to the ages of leucogranites at 442 ± 3 and 441 ± 4 Ma. Within a migmatitic amphibolite, felsic segregations crystallized at 436 ± 2 Ma. Pegmatites, cross-cutting the dominant Caledonian foliation in the Nappe, yield 428 ± 4 and 430 ± 3 Ma ages. The detrital zircon cores in the migmatites and leucogranites provide evidence of Late Palaeoproterozoic, Mesoproterozoic to Early Neoproterozoic source terranes for the metasedimentary rocks. The formation of the ductile and hot Seve migmatites, with their inverted metamorphism and thinning towards the hinterland, can be explained by an extrusion model in which the allochthon stayed ductile for a period of at least 10 million years during cooling from peak-temperature metamorphism early in the Silurian. In our model, Baltica–Laurentia collision occurred in the Late Ordovician–earliest Silurian, with emplacement of the nappes far on to the Baltoscandian platform during the Silurian and early Devonian, Scandian Orogeny lasting until c. 390 Ma.

Cyclic Salinity Model for deposition of large-scale coal deposits

Existing field data for Rangal coals (Late Permian) of the Bowen Basin, Queensland, Australia, are inconsistent with the depositional model generally accepted in the current geological literature to explain coal deposition. Given the apparent unsuitability of the current depositional model to the Bowen Basin coal data, a new depositional model, here named the Cyclic Salinity Model, is proposed and tested in this study.

Eruption centres of the Hamilton area of the Newer Volcanics Province, Victoria, Australia: Pinpointing volcanoes from a multifaceted approach to landform mapping

Volcanic eruption centres of the mostly 4.5 Ma-5000 BP Newer Volcanics Province in the Hamilton area of southeastern Australia were examined in detail using a multifaceted approach, including ground truthing and analysis of ArcGIS Total Magnetic Intensity and seamless geology data, NASA Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) digital elevation models and Google Earth satellite image interpretation. Sixteen eruption centres were recognised in the Hamilton area, including three previously unrecorded volcanoes-one of which, the Cas Maar, constitutes the northernmost maar-cone volcanic complex in the Western Plains subprovince. Seven previously allocated eruption centres were placed into question based on field and laboratory observations. Three phases of volcanic activity have been suggested by other authors and are interpreted to correlate with ages of >4 Ma, ca 2 Ma and <0.5 Ma, which may be further subdivided based on preservation of outcrop. Geochemical compositions of the dominantly basaltic products become increasingly alkaline and enriched in incompatible elements from Phases 1 to 2, with Phase 3 eruptions both covering the entire geochemical range and extending into increasingly enriched compositions. This research highlights the importance of a multifaceted approach to landform mapping and demonstrates that additional volcanic centres may yet be discovered in the Newer Volcanics Province

Criteria for interpreting kimberlite as coherent: Insights from the Muskox and Jericho kimberlites (Nunavut, Canada)

The Jurassic Muskox and Jericho kimberlites (Northern Slave Province, Nunavut, Canada) contain a variety of facies exhibiting different geometries, contact relationships, internal organisation, country rock abundance and olivine shapes, although many have similar matrix/groundmass mineralogies and textures. Five facies are examined that either have characteristics consistent with coherent rocks in general (i.e. intrusive and extrusive non-fragmental rocks) or are mineralogically and texturally similar to kimberlite described as coherent (or apparent coherent). Three facies are interpreted as coherent on the basis of: (1) geological setting, (2) apparent-porphyritic texture, (3) sharp contacts with fragmental kimberlite, (4) relative abundance of elongate and unbroken olivine crystals and (5) paucity of country rock xenoliths, while the remaining two facies are interpreted as fragmental on the basis of: (1) the gradational contacts with demonstrably fragmental kimberlite, (2) relative abundance and range of sizes of country rock lithic clasts and (3) numerous broken olivine crystals. Comparisons are made with coherent and apparent-coherent kimberlite from the literature. Our three coherent facies are similar to literature reported coherent kimberlite dykes hosted in country rock (CKd) in terms of internal organisation, low abundance of country rock xenoliths, and apparent-porphyritic texture. Conversely, our two fragmental facies share attributes with previously described pipe-filling coherent and apparent-coherent kimberlite (CKpf) in terms of geometry, internal organisation and abundance of country rock xenoliths. We conclude that CKd and most CKpf, although similar in matrix/groundmass mineralogy and texture, can be distinguished on the basis of internal organisation, country rock lithic clast abundance, texture (e.g. apparent-porphyritic texture) and possibly olivine crystal shapes and suggest that fragmental kimberlite is more common than reported.

The fossil calibration database - A new resource for divergence dating

Fossils provide the principal basis for temporal calibrations, which are critical to the accuracy of divergence dating analyses. Translating fossil data into minimum and maximum bounds for calibrations is the most important, and often least appreciated, step of divergence dating. Properly justified calibrations require the synthesis of phylogenetic, paleontological, and geological evidence and can be difficult for non-specialists to formulate. The dynamic nature of the fossil record (e.g., new discoveries, taxonomic revisions, updates of global or local stratigraphy) requires that calibration data be updated continually lest they become obsolete. Here, we announce the Fossil Calibration Database (http://fossilcalibrations.org), a new open-access resource providing vetted fossil calibrations to the scientific community. Calibrations accessioned into this database are based on individual fossil specimens and follow best practices for phylogenetic justification and geochronological constraint. The associated Fossil Calibration Series, a calibration-themed publication series at Palaeontologia Electronica, will serve as one key pipeline for peer-reviewed calibrations to enter the database.

A life cycle assessment approach to quantifying greenhouse gas emissions from land-use change for beef production in eastern Australia

In life cycle assessment studies, greenhouse gas (GHG) emissions from direct land-use change have been estimated to make a significant contribution to the global warming potential of agricultural products. However, these estimates have a high uncertainty due to the complexity of data requirements and difficulty in attribution of land-use change. This paper presents estimates of GHG emissions from direct land-use change from native woodland to grazing land for two beef production regions in eastern Australia, which were the subject of a multi-impact life cycle assessment study for premium beef production. Spatially- and temporally consistent datasets were derived for areas of forest cover and biomass carbon stocks using published remotely sensed tree-cover data and regionally applicable allometric equations consistent with Australia's national GHG inventory report. Standard life cycle assessment methodology was used to estimate GHG emissions and removals from direct land-use change attributed to beef production. For the northern-central New South Wales region of Australia estimates ranged from a net emission of 0.03 t CO2-e ha-1 year-1 to net removal of 0.12 t CO2-e ha-1 year-1 using low and high scenarios, respectively, for sequestration in regrowing forests. For the same period (1990-2010), the study region in southern-central Queensland was estimated to have net emissions from land-use change in the range of 0.45-0.25 t CO2-e ha-1 year-1. The difference between regions reflects continuation of higher rates of deforestation in Queensland until strict regulation in 2006 whereas native vegetation protection laws were introduced earlier in New South Wales. On the basis of liveweight produced at the farm-gate, emissions from direct land-use change for 1990-2010 were comparable in magnitude to those from other on-farm sources, which were dominated by enteric methane. However, calculation of land-use change impacts for the Queensland region for a period starting 2006, gave a range from net emissions of 0.11 t CO2-e ha-1 year-1 to net removals of 0.07 t CO2-e ha-1 year-1. This study demonstrated a method for deriving spatially- and temporally consistent datasets to improve estimates for direct land-use change impacts in life cycle assessment. It identified areas of uncertainty, including rates of sequestration in woody regrowth and impacts of land-use change on soil carbon stocks in grazed woodlands, but also showed the potential for direct land-use change to represent a net sink for GHG.

Seasonal and spatial variations in rare earth elements to identify inter-aquifer linkages and recharge processes in an Australian catchment

With the aim of elucidating the seasonal behaviour of rare earth elements (REEs), surface and groundwaters were collected under dry and wet conditions in different hydrological units of the Teviot Brook catchment (Southeast Queensland, Australia). Sampled waters showed a large degree of variability in both REE abundance and normalised patterns. Overall REE abundance ranged over nearly three orders of magnitude, and was consistently lower in the sedimentary bedrock aquifer (18ppt<∑REE<477ppt) than in the other hydrological systems studied. Abundance was greater in springs draining rhyolitic rocks (∑REE=300 and 2054ppt) than in springs draining basalt ranges (∑REE=25 and 83ppt), yet was highly variable in the shallow alluvial groundwater (16ppt<∑REE<5294ppt) and, to a lesser extent, in streamwater (85ppt<∑REE<2198ppt). Generally, waters that interacted with different rock types had different REE patterns. In order to obtain an unbiased characterisation of REE patterns, the ratios between light and middle REEs (R(M/L)) and the ratios between middle and heavy REEs (R(H/M)) were calculated for each sample. The sedimentary bedrock aquifer waters had highly evolved patterns depleted in light REEs and enriched in middle and heavy REEs (0.17geological settings

Low delta18O zircon grains in the Neoarchean Rum Jungle Complex, northern Australia: An indicator of emergent continental crust

The timing of widespread continental emergence is generally considered to have had a dramatic effect on the hydrological cycle, atmospheric conditions, and climate. New secondary ion mass spectrometry (SIMS) oxygen and laser-ablation–multicollector–inductively coupled plasma–mass spectrometry (LA-MC-ICP-MS) Lu-Hf isotopic results from dated zircon grains in the granitic Neoarchean Rum Jungle Complex provide a minimum time constraint on the emergence of continental crust above sea level for the North Australian craton. A 2535 ± 7 Ma monzogranite is characterized by magmatic zircon with slightly elevated δ18O (6.0‰–7.5‰ relative to Vienna standard mean ocean water [VSMOW]), consistent with some contribution to the magma from reworked supracrustal material. A supracrustal contribution to magma genesis is supported by the presence of metasedimentary rock enclaves, a large population of inherited zircon grains, and subchondritic zircon Hf (εHf = −6.6 to −4.1). A separate, distinct crustal source to the same magma is indicated by inherited zircon grains that are dominated by low δ18O values (2.5‰–4.8‰, n = 9 of 15) across a range of ages (3536–2598 Ma; εHf = −18.2 to +0.4). The low δ18O grains may be the product of one of two processes: (1) grain-scale diffusion of oxygen in zircon by exchange with a low δ18O magma or (2) several episodes of magmatic reworking of a Mesoarchean or older low δ18O source. Both scenarios require shallow crustal magmatism in emergent crust, to allow interaction with rocks altered by hydrothermal meteoric water in order to generate the low δ18O zircon. In the first scenario, assimilation of these altered rocks during Neoarchean magmatism generated low δ18O magma with which residual detrital zircons were able to exchange oxygen, while preserving their U-Pb systematics. In the second scenario, wholesale melting of the altered rocks occurred in several distinct events through the Mesoarchean, generating low δ18O magma from which zircon crystallized. Ultimately, in either scenario, the low δ18O zircons were entrained as inherited grains in a Neoarchean granite. The data suggest operation of a modern hydrological cycle by the Neoarchean and add to evidence for the increased emergence of continents by this time

Evolution of a ~2.7 Ga large igneous province: A volcanological, geochemical and geochronological study of the Agnew Greenstone Belt, and new regional correlations for the Kalgoorlie Terrane (Yilgarn Craton, Western Australia)

The thick package of ~2.7 Ga mafic and ultramafic lavas and intrusions preserved among the Neoarchean of the Kalgoorlie Terrene in Western Australia provides valuable insight into geological processes controlling the most prodigious episode of growth and preservation of juvenile continental crust in Earth’s history. Limited exposure of these rocks results in uncertainty about their age, physical and chemical characteristics, and stratigraphic relationships. This in turn prevents confident correlation of regional occurrences of mafic and ultramafic successions (both intrusive and extrusive) and hinders the interpretation of tectonic setting and magmatic evolution. A recent stratigraphic drilling program of the Neoarchean stratigraphy of the Agnew Greenstone Belt in Western Australia has provided continuous exposures through a c. 7 km thick sequence of mafic and ultramafic units. In this study, we present a volcanological, lithogeochemical and chronological study of the Agnew Greenstone Belt, and provide the first pre-2690 Ma regional correlation across the Kalgoorlie Terrane. The Agnew Greenstone Belt records ~30 m.y. of episodic ultramafic-mafic magmatism that includes two cycles, each defined by a komatiite that is overlain by units that become more evolved and contaminated with time. The sequence is divided into nine conformable packages, each consisting of stacked subaqueous lava flows and comagmatic intrusions, as well as two sills without associated extrusions. Lavas, with the exception of intercalations between two units, form a layer-cake stratigraphy and were likely erupted from a system of fissures tapping the same magma source. The komatiites are not contaminated by continental crust ([La/Sm]PM ~0.7) and are of the Al-undepleted Munro-type. Crustal contamination is evident in many units (Songvang Basalt, Never Can Tell Basalt, Redeemer Basalt, and Turrett Dolerite), as judged by [La/Sm]>1, negative Nb and Ti anomalies, and geochemical mixing trends towards felsic contaminants. Crystal fractionation was also significant, with early olivine and chromite (Mg#>65) followed by plagioclase and clinopyroxene removal (Mg<65), and in the most evolved case, titanomagnetite accumulation. Three new TIMS dates on granophyric zones of mafic sills and one ICP-MS date from an interflow felsic tuff are presented and used for regional stratigraphic correlation. Cycle I magmatism began at ~2720 Ma and ended ~2705 Ma, whereas cycle II began ~2705 Ma and ended at 2690.7±1.2 Ma. Regional correlations indicate the western Kalgoorlie Terrane consists of a remarkably similar stratigraphy that can be recognised at Agnew, Ora Banda and Coolgardie, whereas the eastern part of the terrane (e.g., Kambalda Domain) does not include cycle I, but correlates well with cycle II. This research supports an autochthonous model of greenstone formation, in which one large igneous province, represented by two complete cycles, is constructed on sialic crust. New stratigraphic correlations for the Kalgoorlie Terrane indicate that many units can be traced over distances >100 km, which has implications for exploration targeting for stratigraphically hosted ultramafic Ni and VMS deposits.