High-precision tephrostratigraphy : tracking the time-varying eruption pulse of Mt. Taranaki, North Island, New Zealand : a thesis presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Earth Science, Massey University, Palmerston North, New Zealand

In this research it was proposed that a more robust record of volcanic activity for Mt. Taranaki (New Zealand) could be derived from tephras (pyroclastic fall deposits) within cores from several lakes and peatlands across a 120o arc, NE-SE of the volcano, covering a range of prevailing down-wind directions. These data were integrated with previous tephrochronology studies to construct one of the longest and most complete volcanic eruption history records ever developed for an andesitic stratovolcano. Using 44 new radiocarbon dates, electron microprobe analysis of glass shard and titanomagnetite chemical composition, along with whole-rock chemistry, a chrono- and chemostratigraphy was established. The new record identifies at least 272 tephraproducing eruptions over the last 30 cal ka BP. Six chemo-stratigraphic groups were identified: A (0.5 – 3 cal ka BP), B (3 – 4 cal ka BP), C (4 – 9.5 cal ka BP), D (9.5 – 14 cal ka BP), E (14 – 17.5 cal ka BP), and F (23.5 – 30 cal ka BP). These were used to resolve previous stratigraphic uncertainties at upper-flank (proximal) and ring-plain (medial) sites. Several well-known “marker tephras” are now recognized as being ~2000 years older than previously determined (e.g., Waipuku, Tariki, and Mangatoki Tephra units) with the prominent Korito Tephra stratigraphically positioned above the Taupo-derived Stent Tephra. Further, new markers were identified, including the Kokowai Tephra unit (~4.7 cal ka BP), at a beach-cliff exposure, 40-km north-east of the volcano. Once age-models were established for each tephra, units were matched between sites using statistical methods. Initial statistical integration showed that the immediate past high-resolution tephrochronological record suffered from a distinctive “old-carbon” effect on its ages (Lake Rotokare). This had biased the most recent probabilistic forecasting and generated artificially high probability estimates (52-59% eruption chance over the next 50 years). Once the Rotokare record was excluded and chemostratigraphy constraints were applied, a reliable multi-site tephra record could be built only for the last ~14 ka BP. The new data confirms a highly skewed distribution of mainly (98% of cases) short intervals between eruptions (mode of ~9 years and average interval ~65 years). Long intervals (up to 580 years) as seen in earlier records were reduced to 2% of the record, but can now be considered real, rather than missing data. The new data confirm a cyclic pattern of varying eruption frequency (with a five-fold range in annual frequency) on a period of ~1000-1500 years. The new time-varying frequency estimates suggest a lower probability for a new eruption at Mt. Taranaki over the next 50 years of 33-42%. The newly established chemostratigraphy was further used to investigate time-related compositional changes. Whole-lapilli analyses highlighted that a specific very evolved Ca-rich and Fe-poor composition was only found within the easterly and south-easterly depositional sites. This was explained by eruption of a stratified magma reservoir, which holds greater modal proportions of plagioclase and lower proportions of pyroxene within low-density, gas-rich upper conduit regions. During the most explosive phases of eruptions, when plumes reach the stratospheric jetstream, the lowest-density pumice is thus dispersed by high-level stable westerly winds. Further, two distinct evolutional trends were seen in the long and new tephrochronological record; from 17.5 to 3 cal ka BP and <3 cal ka BP; with wholelapilli, glass, and titanomagnetite compositions overall evolving over time. The former compositional trend indicates a crystallising and cooling magma source in the deep crust, with multiple, spatially separated magma source regions forming, each generating magmas (i.e., magma batches) with unique titanomagnetite compositions. This trend is interrupted by a distinct shift towards less-evolved compositions and the initiation of a second parasitic vent (Fanthams Peak at the southern flank of Mt. Taranaki)

Multi-scale rock characterisation and data integration: a case study from the BGS Core Scanning Facility

