Large Silicic Systems Research Articles

High-precision zircon age spectra record the dynamics and evolution of large open-system silicic magma reservoirs

The emplacement history and thermal evolution of subvolcanic magma reservoirs determine their longevity, size, and ability to feed volcanic eruptions. As zircon saturation is dependent on melt temperature and composition, quantitative analysis of zircon age distributions provides insight into the timing and magnitude of intensive parameter variations during the lifetime of a magma reservoir. Here we present chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon crystallization ages and in-situ zircon trace-element geochemistry for a suite of caldera-related, coeval plutonic and volcanic units of the Permian Sesia Magmatic System (northern Italy). This dataset documents the protracted growth and evolution (∼1 Myr) of a voluminous (>1000 km3), upper crustal (3.5–2.0 kbar) silicic magma reservoir. Systematic changes in zircon composition with time reveal episodic intrusions of new magma into a single, progressively differentiating reservoir which was dominantly kept at high crystallinity conditions (>60 vol.%). The volcanic and plutonic units both show dispersed and heterogeneous CA-ID-TIMS age distributions. A stochastic, thermodynamics-based zircon saturation model accounting for fractional crystallization, recharge and thermal rejuvenation effects on zircon growth and stability in a rhyolitic magma body reproduces the observed distributions of zircon U-Pb ages. Model results suggest that zircon compositional variability and crystallization age heterogeneity, characteristic of the Sesia and other large silicic systems, distinguish open- and closed-system magmatic processes. Specifically, an increase of the mass of crystallized zircons over time is the hallmark of thermally and chemically mature magma reservoirs, capable of producing voluminous eruptions upon changes in the thermal and mechanical state of the system. Conversely, a decrease in crystallized zircon mass over time characterizes systems where crystallization dominates over magma mobility, hindering substantial mass release during eruptive events. We suggest that quantitative analyses of dispersed zircon crystallization age distributions obtained with high-precision techniques can be applied to the plutonic and volcanic record to identify mature silicic magma reservoirs whose properties and storage conditions allowed for catastrophic caldera-forming eruptions. This work provides the petrologic community with a new tool to investigate the genesis and evolution of zircon-bearing intermediate to silicic intrusive magmatism and caldera-forming volcanism through time and across tectonic settings.

Mantle source characteristics and magmatic processes during the 2021 La Palma eruption

The 2021 eruption of La Palma (September 19-December 13) was the first subaerial eruption in the Canary Islands in 50 years. Approximately 0.2 km3 of lava erupted from a newly formed, broadly basaltic composite volcanic edifice on the northwestern flank of the Cumbre Vieja volcanic ridge. Comprehensive sampling of the olivine- and clinopyroxene-phyric lavas over the eruption period reveals temporal changes in mineralogy and bulk rock geochemistry from tephrite to basanite. Initial tephrite lavas have low MgO (∼6 wt.%) and elevated TiO2 (∼4 wt.%) and contain amphibole crystals and gabbroic micro-xenoliths. In contrast, lavas with progressively more mafic compositions erupted to approximately day 20 of the eruption and thereafter remained as basanite (∼8 wt.% MgO; 3.7 wt.% TiO2) until eruption termination. Temporal changes in lava chemistry reflect initial eruption of fractionated magmas that crystallized 5-10% olivine and clinopyroxene, as well as minor spinel, sulfide, and magnetite, followed by later eruption of deeper-sourced and more primitive magma. Vanadium-in-olivine oxybarometry indicates parental magmas were oxidized (fO2 = +1.5 to +2 FMQ) with 8.2 ±0.8 wt.% MgO and were generated from between 2.5-3% partial melting of a mantle source potentially containing a pyroxenite component (Xpx = 0.31 ±0.12). Day 1-20 tephrites have more radiogenic 187Os/188Os (0.143-0.148) and lower Pd, Pt, Ir and Os contents than post day 20 basanites (187Os/188Os = 0.141-0.145). Combined with available seismic data, the lavas provide a high-resolution record of eruptive evolution. Initial fractionated tephrite magma was stored in the upper lithosphere up to four years prior to eruption, consistent with pre-cursor seismicity and the presence of partially reacted amphibole and micro-xenoliths. The later lavas of the eruption were fed by more primitive basanitic parental magmas that were likely sourced from the deeper portion of the magma storage system that is underplating the island. Precursor events to the 2021 La Palma eruption involved seismicity and magma emplacement, storage and differentiation, which was followed by mobilisation, eruption, and eventual exhaustion of stored magma and partial melts. This magmatic progression is similar to that documented from the 1949 and 1971 Cumbre Vieja eruptions. Ocean islands with limited basaltic magma supply show similarities to the magmatic evolution observed in large silicic systems, where initial magma emplacement and differentiation is followed by later magma remobilisation that induces volcanic activity.

