Explosive activity on Kīlauea’s Lower East Rift Zone fuelled by a volatile-rich, dacitic melt

Magmas with matrix glass compositions ranging from basalt to dacite erupted from a series of 24 fissures in the first two weeks of the 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano. Eruption styles ranged from low spattering and fountaining to strombolian activity. Major element trajectories in matrix glasses and melt inclusions hosted by olivine, pyroxene and plagioclase are consistent with variable amounts of fractional crystallization, with incompatible elements (e.g., Cl, F, H2O) becoming enriched by 4-5 times as melt MgO contents evolve from 6 to 0.5 wt%. The high viscosity and high H2O contents (~2 wt%) of the dacitic melts erupting at Fissure 17 account for the explosive Strombolian behavior exhibited by this fissure, in contrast to the low fountaining and spattering observed at fissures erupting basaltic to basaltic-andesite melts. Saturation pressures calculated from melt inclusions CO2-H2O contents indicate that the magma reservoir(s) supplying these fissures was located at ~2-3 km depth, which is in agreement with the depth of a dacitic magma body intercepted during drilling in 2005 (~2.5 km) and a seismically-imaged low Vp/Vs anomaly (~2 km depth). Nb/Y ratios in erupted products are similar to lavas erupted between 1955-1960, indicating that melts were stored and underwent variable amounts of crystallization in the LERZ for >60 years before being remobilized by a dike intrusion in 2018. We demonstrate that extensive fractional crystallization generates viscous and volatile-rich magma with potential for hazardous explosive eruptions, which may be lurking undetected at many ocean island volcanoes.

Timing and Mechanisms of Silicate Differentiation and Basalt Magma Generation on the Howardite-Eucrite-Diogenite Asteroid Parent Body

<p>Meteorites provide the only direct record of the chronology and nature of the processes that occurred in the early solar system. In this study, meteorites were examined in order to gain insight into the timing and nature of magmatism and silicate differentiation on asteroidal bodies in the first few million years of the solar system. These bodies are considered the precursors to terrestrial planets, and as such they provide information about conditions in the solar system at the time of planet formation. This study focuses on eucrites, which are basaltic meteorites that are believed to represent the crust of the Howardite-Eucrite-Diogenite (HED) parent body. The processes of silicate differentiation and the relationship between eucrites and the diogenitic mafic cumulate of the HED parent body are poorly understood. The major and trace element chemistry of the minerals in the eucrite suite was measured. There is little variability in mineral major element concentrations in eucrites, however considerable variability was observed in mineral trace element concentrations, particularly with respect to incompatible elements in the mineral phases. Magnesium was separated from digested eucrite samples, and the Mg isotope composition of the eucrites was measured to high precision in order to date the samples using the short-lived ²⁶Al–²⁶Mg chronometer and examine magmatic evolution on the HED parent body. Correlations between incompatible elements in pyroxene and ²⁶Mg anomalies, produced by the decay of ²⁶Al, indicate that the eucrite suite was formed from a single, evolving magma body. Large trace element and Mg isotopic differences between eucrites and diogenites indicate that the two meteorite groups did not, as previously suggested, originate from the same magma body. Instead they may have formed from two large magma bodies, which were spatially or temporally separated on the HED parent body. The application of the short-lived ²⁶Al–²⁶Mg chronometer to this suite of eucrites constrains the onset of eucrite formation to ~3 Myr after the formation of the solar system’s first solids, as a result of rapid accretion and melting of planetesimals due to heating from the decay of ²⁶Al.</p>

Timing and Mechanisms of Silicate Differentiation and Basalt Magma Generation on the Howardite-Eucrite-Diogenite Asteroid Parent Body

<p>Meteorites provide the only direct record of the chronology and nature of the processes that occurred in the early solar system. In this study, meteorites were examined in order to gain insight into the timing and nature of magmatism and silicate differentiation on asteroidal bodies in the first few million years of the solar system. These bodies are considered the precursors to terrestrial planets, and as such they provide information about conditions in the solar system at the time of planet formation. This study focuses on eucrites, which are basaltic meteorites that are believed to represent the crust of the Howardite-Eucrite-Diogenite (HED) parent body. The processes of silicate differentiation and the relationship between eucrites and the diogenitic mafic cumulate of the HED parent body are poorly understood. The major and trace element chemistry of the minerals in the eucrite suite was measured. There is little variability in mineral major element concentrations in eucrites, however considerable variability was observed in mineral trace element concentrations, particularly with respect to incompatible elements in the mineral phases. Magnesium was separated from digested eucrite samples, and the Mg isotope composition of the eucrites was measured to high precision in order to date the samples using the short-lived ²⁶Al–²⁶Mg chronometer and examine magmatic evolution on the HED parent body. Correlations between incompatible elements in pyroxene and ²⁶Mg anomalies, produced by the decay of ²⁶Al, indicate that the eucrite suite was formed from a single, evolving magma body. Large trace element and Mg isotopic differences between eucrites and diogenites indicate that the two meteorite groups did not, as previously suggested, originate from the same magma body. Instead they may have formed from two large magma bodies, which were spatially or temporally separated on the HED parent body. The application of the short-lived ²⁶Al–²⁶Mg chronometer to this suite of eucrites constrains the onset of eucrite formation to ~3 Myr after the formation of the solar system’s first solids, as a result of rapid accretion and melting of planetesimals due to heating from the decay of ²⁶Al.</p>

