Magnetite scavenging and the buoyancy of bubbles in magmas. Part 2: Energetics of crystal-bubble attachment in magmas

The fate of pre-eruptive bubbles depends largely on their buoyancy, which can be strongly modified by the presence of crystals attached to the bubble–melt interface. We define the attachment energy and attachment force as those resulting from the attachment of a crystal to a bubble. The attachment energy is such that (1) attachment of crystals to bubbles is always favored energetically, and (2) oxide minerals attach to bubbles much more strongly than silicates, because the attachment energy is a strong function of the wetting angle. Attaching crystals to bubbles can cause bubble–crystal pairs to become neutrally buoyant. There is a critical bubble radius below which the attachment force will be strong enough to keep the pair together; we show that crystals as large as 1 mm in diameter can form neutrally buoyant pairs. For early erupted Bishop magma, if all magnetite forms neutrally buoyant pairs with gas bubbles, ca. 0.1–0.2 vol% gas can be stored in the magma; 2–3 vol% of gas can be accounted for if all minerals form neutrally buoyant aggregates. These values are an order of magnitude lower than what is inferred from melt inclusions. Hence, both magnetite-free and magnetite-rich bubbles might have existed, but only a very small fraction of them could have been neutrally buoyant. Importantly, an intrinsic association between magnetite crystals and bubbles is expected. However, most magnetite crystals in the early erupted Bishop are free of bubbles; the puzzling conclusion is that nucleation away from crystals is favored over heterogeneous nucleation on crystal substrates. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Compositional Zoning of the Bishop Tuff

Compositional data for >400 pumice clasts, organized accordingto eruptive sequence, crystal content, and texture, providenew perspectives on eruption and pre-eruptive evolution of the>600 km³ of zoned rhyolitic magma ejected as the Bishop Tuffduring formation of Long Valley caldera. Proportions and compositionsof different pumice types are given for each ignimbrite packageand for the intercalated plinian pumice-fall layers that eruptedsynchronously. Although withdrawal of the zoned magma was lesssystematic than previously realized, the overall sequence displaystrends toward greater proportions of less evolved pumice, morecrystals (0·5–24 wt %), and higher FeTi-oxide temperatures(714–818°C). No significant hiatus took place duringthe 6 day eruption of the Bishop Tuff, nearly all of which issuedfrom an integrated, zoned, unitary reservoir. Shortly beforeeruption, however, the zoned melt-dominant portion of the chamberwas invaded by batches of disparate lower-silica rhyolite magma,poorer in crystals than most of the resident magma but slightlyhotter and richer in Ba, Sr, and Ti. Interaction with residentmagma at the deepest levels tapped promoted growth of Ti-richrims on quartz, Ba-rich rims on sanidine, and entrapment ofnear-rim melt inclusions relatively enriched in Ba and CO₂.Varied amounts of mingling, even in higher parts of the chamber,led to the dark gray and swirly crystal-poor pumices sparselypresent in all ash-flow packages. As shown by FeTi-oxide geothermometry,the zoned rhyolitic chamber was hottest where crystal-richest,rendering any model of solidification fronts at the walls orroof unlikely. The main compositional gradient (75–195ppm Rb; 0·8–2·2 ppm Ta; 71–154 ppmZr; 0·40–1·73% FeO*) existed in the melt,prior to crystallization of the phenocryst suite observed, whichincluded zircon as much as 100 kyr older than the eruption.The compositions of crystals, though themselves largely unzoned,generally reflect magma temperature and the bulk compositionalgradient, implying both that few crystals settled or were transportedfar and that the observed crystals contributed little to establishingthat gradient. Upward increases in aqueous gas and dissolvedwater, combined with the adiabatic gradient (for the 5 km depthrange tapped) and the roofward decline in liquidus temperatureof the zoned melt, prevented significant crystallization againstthe roof, consistent with dominance of crystal-poor magma earlyin the eruption and lack of any roof-rind fragments among theBishop ejecta, before or after onset of caldera collapse. Amodel of secular incremental zoning is advanced wherein numerousbatches of crystal-poor melt were released from a mush zone(many kilometers thick) that floored the accumulating rhyoliticmelt-rich body. Each batch rose to its own appropriate levelin the melt-buoyancy gradient, which was self-sustaining againstwholesale convective re-homogenization, while the thick mushzone below buffered it against disruption by the deeper (non-rhyolitic)recharge that augmented the mush zone and thermally sustainedthe whole magma chamber. Crystal–melt fractionation wasthe dominant zoning process, but it took place not principallyin the shallow melt-rich body but mostly in the pluton-scalemush zone before and during batchwise melt extraction. KEY WORDS: Bishop Tuff; ignimbrite; magma zonation; mush model; rhyolite     

