The ~AD1315 Tarawera and Waiotapu eruptions, New Zealand: contemporaneous rhyolite and hydrothermal eruptions driven by an arrested basalt dike system?

A series of large hydrothermal eruptions occurred across the Waiotapu geothermal field at about the same (prehistoric) time as the ~AD1315 Kaharoa rhyolite magmatic eruptions from Tarawera volcano vents, 10–20 km distant. Triggering of the Waiotapu hydrothermal eruptions was previously attributed to displacement of the adjacent Ngapouri Fault. The Kaharoa rhyolite eruptions are now recognised as primed and triggered by multiple basalt intrusions beneath the Tarawera volcano. A ~1000 t/day pulse of CO₂ gas is recorded by alteration mineralogy and fluid inclusions in drill core samples from Waiotapu geothermal wells. This CO₂ pulse is most readily sourced from basalt intruded at depth, and although not precisely dated, it appears to be associated with the Waiotapu hydrothermal eruptions. We infer that the hydrothermal eruptions at Waiotapu were primed by intrusion of the same arrested basalt dike system that drove the rhyolite eruptions at Tarawera. This dike system was likely similar at depth to the dike that generated basalt eruptions from a 17 km-long fissure that formed across the Tarawera region in AD1886. Fault ruptures that occurred in the Waiotapu area in association with both the AD1886 and ~AD1315 eruptions are considered to be a result, rather than a cause, of the dike intrusion processes.Editorial responsibility: J. Donnelly-Nolan     

Crustal structure of the northeastern margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin

The 1000-km-long Darlag–Lanzhou–Jingbian seismic refraction profile is located in the NE margin of the Tibetan plateau. This profile crosses the northern Songpan-Ganzi terrane, the Qinling-Qilian fold system, the Haiyuan arcuate tectonic region, and the stable Ordos basin. The P-wave and S-wave velocity structure and Poisson's ratios reveal many significant characteristics in the profile. The crustal thickness increases from northeast to southwest. The average crustal thickness observed increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane. The crust becomes obviously thicker south of the Haiyuan fault and beneath the West-Qinlin Shan. The crustal velocities have significant variations along the profile. The average P-wave velocities for the crystalline crust vary between 6.3 and 6.4 km/s. Beneath the Songpan-Ganzi terrane, West-Qinling Shan, and Haiyuan arcuate tectonic region P-wave velocities of 6.3 km/s are 0.15 km/s lower than the worldwide average of 6.45 km/s. North of the Kunlun fault, with exclusion of the Haiyuan arcuate tectonic region, the average P-wave velocity is 6.4 km/s and only 0.5 km/s lower than the worldwide average. A combination of the P-wave velocity and Poisson's ratio suggests that the crust is dominantly felsic in composition with an intermediate composition at the base. A mafic lower crust is absent in the NE margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin. There are low velocity zones in the West-Qinling Shan and the Haiyuan arcuate tectonic region. The low velocity zones have low S-wave velocities and high Poisson's ratios, so it is possible these zones are due to partial melting. The crust is divided into two layers, the upper and the lower crust, with crustal thickening mainly in the lower crust as the NE Tibetan plateau is approached. The results in the study show that the thickness of the lower crust increases from 22 to 38 km as the crustal thickness increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane south of the Kunlun fault. Both the Conrad discontinuity and Moho in the West-Qinling Shan and in the Haiyuan arcuate tectonic region are laminated interfaces, implying intense tectonic activity. The arcuate faults and large earthquakes in the Haiyuan arcuate tectonic region are the result of interaction between the Tibetan plateau and the Sino–Korean and Gobi Ala Shan platforms.     

Constraining recycled detritus in quartz‐rich sandstones: Insights from a multi‐proxy provenance study of the Mid‐Carboniferous,Clare Basin,western Ireland