Application of a plate model to reproduce surface deformations observed at Uturuncu volcano, Bolivia

Between 1992 and the mid 2010s, strong uplift has been observed centred on Uturuncu volcano in Bolivia. More recently it was shown that a smaller-amplitude moat of subsidence surrounds the uplift. We propose that this “sombrero” pattern can be the signature of elastic-plate deformation due to a buoyant fluid stored beneath it. We take the base of the “plate” to be the brittle ductile transition which, in this region, is shallow when compared with the crustal thickness. First, we show experimentally that this type of storage can produce the observed deformation pattern. Secondly, we adapt a plate deformation model and use it to invert the surface deformation published by Henderson and Pritchard (2017). A total of 1681 inversions was performed to span the wide range of physical-parameter values. We found that our model can find a source that reproduces the vertical deformation well for a large range of parameters. The viable sources mostly correspond to a thin plate (≈5 km) or a magma with high buoyancy (possibly suggesting the presence of a volatile phase in the storage zone). Interestingly, the petrologically-determined magma storage depth before the last eruption of Uturuncu (≈4±2 km) is approximately the same as the plate thickness deduced. The radial displacement is less correctly reproduced than the vertical. The, at most, 15.6 m magma layer needed to reproduce the total observed uplift since 1992 is smaller than the typical thickness for seismic detection, which can explain why it was not detected. Storage of buoyant silicic magma at a shallow brittle-ductile transition in hot continental crust may be applicable to other large silicic systems.

Active Normal Faulting, Diking, and Doming Above the Rapidly Inflating Laguna del Maule Volcanic Field, Chile, Imaged With CHIRP, Magnetic, and Focal Mechanism Data

AbstractThe Laguna del Maule volcanic field (LdMVF) in Chile, a rapidly inflating silicic volcanic system without historical eruption, is intersected by active regional faults. The LdMVF provides an opportunity to observe how faults influence, accommodate, or are driven by an actively deforming large silicic system. Here we use Compressed High Intensity Radar Pulse (CHIRP) acoustic reflection data to map the fault network in sediments captured within the eponymous lake at the LdMVF and combine our fault maps with the volcanic history, earthquake locations, focal mechanisms, and lacustrine magnetic data to interpret how faults and magmatism interact. Our seismic data image dominantly dip‐slip faults forming grabens within the lake, subparallel to regional faults. No indications exist in the seismic data to suggest that fault patterns were created by the volcanic system, either ring or radial faults. Fault strikes interpreted from seismic and magnetic data are consistent with mapped dike and fault orientations on land. We therefore interpret that active faults at the LdMVF are tectonic rather than volcanic in origin, forming a transtensional zone that hosts the magmatic system. However, vertical motion along a NS‐striking fault near the center of uplift suggests trapdoor‐style faulting above the volcanic center in which tectonic faults are reactivated to accommodate magmatic inflation and overlying deformation. Magnetic anomalies follow regional faults, suggesting that faults also provide migration pathways. Depositional patterns indicate a prior episode of uplift followed by quiescence, indicating that significant magmatically related uplift at the LdMVF can occur without an associated major eruption.

