Fred’s Flow (Canada) and Murphy Well (Australia) : thick komatiitic lava flows with contrasting compositions, emplacement mechanisms and water contents

Two Archaean komatiitic flows, Fred’s Flow in Canada and the Murphy Well Flow in Australia, have similar thicknesses (120 and 160 m) but very different compositions and internal structures. Their contrasting differentiation profiles are keys to determine the cooling and crystallization mechanisms that operated during the eruption of Archaean ultramafic lavas. Fred’s Flow is the type example of a thick komatiitic basalt flow. It is strongly differentiated and consists of a succession of layers with contrasting textures and compositions. The layering is readily explained by the accumulation of olivine and pyroxene in a lower cumulate layer and by evolution of the liquid composition during downward growth of spinifex-textured rocks within the upper crust. The magmas that erupted to form Fred’s Flow had variable compositions, ranging from 12 to 20 wt% MgO, and phenocryst contents from 0 to 20 vol%. The flow was emplaced by two pulses. A first ~20-m-thick pulse was followed by another more voluminous but less magnesian pulse that inflated the flow to its present 120 m thickness. Following the second pulse, the flow crystallized in a closed system and differentiated into cumulates containing 30–38 wt% MgO and a residual gabbroic layer with only 6 wt% MgO. The Murphy Well Flow, in contrast, has a remarkably uniform composition throughout. It comprises a 20-m-thick upper layer of fine-grained dendritic olivine and 2–5 vol% amygdales, a 110–120 m intermediate layer of olivine porphyry and a 20–30 m basal layer of olivine orthocumulate. Throughout the flow, MgO contents vary little, from only 30 to 33 wt%, except for the slightly more magnesian basal layer (38–40 wt%). The uniform composition of the flow and dendritic olivine habits in the upper 20 m point to rapid cooling of a highly magnesian liquid with a composition like that of the bulk of the flow. Under equilibrium conditions, this liquid should have crystallized olivine with the composition Fo94.9, but the most magnesian composition measured by electron microprobe in samples from the flow is Fo92.9. To explain these features, we propose that the parental liquid contained around 32 wt% MgO and 3 wt% H2O. This liquid degassed during the eruption, creating a supercooled liquid that solidified quickly and crystallized olivine with non-equilibrium textures and compositions.

Hydrochemical evolution and groundwater flow processes in the Galilee and Eromanga basins, Great Artesian Basin, Australia: A multivariate statistical approach

The Galilee and Eromanga basins are sub-basins of the Great Artesian Basin (GAB). In this study, a multivariate statistical approach (hierarchical cluster analysis, principal component analysis and factor analysis) is carried out to identify hydrochemical patterns and assess the processes that control hydrochemical evolution within key aquifers of the GAB in these basins. The results of the hydrochemical assessment are integrated into a 3D geological model (previously developed) to support the analysis of spatial patterns of hydrochemistry, and to identify the hydrochemical and hydrological processes that control hydrochemical variability. In this area of the GAB, the hydrochemical evolution of groundwater is dominated by evapotranspiration near the recharge area resulting in a dominance of the Na–Cl water types. This is shown conceptually using two selected cross-sections which represent discrete groundwater flow paths from the recharge areas to the deeper parts of the basins. With increasing distance from the recharge area, a shift towards a dominance of carbonate (e.g. Na–HCO3 water type) has been observed. The assessment of hydrochemical changes along groundwater flow paths highlights how aquifers are separated in some areas, and how mixing between groundwater from different aquifers occurs elsewhere controlled by geological structures, including between GAB aquifers and coal bearing strata of the Galilee Basin. The results of this study suggest that distinct hydrochemical differences can be observed within the previously defined Early Cretaceous–Jurassic aquifer sequence of the GAB. A revision of the two previously recognised hydrochemical sequences is being proposed, resulting in three hydrochemical sequences based on systematic differences in hydrochemistry, salinity and dominant hydrochemical processes. The integrated approach presented in this study which combines different complementary multivariate statistical techniques with a detailed assessment of the geological framework of these sedimentary basins, can be adopted in other complex multi-aquifer systems to assess hydrochemical evolution and its geological controls.