天池火山三期浮岩气孔形态的复杂性及其动力学成因

爆炸式喷发过程中,火山碎屑物气孔记录了挥发分出溶、膨胀和合并等信息,其大小、形态、数量密度、空间分布等局域特征是推断火山喷发动力条件的重要参考。文章基于天池火山三期喷发(50 000年前大喷发的黄色浮岩、千年大喷发的灰白色浮岩和1668年八卦庙期喷发的黑色浮岩)野外地质工作,以非线性火山喷发动力学为指导,开展了火山通道内气泡生长的流体动力学研究,揭示出岩浆流体黏性力和界面张力的共同作用对于岩浆减压和气泡生长过程的约束。在浮岩气孔结构的定量化分析基础上,进一步研究了天池火山三期喷发的浮岩气孔参数,通过气泡生长流体动力学方程得到了千年大喷发灰白色浮岩毛细管数Ca值为253, 明显高于50 000年前大喷发黄色浮岩(Ca值为94)和八卦庙期喷发黑色浮岩(Ca值为111),表明了千年大喷发曾发生过明显的成分变化,推测可能与幔源基性岩浆注入有关;而50 000年前大喷发黄色浮岩气孔不规则形态参数(1-Ω)值为0.098,大于后两期喷发(分别为0.052和0.064),可能意味着天池火山系统动力学平衡的弛豫周期变小或浮岩气泡生长受动力学、流变学改造过程减弱,这可为进一步研究天池火山活动规律提供参考。三期浮岩毛细管数Ca量级为102,气孔不规则形态参数(1-Ω)量级为10-1,从动力学上首次证实了天池火山属于普林尼型或超普林尼型喷发。     

Kilauea East Rift Zone Magmatism: an Episode 54 Perspective

On January 29–30, 1997, prolonged steady-state effusionof lava from Pu'u'O'o was briefly disrupted by shallow extensionbeneath Napau Crater, 1–4 km uprift of the active Kilaueavent. A 23-h-long eruption (episode 54) ensued from fissuresthat were overlapping or en echelon with eruptive fissures formedduring episode 1 in 1983 and those of earlier rift zone eruptionsin 1963 and 1968. Combined geophysical and petrologic data forthe 1994–1999 eruptive interval, including episode 54,reveal a variety of shallow magmatic conditions that persistin association with prolonged rift zone eruption. Near-ventlava samples document a significant range in composition, temperatureand crystallinity of pre-eruptive magma. As supported by phenocryst–liquidrelations and Kilauea mineral thermometers established herein,the rift zone extension that led to episode 54 resulted in mixtureof near-cotectic magma with discrete magma bodies cooled to     

Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

The behavior of chlorine and sulfur during differentiation of the Mt. Somma-Vesuvius magmatic system

Summary Reheated silicate melt inclusions in volcanic rock samples from Mt. Somma-Vesuvius, Italy, have been analyzed for 29 constituents including H₂O, S, Cl, F, B, and P₂O₅. This composite volcano consists of the older Mt. Somma caldera, formed between 14 and 3.55 ka before present, and the younger Vesuvius cone. The melt inclusion compositions provide important constraints on pre-eruptive magma geochemistry, identify relationships that relate to eruption behavior and magma evolution, and provide extensive evidence for magmatic fluid exsolution well before eruption. The melt inclusion data have been categorized by groups that reflect magma compositions, age, and style of eruptions. The data show distinct differences in composition for eruptive products older than 14.0 ka (pre-caldera rocks) versus eruptive products younger than 3.55 ka. Moreover, pre-caldera eruptions were associated with magmas relatively enriched in SiO₂, whereas eruptions younger than 3.55 ka (i.e., the syn- and post-caldera magmas which generated the Somma caldera and the Vesuvius cone) were derived from magmas comparatively enriched in S, Cl, CaO, MgO, P₂O₅, F, and many lithophile trace elements. Melt inclusion data indicate that eruptive behavior at Vesuvius correlates with pre-eruptive volatile enrichments. Most magmas associated with explosive plinian and subplinian events younger than 3.55 ka contained more H₂O, contained significantly more S, and exhibited higher (S/Cl) ratios than syn- and post-caldera magmas which erupted during relatively passive interplinian volcanic phenomena. Received January 10, 2000 Revised version accepted July 17, 2000     

Multiple timescale constraints for high-flux magma chamber assembly prior to the Late Bronze Age eruption of Santorini (Greece)