Chemical heterogeneities reveal early rapid cooling of Apollo Troctolite 76535

The evolution of the lunar interior is constrained by samples of the magnesian suite of rocks returned by the Apollo missions. Reconciling the paradoxical geochemical features of this suite constitutes a feasibility test of lunar differentiation models. Here we present the results of a microanalytical examination of the archetypal specimen, troctolite 76535, previously thought to have cooled slowly from a large magma body. We report a degree of intra-crystalline compositional heterogeneity (phosphorus in olivine and sodium in plagioclase) fundamentally inconsistent with prolonged residence at high temperature. Diffusion chronometry shows these heterogeneities could not have survived magmatic temperatures for >~20 My, i.e., far less than the previous estimated cooling duration of >100 My. Quantitative modeling provides a constraint on the thermal history of the lower lunar crust, and the textural evidence of dissolution and reprecipitation in olivine grains supports reactive melt infiltration as the mechanism by which the magnesian suite formed.

Post-Oruanui supereruption recovery, reconstruction and evolution of Taupo volcano, New Zealand

<p>This thesis research presents geochemical perspectives on the magmatic recovery of Taupo volcano (New Zealand) in the aftermath of the 25.4 ka Oruanui supereruption. Following the Oruanui, and after only ~5 kyr of quiescence, Taupo erupted three small volume (~0.1 km3) dacitic units, followed by another ~5 kyr break, and then the modern sequence from ~12 ka onwards of 25 rhyolitic units organised into 3 geochemically distinct subgroups (SG1-SG3). The eruptive units are stratigraphically constrained over exceptionally short time intervals, providing fine-scale temporal snapshots of the magma system. In this thesis I compare and contrast whole-rock, mineral and glass compositions of Oruanui and post-Oruanui magmas through time to investigate the post-supereruption reconstruction and evolution of Taupo through to the latest eruption. Despite overlapping vent sites and crustal source domains between the Oruanui and post-Oruanui eruptions, U/Th disequilibrium model-ages in zircons from Taupo SG1 rhyolites (erupted 12 ka-10 ka) and SG2 rhyolites (erupted 7 ka-2.6 ka) imply the presence of only minor inheritance of crystals from the Oruanui magma source. Post-Oruanui model-age spectra are instead typically centred close to eruption ages with subordinate older pre-300 ka equiline grains. U-Pb dating of these equiline grains shows that both 300-450 ka plutonic-derived and pre-100 Ma greywacke basement-derived zircons are present. The former largely coincide in age with zircons from the 350 ka Whakamaru eruption products, and are dominant over greywacke in young units which were vented within the published Whakamaru caldera outline. Despite multiple ages and vent sites, trace element compositions are broadly similar in zircons, regardless of their ages. However, a small subset of zircons analysed from SG1 rhyolites have notably high concentrations of U, Th, P, Y+ (REE)3+ and Nb but with only minor changes in Hf and Ti. SG2 zircons typically have higher Sc, reflecting large-scale changes in melt chemistry and crystallising mineral phases with time. The age spectra indicate that most Oruanui zircons were removed by thermally induced dissolution immediately following the supereruption. U-Th ages from individual post-Oruanui eruptions show consistent inheritance of post-Oruanui grains with model ages that centre between the temporally separated but geographically overlapping eruption groups, generating model-age modes. Within the statistical limitations of the isotopic measurements, we interpret these repeated modes to be significant, resulting from incorporation of crystal populations from cyclic post-Oruanui periods of magmatic cooling and crystallisation, acting within a crustal protolith chemically independent of that which built the Oruanui. Cooling periods alternate with times of rejuvenation and eruption, in some cases demonstrably accompanying syn-eruptive regional rifting and mafic injection. Not only were the processes that developed the supersized Oruanui magma body unusually rapid, but this huge magma system was effectively reset and rebuilt on a comparably short timescale. Major and trace element whole rock, glass and mineral chemistry of post-Oruanui eruptive products indicate how the host magma system re-established and evolved. The dacite units show wide variations in melt inclusion compositions and strongly zoned minerals consistent with interaction of less-evolved mafic magmas at a depths of >8 km, overlapping with the inferred base of the old Oruanui mush system. The dacites reflect the first products of the rebuilding silicic magma system, as most of the Oruanui mush was reconfigured or significantly modified in composition following thermal fluxing accompanying post-caldera collapse readjustment. The first (SG1) rhyolites erupted from 12 ka formed through shallow fractionation (4-5 km depth) and cooling of a parental melt similar in composition to the earlier dacite melts, with overlapping melt inclusion and crystal core compositions between the two magma types. For the younger rhyolite units, fine-scale temporal changes in melt chemistry and mineral phase stability occur over time, which are closely linked to the development, stabilisation and maturation of a new and likely unitary rhyolite mush system at Taupo. The new mush system is closely linked to and sometimes physically interacts with the underlying mafic melts, which are similar in composition to those involved in the Oruanui eruption and provide the long-term thermal and chemical driving force for magmatism. We consider that the new mush body has expanded to >250 km3 (and possibly up to 1000 km3) but has not yet been located by geophysical investigations. For the most recent SG3 eruptions, the system once again underwent widespread destabilisation, resulting in increased levels of melt extraction from the silicic mush. Trends in whole-rock chemistry and close links between melt inclusions and mineral zoning with earlier units indicates that the 35 km3 Unit Y (Taupo eruption) melt dominant body formed in response to mafic disruption of the silicic mush pile. Associated Fe-Mg diffusion timescales in orthopyroxene suggest that Taupo is capable of changing behaviour and generating large eruptible melt bodies on timescales as short as decades to centuries. The 232 AD Unit Y eruption culminated from a critical combination of high differential tectonic stress build up, and increased potency in the silicic magma system resulting from elevated levels of mafic magma input, resulting in one of the largest and most violent worldwide Holocene eruptions. The post-Y magma system then responded to further disruption with the eruption of sub-lacustrine dome(s). Taupo is considered to be capable of rapidly recovering in its modern form to continue its hyperactive eruptive behaviour on timescales that are of human interest and concern.</p>