Effects of magnetite on high-frequency ground-penetrating radar

Large concentrations of magnetite in sedimentary deposits and soils with igneous parent material have been reported to affect geophysical sensor performance. We have undertaken the first systematic experimental effort to understand the effects of magnetite for ground-penetrating radar (GPR) characterization of the shallow subsurface. Laboratory experiments were conducted to study how homogeneous magnetite-sand mixtures and magnetite concentrated in layers affect the propagation behavior (velocity, attenuation) of high-frequency GPR waves and the reflection characteristics of a buried target. Important observations were that magnetite had a strong effect on signal velocity and reflection, at magnitudes comparable to what has been observed in small-scale laboratory experiments that measured electromagnetic properties of magnetite-silica mixtures. Magnetite also altered signal attenuation and affected the reflection characteristics of buried targets. Our results indicated important implications for several fields, including land mine detection, Martian exploration, engineering, and moisture mapping using satellite remote sensing and radiometers.

Landform characterization using geophysics—Recent advances, applications, and emerging tools

This paper presents an overview of the strengths and limitations of existing and emerging geophysical tools for landform studies. The objectives are to discuss recent technical developments and to provide a review of relevant recent literature, with a focus on propagating field methods with terrestrial applications. For various methods in this category, including ground-penetrating radar (GPR), electrical resistivity (ER), seismics, and electromagnetic (EM) induction, the technical backgrounds are introduced, followed by section on novel developments relevant to landform characterization. For several decades, GPR has been popular for characterization of the shallow subsurface and in particular sedimentary systems. Novel developments in GPR include the use of multi-offset systems to improve signal-to-noise ratios and data collection efficiency, amongst others, and the increased use of 3D data. Multi-electrode ER systems have become popular in recent years as they allow for relatively fast and detailed mapping. Novel developments include time-lapse monitoring of dynamic processes as well as the use of capacitively-coupled systems for fast, non-invasive surveys. EM induction methods are especially popular for fast mapping of spatial variation, but can also be used to obtain information on the vertical variation in subsurface electrical conductivity. In recent years several examples of the use of plane wave EM for characterization of landforms have been published. Seismic methods for landform characterization include seismic reflection and refraction techniques and the use of surface waves. A recent development is the use of passive sensing approaches. The use of multiple geophysical methods, which can benefit from the sensitivity to different subsurface parameters, is becoming more common. Strategies for coupled and joint inversion of complementary datasets will, once more widely available, benefit the geophysical study of landforms.Three cases studies are presented on the use of electrical and GPR methods for characterization of landforms in the range of meters to 100. s of meters in dimension. In a study of polygonal patterned ground in the Saginaw Lowlands, Michigan, USA, electrical resistivity tomography was used to characterize differences in subsurface texture and water content associated with polygon-swale topography. Also, a sand-filled thermokarst feature was identified using electrical resistivity data. The second example is on the use of constant spread traversing (CST) for characterization of large-scale glaciotectonic deformation in the Ludington Ridge, Michigan. Multiple CST surveys parallel to an ~. 60. m high cliff, where broad (~. 100. m) synclines and narrow clay-rich anticlines are visible, illustrated that at least one of the narrow structures extended inland. A third case study discusses internal structures of an eolian dune on a coastal spit in New Zealand. Both 35 and 200. MHz GPR data, which clearly identified a paleosol and internal sedimentary structures of the dune, were used to improve understanding of the development of the dune, which may shed light on paleo-wind directions.

Variation in parental magmas of Mt Rouse, a complex polymagmatic monogenetic volcano in the basaltic intraplate Newer Volcanics Province, southeast Australia