The dynamics of large silicic systems from satellite remote sensing observations: the intriguing case of Domuyo volcano, Argentina

Silicic magmatic systems are the most dangerous volcanoes on Earth, capable of large and catastrophic eruptions, yet their low eruptive frequency makes it challenging to interpret their short-term unrest. Here we present a decade-plus analysis that integrates, for the first time, time series of satellite interferometric synthetic aperture radar (InSAR) surface deformation and satellite thermal infrared edifice-scale surface warming at a large silicic system: Domuyo volcano, in Argentina. We find that deformation and warming are highly correlated, and depending on the sign and lag between the time series, either shallow sealing or magma influx could drive Domuyo’s ongoing inflation (~ 0.15 m/year; from an InSAR-derived tabular source, ~ 11 × 8 × 1 km; ~ 6.5 km depth; ~ 0.037 km3/year volume-change rate) and warming (0.3–0.4 °C/year). This study shows the potential that combined satellite surface deformation and edifice-scale surface warming time series have on assessing the physical mechanisms of silicic volcanic systems and for constraining deterministic models.

Rapid assembly and rejuvenation of a large silicic magmatic system: Insights from mineral diffusive profiles in the Kidnappers and Rocky Hill deposits, New Zealand

The timescales over which magmas in large silicic systems are reactivated, assembled and stored remains a fundamental question in volcanology. To address this question, we study timescales from Fe–Mg interdiffusion in orthopyroxenes and Ti diffusion in quartz from the caldera-forming 1200 km3 Kidnappers and 200 km3 Rocky Hill eruptions from the Mangakino volcanic centre (Taupo Volcanic Zone, New Zealand). The two eruptions came from the same source area, have indistinguishable 40Ar/39Ar ages (∼1.0 Ma) and zircon U–Pb age spectra, but their respective deposits are separated by a short period of erosion. Compositions of pumice, glass and mineral species in the collective eruption deposits define multiple melt dominant bodies but indicate that these shared a common magmatic mush zone. Diffusion timescales from both eruptions are used to build on chemical and textural crystal signatures and interpret both the crystal growth histories and the timing of magma accumulation. Fe–Mg interdiffusion profiles in orthopyroxenes imply that the three melt-dominant bodies, established through extraction of melt and crystals from the common source, were generated within 600 years and with peak accumulation rates within 100 years of each eruption. In addition, a less-evolved melt interacted with the Kidnappers magma, beginning ∼30 years prior to and peaking within 3 years of the eruption. This interaction did not directly trigger the eruption, but may have primed the magmatic system. Orthopyroxene crystals with the same zoning patterns from the Kidnappers and Rocky Hill pumices yield consistently different diffusion timescales, suggesting a time break between the eruptions of ∼20 years (from core–rim zones) to ∼10 years (outer rim zones). Diffusion of Ti in quartz reveals similarly short timescales and magmatic residence times of <30 years, suggesting quartz is only recording the last period of crystallization within the final eruptible melt. Accumulation of the eruptible magma for these two, closely successive eruptions was accomplished over centuries to decades, in contrast to the gestation time of the magmatic system of ∼200 kyr, as indicated by zircon age patterns. The magmatic system was able to recover after the Kidnappers eruption in only ∼10–20 years to accumulate enough eruptible melt and crystals for a second ∼200 km3 eruption. Our data support concepts of large silicic systems being stored as long-lived crystal mushes, with eruptible melts generated over extraordinarily short timescales prior to eruption.

The role of mantle‐derived magmas in the isotopic evolution of Yellowstone's magmatic system