The dynamics of large-scale silicic magmatic systems: case studies from Mangakino Volcanic Centre, Taupo Volcanic Zone, New Zealand

<p>This thesis research focuses on clast and crystal-specific studies to investigate the pre- and syn-eruptive magmatic processes of two supereruptions in the TVZ: the 1.21 Ma Ongatiti (>500 km3) and the 1.0 Ma Kidnappers (~1200 km3), together with the smaller (~200 km3) 1.0 Ma Rocky Hill eruption from the Mangakino Volcanic Centre (MVC). Crystallisation histories determined through SIMS U-Pb dating of zircons reveal that the paired Kidnappers and Rocky Hill eruptions were products of a common magmatic system, which built over ~200 kyr, in the time break after the Ongatiti eruption. U-Pb age spectra from the Ongatiti show a protracted crystallisation history (over ~250 kyr), in which the majority of zircon crystallised ~100 kyr prior to eruption in a crystal mush. Zircons then ascended with melt during accumulation of the final erupted magma body in the shallow crust. Zircons remained stable in the melt dominant body but underwent little further crystallisation. Zircons from all three systems record common geochemical processes governed by the fractionating assemblage (predominantly plagioclase and amphibole). In particular, the MREE/HREE ratios and Sr concentrations of zircons from the Ongatiti record imply two contrasting source regions governed by different proportions of crystallising amphibole. The in-situ major and trace element chemistry of glass shards and crystals from the Kidnappers fall deposit reveal that magma within the Kidnappers was stored in three discrete bodies, which were systematically tapped during the early stages of eruption. Temperature and pressure (T-P) estimates from amphibole and Fe-Ti oxide equilibria from each magma type are similar and therefore the three magma bodies were adjacent, not vertically stacked, in the crust. Amphibole model T-P estimates range from 770 to 840 °C and 90 to 170 MPa corresponding to pre-eruptive storage depths of ~4.0-6.5 km. The systematic evacuation of the three independent magma bodies implies that there was tectonic triggering and linkage of eruptions. The termination of fall deposition and onset of the overlying ignimbrite emplacement marks the point of widespread caldera collapse and the catastrophic evacuation of a wider variety of melt during the Kidnappers eruption. Pumice compositions from the Kidnappers ignimbrite fall into three groups, two of which (KI-1 and KI-2) can be matched to bodies tapped during the fall phase of the eruption, with the addition of a further discrete batch of lower SiO2 (KI-3) magma. Core-rim textural and chemical variations in major crystal phases (plagioclase, amphibole and orthopyroxene) suggest each compositional group was sourced from a common mush but underwent a unique magmatic history during the development of melt-dominant bodies in the final stages prior to eruption. The field relationships and distinctive appearance of the Rocky Hill ignimbrite (~200 km3 DRE) and the underlying Kidnappers ignimbrite suggests that the two deposits are from distinct eruption events. However, major and trace element chemistry of matrix glass, coupled with the textural and chemical signatures of crystals suggests the magma erupted during the Rocky Hill was generated from the same source or mush zone as the Kidnappers. The two largest melt-dominant bodies (KI-1 and KI-2) within the Kidnappers were renewed, underwent mixing and incorporation of marginal material to form two magma types (RH-1 and RH-2) in the time break prior to the Rocky Hill eruption. Fe-Mg interdiffusion timescales in orthopyroxenes from the Kidnappers and Rocky Hill deposits suggest the establishment of the final melt-dominant bodies, through extraction of melt and crystals from a common mush, occurred within 1000 years, and peaked within centuries of each eruption. In addition, one discrete batch of Kidnappers melt has evidence for interaction with a lesser evolved melt within 50 yrs prior to eruption. This rejuvenation event was not the eruption trigger but may have primed the magma for eruption. The difference in timescales from common zones from both the Kidnappers and Rocky Hill orthopyroxene, recording the same processes reveal the time break between the two eruptions was ~20-40 years. This work highlights the rapidity of rejuvenation and renewal of the melt-dominant bodies within the Kidnappers/Rocky Hill magmatic system. The textural and in-situ compositional signatures of crystals from the Ongatiti ignimbrite imply the final erupted magma body was assembled from a thermally and chemically zoned mush, which extended to the base of the quartzofeldspathic crust (~15km). The mush was close to water saturation and was dominated by amphibole crystallisation. Melt and crystals (including the majority of zircons) were extracted from the mush and ascended to 4-6 km depths during the development of a crystal-rich (20-30%), but melt-dominant body. Significant crystallisation of plagioclase (and lesser proportions of orthopyroxene and amphibole) occurred in an event involving the gradual heating and/or increase of water in the rhyolite, from a broadly andesitic underplated magma. Homogeneous crystal rim and matrix glass compositions imply the final erupted volume of magma was effectively mixed through convection. Eu/Eu* values of whole-rock and matrix glass suggest little crystal-melt separation occurred in the melt-dominant magma body prior to eruption. This work has implications for understanding the generation, storage and eruption of large-scale silicic magma systems. The Ongatiti ignimbrite does not represent either an erupted mush, or a stratified magma chamber, suggesting an alternative model for the development of eruptible magma within large-scale silicic systems. The Kidnappers/Rocky Hill sequence records a complex interplay of multiple melt-dominant bodies, which were established and renewed on rapid timescales. The rapid timescales for the development of melt-dominant bodies and the systematic tapping of magmas in the Kidnappers/Rocky Hill system imply that tectonics may have had a strong external control on the eruptions at Mangakino.</p>

