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.

Sediment concentration prediction and statistical evaluation for annual load estimation

We consider the development of statistical models for prediction of constituent concentration of riverine pollutants, which is a key step in load estimation from frequent flow rate data and less frequently collected concentration data. We consider how to capture the impacts of past flow patterns via the average discounted flow (ADF) which discounts the past flux based on the time lapsed - more recent fluxes are given more weight. However, the effectiveness of ADF depends critically on the choice of the discount factor which reflects the unknown environmental cumulating process of the concentration compounds. We propose to choose the discount factor by maximizing the adjusted R-2 values or the Nash-Sutcliffe model efficiency coefficient. The R2 values are also adjusted to take account of the number of parameters in the model fit. The resulting optimal discount factor can be interpreted as a measure of constituent exhaustion rate during flood events. To evaluate the performance of the proposed regression estimators, we examine two different sampling scenarios by resampling fortnightly and opportunistically from two real daily datasets, which come from two United States Geological Survey (USGS) gaging stations located in Des Plaines River and Illinois River basin. The generalized rating-curve approach produces biased estimates of the total sediment loads by -30% to 83%, whereas the new approaches produce relatively much lower biases, ranging from -24% to 35%. This substantial improvement in the estimates of the total load is due to the fact that predictability of concentration is greatly improved by the additional predictors.

Understanding coal seam gas associated water, regulations and strategies for treatment

The coal seam gas (CSG) industry is globally of potentially great importance economically. This study exemplifies the complex relationship between land use and management, groundwater impact and associated water treatment especially in relation to Queensland where a significant increase in the amount of gas extracted over the past 6 years has occurred. In order to effectively manage the environmental impact of the CSG industry it is necessary to appropriately understand the nature of the gas deposits, methods for gas collection, the physicochemical composition of the by-product associated water and the technologies available for water remediation. Australia is mainly considered arid and semi-arid and thus there is a need to not only beneficially reuse water resources but also protect existing ground water reservoirs such as the Great Artesian Basin (GAB). This paper focussed primarily on the Surat Basin located in Queensland and northern New South Wales. The mechanism for CSG formation, relation to local geological features, extraction approach and the potential impact/benefits of associated water was discussed. An outline of the current legislative requirements on physical and chemical properties of associated water in the Surat Basin was also provided, as well as the current treatment technologies used by the major CSG companies. This review was of significance in relation to the formulation of the most appropriate and cost effective management of associated water, while simultaneously preserving existing water resources and the environment.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

Denitrification losses from an intensively managed sub-tropical pasture: Impact of soil moisture on the partitioning of N2 and N2O emissions

Intensively managed pastures in subtropical Australia under dairy production are nitrogen (N) loaded agro-ecosystems, with an increased pool of N available for denitrification. The magnitude of denitrification losses and N2:N2O partitioning in these agro-ecosystems is largely unknown, representing a major uncertainty when estimating total N loss and replacement. This study investigated the influence of different soil moisture contents on N2 and N2O emissions from a subtropical dairy pasture in Queensland, Australia. Intact soil cores were incubated over 15 days at 80% and 100% water-filled pore space (WFPS), after the application of 15N labelled nitrate, equivalent to 50 kg N ha−1. This setup enabled the direct quantification of N2 and N2O emissions following fertilisation using the 15N gas flux method. The main product of denitrification in both treatments was N2. N2 emissions exceeded N2O emissions by a factor of 8 ± 1 at 80% WFPS and a factor of 17 ± 2 at 100% WFPS. The total amount of N-N2 lost over the incubation period was 21.27 kg ± 2.10 N2-N ha−1 at 80% WFPS and 25.26 kg ± 2.79 kg ha−1 at 100% WFPS respectively. N2 emissions remained high at 100% WFPS, while related N2O emissions decreased. At 80% WFPS, N2 emissions increased constantly over time while N2O fluxes declined. Consequently, N2/(N2 + N2O) product ratios increased over the incubation period in both treatments. N2/(N2 + N2O) product ratios responded significantly to soil moisture, confirming WFPS as a key driver of denitrification. The substantial amount of fertiliser lost as N2 reveals the agronomic significance of denitrification as a major pathway of N loss for sub-tropical pastures at high WFPS and may explain the low fertiliser N use efficiency observed for these agro-ecosystems.