AbstractInjection of mantle‐derived magmas into the Earth's crust provides the heat necessary to develop and maintain large silicic magmatic systems. However, the role of mantle‐derived magmas in controlling the compositional evolution of large silicic systems remains poorly understood. Here we examine the role of mantle‐derived magmas in the postcaldera magmatic system at Yellowstone Plateau, the youngest magmatism associated with the Yellowstone hotspot. Using microbeam techniques, we characterize the age and Hf isotope composition of single zircon crystals hosted in rhyolites from the most recent eruptive episode at Yellowstone Plateau, which produced the Central Plateau Member rhyolites. We place these zircon data into context by comparing them to new solution Hf isotope data for the Central Plateau Member glasses, Yellowstone basalts, and potential local crustal sources. Zircons in the Central Plateau Member rhyolites record a wide range of Hf isotope compositions relative to their host melts and extend from values similar to previously erupted Yellowstone rhyolites to values similar to Yellowstone basalts. Most zircons (∼90%) are in isotopic equilibrium with their host melt, but a significant proportion show εHf values higher than their host melt, thus providing the direct evidence that silicic derivatives of mantle‐derived basalts have recharged Yellowstone's magmatic system. Mixing models confirm that the isotopic characteristics of the youngest Yellowstone rhyolites can be explained by recharge of Yellowstone's magma reservoir with silicic derivatives of underplating, mantle‐derived basalts (∼5–10% material added by mass). This process helps drive the long‐term isotopic evolution of Yellowstone's magmatic system.

Geophysical Evidence for Silicic Crustal Melt in the Continents: Where, What Kind, and How Much?

The accumulation of sizeable volumes of magma in the upper crust may produce plutons and/or result in supereruptions. Geophysical observations provide constraints on the rates, volumes, and melt distributions in magmatic systems, but they suffer from limited resolution and inherent nonuniqueness. Different, yet complementary, geophysical approaches must be combined with petrological, laboratory, and geochemical measurements. We summarize the results from such a combined approach from the central Andes. Taking a global perspective on large silicic systems reveals that several have >10% partial melt over large areas (10s of km2), and there may be localized zones with 50% or more.

Mafic replenishments into floored silicic magma chambers

Commingling between contemporaneous mafic and felsic magmas is now widely recognized in a broad range of intrusions and intrusive complexes. These interactions are important features for two main of reasons: (1) the rapidly chilled margins of mafic magma against silicic magma commonly preserve the compositions of mafic liquids, and (2) because the mafic magma solidifies rapidly, the resulting (final) configurations of mafic and felsic magmas can provide insights into physical processes and changing viscosity contrasts and rheologies of magmas and felsic crystal mush during crystallization of the mafic magma. Mingling of contrasted magmas was first recognized in the 1950s. Wider recognition of interactions between mafic and silicic liquids led to concepts of “net-veining” in the 1960s, “intramagmatic flows” (chilled basaltic layers separated by felsic cumulates) in the 1970s, and in the 1990s to “mafic-silicic layered intrusions” (MASLI), which could be as much as a few kilometers thick and more than 100 km2 in area. It was quickly appreciated that these MASLI preserved stratigraphic records of mafic replenishments into silicic magma chambers floored by felsic crystal mush. Volcanic studies had anticipated the occurrence of this last type of intrusion on the grounds that extensive ponding of basaltic magmas beneath silicic chambers was seen to be essential to keep large silicic systems like Yellowstone active for millions of years. This paper looks at the history of changing perceptions and interpretations of magma mingling and whether or not “sill complexes” are distinct from mafic-silicic layered intrusions. The stratigraphy of mafic-silicic layered intrusions records changing magmatic compositions, events, and processes in a temporal framework comparable to that provided by coeval volcanic rocks. As a result, careful study of MASLI has great potential for linking plutonic and volcanic processes and events.

Catastrophic caldera-forming eruptions II: The subordinate role of magma buoyancy as an eruption trigger