Post-Oruanui supereruption recovery, reconstruction and evolution of Taupo volcano, New Zealand

<p>This thesis research presents geochemical perspectives on the magmatic recovery of Taupo volcano (New Zealand) in the aftermath of the 25.4 ka Oruanui supereruption. Following the Oruanui, and after only ~5 kyr of quiescence, Taupo erupted three small volume (~0.1 km3) dacitic units, followed by another ~5 kyr break, and then the modern sequence from ~12 ka onwards of 25 rhyolitic units organised into 3 geochemically distinct subgroups (SG1-SG3). The eruptive units are stratigraphically constrained over exceptionally short time intervals, providing fine-scale temporal snapshots of the magma system. In this thesis I compare and contrast whole-rock, mineral and glass compositions of Oruanui and post-Oruanui magmas through time to investigate the post-supereruption reconstruction and evolution of Taupo through to the latest eruption. Despite overlapping vent sites and crustal source domains between the Oruanui and post-Oruanui eruptions, U/Th disequilibrium model-ages in zircons from Taupo SG1 rhyolites (erupted 12 ka-10 ka) and SG2 rhyolites (erupted 7 ka-2.6 ka) imply the presence of only minor inheritance of crystals from the Oruanui magma source. Post-Oruanui model-age spectra are instead typically centred close to eruption ages with subordinate older pre-300 ka equiline grains. U-Pb dating of these equiline grains shows that both 300-450 ka plutonic-derived and pre-100 Ma greywacke basement-derived zircons are present. The former largely coincide in age with zircons from the 350 ka Whakamaru eruption products, and are dominant over greywacke in young units which were vented within the published Whakamaru caldera outline. Despite multiple ages and vent sites, trace element compositions are broadly similar in zircons, regardless of their ages. However, a small subset of zircons analysed from SG1 rhyolites have notably high concentrations of U, Th, P, Y+ (REE)3+ and Nb but with only minor changes in Hf and Ti. SG2 zircons typically have higher Sc, reflecting large-scale changes in melt chemistry and crystallising mineral phases with time. The age spectra indicate that most Oruanui zircons were removed by thermally induced dissolution immediately following the supereruption. U-Th ages from individual post-Oruanui eruptions show consistent inheritance of post-Oruanui grains with model ages that centre between the temporally separated but geographically overlapping eruption groups, generating model-age modes. Within the statistical limitations of the isotopic measurements, we interpret these repeated modes to be significant, resulting from incorporation of crystal populations from cyclic post-Oruanui periods of magmatic cooling and crystallisation, acting within a crustal protolith chemically independent of that which built the Oruanui. Cooling periods alternate with times of rejuvenation and eruption, in some cases demonstrably accompanying syn-eruptive regional rifting and mafic injection. Not only were the processes that developed the supersized Oruanui magma body unusually rapid, but this huge magma system was effectively reset and rebuilt on a comparably short timescale. Major and trace element whole rock, glass and mineral chemistry of post-Oruanui eruptive products indicate how the host magma system re-established and evolved. The dacite units show wide variations in melt inclusion compositions and strongly zoned minerals consistent with interaction of less-evolved mafic magmas at a depths of >8 km, overlapping with the inferred base of the old Oruanui mush system. The dacites reflect the first products of the rebuilding silicic magma system, as most of the Oruanui mush was reconfigured or significantly modified in composition following thermal fluxing accompanying post-caldera collapse readjustment. The first (SG1) rhyolites erupted from 12 ka formed through shallow fractionation (4-5 km depth) and cooling of a parental melt similar in composition to the earlier dacite melts, with overlapping melt inclusion and crystal core compositions between the two magma types. For the younger rhyolite units, fine-scale temporal changes in melt chemistry and mineral phase stability occur over time, which are closely linked to the development, stabilisation and maturation of a new and likely unitary rhyolite mush system at Taupo. The new mush system is closely linked to and sometimes physically interacts with the underlying mafic melts, which are similar in composition to those involved in the Oruanui eruption and provide the long-term thermal and chemical driving force for magmatism. We consider that the new mush body has expanded to >250 km3 (and possibly up to 1000 km3) but has not yet been located by geophysical investigations. For the most recent SG3 eruptions, the system once again underwent widespread destabilisation, resulting in increased levels of melt extraction from the silicic mush. Trends in whole-rock chemistry and close links between melt inclusions and mineral zoning with earlier units indicates that the 35 km3 Unit Y (Taupo eruption) melt dominant body formed in response to mafic disruption of the silicic mush pile. Associated Fe-Mg diffusion timescales in orthopyroxene suggest that Taupo is capable of changing behaviour and generating large eruptible melt bodies on timescales as short as decades to centuries. The 232 AD Unit Y eruption culminated from a critical combination of high differential tectonic stress build up, and increased potency in the silicic magma system resulting from elevated levels of mafic magma input, resulting in one of the largest and most violent worldwide Holocene eruptions. The post-Y magma system then responded to further disruption with the eruption of sub-lacustrine dome(s). Taupo is considered to be capable of rapidly recovering in its modern form to continue its hyperactive eruptive behaviour on timescales that are of human interest and concern.</p>

