From mush to eruption in centuries: assembly of the super-sized Oruanui magma body

The magmatic systems that give rise to voluminous crystal-poor rhyolite magma bodies can be considered to operate on two contrasting timescales: Those governed by longer-term processes by which a magma acquires its chemical and isotopic characteristics (e.g., fractional crystallisation and assimilation), and those operating at shorter timescales during the physical accumulation of the melt-dominant magma body that finally erupts. We explore the compositional and textural relationships between amphibole and orthopyroxene crystals from the 25.4 ka, 530 km³ (magma) Oruanui eruption products (Taupo volcano, New Zealand) to investigate how processes related to the physical assembly of the pre-eruptive magma body are represented in the crystal record. Over 90 % of orthopyroxenes from the volumetrically dominant high-SiO₂ (>74 wt%) rhyolite pumices record textural evidence for a significant disequilibrium event (partial dissolution ± resorption of cores and interiors) prior to the growth of 40–500 μm thick rim zones. This dissolution/regrowth history of orthopyroxene is recorded in the chemistry of co-crystallising amphiboles as a prominent inflection in the concentrations of Mn and Zn, two elements notably enriched in orthopyroxene relative to amphibole. Textural and chemical features, linked with in situ thermobarometric estimates, indicate that a major decompression event preceded the formation of the melt-dominant body. The decompression event is inferred to represent the extraction of large volumes of melt plus crystals from the Oruanui crystal mush/source zone at pressures of 140–300 MPa (~6–12 km depth). Orthopyroxene underwent partial dissolution during ascent before reestablishing in the melt-dominant magma body at pressures of 90–140 MPa (~3.5–6 km). We model Fe–Mg diffusion across the core-rim boundaries along the crystallographic a or b-axes to constrain the timing of this decompression event, which marked establishment of the melt-dominant magma body. Maximum modelled ages indicate that this event did not begin until ~1,600 years before eruption, consistent with constraints from zircon model-age spectra. Once extraction began, it underwent runaway acceleration with a peak extraction age of ~230 years, followed by an apparent period of stasis of ~60 years prior to eruption. The rapidity of the extraction and accumulation processes implies the involvement of a dynamic driving force which, in the rifted continental arc setting of the Taupo Volcanic Zone, seems likely to be represented by magma-assisted extensional tectonic processes.     

Textural and micro-analytical insights into mafic–felsic interactions during the Oruanui eruption,Taupo

Quenched juvenile mafic inclusions (enclaves) are an occasional but informative component in the deposits of large felsic eruptions. Typically, the groundmasses of these inclusions rapidly crystallize as the mafic magma is chilled against a more voluminous, cooler felsic host, providing a physical and chemical record of the nature and timing of mafic–felsic interactions. We examine mafic inclusions of two compositional lineages (tholeiitic and calc-alkaline) from deposits of the 25.4 ka Oruanui eruption (Taupo, New Zealand). 2-D quantitative textural data from analysis of back-scattered electron images reveal a marked diversity in the groundmass textures of the inclusions, including median crystal sizes (amphibole: 14–45 µm; plagioclase: 21–75 µm) and aspect ratios (amphibole: 1.7–4.2; plagioclase: 2.1–4.0), area number densities (amphibole: 122–2660 mm^?2; plagioclase: 117–2990 mm^?2), area fractions (?) of minerals (?_plag?=?23–45%, ?_amph?=?0–28%, ?_cpx?=?0–6%, ?_oxides?=?0.6–5.5%), and the relative abundance of plagioclase and amphibole (?_plag/?_amph?=?1.0–4.6). Textural parameters vary more significantly within, rather than between, the two compositional lineages, and in some cases show marked variations across individual clasts, implying that each inclusion’s cooling history, rather than bulk composition, was the dominant control on textural development. Groundmass mineral compositions are also diverse both within and between inclusions (e.g. plagioclase from An_34–92, with typical intra-clast variability of ~?20 mol%), and do not correlate with bulk chemistry. Diverse groundmass textures and mineral and glass chemistries are inferred to reflect complex interplay of a range of factors including the degree and rate of undercooling, bulk composition, water content and, possibly, intensive variables. Our data are inconsistent with breakup of a crystallizing ponded mafic layer at the base of the Oruanui melt-dominant body, instead implying that each inclusion partially crystallized as a discrete body with a unique cooling history. Extensive ingestion of mush-derived macro-crystals suggests that mechanical breakup of mafic feeder dikes occurred within a transition zone between the mush and melt-dominant magma body. In this zone, the mush lacked yield strength, as has been inferred from field studies of narrow (meters to few tens of meters) mush-melt transition zones preserved in composite intrusions. Evidence for plastic deformation of inclusions during eruption and the abundance of fresh residual glass in inclusions from all eruptive phases suggest that the inclusions formed syn-eruptively, and must have been formed recurrently at multiple stages throughout the eruption.     

