Magmatic processes in the Bishop Tuff rhyolitic magma based on trace elements in melt inclusions and pumice matrix glass

To investigate the origin of compositional zonation in the Bishop Tuff magma body, we have analyzed trace elements in the matrix glass of pumice clasts and in quartz-hosted melt inclusions. Our results show contrasting patterns for quartz in different parts of the Bishop Tuff. In all samples from the early part of the eruption, trace element compositions of matrix glasses are similar to but slightly more evolved than quartz-hosted melt inclusions. This indicates a cogenetic relationship between quartz crystals and their surrounding matrix glass, consistent with in situ crystallization. The range of incompatible element concentrations in melt inclusions and matrix glass from single pumice clasts requires 16–20 wt% in situ crystallization. This is greater than the actual crystal content of the pumices (<15 % crystals). In contrast to the pattern for the early pumices, pyroclastic flow samples from the middle part of the eruption show contrasting trends: In some clasts, the matrix is more evolved than the inclusions, whereas in other clasts, the matrix is less evolved. In the late Bishop Tuff, all crystal-rich samples have matrix glasses that are less evolved than the melt inclusions. Trace element abundances indicate that the cores of quartz in the late Bishop Tuff crystallized from more differentiated rhyolitic magma that was similar in many ways, yet distinct from the early-erupted Bishop Tuff. Our results are compatible with a model of secular incremental zoning (Hildreth and Wilson in Compositional zoning of the Bishop Tuff. J Petrol 48(5):951–999, 2007), in which melt batches from underlying crystal mush rise to various levels in a growing magma body according to their buoyancy. Early- and middle-erupted quartz crystallized from highly evolved rhyolitic melt, but then some parts of the middle-erupted magma were invaded by less differentiated rhyolite such that the matrix melt at the time of eruption was less evolved than the melt inclusions. A similar process occurred but to a greater extent in magma that erupted to form the late Bishop Tuff. In addition, there was a final, major magma mixing event in the late magma that formed Ti-rich rims on quartz and Ba-rich rims on sanidine, trapped less evolved rhyolitic melt inclusions, and resulted in dark and swirly crystal-poor pumice that is a rare type throughout much of the Bishop Tuff.     

Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano,Indonesia

The 2010 eruption of Merapi (VEI 4) was the volcano’s largest since 1872. In contrast to the prolonged and effusive dome-forming eruptions typical of Merapi’s recent activity, the 2010 eruption began explosively, before a new dome was rapidly emplaced. This new dome was subsequently destroyed by explosions, generating pyroclastic density currents (PDCs), predominantly consisting of dark coloured, dense blocks of basaltic andesite dome lava. A shift towards open-vent conditions in the later stages of the eruption culminated in multiple explosions and the generation of PDCs with conspicuous grey scoria and white pumice clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxene-hosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H₂O, CO₂) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have re-equilibrated during shallower storage and/or ascent, at depths of ~0.6–9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform “buffered” Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km³ in order to produce the mass of SO₂ emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing.     

Melt Inclusions in Olivine Phenocrysts: Using Diffusive Re-equilibration to Determine the Cooling History of a Crystal, with Implications for the Origin of Olivine-phyric Volcanic Rocks