The Bishop Tuff, California: New Insights into Magmatic Timescales and Processes from Micro-Analytical Approaches

<p>The Bishop Tuff is the product of one of the largest eruptions on Earth in the last 1 Myr. This thesis studies the Bishop Tuff in order to better understand the nature of the pre-eruptive magma body, with an emphasis on the processes that occurred within it and the timescales over which they operated. In situ geochemical analyses of crystals and glass from samples collected throughout the Bishop Tuff stratigraphic succession yields insights into the nature of zoning and mixing within this supervolcanic system. Timescales for zircon growth (inferred to represent longevity of the magma chamber) are investigated using U-Pb dating of zircons. Zircon textural and trace element data obtained by SIMS (SHRIMP-RG) are presented from 15 stratigraphically controlled Bishop Tuff samples and two older Glass Mountain (GM) lava samples. The resulting eruption age estimate derived from the weighted mean of 166 rim ages of 766.6±3.1 ka (95% confidence) is identical within uncertainty to published values from ID-TIMS and 40Ar/39Ar techniques. An eruption age is also derived for GM dome YA (the youngest GM dome) of 862±23 ka (95% confidence), significantly older than the widely used 790±20 ka K-Ar age. The oldest zircon cores from late-erupted Bishop material (including those with GM-type textures) have a weighted mean of 838.5±8.8 ka (95% confidence), implying that the Bishop Tuff system was only active for ~80 kyr, and had effectively no temporal overlap with the GM system. Bishop zircon textures are divided into four suites whose proportions change systematically through the eruptive sequence. Trace element variations in Bishop zircons are influenced strongly by sector zoning for many elements, and thus restrict the value of trace element variations in discerning compositional stratification within the magma chamber. In later-erupted units, bright-rim overgrowths are common, and are inferred to have crystallized from the same „bright-rim‟ magma as generated the contrasting rims seen in CL or BSE imaging on quartz, feldspar and orthopyroxene. From zircon zonation patterns, this less-evolved, slightly hotter magma invaded deeper parts of the chamber represented in the late-erupted northern units possibly up to ~10 kyr prior to eruption. In order to better quantify the timescales of interaction with the „bright-rim‟ magma, two-feldspar thermometry data are presented on multiple Bishop Tuff samples to constrain temperature variations within the pre-eruptive magma body and yield values for diffusion modelling. Two-feldspar thermometry agrees well with published Fe–Ti-oxide thermometry and reveals a ~80 °C uniform thermal gradient between the upper and lower regions of the magma chamber. Using this thermometry, diffusion of Ti in quartz, Ba in sanidine, Sr in sanidine and Fe-Mg interdiffusion in orthopyroxene are modelled to estimate timescales for the formation of overgrowth rims on crystals. Ti in quartz and Fe-Mg in orthopyroxene diffusion both yield timescales of <150 years for the formation of overgrowth rims, although differing by about an order of magnitude in their timing. However, Ba and Sr diffusion modelling in sanidine yields disparate timescales 1-2 orders of magnitude longer than for Ti in quartz. The main cause for this discrepancy is inferred to be an incorrect assumption for the initial profile shape for Ba and Sr diffusion modelling (i.e. the profile is influenced by growth zoning). Using the comparison with Sr, constraints are placed on the initial width of the core-rim interface and the initial conditions can be refined, bringing Ba and Sr diffusion timescales into mutual alignment and closer to the values from Ti in quartz. This modelling shows that piecemeal rejuvenation of lower Bishop Tuff magma chamber occurred over a period of ~500 years leading up to eruption. In situ major and trace element analyses of sanidine, plagioclase, biotite, orthopyroxene, clinopyroxene, zircon and matrix glass from the Bishop Tuff and two GM lavas are presented to investigate the pre-eruptive stratification of the Bishop magma chamber and its chemical relationship to the GM system. Analyses of samples from the entire Bishop stratigraphy confirm that the magma chamber was thermally and compositionally zoned prior to growth of crystals and the intrusion of the „bright-rim‟ forming magma. Study of rare mixed swirly and dacitic pumice samples shows enrichments in Ba, Sr and Ti (the elements responsible for bright-rim overgrowths in phenocryst phases) and identifies these pumices as possible representatives of the „bright-rim‟ magma. This integrated study of phenocrysts and glass from the Bishop Tuff leads to development of a revised magma chamber model, in which there is a unitary chamber with a stepped or sloping roof. The chamber has an upper, volumetrically dominant (~2/3) part showing no evidence for convection and with unzoned crystals, and a lower part which had experienced mixing of crystals and interaction with the „bright-rim‟ magma. Intrusion of the „bright-rim‟ magma introduced orthopyroxene and dominantly bright zircon crystals, and caused overgrowth of bright rims enriched in Ti, Sr and Ba on sanidine and quartz phenocrysts. Chemical compositions of GM and Bishop Tuff materials show a shared consanguinity, implying common modes of magma generation, yet the generation of GM and Bishop eruptible magma bodies were physically and temporally separate events.</p>