Petrological imaging of an active pluton beneath Cerro Uturuncu,Bolivia

Uturuncu is a dormant volcano in the Altiplano of SW Bolivia. A present day ~70 km diameter interferometric synthetic aperture radar (InSAR) anomaly roughly centred on Uturuncu’s edifice is believed to be a result of magma intrusion into an active crustal pluton. Past activity at the volcano, spanning 0.89 to 0.27 Ma, is exclusively effusive and almost all lavas and domes are dacitic with phenocrysts of plagioclase, orthopyroxene, biotite, ilmenite and Ti-magnetite plus or minus quartz, and microlites of plagioclase and orthopyroxene set in rhyolitic groundmass glass. Plagioclase-hosted melt inclusions (MI) are rhyolitic with major element compositions that are similar to groundmass glasses. H₂O concentrations plotted versus incompatible elements for individual samples describe a trend typical of near-isobaric, volatile-saturated crystallisation. At 870 °C, the average magma temperature calculated from Fe–Ti oxides, the average H₂O of 3.2 ± 0.7 wt% and CO₂ typically <160 ppm equate to MI trapping pressures of 50–120 MPa, approximately 2–4.5 km below surface. Such shallow storage precludes the role of dacite magma emplacement into pre-eruptive storage regions as being the cause of the observed InSAR anomaly. Storage pressures, whole-rock (WR) chemistry and phase assemblage are remarkably consistent across the eruptive history of the volcano, although magmatic temperatures calculated from Fe–Ti oxide geothermometry, zircon saturation thermometry using MI and orthopyroxene-melt thermometry range from 760 to 925 °C at NNO ± 1 log. This large temperature range is similar to that of saturation temperatures of observed phases in experimental data on Uturuncu dacites. The variation in calculated temperatures is attributed to piecemeal construction of the active pluton by successive inputs of new magma into a growing volume of plutonic mush. Fluctuating temperatures within the mush can account for sieve-textured cores and complex zoning in plagioclase phenocrysts, resorption of quartz and biotite phenocrysts and apatite microlites. That Fe–Ti oxide temperatures vary by ~50–100 °C in a single thin section indicates that magmas were not homogenised effectively prior to eruption. Phenocryst contents do not correlate with calculated magmatic temperatures, consistent with crystal entrainment from the mush during magma ascent and eruption. Microlites grew during ascent from the magma storage region. Variability in the proportion of microlites is attributed to differing ascent and effusion rates with faster rates in general for lavas >0.5 Ma compared to those <0.5 Ma. High microlite contents of domes indicate that effusion rates were probably slowest in dome-forming eruptions. Linear trends in WR major and trace element chemistries, highly variable, bimodal mineral compositions, and the presence of mafic enclaves in lavas demonstrate that intrusion of more mafic magmas into the evolving, shallow plutonic mush also occurred further amplifying local temperature fluctuations. Crystallisation and resorption of accessory phases, particularly ilmenite and apatite, can be detected in MI and groundmass glass trace element covariation trends, which are oblique to WRs. Marked variability of Ba, Sr and La in MI can be attributed to temperature-controlled, localised crystallisation of plagioclase, orthopyroxene and biotite within the evolving mush.     

The 26{middle dot}5 ka Oruanui Eruption, Taupo Volcano, New Zealand: Development, Characteristics and Evacuation of a Large Rhyolitic Magma Body

The caldera-forming 26·5 ka Oruanui eruption (Taupo,New Zealand) erupted 530 km³ of magma, >99% rhyolitic, <1%mafic. The rhyolite varies from 71·8 to 76·7 wt% SiO₂ and 76 to 112 ppm Rb but is dominantly 74–76 wt% SiO₂. Average rhyolite compositions at each stratigraphiclevel do not change significantly through the eruption sequence.Oxide geothermometry, phase equilibria and volatile contentsimply magma storage at 830–760°C, and 100–200MPa. Most rhyolite compositional variations are explicable by28% crystal fractionation involving the phenocryst and accessoryphases (plagioclase, orthopyroxene, hornblende, quartz, magnetite,ilmenite, apatite and zircon). However, scatter in some elementconcentrations and ⁸⁷Sr/⁸⁶Sr ratios, and the presence of non-equilibriumcrystal compositions imply that mixing of liquids, phenocrystsand inherited crystals was also important in assembling thecompositional spectrum of rhyolite. Mafic compositions comprisea tholeiitic group (52·3–63·3 wt % SiO₂)formed by fractionation and crustal contamination of a contaminatedtholeiitic basalt, and a calc-alkaline group (56·7–60·5wt % SiO₂) formed by mixing of a primitive olivine–plagioclasebasalt with rhyolitic and tholeiitic mafic magmas. Both maficgroups are distinct from other Taupo Volcanic Zone eruptivesof comparable SiO₂ content. Development and destruction by eruptionof the Oruanui magma body occurred within 40 kyr and Oruanuicompositions have not been replicated in vigorous younger activity.The Oruanui rhyolite did not form in a single stage of evolutionfrom a more primitive forerunner but by rapid rejuvenation ofa longer-lived polygenetic, multi-age ‘stockpile’of silicic plutonic components in the Taupo magmatic system. KEY WORDS: Taupo Volcanic Zone; Taupo volcano; Oruanui eruption; rhyolite, zoned magma chamber; juvenile mafic compositions; eruption withdrawal systematics     