A technique is described for determining the cooling historyof olivine phenocrysts. The technique is based on the analysisof the diffusive re-equilibration of melt inclusions trappedby olivine phenocrysts during crystallization. The mechanismof re-equilibration involves diffusion of Fe from and Mg intothe initial volume of the inclusion. The technique applies toa single crystal, and thus the cooling history of differentphenocrysts in a single erupted magma can be established. Weshow that melt inclusions in high-Fo olivine phenocrysts frommantle-derived magmas are typically partially re-equilibratedwith their hosts at temperatures below trapping. Our analysisdemonstrates that at a reasonable combination of factors suchas (1) cooling interval before eruption (<350°C), (2)eruption temperatures (>1000°C), and (3) inclusion size(<70 µm in radius), partial re-equilibration of upto 85% occurs within 3–5 months, corresponding to coolingrates faster than 1–2°/day. Short residence timesof high-Fo phenocrysts suggest that if eruption does not happenwithin a few months after a primitive magma begins cooling andcrystallization, olivines that crystallize from it are unlikelyto be erupted as phenocrysts. This can be explained by efficientseparation of olivine crystals from the melt, and their rapidincorporation into the cumulate layer of the chamber. Theseresults also suggest that in most cases erupted high-Fo olivinephenocrysts retain their original composition, and thus compositionsof melt inclusions in erupted high-Fo olivine phenocrysts donot suffer changes that cannot be reversed. Short residencetimes also imply that large unzoned cores of high-Fo phenocrystscannot reflect diffusive re-equilibration of originally zonedphenocrysts. The unzoned cores are a result of fast efficientaccumulation of olivines from the crystallizing magma, i.e.olivines are separated from the magma faster than melt changesits composition. Thus, the main source of high-Fo crystals inthe erupted magmas is the cumulate layers of the magmatic system.In other words, olivine-phyric rocks represent mixtures of anevolved transporting magma (which forms the groundmass of therock) with crystals that were formed during crystallizationof more primitive melt(s). Unlike high-Fo olivine phenocrysts,the evolved magma may reside in the magmatic system for a longtime. This reconciles long magma residence times estimated fromthe compositions of rocks with short residence times of high-Foolivine phenocrysts. KEY WORDS: melt inclusions; olivine; picrites; residence time; diffusion     

Melt and fluid inclusion evidence for a genetic relationship between magmatism in the Shuangliao volcanic field and inorganic CO2 gas reservoirs in the southern Songliao Basin,China

We have studied melt and fluid inclusions in minerals from alkali basalts, mantle xenoliths, and dawsonite-bearing sandstones from the Shuangliao volcanic field in southern Songliao Basin, Northeast China. The inclusions have been investigated using petrographic, geochemical, and laser Raman spectroscopic techniques. Volcanic rocks of the Shuangliao field are predominantly alkali olivine basalts that contain rare mantle xenoliths. Silicate melt and fluid inclusions are common in both olivine phenocrysts and the mantle xenoliths. The fluid inclusions are mainly composed of CO₂ with small amounts of CO, CH₄, N₂, and H₂O, which is consistent with an upper mantle origin. CO₂ gas reservoirs in the southern Songliao Basin are mostly derived from a mantle–magmatic source. Coeval fluid-inclusion homogenization temperatures, coupled with the thermal burial history, show that the CO₂ gas reservoirs in the southern Songliao Basin are Cenozoic (40–63 Ma) and coeval with the magmatism in the Shuangliao volcanic field. Despite the relatively small scale of this volcanic activity, it released large amounts of CO₂. Much of the magma was not erupted, and CO₂- and H₂O-rich magma was probably intruded into the basin along deep faults, acting as a major source of inorganic CO₂ gas in the southern Songliao Basin.     

The 1783 Lakagigar eruption in Iceland: geochemistry,CO2 and sulfur degassing