Monogenetic volcanoes have long been regarded as simple in nature, involving single magma batches and uncomplicated evolutions; however, recent detailed research into individual centres is challenging that assumption. Mt Rouse (Kolor) is the volumetrically largest volcano in the monogenetic Newer Volcanics Province of southeast Australia. This study presents new major, trace and Sr–Nd–Pb isotope data for samples selected on the basis of a detailed stratigraphic framework analysis of the volcanic products from Mt Rouse. The volcano is the product of three magma batches geochemically similar to Ocean–Island basalts, featuring increasing LREE enrichment with each magma batch (batches A, B and C) but no evidence of crustal contamination; the Sr–Nd–Pb isotopes define two groupings. Modelling suggests that the magmas were sourced from a zone of partial melting crossing the lithosphere–asthenosphere boundary, with batch A forming a large volume partial melt in the deep lithosphere (1.7 GPa/55.5 km); and batches B and C from similar areas within the shallow asthenosphere (1.88 GPa/61 km and 1.94 GPa/63 km, respectively). The formation and extraction of these magmas may have been due to high deformation rates in the mantle caused by edge-driven convection and asthenospheric upwelling. The lithosphere– asthenosphere boundary is important with respect to NVP volcanism. An eruption chronology involves sequential eruption of magma batches A, C and B, followed by simultaneous eruption of batches A and B. Mt Rouse is a complex polymagmatic monogenetic volcano that illustrates the complexity of monogenetic volcanism and demonstrates the importance of combining detailed stratigraphic analysis alongside systematic geochemical sampling.

Lower mantle diamonds from Rio Soriso (Juina area, Mato Grosso, Brazil)