The petrologic evolution and pre-eruptive conditions of the rhyolitic Kos Plateau Tuff (Aegean arc)

The Kos Plateau Tuff is a large (>60 km³) and young (160 k.y.) calc-alkaline, high-SiO₂ rhyolitic ignimbrite from the active Kos-Nisyros volcanic center in the Aegean arc (Greece). Combined textural, petrological and geochemical information suggest that (1) the system evolved dominantly by crystal fractionation from (mostly unerupted) more mafic parents, (2) the magma chamber grew over ≥ 250 000 years at shallow depth (～1.5-2.5 kb) and was stored as a H₂O-rich crystalline mush close to its solidus (～670-750°C), (3) the eruption occurred after a reheating event triggered by the intrusion of hydrous mafic magma at the base of the rhyolitic mush. Rare banded pumices indicate that the mafic magma only mingled with a trivial portion of resident crystal-rich rhyolite; most of the mush was remobilized following partial melting of quartz and feldspars induced by advection of heat and volatiles from the underplated, hotter mafic influx.     

The role of magma mixing/mingling and cumulate melting in the Neapolitan Yellow Tuff caldera-forming eruption (Campi Flegrei,Southern Italy)

Understanding the mechanisms responsible for the generation of chemical gradients in high-volume ignimbrites is key to retrieve information on the processes that control the maturation and eruption of large silicic magmatic reservoirs. Over the last 60 ky, two large ignimbrites showing remarkable zoning were emplaced during caldera-forming eruptions at Campi Flegrei (i.e., Campanian Ignimbrite, CI, ~?39 ka and Neapolitan Yellow Tuff, NYT, ~?15 ka). While the CI displays linear compositional, thermal and crystallinity gradients, the NYT is a more complex ignimbrite characterized by crystal-poor magmas ranging in composition from trachy-andesites to phonolites. By combining major and trace element compositions of matrix glasses and mineral phases from juvenile clasts located at different stratigraphic heights along the NYT pyroclastic sequence, we interpret such compositional gradients as the result of mixing/mingling between three different magmas: (1) a resident evolved magma showing geochemical characteristics of a melt extracted from a cumulate mush dominated by clinopyroxene, plagioclase and oxides with minor sanidine and biotite; (2) a hotter and more mafic magma from recharge providing high-An plagioclase and high-Mg clinopyroxene crystals and (3) a compositionally intermediate magma derived from remelting of low temperature mineral phases (i.e., sanidine and biotite) within the cumulate crystal mush. We suggest that the presence of a refractory crystal mush, as documented by the occurrence of abundant crystal clots containing clinopyroxene, plagioclase and oxides, is the main reason for the lack of erupted crystal-rich material in the NYT. A comparison between the NYT and the CI, characterized by both crystal-poor extracted melts and crystal-rich magmas representing remobilized portions of a “mature” (i.e., sanidine dominated) cumulate residue, allows evaluation of the capability of crystal mushes of becoming eruptible upon recharge.     

Evolution of the Taupo Volcanic Center, New Zealand: petrological and thermal constraints from the Omega dacite