Physical and chemical characterisation of sedimentary successions from core underpin a number of areas of geoscience, including source rock and reservoir evaluation, correlation, and, palaeoenvironmental analyses. Recent technical advances in the capabilities and analytical precision of X-ray florescence (XRF), computed tomography (CT), and multi-sensor (MSCLS) core scanners mean that they are increasingly applicable to such studies. In contrast to traditional semi-destructive sampling techniques, they offer the ability to rapidly acquire non-destructive and nearcontinuous records. These measurements can be integrated with other datasets at different scale, thereby increasing analytical resolution and the value of core material at the same time as reducing costs. The new Core-Scanning Facility (CSF) at the headquarters of the British Geological Survey (BGS) in Keyworth houses four high-resolution core scanners – two XRF’s, a CT, and a MSCLS – enabling geophysical, mineralogical, and geochemical characterisation of core and allowing high-definition optical, near-infrared (NIR), ultraviolet (UV), and Xradiographic images to be collected. The facility was commissioned to provide fundamental support to the NERC/UKRI funded UK Geo-Energy Observatories (UKGEOS; https://www.ukgeos.ac.uk/), an initiative to characterise the subsurface environment at 2 localities (Glasgow and Cheshire) in order to better understand the impact of subterranean infrastructure and the role rocks may play in decarbonising energy supplies. Here, we use the preliminary results from the Glasgow Observatory as a case study to demonstrate how core scanner datasets are generated, evaluated, and visualised. This involves integrating core scanner outputs with other data generated at widely different scales (e.g. wireline logs, core plug analysis). We highlight the value of core scanning as a mechanism for correctly reconstructing the vertical position of core, for identifying specific intervals of interest, and for resolving new relationships between physical and chemical properties. The core scanner data can also highlight changes in the rock properties that would not be easily identified by conventional downhole geophysical logging and point geochemical sampling. Our ultimate goal is to generate a multi-scale “digital twin” for each core, freely available online, that will allow a more complete understanding of the variability of rock properties

CoreScore: an ML approach to assess legacy core condition

Today's geoscience challenges often require repurposing of data and samples from legacy boreholes. Collection of new deep core is expensive; maximising this investment is vital. However, condition of legacy cores varies due to factors including recovery, sampling, lithology, and storage. Rock Quality Designation analysis is often undertaken on new core but this only provides a snapshot of core condition and will not be indicative of subsequent condition. Poor core condition can make destructive analytical techniques impossible and also impacts non-destructive techniques including core scanning. Since 2011, BGS have systematically collected 125,000 core images. This study investigates if core condition of this archive can be assessed using automated analysis by machine learning. A neural network-based approach was used to segment these images. By differentiating imaged core from their background, properties such as number of fragments and total rock area were determined and used to assess core condition. Analysis of outputs demonstrate that with minimal input data, core condition can be rapidly assessed. This allows users to better understand and visualise core. This can be used to qualitatively assess non-destructive data, improve success of destructive sampling through targeted sampling and reduce the time and effort spent interacting with physical material

Environmental changes during the onset of the Late Pliensbachian Event (Early Jurassic) in the Cardigan Bay Basin, Wales

The Late Pliensbachian Event (LPE), in the Early Jurassic, is associated with a perturbation in the global carbon cycle (positive carbon isotope excursion (CIE) of ∼2 ‰), cooling of ∼5 ∘C, and the deposition of widespread regressive facies. Cooling during the late Pliensbachian has been linked to enhanced organic matter burial and/or disruption of thermohaline ocean circulation due to a sea level lowstand of at least regional extent. Orbital forcing had a strong influence on the Pliensbachian environments and recent studies show that the terrestrial realm and the marine realm in and around the Cardigan Bay Basin, UK, were strongly influenced by orbital climate forcing. In the present study we build on the previously published data for long eccentricity cycle E459 ± 1 and extend the palaeoenvironmental record to include E458 ± 1. We explore the environmental and depositional changes on orbital timescales for the Llanbedr (Mochras Farm) core during the onset of the LPE. Clay mineralogy, X-ray fluorescence (XRF) elemental analysis, isotope ratio mass spectrometry, and palynology are combined to resolve systematic changes in erosion, weathering, fire, grain size, and riverine influx. Our results indicate distinctively different environments before and after the onset of the LPE positive CIE and show increased physical erosion relative to chemical weathering. We also identify five swings in the climate, in tandem with the 405 kyr eccentricity minima and maxima. Eccentricity maxima are linked to precessionally repeated occurrences of a semi-arid monsoonal climate with high fire activity and relatively coarser sediment from terrestrial runoff. In contrast, 405 kyr minima in the Mochras core are linked to a more persistent, annually wet climate, low fire activity, and relatively finer-grained deposits across multiple precession cycles. The onset of the LPE positive CIE did not impact the expression of the 405 kyr cycle in the proxy records; however, during the second pulse of heavier carbon (13C) enrichment, the clay minerals record a change from dominant chemical weathering to dominant physical erosion

Petrological and geochemical characterisation of the sarsen stones at Stonehenge.

Little is known of the properties of the sarsen stones (or silcretes) that comprise the main architecture of Stonehenge. The only studies of rock struck from the monument date from the 19th century, while 20th century investigations have focussed on excavated debris without demonstrating a link to specific megaliths. Here, we present the first comprehensive analysis of sarsen samples taken directly from a Stonehenge megalith (Stone 58, in the centrally placed trilithon horseshoe). We apply state-of-the-art petrographic, mineralogical and geochemical techniques to two cores drilled from the stone during conservation work in 1958. Petrographic analyses demonstrate that Stone 58 is a highly indurated, grain-supported, structureless and texturally mature groundwater silcrete, comprising fine-to-medium grained quartz sand cemented by optically-continuous syntaxial quartz overgrowths. In addition to detrital quartz, trace quantities of silica-rich rock fragments, Fe-oxides/hydroxides and other minerals are present. Cathodoluminescence analyses show that the quartz cement developed as an initial <10 μm thick zone of non-luminescing quartz followed by ~16 separate quartz cement growth zones. Late-stage Fe-oxides/hydroxides and Ti-oxides line and/or infill some pores. Automated mineralogical analyses indicate that the sarsen preserves 7.2 to 9.2 area % porosity as a moderately-connected intergranular network. Geochemical data show that the sarsen is chemically pure, comprising 99.7 wt. % SiO2. The major and trace element chemistry is highly consistent within the stone, with the only magnitude variations being observed in Fe content. Non-quartz accessory minerals within the silcrete host sediments impart a trace element signature distinct from standard sedimentary and other crustal materials. 143Nd/144Nd isotope analyses suggest that these host sediments were likely derived from eroded Mesozoic rocks, and that these Mesozoic rocks incorporated much older Mesoproterozoic material. The chemistry of Stone 58 has been identified recently as representative of 50 of the 52 remaining sarsens at Stonehenge. These results are therefore representative of the main stone type used to build what is arguably the most important Late Neolithic monument in Europe

Petrography, Geochemistry and Mineralogy of the Stonehenge Sarsens: Digital Data Collection

This collection includes a suite of digital materials that, in combination, characterise the petrography, mineralogy and geochemistry of a sarsen upright (Stone 58) from the central trilithon horseshoe at Stonehenge. The collection arises from work undertaken during the British Academy/Leverhulme Trust project "Geochemical fingerprinting the sarsen stones at Stonehenge" (Small Research Grant SG-170610), led by the University of Brighton. The data accompanies the publication: Nash, D. J., Ciborowski, T. J. R., Darvill, T., Parker Pearson, M., Ullyott, J. S., Damaschke, M., Evans, J. A., Goderis, S., Greaney, S., Huggett, J. M., Ixer, R. A., Pirrie, D., Power, M. R., Salge, T. & Whitaker, N. (2021, in review) Petrological and geochemical characterisation of the sarsen stones at Stonehenge. PLoS ONE

Sedimentology on profiles from Lake Prespa, Greece

Hydro-acoustic surveys and coring campaigns at Lake Prespa were carried out between 2007 and 2009. This paper presents hydro-acoustic profiles and provide lithological and chronostratigraphical information from three up to 15.75 m long sediment sequences from the Macedonian side of the lake. The sediment sequences comprise glacial and interglacial sediments likely deposited from the end of Marine Isotope Stage (MIS) 5 to present day. The information implies a distinct change of sedimentation patterns at the Pleistocene/Holocene transition and the establishment of a relatively strong Holocene current system and deposition of channel-related contourite drift in Lake Prespa. Potential causes for the establishment of this current during the Holocene include significant lake level change, reduced winter ice cover, and/or higher aeolian activity