The rhyodacitic magma discharged during the 30–80 km³ DRE (dense rock equivalent) Late Bronze Age (LBA; also called ‘Minoan’) eruption of Santorini caldera is known from previous studies to have had a complex history of polybaric ascent and storage prior to eruption. We refine the timescales of these processes by modelling Mg–Fe diffusion profiles in orthopyroxene and clinopyroxene crystals. The data are integrated with previously published information on the LBA eruption (phase equilibria studies, melt inclusion volatile barometry, Mg-in-plagioclase diffusion chronometry), as well as new plagioclase crystal size distributions and the established pre-LBA history of the volcano, to reconstruct the events that led up to the assembly and discharge of the LBA magma chamber. Orthopyroxene, clinopyroxene and plagioclase crystals in the rhyodacite have compositionally distinct rims, overgrowing relict, probably source-derived, more magnesian (or calcic) cores, and record one or more crystallization (plag???opx?>?cpx) events during the few centuries to years prior to eruption. The crystallization event(s) can be explained by the rapid transfer of rhyodacitic melt from a dioritic/gabbroic region of the subcaldera pluton (mostly in the 8–12 km depth range), followed by injection, cooling and mixing in a large melt lens at 4–6 km depth (the pre-eruptive magma chamber). Since crystals from all eruptive phases yield similar timescales, the melt transfer event(s), the last of which took place less than 2 years before the eruption, must have involved most of the magma that subsequently erupted. The data are consistent with a model in which prolonged generation, storage and segregation of silicic melts were followed by gravitational instability in the subcaldera pluton, causing the rapid interconnection and amalgamation of melt-rich domains. The melts then drained to the top of the pluton, at fluxes of up to 0.1–1 km³ year^??1, where steep vertical gradients of density and rheology probably caused them to inject laterally, forming a short-lived holding chamber prior to eruption. This interpretation is consistent with growing evidence that some large silicic magma chambers are transient features on geological timescales. A similar process preceded at least one earlier caldera-forming eruption on Santorini, suggesting that it may be a general feature of this rift-hosted magmatic system.     

Development of the Long Valley,California, magma chamber recorded in precaldera rhyolite lavas of Glass Mountain

Glass Mountain, California, consists of >50 km³ of high-silica rhyolite lavas and associated pyroclastic deposits that erupted over a period of >1 my preceding explosive eruption of the Bishop Tuff and formation of the Long Valley caldera at 0.73 Ma. These “minimum-melt” rhyolites yield Fe-Ti-oxide temperatures of 695–718°C and contain sparse phenocrysts of plagioclase+quartz+magnetite+apatite±sanidine, biotite, ilmenite, allanite, and zircon. Incompatible trace elements show similar or larger ranges within the Glass Mountain suite than within the Bishop Tuff, despite a much smaller range of major-element concentrations, largely due to variability among the older lavas (erupted between 2.1 and 1.2 Ma). Ratios of the most incompatible elements have larger ranges in the older lavas than in the younger lavas (1.2–0.79 Ma), and concentrations of incompatible elements span wide ranges at nearly constant Ce/Yb, suggesting that the highest concentrations of these elements are not the result of extensive fractional crystallization alone; rather, they are inherited from parental magmas with a larger proportion of crustal partial melt. Evidence for the nature of this crustal component comes from the presence of scarce, tiny xenocrysts derived from granitic and greenschist-grade metamorphic rocks. The wider range of chemical and isotopic compositions in the older lavas, the larger range in phenocryst modes, the eruption of magmas with different compositions at nearly the same time in different parts of the field, and the smaller volume of individual lavas suggest either that more than one magma body was tapped during eruption of the older lavas or that a single chamber tapped by all lavas was small enough that the composition of its upper reaches easily affected by new additions of crustal melts. We interpret the relative chemical, mineralogical, and isotopic homogeneity of the younger Glass Mountain lavas as reflecting eruptions from a large, integrated magma chamber. The small number of cruptions between 1.4 and 1.2 ma may have allowed time for a large magma body to coalesce, and, as the chamber grew, its upper reaches became less affected by new inputs of crustal melts, so that trace-element trends in magmas erupted after 1.2 Ma are largely controlled by fractional crystallization. The extremely low Sr concentrations of Glass Mountain lavas imply extensive crystallization in chambers at least hundreds of cubic kilometers in volume. The close similarity in Sr, Nd, and Pb isotopic ratios between the younger Glass Mountain lavas and unaltered Bishop Tuff indicates that they tapped the same body of magma, which had become isotopically homogenous by 1.2 Ma but continued to differentiate after that time. From 1.2 to 0.79 Ma, volumetric eruptive rates may have exceeded rates of differentiation, as younger Glass Mountain lavas become slightly less evolved with time. Early-erupted Bishop Tuff is more evolved than the youngest of the Glass Mountain lavas and is characterized by slightly different trace element ratios. This suggests that although magma had been present for 0.5 my, the composiional gradient exhibited by the Bishop Tuff had not been a long-term, steady-state condition in the Long Valley magma chamber, but developed at least in part during the 0.06-my hiatus between extrusion of the last Glass Mountain lava and the climactic eruption.     