The dynamics of large-scale silicic magmatic systems: case studies from Mangakino Volcanic Centre, Taupo Volcanic Zone, New Zealand

<p>This thesis research focuses on clast and crystal-specific studies to investigate the pre- and syn-eruptive magmatic processes of two supereruptions in the TVZ: the 1.21 Ma Ongatiti (>500 km3) and the 1.0 Ma Kidnappers (~1200 km3), together with the smaller (~200 km3) 1.0 Ma Rocky Hill eruption from the Mangakino Volcanic Centre (MVC). Crystallisation histories determined through SIMS U-Pb dating of zircons reveal that the paired Kidnappers and Rocky Hill eruptions were products of a common magmatic system, which built over ~200 kyr, in the time break after the Ongatiti eruption. U-Pb age spectra from the Ongatiti show a protracted crystallisation history (over ~250 kyr), in which the majority of zircon crystallised ~100 kyr prior to eruption in a crystal mush. Zircons then ascended with melt during accumulation of the final erupted magma body in the shallow crust. Zircons remained stable in the melt dominant body but underwent little further crystallisation. Zircons from all three systems record common geochemical processes governed by the fractionating assemblage (predominantly plagioclase and amphibole). In particular, the MREE/HREE ratios and Sr concentrations of zircons from the Ongatiti record imply two contrasting source regions governed by different proportions of crystallising amphibole. The in-situ major and trace element chemistry of glass shards and crystals from the Kidnappers fall deposit reveal that magma within the Kidnappers was stored in three discrete bodies, which were systematically tapped during the early stages of eruption. Temperature and pressure (T-P) estimates from amphibole and Fe-Ti oxide equilibria from each magma type are similar and therefore the three magma bodies were adjacent, not vertically stacked, in the crust. Amphibole model T-P estimates range from 770 to 840 °C and 90 to 170 MPa corresponding to pre-eruptive storage depths of ~4.0-6.5 km. The systematic evacuation of the three independent magma bodies implies that there was tectonic triggering and linkage of eruptions. The termination of fall deposition and onset of the overlying ignimbrite emplacement marks the point of widespread caldera collapse and the catastrophic evacuation of a wider variety of melt during the Kidnappers eruption. Pumice compositions from the Kidnappers ignimbrite fall into three groups, two of which (KI-1 and KI-2) can be matched to bodies tapped during the fall phase of the eruption, with the addition of a further discrete batch of lower SiO2 (KI-3) magma. Core-rim textural and chemical variations in major crystal phases (plagioclase, amphibole and orthopyroxene) suggest each compositional group was sourced from a common mush but underwent a unique magmatic history during the development of melt-dominant bodies in the final stages prior to eruption. The field relationships and distinctive appearance of the Rocky Hill ignimbrite (~200 km3 DRE) and the underlying Kidnappers ignimbrite suggests that the two deposits are from distinct eruption events. However, major and trace element chemistry of matrix glass, coupled with the textural and chemical signatures of crystals suggests the magma erupted during the Rocky Hill was generated from the same source or mush zone as the Kidnappers. The two largest melt-dominant bodies (KI-1 and KI-2) within the Kidnappers were renewed, underwent mixing and incorporation of marginal material to form two magma types (RH-1 and RH-2) in the time break prior to the Rocky Hill eruption. Fe-Mg interdiffusion timescales in orthopyroxenes from the Kidnappers and Rocky Hill deposits suggest the establishment of the final melt-dominant bodies, through extraction of melt and crystals from a common mush, occurred within 1000 years, and peaked within centuries of each eruption. In addition, one discrete batch of Kidnappers melt has evidence for interaction with a lesser evolved melt within 50 yrs prior to eruption. This rejuvenation event was not the eruption trigger but may have primed the magma for eruption. The difference in timescales from common zones from both the Kidnappers and Rocky Hill orthopyroxene, recording the same processes reveal the time break between the two eruptions was ~20-40 years. This work highlights the rapidity of rejuvenation and renewal of the melt-dominant bodies within the Kidnappers/Rocky Hill magmatic system. The textural and in-situ compositional signatures of crystals from the Ongatiti ignimbrite imply the final erupted magma body was assembled from a thermally and chemically zoned mush, which extended to the base of the quartzofeldspathic crust (~15km). The mush was close to water saturation and was dominated by amphibole crystallisation. Melt and crystals (including the majority of zircons) were extracted from the mush and ascended to 4-6 km depths during the development of a crystal-rich (20-30%), but melt-dominant body. Significant crystallisation of plagioclase (and lesser proportions of orthopyroxene and amphibole) occurred in an event involving the gradual heating and/or increase of water in the rhyolite, from a broadly andesitic underplated magma. Homogeneous crystal rim and matrix glass compositions imply the final erupted volume of magma was effectively mixed through convection. Eu/Eu* values of whole-rock and matrix glass suggest little crystal-melt separation occurred in the melt-dominant magma body prior to eruption. This work has implications for understanding the generation, storage and eruption of large-scale silicic magma systems. The Ongatiti ignimbrite does not represent either an erupted mush, or a stratified magma chamber, suggesting an alternative model for the development of eruptible magma within large-scale silicic systems. The Kidnappers/Rocky Hill sequence records a complex interplay of multiple melt-dominant bodies, which were established and renewed on rapid timescales. The rapid timescales for the development of melt-dominant bodies and the systematic tapping of magmas in the Kidnappers/Rocky Hill system imply that tectonics may have had a strong external control on the eruptions at Mangakino.</p>