Recent analytical investigations have suggested that magma buoyancy is critical for triggering catastrophic caldera forming eruptions. Through detailed assessment of these approaches, we illustrate how analytical models have been misapplied for investigating buoyancy and are, therefore, incorrect and inconclusive. Nevertheless, the hypothesis that buoyancy is the critical trigger for larger eruptions warrants further investigation. As such, we utilize viscoelastic finite element models that incorporate buoyancy to test overpressure evolution and mechanical failure in the roof due to the coalescence of large buoyant magma bodies for two model cases. In the first case, we mimic empirical approaches and include buoyancy as an explicit boundary condition. In the second set of models, buoyancy is calculated implicitly due to the density contrast between the magma in the reservoir and the host rock. Results from these numerical experiments indicate that buoyancy promotes only minimal overpressurization of large silicic magma reservoirs (<0.1MPa). Furthermore, no Mode-1 tensile failure is predicted along the magma chamber boundary due to buoyancy in large reservoirs. Rather, compressional stresses are observed due to buoyant magma focusing away from the edges of the reservoir and toward the center. Given the shortcomings of the analytical implementations and the results from the numerical experiments, we conclude that buoyancy does not provide an eruption triggering mechanism for large silicic systems. Therefore, correlations of buoyancy with magma residence times, the eruption frequency–volume relationship, and the dimensions of calderas are re-assessed. We find a causal relationship with magma reservoir volume that implicates the mechanical conditions of the host rock as a primary control on eruption frequency. As magma reservoirs grow in size (>100km3) they surpass a rheological threshold where their subsequent evolution is controlled by host rock mechanics. Consequently, this results in a thermomechanical division between small systems that are triggered “internally” by magmatic processes and large systems that are triggered “externally” by faulting related to roof uplift or tectonism. Finally, critical assessment of recent analytical approaches illustrates that care must be used when applying previously derived analytical solutions to ensure that assumptions used in the original formulation are not violated during application to new geologic problems.

Timescales of Magma Recharge and Reactivation of Large Silicic Systems from Ti Diffusion in Quartz

Timescales of magma chamber assembly and recharge are investigated here by applying 1D and 2D diffusion modeling techniques to high-resolution maps of titanium in quartz from a large-volume ignimbrite eruption in the Taupo Volcanic Zone, New Zealand. We compare quartz zonation patterns and associated diffusion timescales from the 340 kaWhakamaru super-eruption (magma volume 1000 km) with the Younger Toba Tuff super-eruption, 74 ka (2000 km), Sumatra, and the smaller volume 50 ka Earthquake Flat eruption (10 km), Okataina Caldera Complex, New Zealand. Two principal timescales are presented: that of chamber recharge and eruption triggering events, and that of magma generation (involving long-term assembly, stirring and reactivation). Synchrotron micro-X-ray fluorescence maps of core^rim quartz transects provide a high-resolution record of magma chamber conditions throughout quartz crystallization. Quartz crystals from the Whakamaru magma display complex zonation patterns indicating fluctuating pressure^temperature conditions throughout the crystallization history. Toba and Earthquake Flat, in contrast, display simple quartz-zoning patterns and record slightly longer periods of crystal residence in the chamber that fed the eruption. We apply Lattice Boltzmann 2D diffusion modeling to reconstruct the timescales of quartz crystal zonation, accounting for crystal boundary complexities. Quartz crystal orientation is also accounted for by using geometry constraints from the synchrotron data. Our calculations suggest that crystal-mush reactivation for the main Whakamaru magma reservoir occured over a period of the order of 10^10 years. Both the Earthquake Flat and Toba eruptions experienced a significant recharge event (causing a temperature and pressure change), which occurred within 100 years of eruption. In comparison, the complex Whakamaru quartz zoning patterns suggest that the magma body experienced numerous thermal and compositional fluctuations in the lead-up to eruption. The final magma recharge event, which most probably triggered the eruption, occurred within 10^60 years of the eruption. Even though the volume of these systems spans two orders of magnitude, there does not appear to be a relationship between magma volume and diffusion timescale, suggesting similar histories before eruption.

Decadal to monthly timescales of magma transfer and reservoir growth at a caldera volcano