The Cougar Point Tuff: Implications for Thermochemical Zonation and Longevity of High-Temperature, Large-Volume Silicic Magmas of the Miocene Yellowstone Hotspot

The 12·7–10·5 Ma Cougar Point Tuff in southernIdaho, USA, consists of 10 large-volume (>10²–10³ km³each), high-temperature (800–1000°C), rhyolitic ash-flowtuffs erupted from the Bruneau–Jarbidge volcanic centerof the Yellowstone hotspot. These tuffs provide evidence forcompositional and thermal zonation in pre-eruptive rhyolitemagma, and suggest the presence of a long-lived reservoir thatwas tapped by numerous large explosive eruptions. Pyroxene compositionsexhibit discrete compositional modes with respect to Fe andMg that define a linear spectrum punctuated by conspicuous gaps.Airfall glass compositions also cluster into modes, and thepresence of multiple modes indicates tapping of different magmavolumes during early phases of eruption. Equilibrium assemblagesof pigeonite and augite are used to reconstruct compositionaland thermal gradients in the pre-eruptive reservoir. The recurrenceof identical compositional modes and of mineral pairs equilibratedat high temperatures in successive eruptive units is consistentwith the persistence of their respective liquids in the magmareservoir. Recurrence intervals of identical modes range from0·3 to 0·9 Myr and suggest possible magma residencetimes of similar duration. Eruption ages, magma temperatures,Nd isotopes, and pyroxene and glass compositions are consistentwith a long-lived, dynamically evolving magma reservoir thatwas chemically and thermally zoned and composed of multiplediscrete magma volumes. KEY WORDS: ash-flow tuff; Bruneau–Jarbidge; rhyolite; Yellowstone hotspot; residence time     

Magma Evolution and Open-System Processes at Shiveluch Volcano: Insights from Phenocryst Zoning

Phenocryst zoning patterns are used to identify open-systemmagmatic processes in the products of the 2001 eruption of ShiveluchVolcano, Kamchatka. The lavas and pumices studied are hornblende–plagioclaseandesites with average pre-eruptive temperatures of 840°Cand fO₂ of 1·5–2·1 log units above nickel–nickeloxide (NNO). Plagioclase zoning includes oscillatory and patchyzonation and sieve textures. Hornblendes are commonly unzoned,but some show simple, multiple or patchy zoning. Apatite microphenocrystsdisplay normal and reverse zoning of sulphur. The textural similarityof patchy hornblende and plagioclase, together with Ba–Srsystematics in patchy plagioclase, indicate that the cores ofthese crystals were derived from cumulate material. Plagioclase–liquidequilibria suggest that the patchy texture develops by resorptionduring H₂O-undersaturated decompression. When H₂O-saturatedcrystallization recommences at lower pressure, reduced pH₂Oresults in lower X_An in plagioclase, causing more Al-rich hornblendeto crystallize. Plagioclase cores with diffuse oscillatory zoning,and unzoned hornblende crystals, probably represent a populationof crystals resident in the magma chamber for long periods oftime. In contrast, oscillatory zoning in the rims of plagioclasephenocrysts may reflect eruption dynamics during decompressioncrystallization. Increasing Fe/Al in oscillatory zoned rimssuggests oxidation as a result of degassing of H₂O during decompression.A general lack of textural overlap between phenocryst typessuggests that different phenocryst populations were spatiallyor temporally isolated during crystallization. We present evidencethat the host andesite has mixed with both more felsic and moremafic magmas. Olivine and orthopyroxene xenocrysts with reactionor overgrowth rims and strong normal zoning indicate mixingwith basalt. Sieve-textured plagioclase resulted from mixingof a more felsic magma with the host andesite. The mineralogyand mineral compositions of a mafic andesite enclave are identicalto those of the host magma, which implies efficient thermalquenching, and thus small volumes of intruding magma. Mixingof this magma with the host andesite results in phenocryst zoningbecause of differences in dissolved volatile contents. We suggestthat small magma pulses differentiated at depth and ascendedintermittently into the growing magma chamber, producing incrementalvariations in whole-rock compositions. KEY WORDS: patchy zoning; magma mixing; Shiveluch     

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.     