About 12.3 km³ of basaltic magma were erupted from the Lakagigar fissure in Iceland in 1783, which may have been derived from the high-level reservoir of Grimsvotn central volcano, by lateral flow within the rifted crust. We have studied the petrology of quenched, glassy tephra from sections through pyroclastic cones along the fissure. The chemical composition of matrix glass of the 1783 tephra is heterogeneous and ranges from olivine tholeiite to Fe–Ti rich basalt, but the most common magma erupted is quartz tholeiite (Mg#43.6 to 37.2). The tephra are characterized by low crystal content (5 to 9 vol%). Glass inclusions trapped in plagioclase and Fo₈₆ to Fo₇₅ olivine phenocrysts show a large range of compositions, from primitive olivine tholeiite (Mg#64.3), quartz tholeiite (Mg#43–37), to Fe–Ti basalts (Mg#33.5) which represent the most differentiated liquids and are trapped as rare melt inclusions in clinopyroxene. Both matrix glass and melt inclusion data indicate a chemically heterogeneous magma reservoir, with quartz tholeiite dominant. LREE-depleted olivine-tholeiite melt-inclusions in Mg-rich olivine and anorthitic-plagioclase phenocrysts may represent primitive magma batches ascending into the reservoir at the time of the eruption. Vesicularity of matrix glasses correlates with differentiation, ranging from 10 to 60 vol.% in evolved quartz-tholeiite glasses, whereas olivine-tholeiite glasses contain less than 10 vol.% vesicles. FTIR analyses of olivine-tholeiite melt-inclusions indicate concentrations of 0.47 wt% H₂O and 430 to 510 ppm for CO₂. Chlorine in glass inclusions and matrix glasses increases from 50 ppm in primitive tholeiite to 230 ppm in Fe–Ti basalts, without clear evidence of degassing. Melt inclusion analyses show that sulfur varies from 915 ppm to 1970 ppm, as total FeO^* increases from 9 to 13.5 wt%. Sulfur degassing correlates both with vesicularity and magma composition. Thus sulfur in matrix glasses decreases from 1490 ppm to 500 ppm, as Mg # decreases from 47 to 37 and vesicularity of the magma strongly increases. These results indicate loss of at least 75% of sulfur during the eruption. The correlation of low sulfur content in matrix glasses with high vesicularity is regarded as evidence of the control of a major exsolving volatile phase on the degassing efficiency of the magma. Our model is consistent a quasi-permanent CO₂ flux through the shallow-level magmatic reservoir of Grimsvotn. Following magma withdrawal from the reservoir and during eruption from the Lakagigar fissure, sulfur degassing was controlled by inherent CO₂-induced vesicularity of the magma.     

Evolution of Bishop Tuff Rhyolitic Magma Based on Melt and Magnetite Inclusions and Zoned Phenocrysts

The evolution of large bodies of silicic magma is an importantaspect of planetary differentiation. Melt and mineral inclusionsin phenocrysts and zoned phenocrysts can help reveal the processesof differentiation such as magma mixing and crystal settling,because they record a history of changing environmental conditions.Similar major element compositions and unusually low concentrationsof compatible elements (e.g. 0·45–4·6 ppmBa) in early-erupted melt inclusions, matrix glasses and bulkpumice from the Bishop Tuff, California, USA, suggest eutectoidfractional crystallization. On the other hand, late-eruptedsanidine phenocrysts have rims rich in Ba, and late-eruptedquartz phenocrysts have CO₂-rich melt inclusions closest tocrystal rims. Both features are the reverse of in situ crystallizationdifferentiation, and they might be explained by magma mixingor crystal sinking. Log(Ba/Rb) correlates linearly with log(Sr/Rb)in melt inclusions, and this is inconsistent with magma mixing.Melt inclusion gas-saturation pressure increases with CO₂ fromphenocryst core to rim and suggests crystal sinking. Some inclusionsof magnetite in late-erupted quartz are similar to early-eruptedmagnetite phenocrysts, and this too is consistent with crystalsinking. We argue that some large phenocrysts of late-eruptedquartz and sanidine continued to crystallize as they sank severalkilometers through progressively less differentiated melts.Probable diffusive modification of Sr in sanidine phenocrystsand the duration of crystal sinking are consistent with an evolutionaryinterval of some 100 ky or more. Crystal sinking enhanced thedegree of differentiation of the early-erupted magma and pointsto the importance of H₂O (to diminish viscosity and enhancethe rate of crystal sinking) in the evolution of silicic magmas. KEY WORDS: crystal settling; differentiation; melt inclusions; rhyolite; trace elements     

Melt inclusions and crystal-liquid separation in rhyolitic magma of the Bishop Tuff

Melt inclusions in quartz phenocrysts from a single clast of pumice near the base of the plinian pumice fall of the Bishop Tuff were studied to test ideas concerning separation of melt and crystals in silicic magmas. Ten analyzed inclusions from the pumice clast are of high silica rhyolite composition with very low contents of the highly compatible elements Ba, Sr, and Eu, consistent with extensive fractionation. The concentrations of U, La, Ce, Mg, and Ca of these ten melt inclusions vary considerably as determined by ion microprobe. Petrologic considerations indicate that uranium is an incompatible element with a maximum bulk partition coefficient D of about 0.2 and that the evolution of the uranium content of the melt was controlled by crystallization of the magma. A minimum of 33 wt% perfect fractional crystallization is required to explain the observed range in uranium. However, only 17 wt% crystals occurred in the pumice clast. The greater calculated fraction of crystals requires significant separation of crystals and melt before the eruption of the plinian pumice fall in spite of the fact that crystal mixing (settling, etc.) did not occur in the Bishop magma.     

Rapid open-system assembly of a large silicic magma body: time-resolved evidence from cored plagioclase crystals in the Oruanui eruption deposits, New Zealand

Rhyolitic pumices in the 26.5 ka Oruanui eruption (Taupo volcano, New Zealand) contain an average of 10 wt% crystals. About 2 wt% of the crystal population is feldspar crystals that display bluish–grey cloudy cores, the colour being imparted by exsolved needles of rutile. The volume of cloudy-cored feldspars thus amounts to ~1.0 km³ in a total magma volume of ~530 km³. The cored feldspars show great variability in detail, but in general have a rounded cloudy core bounded by a zone rich in glass and mineral inclusions, that was then overgrown by a euhedral clear rim. Sr-isotopic variations in eight representative crystals were measured on micromilled samples of selected growth zones in the cores and rims, and linked to feldspar compositions through microprobe traverses. The cloudy cores range from ⁸⁷Sr/⁸⁶Sr = 0.70547 to 0.70657, with compositions of An₄₃ to An₇₈. The overgrowth rims display wider variations: inner parts show extreme ranges in composition (maxima ⁸⁷Sr/⁸⁶Sr = 0.70764 and An₇₈), while outer parts in seven of eight crystals are zoned, with outward-decreasing Sr-isotopic and An values to figures that are in accord with the bulk pumice and other, clear-feldspar values, respectively. The three parts of the crystals represent distinct regimes. The cloudy cores are inherited from an intermediate plutonic protolith that has been subjected to melting. The inner overgrowth rims were crystallised from a high temperature, relatively radiogenic melt derived from Mesozoic-Palaeozoic metasedimentary rocks (“greywacke”). The outer euhedral rims reflect mixing into and continued growth within the growing Oruanui magma body. The cloudy-cored feldspars also contain rare zircon inclusions. Twenty one zircons were recovered by HF digestion of a bulk sample of cloudy feldspars and analysed by SHRIMP for U–Th isotopes with which to calculate model ages. Eighteen of 21 crystals returned finite ages, the model-age spectrum of which is similar to the age spectra from free zircons in Oruanui pumices. Assembly of the Oruanui magma body was not only rapid (over ~40 kyr, as shown by other data) but involved a wide open system, with significant contributions from partly-melted intermediate-composition igneous intrusions (cloudy cores) and greywacke melts (inner overgrowths) being introduced into the magma body up to the point of eruption. Such open system behaviour contrasts with that proposed in models for comparably voluminous silicic magmas derived dominantly by fractionation (such as the Bishop Tuff) where the magma and its crystal cargo were better insulated thermally and chemically from country-rock interaction.     

Remobilization of Highly Crystalline Felsic Magma by Injection of Mafic Magma: Constraints from the Middle Sixth Century Eruption at Haruna Volcano, Honshu, Japan

The latest eruption of Haruna volcano at Futatsudake took placein the middle of the sixth century, starting with a Plinianfall, followed by pyroclastic flows, and ending with lava domeformation. Gray pumices found in the first Plinian phase (lowerfall) and the dome lavas are the products of mixing betweenfelsic (andesitic) magma having 50 vol. % phenocrysts and maficmagma. The mafic magma was aphyric in the initial phase, whereasit was relatively phyric during the final phase. The aphyricmagma is chemically equivalent to the melt part of the phyricmafic magma and probably resulted from the separation of phenocrystsat their storage depth of 15 km. The major part of the felsicmagma erupted as white pumice, without mixing and heating priorto the eruption, after the mixed magma (gray pumice) and heatedfelsic magma (white pumice) of the lower fall deposit. Althoughthe mafic magma was injected into the felsic magma reservoir(at 7 km depth), part of the product (lower fall ejecta) precedederuption of the felsic reservoir magma, as a consequence ofupward dragging by the convecting reservoir of felsic magma.The mafic magma injection made the nearly rigid felsic magmaerupt, letting low-viscosity mixed and heated magmas open theconduit and vent. Indeed the lower fall white pumices preservea record of syneruptive slow ascent of magma to 2 km depth,probably associated with conduit formation. KEY WORDS: high-crystallinity felsic magma; magma plumbing system; multistage magma mixing; upward dragging of injected magma; vent opening by low-viscosity magma     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

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.     

A model for mixing basaltic and dacitic magmas as deduced from experimental data

It has been demonstrated experimentally that basaltic and dacitic magmas can be easily mixed to form both banded dacite and homogeneous andesite in less than a few hours. The presence of phenocrysts larger than 0.5 mm increased considerably the efficiency of mixing. Flow patterns in the experimental system were visualized using Pt spheres, which indicated that convection occurs in basalt melt, but not in dacite melt. The Reynolds numbers of the basaltic and dacitic melts in the experimental system were calculated to be about 10^–3 and 10^–6, respectively. Mixing proceeds initially by mechanical mixing of the two magmas in a large scale, but later by coupling interfacial convection and mutual diffusion. Thus, depending on the depth where vesiculation and following disruption of the magma occurs, banded pumice, homogeneous pumice and homogeneous andesite lava are erupted. The observed textures of mixed rocks of Plinian type eruption and the limiting occurrence of banded pumice are satisfactorily accounted for on this model.     

Experimental Petrology of the 1991-1995 Unzen Dacite, Japan. Part I: Phase Relations, Phase Composition and Pre-eruptive Conditions

Crystallization experiments were conducted on dry glasses fromthe Unzen 1992 dacite at 100–300 MPa, 775–875°C,various water activities, and fO₂ buffered by the Ni–NiObuffer. The compositions of the experimental products and naturalphases are used to constrain the temperature and water contentsof the low-temperature and high-temperature magmas prior tothe magma mixing event leading to the 1991–1995 eruption.A temperature of 1050 ± 75°C is determined for thehigh-temperature magma based on two-pyroxene thermometry. Theinvestigation of glass inclusions suggests that the water contentof the rhyolitic low-temperature magma could be as high as 8wt % H₂O. The phase relations at 300 MPa and in the temperaturerange 870–900°C, which are conditions assumed to berepresentative of the main magma chamber after mixing, showthat the main phenocrysts (orthopyroxene, plagioclase, hornblende)coexist only at reduced water activity; the water content ofthe post-mixing dacitic melt is estimated to be 6 ± 1wt % H₂O. Quartz and biotite, also present as phenocrysts inthe dacite, are observed only at low temperature (below 800–775°C).It is concluded that the erupted dacitic magma resulted fromthe mixing of c. 35 wt % of an almost aphyric pyroxene-bearingandesitic magma (1050 ± 75°C; 4 ± 1 wt % H₂Oin the melt) with 65 wt % of a phenocryst-rich low-temperaturemagma (760–780°C) in which the melt phase was rhyolitic,containing up to 8 ± 1 wt % H₂O. The proportions of rhyoliticmelt and phenocrysts in the low-temperature magma are estimatedto be 65% and 35%, respectively. It is emphasized that the strongvariations of phenocryst compositions, especially plagioclase,can be explained only if there were variations of temperatureand/or water activity (in time and/or space) in the low-temperaturemagma. KEY WORDS: Unzen volcano; magma mixing; experimental study     

The Bishop Tuff giant magma body: an alternative to the Standard Model

The Bishop Tuff, one of the most extensively studied high-silica rhyolite bodies in the world, is usually considered as the archetypical example of a deposit formed from a magma body characterized by thermal and compositional vertical stratification—what we call the Standard Model for the Bishop magma body. We present here new geothermometry and geobarometry results derived using a large database of previously published quartz-hosted glass inclusion compositions. Assuming equilibrium between melt and an assemblage composed of quartz, ±plagioclase, ±sanidine, +zircon, ±fluid, we use Zr contents in glass inclusions to derive quartz crystallization temperatures, and we use (1) silica contents in glass, (2) projection of glass compositions onto the haplogranitic (quartz-albite-orthoclase) ternary, and (3) phase equilibria calculations using rhyolite-MELTS, to constrain crystallization pressures. We find crystallization temperatures of ~740–750 °C for all inclusions from both early- and late-erupted pumice. Crystallization pressures for both early- and late-erupted inclusions are also very similar to each other, with averages of ~175–200 MPa. We find no evidence of late-erupted inclusions having been entrapped at higher temperatures or pressures than early-erupted inclusions, as would be expected by the Standard Model. We argue that the thermal gradient inferred from Fe–Ti oxides—the backbone of the Standard Model—does not reflect equilibrium pre-eruptive conditions; we also note that H₂O–CO₂ systematics of glass inclusions yields overlapping pressure ranges for early- and late-erupted inclusions, similar to the results presented here; and we show that glass inclusion and phenocryst compositions show bimodal distributions, suggestive of compositional separation between early- and late-erupted populations. These findings are inconsistent with the Standard Model. The similarity in crystallization conditions and the compositional separation between early- and late-erupted magmas suggest that two laterally juxtaposed independent magma reservoirs existed in the same region at the same time and co-erupted to form the Long Valley Caldera and the Bishop Tuff. This hypothesis would explain the lack of mixing between early- and late-erupted crystal populations in pumice clasts; it could also explain the inferred eruption pattern—which resulted in early-erupted magmas being deposited only to the south of the caldera—if the early-erupted magma body resided to the south and the late-erupted magma body was located to the north. Our alternative model is consistent with the patchy distribution of thermal anomalies and the inference of co-eruption of distinct magma types in active volcanic areas such as the central Taupo Volcanic Zone.     

Volatiles H2O,CO2, and Cl in a subduction related basalt

The products of the 1974 eruption of Fuego, a subduction zone volcano in Guatemala, have been investigated through study of silicate melt inclusions in olivine. The melt inclusions sampled liquids in regions where olivine, plagioclase, magnetite, and augite were precipitating. Comparisons of the erupted ash, groundmass, and melt inclusion compositions suggest that the inclusions represent samples of liquids present in a thermal boundary layer of the magma body. The concentrations of H₂O and CO₂ in glass inclusions were determined by a vacuum fusion manometric technique using individual olivine crystals (Fo77 to Fo71) with glass inclusion compositions that ranged from high-alumina basalt to basaltic andesite. Water, Cl, and K₂O concentrations increased by a factor of two as the olivine crystals became more iron-rich (Fo77 to Fo71) and as the glass inclusions increased in SiO₂ from 51 to 54 wt.% SiO₂. The concentration of H₂O in the melt increased from 1.6 wt.% in the least differentiated liquid to about 3.5% in a more differentiated liquid. Carbon dioxide is about an order of magnitude less abundant than H₂O in these inclusions. The gas saturation pressures for pure H₂O in equilibrium with the melt inclusions, which were calculated from the glass inclusion compositions using the solubility model of Burnham (1979), are given approximately by P(H₂O)(Pa)=(SiO₂−48.5 wt.%) × 1.45 × 10⁷. The concentrations of water in the melt and the gas saturation pressures increased from about 1.5% to 3.5% and from 300 to 850 bars, respectively, during pre-eruption crystallization.     

Millennial timescale resolution of rhyolite magma recharge at Tarawera volcano: insights from quartz chemistry and melt inclusions

Most rhyolite eruption episodes of Tarawera volcano have emitted several physiochemically distinct magma batches (∼1–10 km³). These episodes were separated on a millennial timescale. The magma batches were relatively homogeneous in temperature and composition at pumice scale (>4 cm), but experienced isolated crystallisation histories. At the sub-cm scale, matrix glasses have trace element compositions (Sr, Ba, Rb) that vary by factors up to 2.5, indicating incomplete mixing of separate melts. Some quartz-hosted melt inclusions are depleted in compatible trace elements (Sr, Ti, Ba) compared to enclosing matrix glasses. This could reflect re-melting of felsic crystals deeper in the crystal pile. Individual quartz crystals display a variety of cathodoluminescence brightness and Ti zoning patterns including rapid changes in melt chemistry and/or temperature (∼50–100°C), and point to multi-cycle crystallisation histories. The Tarawera magma system consisted of a crystal-rich mass containing waxing and waning melt pockets that were periodically recharged by silicic melts driven by basaltic intrusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Evolution of the trachydacite and pantellerite magmas of the bimodal volcanic association of Dzarta-Khuduk, Central Mongolia: Investigation of inclusions in minerals

Using various methods of melt inclusion investigation, including electron and ion microprobe techniques, we estimated the composition, evolution, and formation conditions of melts producing the trachydacites and pantellerites of the Late Paleozoic bimodal volcanic association of Dzarta-Khuduk, Central Mongolia. Primary crystalline and melt inclusions were detected in anorthoclase from trachydacites and quartz from pantellerites and pantelleritic tuffs. Among the crystalline inclusions, we identified hedenbergite, fluorapatite, and pyrrhotite in the trachydacites and F-arfvedsonite, fluorite, ilmenite, and the rare REE diorthosilicate chevkinite in the pantellerites. Melt inclusions in anorthoclase from the trachydacites are composed of glass, a gas phase, and daughter minerals (F-arfvedsonite, fluorite, villiaumite, and anorthoclase rim on the inclusion wall). Melt inclusions in quartz from the pantellerites are composed of glass, a gas phase, and a fine-grained salt aggregate consisting of Li, Na, and Ca fluorides (griceite, villiaumite, and fluorite). Melt inclusions in quartz crystalloclasts from the pantelleritic tuffs are composed of homogeneous silicate glasses. The phenocrysts of the trachydacites and pantellerites crystallized at temperatures of 1060–1000°C. During thermometric experiments with quartz-hosted melt inclusions from the pantellerites, the formation of immiscible silicate and salt (fluoride) melts was observed at a temperature of 800°C. Homogeneous melt inclusions in anorthoclase from the trachydacites have both trachydacite and rhyolite compositions (wt %): 68–70 SiO₂, 12–13 Al₂O₃, 0.34–0.74 TiO₂, 5–7 FeO, 0.4–0.9 CaO, and 9–12 Na₂O + K₂O. The agpaitic index ranges from 0.92 to 1.24. The glasses of homogenized melt inclusions in quartz from the pantellerites and pantelleritic tuffs have rhyolitic compositions. Compared with the homogeneous glasses trapped in anorthoclase of the trachydacites, quartz-hosted inclusions from the pantellerites show higher SiO₂ (72–78 wt %) and lower Al₂O₃ contents (7.8–10.0 wt %). They also contain 0.14–0.26 wt % TiO₂, 2.5–4.9 wt % FeO, 9–11 wt % Na₂O + K₂O, and 0.9–0.15 wt % CaO and show an agpaitic index of 1.2–2.05. Homogeneous melt inclusions in quartz from the pantelleritic tuffs contain 69–72 wt % SiO₂. The contents of other major components, including TiO₂, Al₂O₃, FeO, and CaO, are close to those in the homogeneous glasses of quartzhosted melt inclusions in the pantellerites. The contents of Na₂O + K₂O are 4–10 wt %, and the agpaitic index is 1.0–1.6. The glasses of melt inclusions from each rock group show distinctive volatile compositions. The H₂O content is up to 0.08 wt % in anorthoclase of the trachydacites, 0.4–1.4 wt % in quartz of the pantellerites, and up to 5 wt % in quartz of the pantelleritic tuffs. The content of F in the glasses of melt inclusions in the phenocrysts of the trachydacites is no higher than 0.67 wt %, and up to 1.4–2.8 wt % in quartz from the pantellerites. The Cl content is up to 0.2 wt % in the glasses of melt inclusions in the minerals of the trachydacites and up to 0.5 wt % in the glasses of quartz-hosted melt inclusions from the pantellerites. The investigation of trace elements in the homogenized glasses of melt inclusions in minerals showed that the trachydacites and pantellerites were formed from strongly evolved rare-metal alkaline silicate melts with high contents of Li, Zr, Rb, Y, Hf, Th, U, and REE. The analysis of the composition of homogeneous melt inclusions in the minerals of the above rocks allowed us to distinguish magmatic processes resulting in the enrichment of these rocks in trace and rare earth elements. The most important processes are the crystallization differentiation and immiscible separation of silicate and fluoride salt melts. It was also shown that all the melts studied evolved in spatially separated magma chambers. This caused the differences in the character of melt evolution between the trachydacites and pantellerites. During the final stages of differentiation, when the magmatic system was saturated with respect to ore elements, Na-Ca fluoride melts were separated and extracted considerable amounts of Li.     

长白山天池火山千年大喷发岩浆含水量研究——熔融包裹体含水量的红外光谱测试

长白山天池火山在公元一千年左右曾发生过大规模喷发，其产物为大面积分布的灰白色碱流质浮岩和碎屑流，在其斑晶矿物橄榄石、钙铁辉石和碱性长石中均可见到熔融包裹体。在最主要的斑晶矿物——碱性长石中含有数量众多且个体也较大的熔融包裹体，多数含有一个以上的气泡，其中部分含有子晶，根据形貌特征的不同可分为截然不同的两组包裹体。这些熔融包裹体带有大量喷发前地下岩浆的信息，成为研究地下深部的岩浆在复杂的溢流-爆炸喷发中所发生变化的最好媒介。也是本文的研究对象，通过其中挥发份尤其是水的含量，可以推知天池火山发生大喷发的原因。经Nicolet Magna-IR 550红外光谱仪测定，这些熔融包裹体的含水量较高，达1．6％-3．6％，为当时天池火山发生了巨大规模爆炸喷发的原因提供了强有力的证据。但目前红外光谱仪的应用范围还比较有限，有待今后拓宽其应用领域。     

Reconstruction of an Evolving Magma Chamber beneath Usu Volcano since the 1663 Eruption

Usu volcano, located in northern Japan, has erupted seven timessince AD 1663. Before these seven eruptions, the volcano hada long repose period ( 5000 yr). The 1663 eruption was thefirstand by far the largest among the seven, producing nearlyaphyric rhyolitic pumice. Small mafic inclusions (‘micro-clots’J,consisting of glass, quenched crystals and abundant vesiclesoccur in the pumice. On the basis of petrological studies ofthe microclots, it is concluded that these are quenched meltsof a mafic magma injected into the rhyolite. The products ofthe 1769 eruption (and those of the following five eruptions)were dacites with abundant (10–15 vol %) microphenocrysts.According to crystal size distribution (CSD) analysis, the newmicrophenocrysts appear to have crystallized at a considerablyhigher cooling rate ( 300 times) than the phenocrysts in the1663 eruptive products. The contrasting petrologic featuresof the aphyric rhyolite and the following microphenocryst-richdacites can be explained by mixing and rapid cooling of a maficmagma injected during the 1663 eruption. We estimate the sizeof the magma chamber beneath Usu volcano just after the 1663eruption, using numerical calculations for a cooling magma chamber.If the magma chamber was sill-like, its thickness is estimatedto have been several hundreds of meters. KEY WORDS: Usu volcano; Japan; magma chamber evolution ^*Corresponding author. Present address: Geomechanics, Earthquake Research Institute, The University of Tokyo, I-I-I Yayoi, Bunkyo-ku, Tokyo 113, Japan