The 20 ka ~0.1 km³ Omega dacite, which erupted shortly after the 26.5 ka Oruanui super-eruption, compositionally stands out among Taupo Volcanic Zone (TVZ) magmas, which are overwhelmingly characterized by rhyolites (>90 % by volume). The previously reported presence of inherited zircons in this zircon-undersaturated magma has provided unequivocal evidence for the involvement of upper-crustal material in a 1–10 year timescale prior to the Omega eruption. However, whether this crustal involvement is characterized by wholesale, melting of preexisting crust or subordinate bulk assimilation into an already differentiated magma body remains unclear. To disentangle these processes, we describe the mineral chemistry of the major phases present in the Omega dacite and determine intensive parameters describing magma chamber conditions. Dominantly unimodal populations of plagioclase (An_50–60), orthopyroxene (Mg# from 58 to 68), and clinopyroxene (Mg# from 65 to 73), along with coexisting equilibrium pairs of Fe–Ti oxides, constrain pre-eruptive temperatures to 850–950 °C, a pressure between ~3 and 7 kbars, and an oxygen fugacity of ~NNO. MELTS thermodynamic modeling suggests that this phase assemblage is in equilibrium with the bulk rock and glass compositions of the Omega dacite at these estimated P–T–fO₂ pre-eruptive conditions. Combining these petrological observations with insights into conductive thermal models of magma–crust interactions, we argue that the Omega dacite more likely formed in the mid-to-lower crust via protracted processing through fractional crystallization coupled with some assimilation (AFC). Incorporation of crustal material is likely to have occurred at various stages, with the inherited zircons (and potentially parts of glomerocrysts) representing late and subordinate upper-crustal assimilants. This petrogenetic model is consistent with the presence of a differentiating crustal column, consisting of a polybaric fractional crystallization and assimilation history. On the basis of petrological, thermal, and geophysical considerations, upper-crustal reservoirs, which feed large-scale rhyolitic volcanism in the TVZ, most likely take the form of large, long-lived crystal mush zones. Following large eruptions, such as the Oruanui event, this mush is expected to crystallize significantly (up to 70–80 vol% crystals) due to syn-eruptive decompression. Hence, the Omega dacite, immediately post-dating the Oruanui event, potentially represents incoming deeper recharge of less-evolved magma that was able to penetrate the nearly solidified upper-crustal mush. Over the past 20,000 years, similar intermediate recharge magmas have incrementally reheated, reconstructed, and reactivated the upper-crustal mush zone, allowing a gradual return to rhyolitic volcanism at the Taupo Volcanic Center.     

Building zoned ignimbrites by recycling silicic cumulates: insight from the 1,000 km<Superscript>3</Superscript> Carpenter Ridge Tuff,CO

The ~1,000 km³ Carpenter Ridge Tuff (CRT), erupted at 27.55 Ma during the mid-tertiary ignimbrite flare-up in the western USA, is among the largest known strongly zoned ash-flow tuffs. It consists primarily of densely welded crystal-poor rhyolite with a pronounced, highly evolved chemical signature (high Rb/Sr, low Ba, Zr, Eu), but thickly ponded intracaldera CRT is capped by a more crystal-rich, less silicic facies. In the outflow ignimbrite, this upper zone is defined mainly by densely welded crystal-rich juvenile clasts of trachydacite composition, with higher Fe–Ti oxide temperatures, and is characterized by extremely high Ba (to 7,500 ppm), Zr, Sr, and positive Eu anomalies. Rare mafic clasts (51–53 wt% SiO₂) with Ba contents to 4,000–5,000 ppm and positive Eu anomalies are also present. Much of the major and trace-element variations in the CRT juvenile clasts can be reproduced via in situ differentiation by interstitial melt extraction from a crystal-rich, upper-crustal mush zone, with the trachydacite, crystal-rich clasts representing the remobilized crystal cumulate left behind by the melt extraction process. Late recharge events, represented by the rare mafic clasts and high-Al amphiboles in some samples, mixed in with parts of the crystal cumulate and generated additional scatter in the whole-rock data. Recharge was important in thermally remobilizing the silicic crystal cumulate by partially melting the near-solidus phases, as supported by: (1) ubiquitous wormy/sieve textures and reverse zoning patterns in feldspars and biotites, (2) absence of quartz in this very silicic unit stored at depths of >4–5 km, and (3) heterogeneous melt compositions in the trachydacite fiamme and mafic clasts, particularly in Ba, indicating local enrichment of this element due mostly to sanidine and biotite melting. The injection of hot, juvenile magma into the upper-crustal cumulate also imparted the observed thermal gradient to the deposits and the mixing overprint that partly masks the in situ differentiation process. The CRT provides a particularly clear perspective on processes of in situ crystal-liquid separation into a lower crystal-rich zone and an upper eruptible cap, which appears common in incrementally built upper-crustal magma reservoirs of high-flux magmatic provinces.     

Quartz zoning and the pre-eruptive evolution of the ~340-ka Whakamaru magma systems,New Zealand

Cathodoluminescence (CL) zoning in quartz crystals from rhyolitic pumices in two ignimbrite members of the ~340-ka Whakamaru super-eruption deposits, Taupo Volcanic Zone, New Zealand, is investigated in conjunction with the analysis of Ti concentration in quartz to reconstruct the history of changing magma chamber conditions and to elucidate the eruption-triggering processes. CL intensity images are taken as a proxy for Ti concentration and thus temperature and/or pressure and/or compositional variations during crystal growth history. Estimates of the maximum temperature changes (i.e., assuming other factors influencing Ti uptake remain constant) are made using the TitaniQ geothermometer based on the Ti concentration in quartz. These results are reviewed in comparison with Fe–Ti oxide, feldspar-melt and amphibole geothermometry. Core-to-rim quartz Ti profiles record a marked change in conditions (temperature increase and/or pressure decrease and/or change in melt composition) causing and then following a significant resorption horizon in the outer parts of the crystals. Two alternative models that could explain the quartz Ti zonation invoke a temperature increase caused by mafic recharge and/or a pressure decrease due to magma ponding and re-equilibration at shallow crustal levels. Concomitant changes in melt composition and Ti activity may, however, also have strongly influenced Ti uptake into the quartz. Some crystals also show other marked increases in CL brightness internally, but any accompanying magmatic changes did not result in eruption. Diffusion modelling indicates that this significant change in conditions occurred over ~10–85 years prior to caldera-forming eruption. This rapid thermal pulse or pressure change is interpreted as evidence for open-system processes, and appears to record a magma chamber recharge event that rejuvenated the Whakamaru magma system (melt-dominant magma plus crystal mush), and potentially acted as a trigger for processes that led to eruption.     

Evacuation of multiple magma bodies and the onset of caldera collapse in a supereruption,captured in glass and mineral compositions

Complexities in the nature of large-scale silicic eruptions and their magmatic systems can be discerned through micro-analytical geochemical studies. We present high-resolution, stratigraphically constrained compositional data on glassy matrix material and feldspar crystals from the initial fall deposits and earliest ignimbrite (base of member A) of the 2.08 Ma, ~?2500 km³ Huckleberry Ridge Tuff (HRT), Yellowstone. We use these data to document the nature of the magmatic system and compositional changes related to the transition from fall to widespread ignimbrite deposition, inferred to mark the onset of caldera collapse. Although major element glass compositions are relatively uniform, trace elements span a large range (e.g. Ba 10–900 ppm, Sr/Rb?=?0.005–0.09), with highly evolved glasses dominating in the fall deposits. Several trace elements (e.g. Ba and light rare earth elements) in the glass samples serve to define statistically significant compositional clustering in the fall deposits and basal ignimbrite. These clusters are inferred to reflect melt compositions controlled by fractional crystallisation processes and are interpreted to represent multiple, discrete melt-dominant domains that were tapped by multiple vents. The onset of widespread ignimbrite deposition is marked by an increase in the number of erupted melt compositional clusters from four in the fall deposits to eight, representing nine melt-dominant domains. There is an absence of geographical variation of glass compositions within the basal ignimbrite, with samples from proximal to distal localities north, west and south of the HRT caldera exhibiting similar variability. Pairing of glass analyses with sanidine major and minor element compositional data suggests that the nine melt compositional domains converged at depth into two compositionally distinct upper-crustal magmatic lineages that were both active during these early stages of the eruption. Our data collectively indicate the evacuation of an exceptionally complex and heterogeneous magma system. The onset of widespread ignimbrite deposition, inferred to relate to caldera collapse, occurred after ~ 50 km³ of magma had been discharged. Although external controls were important as an eruption trigger, depressurisation of the system led to caldera collapse with the eruption of numerous discrete melt-dominant domains.     

Large-scale magmatic layering in the Main Zone of the Bushveld Complex and episodic downward magma infiltration

The Bellevue drillcore intersects ~3 km of Main and Upper Zone cumulates in the Northern Limb of the Bushveld Complex. Main Zone cumulates are predominately gabbronorites, with localized layers of pyroxenite and anorthosite. Some previous workers, using bulk rock major, trace and isotopic compositions, have suggested that the Main Zone crystallized predominantly from a single pulse of magma. However, density measurements throughout the Bellevue drillcore reveal intervals that show up-section increases in bulk rock density, which are difficult to explain by crystallization from a single batch of magma. Wavelet analysis of the density data suggests that these intervals occur on length-scales of ~40 to ~170 m, thus defining a scale of layering not previously described in the Bushveld Complex. Upward increases in density in the Main Zone correspond to upward increases in modal pyroxene, producing intervals that grade from a basal anorthosite (with 5% pyroxene) to gabbronorite (with 30–40% pyroxene). We examined the textures and mineral compositions of a ~40 m thick interval showing upwardly increasing density to establish how this type of layering formed. Plagioclase generally forms euhedral laths, while orthopyroxene is interstitial in texture and commonly envelops finer-grained and embayed plagioclase grains. Minor interstitial clinopyroxene was the final phase to crystallize from the magma. Plagioclase compositions show negligible change up-section (average An₆₂), with local reverse zoning at the rims of cumulus laths (average increase of 2 mol%). In contrast, interstitial orthopyroxene compositions become more primitive up-section, from Mg# 57 to Mg# 63. Clinopyroxene similarly shows an up-section increase in Mg#. Pyroxene compositions record the primary magmatic signature of the melt at the time of crystallization and are not an artefact of the trapped liquid shift effect. Combined, the textures and decoupled mineral compositions indicate that the upward density increase is produced by the downward infiltration of noritic magma into a previously emplaced plagioclase-rich crystal mush. Fresh noritic magma soaked down into the crystallizing anorthositic mush, partially dissolving plagioclase laths and assimilating Fe-enriched pore melt. The presence of multiple cycles showing upward increases in density in the Bellevue drillcore suggests that downward magma infiltration occurred episodically during crystallization of the Main Zone.     

Small volume andesite magmas and melt–mush interactions at Ruapehu, New Zealand: evidence from melt inclusions

Historical eruptions from Mt. Ruapehu (New Zealand) have been small (<0.001 km³ of juvenile magma) and have often occurred without significant warning. Developing better modelling tools requires an improved understanding of the magma storage and transport system beneath the volcano. Towards that end, we have analysed the volatile content and major element chemistry of groundmass glass and phenocryst-hosted melt inclusions in erupted samples from 1945 to 1996. We find that during this time period, magma has been stored at depths of ~2–9 km, consistent with inferences from geophysical data. Our data also show that Ruapehu magmas are relatively H₂O-poor (<2 wt%) and CO₂-rich (≤1,000 ppm) compared to typical arc andesites. Surprisingly, we find that melt inclusions are often more evolved than their transporting melt (as inferred from groundmass glass compositions). Furthermore, even eruptions that are separated by less than 2 years exhibit distinct major element chemistry, which suggests that each eruption involved magma with a unique ascent history. From these data, we infer that individual melt batches rise through, and interact with, crystal mush zones formed by antecedent magmas. From this perspective, we envision the magmatic system at Ruapehu as frequently recharged by small magma inputs that, in turn, cool and crystallise to varying degrees. Melts that are able to erupt through this network of crystal mush entrain (to a greater or lesser extent) exotic crystals. In the extreme case (such as the 1996 eruption), the resulting scoria contain melt inclusion-bearing crystals that are exotic to the transporting magma. Finally, we suggest that complex interactions between recharge and antecedent magmas are probably common, but that the small volumes and short time scales of recharge at Ruapehu provide a unique window into these processes.     

Decoupled crystallization and eruption histories of the rhyolite magmatic system at Tarawera volcano revealed by zircon ages and growth rates

We obtained U–Th disequilibrium age data on zircons from each of the four rhyolite eruptions that built Tarawera volcano in the last 22 ka within the Okataina Volcanic Center (OVC), caldera, New Zealand. Secondary ion mass spectrometry analyses on unpolished euhedral crystal faces that lack resorption features show that crystal growth variously terminated from near-eruption age to ~100 ka prior to eruption. Age-depth profiling of crystals reveals long periods of continuous (~34 ka) and discontinuous growth (~90 ka). Growth hiatuses of up to ~40 ka duration occur, but do not all relate to obvious resorption surfaces. Age differences up to similar magnitude are found on opposing faces of some crystals suggesting episodes of partial exposure to melts. These features are best explained by periodic, complete, or partial, sub-solidus storage and/or inclusion in larger crystal phases, followed by rapid liberation prior to eruption. This is supported by high abundances of U and Th (~500 − >2,000 ppm) in some zircons consistent with periods of high crystallinity (>70%) in the magmatic system, based on crystal/melt partitioning. Contemporaneous but contrasting rim-ward trends of these elements within crystals, even in the same lava hand sample, require synchronous growth in separate melt bodies and little connectivity within the system, but also significant crystal transport and mixing prior to eruption. Many crystals record continuity of growth through the preceding ~60 ka OVC caldera-collapse and subsequent eruptions from Tarawera. This demonstrates a decoupling between eruption triggers, such as shallow crustal extension and mafic intrusion, and the crystallization state of the OVC silicic magmatic system. The data highlights the need to distinguish between the time for accumulation of eruptible magma and the long-term magma residence time based on the age of crystals with high closure temperatures, when assessing the potential for catastrophic eruptions.     

Amphibole perspective to unravel pre-eruptive processes and conditions in volcanic plumbing systems beneath intermediate arc volcanoes: a case study from Ciomadul volcano (SE Carpathians)

Ciomadul is the youngest volcano in the Carpathian–Pannonian region produced crystal-rich high-K dacites that contain abundant amphibole phenocrysts. The amphiboles in the studied dacites are characterized by large variety of zoning patterns, textures, and a wide range of compositions (e.g., 6.4–15 wt% Al₂O₃, 79–821 ppm Sr) often in thin-section scale and even in single crystals. Two amphibole populations were observed in the dacite: low-Al hornblendes represent a cold (<800 °C) silicic crystal mush, whereas the high-Al pargasites crystallized in a hot (>900 °C) mafic magma. Amphibole thermobarometry suggests that the silicic crystal mush was stored in an upper crustal storage (~8–12 km). This was also the place where the erupted dacitic magma was formed during the remobilization of upper crustal silicic crystal mush body by hot mafic magma indicated by simple-zoned and composite amphiboles. This includes reheating (by ~200 °C) and partial remelting of different parts of the crystal mush followed by intensive crystallization of the second mineral population (including pargasites). Breakdown textures of amphiboles imply that they were formed by reheating in case of hornblendes, suggesting that pre-eruptive heating and mixing could take place within days or weeks before the eruption. The decompression rim of pargasites suggests around 12 days of magma ascent in the conduit. Several arc volcanoes produce mixed intermediate magmas with similar bimodal amphibole cargo as the Ciomadul, but in our dacite the two amphibole population can be found even in a single crystal (composite amphiboles). Our study indicates that high-Al pargasites form as a second generation in these magmas after the mafic replenishment into a silicic capture zone; thus, they cannot unambiguously indicate a deeper mafic storage zone beneath these volcanoes. The simple-zoned and composite amphiboles provide direct evidence that significant compositional variations of amphiboles do not necessarily mean variation in the pressure of crystallization even if the Al-tschermak substitution can be recognized, suggesting that amphibole barometers that consider only amphibole composition may often yield unrealistic pressure variation.     

Geochemical and isotopic insights into the assembly,evolution and disruption of a magmatic plumbing system before and after a cataclysmic caldera-collapse eruption at Ischia volcano (Italy)

New geochemical and isotopic data on volcanic rocks spanning the period ~75–50 ka BP on Ischia volcano, Italy, shed light on the evolution of the magmatic system before and after the catastrophic, caldera-forming Monte Epomeo Green Tuff (MEGT) eruption. Volcanic activity during this period was influenced by a large, composite and differentiating magmatic system, replenished several times with isotopically distinct magmas of deep provenance. Chemical and isotopic variations highlight that the pre-MEGT eruptions were fed by trachytic/phonolitic magmas from an isotopically zoned reservoir that were poorly enriched in radiogenic Sr and became progressively less radiogenic with time. Just prior to the MEGT eruption, the magmatic system was recharged by an isotopically distinct magma, relatively more enriched in radiogenic Sr with respect to the previously erupted magmas. This second magma initially fed several SubPlinian explosive eruptions and later supplied the climactic, phonolitic-to-trachytic MEGT eruption(s). Isotopic data, together with erupted volume estimations obtained for MEGT eruption(s), indicate that >5–10 km³ of this relatively enriched magma had accumulated in the Ischia plumbing system. Geochemical modelling indicates that it accumulated at shallow depths (4–6 km), over a period of ca. 20 ka. After the MEGT eruption, volcanic activity was fed by a new batch of less differentiated (trachyte-latite) magma that was slightly less enriched in radiogenic Sr. The geochemical and Sr–Nd-isotopic variations through time reflect the upward flux of isotopically distinct magma batches, variably contaminated by Hercynian crust at 8–12 km depth. The deep-sourced latitic to trachytic magmas stalled at shallow depths (4–6 km depth), differentiated to phonolite through crystal fractionation and assimilation of a feldspar-rich mush, or ascended directly to the surface and erupted.     

Evolution of silicic magmas in the Kos-Nisyros volcanic center,Greece: a petrological cycle associated with caldera collapse

Multiple eruptions of silicic magma (dacite and rhyolites) occurred over the last ~3 My in the Kos-Nisyros volcanic center (eastern Aegean sea). During this period, magmas have changed from hornblende-biotite-rich units with low eruption temperatures (≤750–800°C; Kefalos and Kos dacites and rhyolites) to hotter, pyroxene-bearing units (>800–850°C; Nisyros rhyodacites) and are transitioning back to cooler magmas (Yali rhyolites). New whole-rock compositions, mineral chemistry, and zircon Hf isotopes show that these three types of silicic magmas followed the same differentiation trend: they all evolved by crystal fractionation and minor crustal assimilation (AFC) from parents with intermediate compositions characterized by high Sr/Y and low Nb content, following a wet, high oxygen fugacity liquid line of descent typical of subduction zones. As the transition between the Kos-Kefalos and Nisyros-type magmas occurred immediately and abruptly after the major caldera collapse in the area (the 161 ka Kos Plateau Tuff; KPT), we suggest that the efficient emptying of the magma chamber during the KPT drew out most of the eruptible, volatile-charged magma and partly solidified the unerupted mush zone in the upper crust due to rapid unloading, decompression, and coincident crystallization. Subsequently, the system reestablished a shallow silicic production zone from more mafic parents, recharged from the mid to lower crust. The first silicic eruptions evolving from these parents after the caldera collapse (Nisyros units) were hotter (up to >100°C) than the caldera-forming event and erupted from reservoirs characterized by different mineral proportions (more plagioclase and less amphibole). We interpret such a change as a reflection of slightly drier conditions in the magmatic column after the caldera collapse due to the decompression event. With time, the upper crustal intermediate mush progressively transitioned into the cold-wet state that prevailed during the Kefalos-Kos stage. The recent eruptions of the high-SiO₂ rhyolite on Yali Island, which are low temperature and hydrous phases (sanidine, quartz, biotite), suggest that another large, potentially explosive magma chamber is presently building under the Kos-Nisyros volcanic center.     

Quantitative field constraints on the dynamics of silicic magma chamber rejuvenation and overturn

The catastrophic eruption of large-volume, crystal-rich silicic magmas is often proposed to be a consequence of reheating, melting and overturn of partially molten, buoyant silicic material following repeated injection of dense, hot mafic magma. To test this “rejuvenation hypothesis”, we analyze at high spatial resolution 33 examples of deformed interfaces between intrusive mafic and silicic layers in two plutons of the Coastal Maine Magmatic Province, USA. These deformed interfaces are thought to record the buoyant overturn of silicic crystal mush layers, apparently in response to the injection and cooling of hot, dense mafic magmas. We use spectral analysis and scaling theory along with petrologic and textural data to identify, characterize, and understand periodic deformations from the scale of individual crystals (≈1 cm) to the thicknesses of mafic and silicic layers. Deformations at the largest scale lengths (>100 m) are at wavelengths comparable to, or greater than, silicic layer thicknesses and support a conjecture that mafic recharge can cause large-scale Rayleigh–Taylor-type overturning of silicic mushy layers. By contrast, the smallest scales of individual crystals probably record effects related to production and buoyancy-driven rise of melt from the tops of silicic mushes in contact with overlying hot basalt, whereas intermediate scales are explained by compaction. Our results constrain the evolution of a thermal rejuvenation event and potentially identify a condition for a large-scale overturn of the magma chamber that may lead to eruption. This work provides the first quantitative field-based constraints on some of the key physical processes inherent to the rejuvenation hypothesis.     

The 2011 eruption of Nabro volcano,Eritrea: perspectives on magmatic processes from melt inclusions

The 2011 eruption of Nabro volcano, Eritrea, produced one of the largest volcanic sulphur inputs to the atmosphere since the 1991 eruption of Mt. Pinatubo, yet has received comparatively little scientific attention. Nabro forms part of an off-axis alignment, broadly perpendicular to the Afar Rift, and has a history of large-magnitude explosive silicic eruptions, as well as smaller more mafic ones. Here, we present and analyse extensive petrological data obtained from samples of trachybasaltic tephra erupted during the 2011 eruption to assess the pre-eruptive magma storage system and explain the large sulphur emission. We show that the eruption involved two texturally distinct batches of magma, one of which was more primitive and richer in sulphur than the other, which was higher in water (up to 2.5 wt%). Modelling of the degassing and crystallisation histories demonstrates that the more primitive magma rose rapidly from depth and experienced degassing crystallisation, while the other experienced isobaric cooling in the crust at around 5 km depth. Interaction between the two batches occurred shortly before the eruption. The eruption itself was likely triggered by recharge-induced destabilisation of vertically extensive mush zone under the volcano. This could potentially account for the large volume of sulphur released. Some of the melt inclusions are volatile undersaturated, and suggest that the original water content of the magma was around 1.3 wt%, which is relatively high for an intraplate setting, but consistent with seismic studies of the Afar plume. This eruption was smaller than some geological eruptions at Nabro, but provides important insights into the plumbing systems and dynamics of off-axis volcanoes in Afar.     

Petrogenesis of Mafic Inclusions in Rhyolitic Lavas from Narugo Volcano, Northeastern Japan

Mafic inclusions present in the rhyolitic lavas of Narugo volcano,Japan, are vesiculated andesites with diktytaxitic texturesmainly composed of quenched acicular plagioclase, pyroxenes,and interstitial glass. When the mafic magma was incorporatedinto the silica-rich host magma, the cores of pyroxenes andplagioclase began to crystallize (>1000°C) in a boundarylayer between the mafic and felsic magmas. Phenocryst rim compositionsand interstitial glass compositions (average 78 wt % SiO₂) inthe mafic inclusions are the same as those of the phenocrystsand groundmass glass in the host rhyolite. This suggests thatthe host felsic melt infiltrated into the incompletely solidifiedmafic inclusion, and that the interstitial melt compositionin the inclusions became close to that of the host melt (c.850°C). Infiltration was enhanced by the vesiculation ofthe mafic magma. Finally, hybridized and density-reduced portionsof the mafic magma floated up from the boundary layer into thehost rhyolite. We conclude that the ascent of mafic magma triggeredthe eruption of the host rhyolitic magma. KEY WORDS: mafic inclusion; stratified magma chamber; magma mixing; mingling; Narugo volcano; Japan     

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.