Oxygen isotope study of the Long Valley magma system, California: isotope thermometry and convection in large silicic magma bodies

Products of voluminous pyroclastic eruptions with eruptive draw-down of several kilometers provide a snap-shot view of batholith-scale magma chambers, and quench pre-eruptive isotopic fractionations (i.e., temperatures) between minerals. We report analyses of oxygen isotope ratio in individual quartz phenocrysts and concentrates of magnetite, pyroxene, and zircon from individual pumice clasts of ignimbrite and fall units of caldera-forming 0.76 Ma Bishop Tuff (BT), pre-caldera Glass Mountain (2.1-0.78 Ma), and post-caldera rhyolites (0.65-0.04 Ma) to characterize the long-lived, batholith-scale magma chamber beneath Long Valley Caldera in California. Values of '¹⁸O show a subtle 1‰ decrease from the oldest Glass Mountain lavas to the youngest post-caldera rhyolites. Older Glass Mountain lavas exhibit larger (~1‰) variability of '¹⁸O(quartz). The youngest domes of Glass Mountain are similar to BT in '¹⁸O(quartz) values and reflect convective homogenization during formation of BT magma chamber surrounded by extremely heterogeneous country rocks (ranging from 2 to +29‰). Oxygen isotope thermometry of BT confirms a temperature gradient between "Late" (815 °C) and "Early" (715 °C) BT. The '¹⁸O(quartz) values of "Early" and "Late" BT are +8.33 and 8.21‰, consistent with a constant '¹⁸O(melt)=7.8ǂ.1‰ and 100 °C temperature difference. Zircon-melt saturation equilibria gives a similar temperature range. Values of '¹⁸O(quartz) for different stratigraphic units of BT, and in pumice clasts ranging in pre-eruptive depths from 6 to 11 km (based on melt inclusions), and document vertical and lateral homogeneity of '¹⁸O(melt). Worldwide, five other large-volume rhyolites, Lava Creek, Lower Bandelier, Fish Canyon, Cerro Galan, and Toba, exhibit equal '¹⁸O(melt) values of earlier and later erupted portions in each of the these climactic caldera-forming eruptions. We interpret the large-scale '¹⁸O homogeneity of BT and other large magma chambers as evidence of their longevity (>10⁵ years) and convection. However, remaining isotopic zoning in some quartz phenocrysts, trace element gradients in feldspars, and quartz and zircon crystal size distributions are more consistent with far shorter timescales (10²-10⁴ years). We propose a sidewall-crystallization model that promotes convective homogenization, roofward accumulation of more evolved and stagnant, volatile-rich liquid, and develops compositional and temperature gradients in pre-climactic magma chamber. Crystal + melt + gas bubbles mush near chamber walls of variable '¹⁸O gets periodically remobilized in response to chamber refill by new hotter magmas. One such episode of chamber refill by high-Ti, Sr, Ba, Zr, and volatile-richer magma happened 10³-10⁴ years prior to the 0.76-Ma caldera collapse that caused magma mixing at the base, mush thawing near the roof and walls, and downward settling of phenocrysts into this hybrid melt.     

Experimental constraints on permeable gas transport in crystalline silicic magmas

The gas and fluid transport in magmas via permeable flow through interconnected bubble networks controls the rate of outgassing from magmas ascending in volcanic conduits and the fluid transport in the mushy boundary layer of magma reservoirs. Hence, clarifying its mechanism and rate is crucial to understanding the explosivity of volcanic eruptions and the evolution and dynamics of a magma reservoir. Recent experimental studies have determined the gas permeabilities in crystal-free rhyolite and basalt. However, no experimental study has investigated the effect of the crystal contents on the permeable gas transport in magmas. In this study, we performed decompression experiments for hydrous rhyolitic melts having crystallinities of 30 and 50 vol% to examine the effect of crystals on the bubble microstructure and gas permeability during magma vesiculation. Size-controlled (100-meshed) corundum crystals were used as an analog of the phenocrysts in silicic magmas. Microstructural analyses using X-ray CT showed that bubbles coalesce and their connectivity increases with a decrease in the final pressure after the decompression, that is, an increase in the vesicularity. As long as the vesicularities of melt part in the crystal-free basis (melt vesicularity) were similar, no clear effect of the crystallinity on the degree of bubble coalescence and connectivity was observed at melt vesicularities <68 vol%. The corundum showed a large contact angle with aqueous fluid as well as plagioclase and alkaline feldspar; this failed to induce the efficient heterogeneous nucleation and coalescence of bubbles on its surface. The gas permeabilities of all the run products were lower than the detection limits of the present analysis (the order of 10⁻¹⁶ m²) at melt vesicularities <68 vol%. These results show that silicic magmas containing 30 and 50 vol% phenocrysts with a large contact angle have low gas permeabilities until the vesicularity becomes large (at least >68 vol%). This result indicates that the permeable fluid transport through a deep volcanic conduit, which has been proposed on the basis of the observations of volcanic gases and natural products, is so slow that other processes, like shear deformation or magma convection, may be needed to explain the observations.     