The morphology, colour, fluorescence, cathodoluminescence, nitrogen content and aggregation state, internal structure and mineral inclusions have been studied for 69 alluvial diamonds from the Rio Soriso (Juina area, Mato Grosso State, Brazil). Nitrogen in most diamonds (53%) is fully aggregated as B centres, but there is also a large proportion of N-free stones (38%). A strong positive correlation between nitrogen and IR-active hydrogen concentrations is observed. The diamonds contain (in order of decreasing abundance) ferropericlase, CaSi-perovskite, magnetite, MgSi-perovskite, pyrrhotite, 'olivine', SiO2, perovskite, tetragonal almandine-pyrope phase and some other minerals represented by single grains. The Rio Soriso diamond suite is subdivided into several subpopulations that originated in upper and lower mantle of ultramafic and mafic compositions, with the largest subgroup forming in the ultramafic lower mantle. Analysed ferropericlase grains are enriched in Fe (Mg#=0.43-0.89), which is ascribed to their origin in the lowermost mantle. The Juina kimberlites may be unique in sampling the material from depths below 1,700 km that ascended in a plume formed at the core-mantle boundary.

The Cougar Point tuff: Implications for thermochemical zonation and longevity of high-temperature, large-volume silicic magmas of the Miocene Yellowstone hotspot

The 12.7-10.5 Ma Cougar Point Tuff in southern Idaho, USA, consists of 10 large-volume (>10²-10³ km³ each), high-temperature (800-1000 °C), rhyolitic ash-flow tuffs erupted from the Bruneau-Jarbidge volcanic center of the Yellowstone hotspot. These tuffs provide evidence for compositional and thermal zonation in pre-eruptive rhyolite magma, and suggest the presence of a long-lived reservoir that was tapped by numerous large explosive eruptions. Pyroxene compositions exhibit discrete compositional modes with respect to Fe and Mg that define a linear spectrum punctuated by conspicuous gaps. Airfall glass compositions also cluster into modes, and the presence of multiple modes indicates tapping of different magma volumes during early phases of eruption. Equilibrium assemblages of pigeonite and augite are used to reconstruct compositional and thermal gradients in the pre-eruptive reservoir. The recurrence of identical compositional modes and of mineral pairs equilibrated at high temperatures in successive eruptive units is consistent with the persistence of their respective liquids in the magma reservoir. Recurrence intervals of identical modes range from 0.3 to 0.9 Myr and suggest possible magma residence times of similar duration. Eruption ages, magma temperatures, Nd isotopes, and pyroxene and glass compositions are consistent with a long-lived, dynamically evolving magma reservoir that was chemically and thermally zoned and composed of multiple discrete magma volumes.

High-temperature, low-H₂O silicic magmas of the Yellowstone Hotspot: an experimental study of Rhyolite from the Bruneau-Jarbidge eruptive center, Central Snake River Plain, USA

The phase relations have been investigated experimentally at 200 and 500 MPa as a function of water activity for one of the least evolved (Indian Batt Rhyolite) and of a more evolved rhyolite composition (Cougar Point Tuff XV) from the 12·8-8·1 Ma Bruneau-Jarbidge eruptive center of the Yellowstone hotspot. Particular priority was given to accurate determination of the water content of the quenched glasses using infrared spectroscopic techniques. Comparison of the composition of natural and experimentally synthesized phases confirms that high temperatures (>900°C) and extremely low melt water contents (<1·5 wt % H₂O) are required to reproduce the natural mineral assemblages. In melts containing 0·5-1·5 wt % H₂O, the liquidus phase is clinopyroxene (excluding Fe-Ti oxides, which are strongly dependent on fO₂), and the liquidus temperature of the more evolved Cougar Point Tuff sample (BJR; 940-1000°C) is at least 30°C lower than that of the Indian Batt Rhyolite lava sample (IBR2; 970-1030°C). For the composition BJR, the comparison of the compositions of the natural and experimental glasses indicates a pre-eruptive temperature of at least 900°C. The composition of clinopyroxene and pigeonite pairs can be reproduced only for water contents below 1·5 wt % H₂O at 900°C, or lower water contents if the temperature is higher. For the composition IBR2, a minimum temperature of 920°C is necessary to reproduce the main phases at 200 and 500 MPa. At 200 MPa, the pre-eruptive water content of the melt is constrained in the range 0·7-1·3 wt % at 950°C and 0·3-1·0 wt % at 1000°C. At 500 MPa, the pre-eruptive temperatures are slightly higher (by 30-50°C) for the same ranges of water concentration. The experimental results are used to explore possible proxies to constrain the depth of magma storage. The crystallization sequence of tectosilicates is strongly dependent on pressure between 200 and 500 MPa. In addition, the normative Qtz-Ab-Or contents of glasses quenched from melts coexisting with quartz, sanidine and plagioclase depend on pressure and melt water content, assuming that the normative Qtz and Ab/Or content of such melts is mainly dependent on pressure and water activity, respectively. The combination of results from the phase equilibria and from the composition of glasses indicates that the depth of magma storage for the IBR2 and BJR compositions may be in the range 300-400 MPa (13 km) and 200-300 MPa (10 km), respectively.

Multiple sulfur isotope analyses support a magmatic model for the volcanogenic massive sulfide deposits of the Teutonic Bore volcanic complex, Yilgarn Craton, Western Australia

We report sensitive high mass resolution ion microprobe, stable isotopes (SHRIMP SI) multiple sulfur isotope analyses (32S, 33S, 34S) to constrain the sources of sulfur in three Archean VMS deposits—Teutonic Bore, Bentley, and Jaguar—from the Teutonic Bore volcanic complex of the Yilgarn Craton, Western Australia, together with sedimentary pyrites from associated black shales and interpillow pyrites. The pyrites from VMS mineralization are dominated by mantle sulfur but include a small amount of slightly negative mass-independent fractionation (MIF) anomalies, whereas sulfur from the pyrites in the sedimentary rocks has pronounced positive MIF, with ∆33S values that lie between 0.19 and 6.20‰ (with one outlier at −1.62‰). The wall rocks to the mineralization include sedimentary rocks that have contributed no detectable positive MIF sulfur to the VMS deposits, which is difficult to reconcile with the leaching model for the formation of these deposits. The sulfur isotope data are best explained by mixing between sulfur derived from a magmatic-hydrothermal fluid and seawater sulfur as represented by the interpillow pyrites. The massive sulfide lens pyrites have a weighted mean ∆33S value of −0.27 ± 0.05‰ (MSWD = 1.6) nearly identical with −0.31 ± 0.08‰ (MSWD = 2.4) for pyrites from the stringer zone, which requires mixing to have occurred below the sea floor. We employed a two-component mixing model to estimate the contribution of seawater sulfur to the total sulfur budget of the two Teutonic Bore volcanic complex VMS deposits. The results are 15 to 18% for both Teutonic Bore and Bentley, much higher than the 3% obtained by Jamieson et al. (2013) for the giant Kidd Creek deposit. Similar calculations, carried out for other Neoarchean VMS deposits give value between 2% and 30%, which are similar to modern hydrothermal VMS deposits. We suggest that multiple sulfur isotope analyses may be used to predict the size of Archean VMS deposits and to provide a vector to ore deposit but further studies are needed to test these suggestions.

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.17sedimentary bedrock from upgradient to downgradient; typically they showed flat patterns in the upstream section of the alluvium (median R(M/L)=0.21 and median R(H/M)=-0.06) gradually evolving towards patterns depleted in light REEs and enriched in middle and heavy REEs downgradient (median R(M/L)=0.48 and median R(H/M)=0.38). To document the seasonal variations in REE patterns, the difference in ratios between dry and wet sampling campaigns was determined for each repeated sampling location. Contributions from the sedimentary bedrock water to the alluvium during the wet season were identified at two locations (increase from R(H/M)=0.03 and 0.35 to R(H/M)=0.62 and 0.89). The effect of recharge through fractured igneous rocks was also observed in two boreholes intercepting the sedimentary bedrock, where the freshly recharged waters likely contributed to the deeper groundwater flow during the wet season (decrease from R(M/L)=0.81 and 0.56 to R(M/L)=0.46 and 0.17). Results from this study suggest that REEs may be usefully applied as indicators of recharge processes and inter-aquifer mixing. They also underline the importance of conducting seasonal sampling campaigns to capture possible short-term variations in REE patterns and abundance, which is essential to enable a better understanding of hydrological and hydrochemical processes in complex geological settings

When using small samples to evaluate hydrocarbon reservoirs, proceed with caution

Small-angle neutron scattering (SANS) and ultra-small angle neutron scattering (USANS) with contrast matching techniques (Melnichenko and others, 2012) were used to investigate size distribution and gas accessibility in pores in an approximately 10.6 cm long Mississippian Barnett Shale butt core from the Fort Worth Basin, Texas, USA. SANS and USANS measurements record scattering from all pores, both open and closed, in the size range 10nm - ~10 μ. The techniques can also be used to determine the material that contains pores and the number of pores as a function of size. By injecting deuterated methane gas (CD4) at contrast matching pressure it is possible to distinguish which pores are accessible, or open, to fluids and which ones are not.

The tectonic significance of lower Permain successions in the Texas Orocline (Eastern Australia)

The Texas Orocline is a prominent orogenic curvature that developed during the early Permian in the southern New England Orogen. Outliers preserving lower Permian sedimentary successions (Bondonga, Silver Spur, Pikedale, Terrica, Alum Rock and Ashford beds) approximately outline the oroclinal structure, but the tectonic processes responsible for the development of these basinal successions, and their relationships to the Texas Orocline, are unclear. Here we address this shortcoming by providing new U–Pb detrital and primary zircon ages from these successions, as well as detailed stratigraphic and structural data from the largest exposed succession (Bondonga beds). Field observations and U–Pb geochronological data suggest that the lower Permian successions in the Texas Orocline are remnants of a single, formerly larger basin that was deposited after ca 302 Ma. Time constraints for formation of this basin are correlative with constraints from the lower Permian Nambucca Block, which was likely deposited in response to regional back-arc extension during and/or after the development of the Texas Orocline. The conclusion that the lower Permian sedimentary basins in the Texas Orocline belong to this back-arc extensional system supports the suggestion that oroclinal bending in the New England Orogen was primarily controlled by trench retreat and associated overriding-plate extension.

Structural-stratigraphic reconstruction of the lower Whakataki formation, north island, New Zealand

This thesis details a Miocene aged sedimentary rock formation located in north island New Zealand. Mapping, stratigraphic logging and petrographic analysis of the rock formation ascertained that it was deposited in a deep-marine, tectonically active region. The work details the make-up of the sedimentary rocks using geochemistry and microscopy to define their origin. This definition was used to interpret the depositional model of the sediments detailing how they were transported and how they were emplaced.