Caldera-forming volcanic eruptions are low-frequency, high-impact events capable of discharging tens to thousands of cubic kilometres of magma explosively on timescales of hours to days, with devastating effects on local and global scales. Because no such eruption has been monitored during its long build-up phase, the precursor phenomena are not well understood. Geophysical signals obtained during recent episodes of unrest at calderas such as Yellowstone, USA, and Campi Flegrei, Italy, are difficult to interpret, and the conditions necessary for large eruptions are poorly constrained. Here we present a study of pre-eruptive magmatic processes and their timescales using chemically zoned crystals from the 'Minoan' caldera-forming eruption of Santorini volcano, Greece, which occurred in the late 1600s BC. The results provide insights into how rapidly large silicic systems may pass from a quiescent state to one on the edge of eruption. Despite the large volume of erupted magma (40-60 cubic kilometres), and the 18,000-year gestation period between the Minoan eruption and the previous major eruption, most crystals in the Minoan magma record processes that occurred less than about 100 years before the eruption. Recharge of the magma reservoir by large volumes of silicic magma (and some mafic magma) occurred during the century before eruption, and mixing between different silicic magma batches was still taking place during the final months. Final assembly of large silicic magma reservoirs may occur on timescales that are geologically very short by comparison with the preceding repose period, with major growth phases immediately before eruption. These observations have implications for the monitoring of long-dormant, but potentially active, caldera systems.

Geochemical homogeneity of a long-lived, large silicic system; evidence from the Cerro Galán caldera, NW Argentina

By applying a number of analytical techniques across a spectrum of spatial scales (centimeter to micrometer) in juvenile components, we show that the Cerro Galan volcanic system has repeatedly erupted magmas with nearly identical geochemistries over >3.5 Myr. The Cerro Galan system produced nine ignimbrites (∼5.6 to 2 Ma) with a cumulative volume of >1,200 km3 (DRE; dense rock equivalent) of calc-alkaline, high-K rhyodacitic magmas (68–71 wt.% SiO2). The mineralogy is broadly constant throughout the eruptive sequence, comprising plagioclase, quartz, biotite, Fe–Ti oxides, apatite, and titanite. Early ignimbrite magmas also contained amphibole, while the final eruption, the most voluminous Cerro Galan ignimbrite (CGI; 2.08 ± 0.02 Ma) erupted a magma containing rare amphibole, but significant sanidine. Each ignimbrite contains two main juvenile clast types; dominant “white” pumice and ubiquitous but subordinate “grey” pumice. Fe–Ti oxide and amphibole-plagioclase thermometry coupled with amphibole barometry suggest that the grey pumice originated from potentially hotter and deeper magmas (800–840°C, 3–5 kbar) than the more voluminous white pumice (770–810°C, 1.5–2.5 kbar). The grey pumice is interpreted to represent the parental magmas to the Galan system emplaced into the upper crust from a deeper storage zone. Most inter-ignimbrite variations can be accounted for by differences in modal mineralogy and crystal contents that vary from 40 to 55 vol.% on a vesicle-free basis. Geochemical modeling shows that subtle bulk-rock variations in Ta, Y, Nb, Dy, and Yb between the Galan ignimbrites can be reconciled with differences in amounts of crystal fractionation from the “grey” parent magma. The amount of fractionation is inversely correlated with volume; the CGI (∼630 km3) and Real Grande Ignimbrite (∼390 km3) return higher F values (proportion of liquid remaining) than the older Toconquis Group ignimbrites (<50 km3), implying less crystal fractionation took place during the upper-crustal evolution of these larger volume magmas. We attribute this relationship to variations in magma chamber geometry; the younger, largest volume ignimbrites came from flat sill-like magma chambers, reducing the relative proportion of sidewall crystallization and fractionation compared to the older, smaller-volume ignimbrite eruptions. The grey pumice clasts also show evidence of silicic recharge throughout the history of the Cerro Galan system, and recharge days prior to eruption has previously been suggested based on reversely zoned (OH and Cl) apatite phenocrysts. A rare population of plagioclase phenocrysts with thin An-rich rims in juvenile clasts in many ignimbrites supports the importance of recharge in the evolution and potential triggering of eruptions. This study extends the notion that large volumes of nearly identical silicic magmas can be generated repeatedly, producing prolonged geochemical homogeneity from a long-lived magma source in a subduction zone volcanic setting. At Cerro Galan, we propose that there is a zone between mantle magma input and upper crustal chambers, where magmas are geochemically “buffered”, producing the underlying geochemical and isotopic signatures. This produces the same parental magmas that are delivered repeatedly to the upper crust. A lower-crustal MASH (melting, assimilation, storage, and homogenization) zone is proposed to act as this buffer zone. Subsequent upper crustal magmatic processes serve only to slightly modify the geochemistry of the magmas.