Magma storage and ascent during the 1995 eruption of Fogo,Cape Verde Archipelago

The 1995 eruption of Fogo (Cape Verde Islands) differed from previous eruptions by the occurrence of evolved lavas, the SW-orientation of vents, and pre-eruptive seismicity between Fogo and the adjacent (~20 km) island of Brava. We have conducted a thermobarometric and chemical study of this eruption in order to reconstruct its magma plumbing system and to test for possible connections to Brava. The bimodal eruption produced basanites (5.2–6.7 wt% MgO) and phonotephrites (2.4–2.8 wt% MgO) that are related by fractional crystallization. Clinopyroxene-melt-barometry of phenocrysts yields pressure ranges of 460–680 MPa for the basanites and 460–520 MPa for the phonotephrites. Microthermometry of CO₂-dominated fluid inclusions in olivine and clinopyroxene phenocrysts yields systematically lower pressure ranges of 200–310 MPa for basanites and 270–470 MPa for phonotephrites. The combined data indicate pre-eruptive storage of the 1995 magmas within the lithospheric mantle between 16 and 24 km depth. During eruption, the ascending magmas stalled temporarily at 8–11 km depth, within the lower crust, before they ascended to the surface in a few hours as indicated by zonations of olivine phenocrysts. Our data provide no evidence for magma storage at shallow levels (<200 MPa) or lateral magma movements beneath the Fogo-Brava platform. Sr–Nd–Pb isotope ratios of samples from Brava differ significantly from those of the 1995 and older Fogo lavas, which rules out contamination of the 1995 magmas by Brava material and indicates different mantle sources and magma plumbing systems for both islands.     

Petrology and Petrogenesis of the Quaternary Pumice Ash in the Taupo Area, New Zealand

Systematic mineralogical and chemical study of the younger Taupovolcanic ash showers, collected from the Terraces Quarry, Taupo,has enabled the showers to be grouped into ten eruptive sequences(numbers refer to the stratigraphical position of the deposits,after Baumgart (1954)): 3; 4; 5–8; 9–13; 14–15;16; 17–18; 19–22; 23–25; 26. The showers arecomposed of rhyolitic pumice and finer glass fragments; accessoryfragments dominantly of rhyolite, with minor dacite, ignimbrite,and andesite; crystals of plagioclase, hypersthene, and magnetite.The proportions of these components provide a useful guide todistinguishing the showers for correlation purposes. Refractiveindices of the glasses show wide variations, even in one pumicefragment, and increase markedly with increased weathering. Correlationof natural glasses by this method must, therefore, be done withextreme care. Vesiculation of pumice from most of the sequences has been studiedquantitatively by measurements of density and porosity (by modalanalysis). These give a measure of the extent of vesiculation,and in the sequences studied it is shown that there is a progressivedecrease in vesiculation with time in each sequence, here attributedto progressive volatile loss. In the large Taupo (3) and Waimihia(15) deposits, there was a rapid initial rise of intensity ofvesiculation. The last phase of several sequences appears tohave been relatively hotter than the earlier magma, as shownby more calcic plagioclase and more magnesian-rich pyroxenes. As a result of differences in nucleation of the various pumices,three general textures can be recognized: (i) finely cellularwith thick intercellular walls; (ii) finely cellular with thinintercellular walls; (iii) coarsely cellular pumice, usuallywith widely varying vesicle size. Glasses from seven members have been chemically analysed, andshow a close similarity. However, they fall into two distinctgroups, comprising the older showers (19, 24, 25) and the youngershowers (3, 5, 8, 15). The former have higher normative quartz,orthoclase, and anorthite. The compositions are rhyolitic. Plagioclase (An₃₆-An₄₈), hypersthene (Mg₄₈-Mg₆₄), and magnetiteform the characteristic mineralogical assemblage. The plagioclasehas well-developed normal oscillatory zoning, which is believedto be due to movement of crystals within the magma into zonesof varying vapour pressure. Many crystals have sharp discontinuitiesof some zones, which are attributed to sharp drops in vapourpressure, possibly due to eruption of preceding members of thesequence. Three chemical analyses of plagioclase are presented. From vesiculation and chemical results, it is believed thatsequences 3, 5–8, 14–15, and 17–18 crystallizedunder a very high vapour pressure, and were probably eruptedas a result of this. Sequences 19–22 and 23–25 crystallizedunder much lower vapour pressure. From the Ab-Or-Q-H₂O diagram,a vapour pressure of 2,000–3,000 kg/cm² is suggested forthe former sequences, equivalent to about 4 miles lithostaticload. Comparison is made between the mineralogy of the Tauporhyolitic glasses and recently published work on some Britishand Icelandic Tertiary rhyolitic glasses.     