Quartz‐rich sandstones can be produced through multiple sedimentary processes, potentially acting in combination, such as extensive sedimentary recycling or intense chemical weathering. Determining the provenance of such sedimentary rocks can be challenging due to low amounts of accessory minerals, the fact that the primary mineralogy may have been altered during transport, storage or burial and difficulties in the recognition of polycyclic components. This study uses zircon and apatite U‐Pb geochronology, apatite trace elements, zircon‐tourmaline‐rutile indices and petrographic observations to investigate the sedimentary history of mineralogically mature mid‐Carboniferous sandstones of the Tullig Cyclothem, Clare Basin, western Ireland. The provenance data show that the sandstones have been dominantly and ultimately sourced from three basement terranes: older Laurentian‐ associated rocks (ca. 900–2500 Ma) which lay to the north of the basin, peri‐Gondwanan terranes (ca. 500–700 Ma) to the south and igneous intrusive rocks associated with the Caledonian Orogenic Cycle (ca. 380–500 Ma). However, the multi‐proxy approach also helps constrain the sedimentary history and suggests that not all grain populations were derived directly from their original source. Grains with a Laurentian or a Caledonian affinity have likely been recycled through Devonian basins to the south. Grains with a peri‐Gondwanan affinity appear to be first cycle and are potentially derived from south/southwest of the basin. Taken as a whole, these data are consistent with input into the basin from the south and southwest, with the reworking of older sedimentary rocks, rather than intensive first‐cycle chemical weathering, likely explaining the compositional maturity of the sandstones. This study highlights the need for a multi‐proxy provenance approach to constrain sedimentary recycling, particularly in compositionally mature sandstones, as the use of zircon geochronology alone would have led to erroneous provenance interpretations. Zircon, together with U‐Pb geochronology from more labile phases such as apatite, can help distinguish first‐cycle versus polycyclic detritus.     

Automated statistical matching of multiple tephra records exemplified using five long maar sequences younger than 75 ka,Auckland, New Zealand

Detailed tephrochronologies are built to underpin probabilistic volcanic hazard forecasting, and to understand the dynamics and history of diverse geomorphic, climatic, soil-forming and environmental processes. Complicating factors include highly variable tephra distribution over time; difficulty in correlating tephras from site to site based on physical and chemical properties; and uncertain age determinations. Multiple sites permit construction of more accurate composite tephra records, but correctly merging individual site records by recognizing common events and site-specific gaps is complex. We present an automated procedure for matching tephra sequences between multiple deposition sites using stochastic local optimization techniques. If individual tephra age determinations are not significantly different between sites, they are matched and a more precise age is assigned. Known stratigraphy and mineralogical or geochemical compositions are used to constrain tephra matches. We apply this method to match tephra records from five long sediment cores (≤ 75 cal ka BP) in Auckland, New Zealand. Sediments at these sites preserve basaltic tephras from local eruptions of the Auckland Volcanic Field as well as distal rhyolitic and andesitic tephras from Okataina, Taupo, Egmont, Tongariro, and Tuhua (Mayor Island) volcanic centers. The new correlated record compiled is statistically more likely than previously published arrangements from this area.     

A cascade of magmatic events during the assembly and eruption of a super-sized magma body