Quartz crystal size distribution of the Ordovician Millbrig K-bentonite

Crystal size distribution (CSD) analysis has been applied to quartz crystals of the Ordovician Millbrig K-bentonite, which represents one of the largest known fallout ash deposits in the Phanerozoic Era, to establish crystal growth histories and conditions in the magma chamber prior to eruption. Specific CSDs were examined for crystal growth conditions. All Millbrig samples exhibit lognormal concave-down shapes in agreement with previously reported CSDs on large silicic systems but in contrast to more mafic systems characterized by linear CSDs. Concave-down CSDs are possibly an intrinsic property of phenocryst growth in large volume silicic magmas as documented in the literature [Bindeman, I.N., 2003. Crystal sizes in evolving silicic magma chambers. Geology, 31: 367–370.]. Quartz crystals of the Millbrig must have been produced by surface-controlled, size-dependent growth with a minor Ostwald ripening component, which could have easily changed to supply-controlled growth due to a limited supply of nutrients, as the crystals grew larger. CSDs were also used to examine the possible origin of multiple ash beds in the Millbrig. These ash beds appear to be a product of a series of separate eruptions that represent separate magma layers or batches, each with slightly different crystal growth conditions.

Gas percolation in upper-crustal silicic crystal mushes as a mechanism for upward heat advection and rejuvenation of near-solidus magma bodies

Several crystal-rich, intermediate to silicic magmas erupted at arc volcanoes record a reheating event shortly prior to eruption: they provide evidence for remobilization of crystal mushes by mafic magmas. As hybridization between the mush and the mafic magma is often limited, bulk mixing could not be the dominant process in transferring heat. Conductive heating from a basaltic underplate plays a role, but a few characteristics of these rejuvenated mushes suggest that reheating occurs faster than predicted by conduction. In the upper crust, a process that can transport heat faster than conduction, and still remain chemically nearly imperceptible, is the upward migration of a hot volatile phase (“gas sparging”) that originates in underplated mafic magmas. Using numerical simulations, we quantified the thermal effects of two-phase flow (a silicic melt phase and a H 2O–CO 2 fluid phase) in the pore space of shallow silicic mushes that have reached their rheological lock-up point (i.e., rigid porous medium, crystallinity ≥ 50 vol.%). Results show that the reheating rates are significantly faster than conduction for volatile fluxes > 0.1 m 3/m 2 yr. Considering that volatiles can be rapidly exsolved from the underplated mafic magma, these high fluxes can be promptly reached, leading to fast reheating; sill-like batches of mushes with volumes similar to the 1995–present eruption of the Soufrière Hills (Montserrat, W.I.) can be reheated by a few tens of degrees and remobilized within days to weeks. At these high fluxes, a considerable volume of volatiles is needed (similar to the volume of mush being reheated). Large silicic systems (> 100–1000 km 3) require unrealistic amounts of volatiles to be reheated in a continuous, high-flux sparging event. Rejuvenation of batholithic mushes therefore requires multiple sparging episodes separated by periods dominated by near-conductive heat transfer at low-flux sparging (< 0.1 m 3/m 2 yr) and may take up to 100–200 ky.

Time scales of large volume silicic magma systems: Sr isotopic systematics of phenocrysts and glass from the Bishop Tuff, Long Valley, California