Sixty thousand years of magmatic volatile history before the caldera-forming eruption of Mount Mazama, Crater Lake, Oregon

The well-documented eruptive history of Mount Mazama, Oregon, provides an excellent opportunity to use pre-eruptive volatile concentrations to study the growth of an explosive silicic magmatic system. Melt inclusions (MI) hosted in pyroxene and plagioclase crystals from eight dacitic–rhyodacitic eruptive deposits (71–7.7?ka) were analyzed to determine variations in volatile-element concentrations and changes in magma storage conditions leading up to and including the climactic eruption of Crater Lake caldera. Temperatures (Fe–Ti oxides) increased through the series of dacites, then decreased, and increased again through the rhyodacites (918–968 to ~950 to 845–895?°C). Oxygen fugacity began at nickel–nickel-oxide buffer (NNO) +0.8 (71?ka), dropped slightly to NNO +0.3, and then climbed to its highest value with the climactic eruption (7.7?ka) at NNO +1.1 log units. In parallel with oxidation state, maximum MI sulfur concentrations were high early in the eruptive sequence (~500?ppm), decreased (to ~200?ppm), and then increased again with the climactic eruption (~500?ppm). Maximum MI sulfur correlates with the Sr content (as a proxy for LREE, Ba, Rb, P₂O₅) of recharge magmas, represented by basaltic andesitic to andesitic enclaves and similar-aged lavas. These results suggest that oxidized Sr-rich recharge magmas dominated early and late in the development of the pre-climactic dacite–rhyodacite system. Dissolved H₂O concentrations in MI do not, however, correlate with these changes in dominant recharge magma, instead recording vapor solubility relations in the developing shallow magma storage and conduit region. Dissolved H₂O concentrations form two populations through time: the first at 3–4.6 wt% (with a few extreme values up to 6.1 wt%) and the second at ≤2.4 wt%. CO₂ concentrations measured in a subset of these inclusions reach up to 240?ppm in early-erupted deposits (71?ka) and are below detection in climactic deposits (7.7?ka). Combined H₂O and CO₂ concentrations and solubility models indicate a dominant storage region at 4–7?km (up to 12?km), with drier inclusions that diffusively re-equilibrated and/or were trapped at shallower depths. Boron and Cl (except in the climactic deposit) largely remained in the melt, suggesting vapor–melt partition coefficients and gas fractions were low. Modeled Li, F, and S vapor–melt partition coefficients are higher than those of B and Cl. The decrease in maximum MI CO₂ concentration following the earliest dacitic eruptions is interpreted to result from a broadening of the shallow storage region to greater than the diameter of subjacent feeders, so that greater proportions of reservoir magma were to the side of CO₂-bearing vapor bubbles ascending vertically from the locus of recharge magma injection, thereby escaping recarbonation by streaming vapor bubbles. The Mazama melt inclusions provide a picture of a growing magma storage region, where chemical variations in melt and magma occur due to changes in the nature and supply rate of magma recharge, the timing of degassing, and the possible degree of equilibration with gases from below.     

Timescales of mixing and mobilisation in the Bishop Tuff magma body: perspectives from diffusion chronometry