The Bishop Tuff, California: New Insights into Magmatic Timescales and Processes from Micro-Analytical Approaches

<p>The Bishop Tuff is the product of one of the largest eruptions on Earth in the last 1 Myr. This thesis studies the Bishop Tuff in order to better understand the nature of the pre-eruptive magma body, with an emphasis on the processes that occurred within it and the timescales over which they operated. In situ geochemical analyses of crystals and glass from samples collected throughout the Bishop Tuff stratigraphic succession yields insights into the nature of zoning and mixing within this supervolcanic system. Timescales for zircon growth (inferred to represent longevity of the magma chamber) are investigated using U-Pb dating of zircons. Zircon textural and trace element data obtained by SIMS (SHRIMP-RG) are presented from 15 stratigraphically controlled Bishop Tuff samples and two older Glass Mountain (GM) lava samples. The resulting eruption age estimate derived from the weighted mean of 166 rim ages of 766.6±3.1 ka (95% confidence) is identical within uncertainty to published values from ID-TIMS and 40Ar/39Ar techniques. An eruption age is also derived for GM dome YA (the youngest GM dome) of 862±23 ka (95% confidence), significantly older than the widely used 790±20 ka K-Ar age. The oldest zircon cores from late-erupted Bishop material (including those with GM-type textures) have a weighted mean of 838.5±8.8 ka (95% confidence), implying that the Bishop Tuff system was only active for ~80 kyr, and had effectively no temporal overlap with the GM system. Bishop zircon textures are divided into four suites whose proportions change systematically through the eruptive sequence. Trace element variations in Bishop zircons are influenced strongly by sector zoning for many elements, and thus restrict the value of trace element variations in discerning compositional stratification within the magma chamber. In later-erupted units, bright-rim overgrowths are common, and are inferred to have crystallized from the same „bright-rim‟ magma as generated the contrasting rims seen in CL or BSE imaging on quartz, feldspar and orthopyroxene. From zircon zonation patterns, this less-evolved, slightly hotter magma invaded deeper parts of the chamber represented in the late-erupted northern units possibly up to ~10 kyr prior to eruption. In order to better quantify the timescales of interaction with the „bright-rim‟ magma, two-feldspar thermometry data are presented on multiple Bishop Tuff samples to constrain temperature variations within the pre-eruptive magma body and yield values for diffusion modelling. Two-feldspar thermometry agrees well with published Fe–Ti-oxide thermometry and reveals a ~80 °C uniform thermal gradient between the upper and lower regions of the magma chamber. Using this thermometry, diffusion of Ti in quartz, Ba in sanidine, Sr in sanidine and Fe-Mg interdiffusion in orthopyroxene are modelled to estimate timescales for the formation of overgrowth rims on crystals. Ti in quartz and Fe-Mg in orthopyroxene diffusion both yield timescales of <150 years for the formation of overgrowth rims, although differing by about an order of magnitude in their timing. However, Ba and Sr diffusion modelling in sanidine yields disparate timescales 1-2 orders of magnitude longer than for Ti in quartz. The main cause for this discrepancy is inferred to be an incorrect assumption for the initial profile shape for Ba and Sr diffusion modelling (i.e. the profile is influenced by growth zoning). Using the comparison with Sr, constraints are placed on the initial width of the core-rim interface and the initial conditions can be refined, bringing Ba and Sr diffusion timescales into mutual alignment and closer to the values from Ti in quartz. This modelling shows that piecemeal rejuvenation of lower Bishop Tuff magma chamber occurred over a period of ~500 years leading up to eruption. In situ major and trace element analyses of sanidine, plagioclase, biotite, orthopyroxene, clinopyroxene, zircon and matrix glass from the Bishop Tuff and two GM lavas are presented to investigate the pre-eruptive stratification of the Bishop magma chamber and its chemical relationship to the GM system. Analyses of samples from the entire Bishop stratigraphy confirm that the magma chamber was thermally and compositionally zoned prior to growth of crystals and the intrusion of the „bright-rim‟ forming magma. Study of rare mixed swirly and dacitic pumice samples shows enrichments in Ba, Sr and Ti (the elements responsible for bright-rim overgrowths in phenocryst phases) and identifies these pumices as possible representatives of the „bright-rim‟ magma. This integrated study of phenocrysts and glass from the Bishop Tuff leads to development of a revised magma chamber model, in which there is a unitary chamber with a stepped or sloping roof. The chamber has an upper, volumetrically dominant (~2/3) part showing no evidence for convection and with unzoned crystals, and a lower part which had experienced mixing of crystals and interaction with the „bright-rim‟ magma. Intrusion of the „bright-rim‟ magma introduced orthopyroxene and dominantly bright zircon crystals, and caused overgrowth of bright rims enriched in Ti, Sr and Ba on sanidine and quartz phenocrysts. Chemical compositions of GM and Bishop Tuff materials show a shared consanguinity, implying common modes of magma generation, yet the generation of GM and Bishop eruptible magma bodies were physically and temporally separate events.</p>