We use comprehensive geochemical and petrological records from whole-rock samples, crystals, matrix glasses and melt inclusions to derive an integrated picture of the generation, accumulation and evacuation of 530 km³ of crystal-poor rhyolite in the 25.4 ka Oruanui supereruption (New Zealand). New data from plagioclase, orthopyroxene, amphibole, quartz, Fe–Ti oxides, matrix glasses, and plagioclase- and quartz-hosted melt inclusions, in samples spanning different phases of the eruption, are integrated with existing data to build a history of the magma system prior to and during eruption. A thermally and compositionally zoned, parental crystal-rich (mush) body was developed during two periods of intensive crystallisation, 70 and 10–15 kyr before the eruption. The mush top was quartz-bearing and as shallow as ~3.5 km deep, and the roots quartz-free and extending to >10 km depth. Less than 600 year prior to the eruption, extraction of large volumes of ~840 °C low-silica rhyolite melt with some crystal cargo (between 1 and 10%), began from this mush to form a melt-dominant (eruptible) body that eventually extended from 3.5 to 6 km depth. Crystals from all levels of the mush were entrained into the eruptible magma, as seen in mineral zonation and amphibole model pressures. Rapid translation of crystals from the mush to the eruptible magma is reflected in textural and compositional diversity in crystal cores and melt inclusion compositions, versus uniformity in the outermost rims. Prior to eruption the assembled eruptible magma body was not thermally or compositionally zoned and at temperatures of ~790 °C, reflecting rapid cooling from the ~840 °C low-silica rhyolite feedstock magma. A subordinate but significant volume (3–5 km³) of contrasting tholeiitic and calc-alkaline mafic material was co-erupted with the dominant rhyolite. These mafic clasts host crystals with compositions which demonstrate that there was some limited pre-eruptive physical interaction of mafic magmas with the mush and melt-dominant body. However, the mafic magmas do not appear to have triggered the eruption or controlled magmatic temperatures in the erupted rhyolite. Integration of textural and compositional data from all available crystal types, across all dominant and subordinate magmatic components, allow the history of the Oruanui magma body to be reconstructed over a wide range of temporal scales using multiple techniques. This history spans the tens of millennia required to grow the parental magma system (U–Th disequilibrium dating in zircon), through the centuries and decades required to assemble the eruptible magma body (textural and diffusion modelling in orthopyroxene), to the months, days, hours and minutes over which individual phases of the eruption occurred, identified through field observations tied to diffusion modelling in magnetite, olivine, quartz and feldspar. Tectonic processes, rather than any inherent characteristics of the magmatic system, were a principal factor acting to drive the rapid accumulation of magma and control its release episodically during the eruption. This work highlights the richness of information that can be gained by integrating multiple lines of petrologic evidence into a holistic timeline of field-verifiable processes.     

Zircon trace element chemistry at sub-micrometer resolution for Tarawera volcano,New Zealand,and implications for rhyolite magma evolution

Zoned crystals can be important recorders of magmatic processes in space and time. However, in most situations, the temporal dimension is difficult to quantify. Here, we have employed secondary ion mass spectrometry depth profiling to excavate parallel pits into non-polished crystal faces of zircon to obtain ~5 μm resolution U–Th disequilibrium ages (one pit) that can be correlated with trace element zoning at sub-μm resolution derived from a second pit. Data from 17 crystals representing each of the four rhyolite eruptions of Tarawera volcano, an intra-caldera edifice within the Okataina Volcanic Centre, reveal diverse zircon growth conditions over time. Most crystals display rimward depletions in Zr/Hf and Ti, broadly consistent with cooling and crystallization. However, a significant fraction of crystals lacks these patterns and displays rimward trace element variations consistent with isothermal or prograde crystallization. Oscillatory zonation patterns in Y, Th, and U are superimposed on the Zr/Hf and Ti trends. Despite the limited number of crystals analyzed in this way, the striking lack of ubiquitous trace element zoning patterns in crystals from the same hand sample implies that fractional crystallization upon cooling was punctuated by magma recharge and crystal mixing affecting different parts of the magma reservoir. By combining data from all crystals, a systematic change to more heterogeneous trace element abundances is revealed by zircon crystal domains <45 ka following the Rotoiti caldera-forming eruption. This contrasts with the more uniform conditions of zircon crystallization lasting >100 ka prior to caldera formation and is best explained by the post-caldera system consisting of small, isolated melt pockets that evolved independently. An important conclusion is that the zircon ‘cargo’ in volcanic rocks reflects thermally and compositionally divergent processes that act near simultaneously in a magma storage region and not exclusively the conditions in the eruptible magma.