We present two-feldspar thermometry and diffusion chronometry from sanidine, orthopyroxene and quartz from multiple samples of the Bishop Tuff, California, to constrain the temperature stratification within the pre-eruptive magma body and the timescales of magma mixing prior to its evacuation. Two-feldspar thermometry yields estimates that agree well with previous Fe–Ti oxide thermometry and gives a ~80 °C temperature difference between the earlier- and later-erupted regions of the magma chamber. Using the thermometry results, we model diffusion of Ti in quartz, and Ba and Sr in sanidine as well as Fe–Mg interdiffusion in orthopyroxene to yield timescales for the formation of overgrowth rims on these crystal phases. Diffusion profiles of Ti in quartz and Fe–Mg in orthopyroxene both yield timescales of <150 years for the formation of overgrowth rims. In contrast, both Ba and Sr diffusion in sanidine yield nominal timescales 1–2 orders of magnitude longer than these two methods. 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. growth zoning exists). Utilising the divergent diffusion behaviour of Ba and Sr, we place constraints on the initial width of the interface and can refine our initial conditions considerably, bringing Ba and Sr data into alignment, and yielding timescales closer to 500 years, the majority of which are then within uncertainty of timescales modelled from Ti diffusion in quartz. Care must be thus taken when using Ba in sanidine geospeedometry in evolved magmatic systems where no other phases or elements are available for comparative diffusion profiling. Our diffusion modelling reveals piecemeal rejuvenation of the lower parts of the Bishop Tuff magma chamber at least 500 years prior to eruption. Timescales from our mineral profiling imply either that diffusion coefficients currently used are uncertain by 1–2 orders of magnitude, or that the minerals concerned did not experience a common history, despite being extracted from the same single pumice clasts. Introduction of the magma initiating crystallisation of the contrasting rims on sanidine, quartz, orthopyroxene and zircon was prolonged, and may be a marker of other processes that initiated the Bishop Tuff eruption rather than the trigger itself.     

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

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

The Campanian Ignimbrite (southern Italy) geochemical zoning: insight on the generation of a super-eruption from catastrophic differentiation and fast withdrawal

More than ca 100 km³ of nearly homogeneous crystal-poor phonolite and ca 100 km³ of slightly zoned trachyte were erupted 39 ka during the Campanian Ignimbrite super eruption, the most powerful in the Neapolitan area. Partition coefficient calculations, equilibrium mineral assemblages, glass compositions and texture were used to reconstruct compositional, thermal and pressure gradients in the pre-eruptive reservoir as well as timing and mechanisms of evolution towards magma chamber overpressure and eruption. Our petrologic data indicate that a wide sill-like trachytic magma chamber was active under the Campanian Plain at 2.5 kbar before CI eruption. Thermal exchange between high liquidus (1199°C) trachytic sill and cool country rocks caused intense undercooling, driving a catastrophic and fast (10² years) in situ fractional crystallization and crustal assimilation that produced a water oversaturated phonolitic cap and an overpressure in the chamber that triggered the super eruption. This process culminated in an abrupt reservoir opening and in a fast single-step high decompression. Sanidine phenocrysts crystal size distributions reveal high differentiation rate, thus suggesting that such a sill-like magmatic system is capable of evolving in a very short time and erupting suddenly with only short-term warning.     

Chemical zoning and crystallization mechanisms in the magma chamber of the Pomici di Base plinian eruption of Somma-Vesuvius (Italy)

Products of the Pomici di Base plinian eruption of Somma-Vesuvius consist of pumice and scoria fall deposits overlain by lithic-rich phreatomagmatic deposits. The plinian fall, which represents most of the magma volume involved in the eruption, ranges in composition from trachyte (SiO₂ = 62.5 wt%) to latite (SiO₂ ≈ 58 wt%) in the lower one-third of the deposit, whereas the upper two-thirds of the total thickness consists of latitic scoriae with fairly uniform composition (SiO₂ ≈ 55–56 wt%). All the products have very low content of phenocrysts (from 4 wt% in trachyte pumice to 1 wt% in the latite scoriae), most of which are not in equilibrium with the host rock. Minerals not in equilibrium, both in trachytic and latitic rocks, consist of discrete crystals of sanidine and plagioclase wetted by trachytic glass and felsic aggregates with interstital trachytic glass. Trends of major and trace elements are consistent with crystal-liquid fractionation processes and rule out syn-eruptive mixing processes between latitic and trachytic magmas. We suggest that discrete crystals and crystal aggregates not in equilibrium with the host rock represent fragments of the crystallising boundary layer at the upper walls of the magma chamber, which were wrenched and admixed into the magma during the ascent. This process diversifies the mineral assemblage and increases the crystal content of the rocks. We propose that diffusive crystallization processes operating at the wall of the chamber allowed the formation of a two-fold layered reservoir with a more mafic, homogeneous lower body and a more evolved, compositionally graded upper body. Around one-quarter of crystals adhering to the upper part of the magma chamber were admixed into the magma during the eruption. The absence of significant syn-eruptive mixing processes and the major role played by diffusive crystallization are consistent with a low aspect ratio magma chamber (width/height <1). Received: 23 March 1998 / Accepted: 11 December 1998