The initial Sr isotopic compositions of glass and mineral separates from the 0.74 Ma Bishop Tuff ashflow in eastern California were determined to investigate the time scales of magmatic processes in a large silicic system. It was found that there is substantial isotopic heterogeneity, both between eruptive units and between glass and phenocryst phases of individual units. The frist-erupted, lower temperature units generally have higher initial 87Sr/86Sr than later crupted, higher temperature units. Within each unit, feldspar phenocrysts have the lowest 87Sr/86Sr, associated glass has higher 87Sr/86Sr, and biotite phenocrysts have the highest 87Sr/86Sr. These isotopic differences were produced by processes in the magma chamber and not by post-eruptive alteration. Two samples were similar Fe−Ti oxide temperatures but from widely separated localities have nearly identifical Sr isotopic characteristics, indicating the existence of compositionally uniform layers of substantial volume within the chamber. Trace element data indicate that the feldspars crystallized from a liquid represented by the associated glass, and that the feldspar-glass pairs are not accidental. The rhyolitic liquids of the Bishoptuff magma chamber apparently experienced increasing 87Sr/86Sr at a rate too fast for feldspar phenocrysts to remain in isotopic equilibrium. The increasing 87Sr/86Sr is caused primarily by radioactive decay of 87Rb in the high-Rb/Sr liquids and not primarily by assimilation of radiogenic wall-rock material. A self-consistent model can be constructed to account for all of the isotopic data except for those on biotite phenocrysts. The time scale for evolution of the system is bounded on the high side at about 500 ky by observations made on precaldera lavas, and on the low side at approximately 300 ky by the time necessary to establish homogeneous layers in an actively differentiating chamber. The deduced time scale is consistent with model Rb−Sr ages, which date the differentiation of low temperature liquids from higher temperature liquids, and is compatible with the observed isotopic disequilibrium between feldspars and glass because of the low diffusivity of Sr in fieldspars ( 3x105 years. The silicic-magma production rate within the Bishop Tuff magma chamber is estimated to be 10-3km3/y. The growth rate of alkali feldspar is estimated to be about 10-14 cm/s based on the Sr isotopic difference between sanidine and glass of the lower Bishop Tuff. The biotite population is inferred to be partially (>50 ppm) xenocrystic, the xenocrysts being introduced to the chamber less than one year prior to eruption.

Gradients in physical parameters in zoned felsic magma bodies: Implications for evolution and eruptive withdrawal

Abstract Five diverse, well documented, chemically zoned magmas have been chosen from the literature to demonstrate the extent and patterns of density and viscosity gradients in zoned magma chambers. The patterns are used to assess implications for development of zonation, and withdrawal dynamics and preservation of systematic chemical variations in the final pyroclastic deposit. These examples are: Bishop Tuff, California (high-silica rhyolite); Los Humeros, Mexico (calc-alkaline rhyolite to andesite); Fogo A, Azores (trachyte); Laacher See, Eifel (phonolite) and Tenerife, Canary Islands (phonolite). It was necessary to make several simplifying assumptions in order to calculate viscosity and density profiles through each system; results are particularly sensitive to magmatic water and crystal contents. Nevertheless, the following conclusions can be drawn: 1. (1) Small, strongly zoned, alkaline magma systems which evolved through fractional crystallisation of a basaltic parent (Fogo A, Laacher See) have suffered a partial time-integrated volatile depletion prior to eruption. The most likely mechanism of volatile loss is degassing of the uppermost, highly differentiated, “cupola” magma layer. 2. (2) Eruption withdrawal dynamics are critically dependent on density gradients (and therefore on volatile content and phenocryst abundance), while viscosity variations play a subordinate role in the chosen examples. 3. (3) Formation of a chemically zoned tephra sequence by eruption of chemically zoned felsic magma requires a pre-eruptive volatile gradient in the magma chamber. 4. (4) Withdrawal-layer thicknesses during eruptions from naturally zoned magma chambers are of the order of 100 m. 5. (5) The quantitative treatment of gravitational liquid segregation processes by Nilson et al. (1985) successfully predicts times required for zonation of magma bodies: typically 103–104 years for small alkaline systems, and > 105 years for large silicic systems.