Petrology of the zoned calcalkaline magma chamber of Mount Mazama,Crater Lake,Oregon

Evolution of the magma chamber at Mount Mazama involved repeated recharge by two types of andesite (high-Sr and low-Sr), crystal fractionation, crystal accumulation, assimilation, and magma mixing (Bacon and Druitt 1988). This paper addresses the modal compositions, textures, mineral chemistry and magmatic temperatures of (i) products of the 6845±50 BP climactic eruption, (ii) blocks of partially fused granitoid wallrock found in the ejecta, and (iii) preclimactic rhyodacitic lavas leaked from the chamber in late Pleistocene and early Holocene time. Immediately prior to the climactic eruption the chamber contained ≳ 40 km³ of rhyodacite (10 vol% plag + opx + aug + hb + mt + ilm, ∼880° C) overlying high-Sr andesite and cumulus-crystal mush (28–51 vol% plag + hb ± opx ± aug + mt ± ilm, 880° to ≥950° C), which in turn overlay low-Sr crystal mush (50–66 vol% plag + opx + aug ± hb ± ol + mt + ilm, 890° to ≥950† C). Despite the well known compositional gap in the ejecta, no thermal discontinuity existed in the chamber. Pre-eruptive water contents of pore liquids in most high-Sr and low-Sr mushes were 4–6 wt%, but on average the high-Sr mushes were slightly richer in water. Although parental magmas of the crystal mushes were andesitic, xenocrysts of bytownite and Ni-rich magnesian olivine in some scoriae record the one-time injection of basalt into the chamber. Textures in ol-bearing scoriae preserve evidence for the reactions ol + liq = opx and ol + aug + liq(+ plag?) = hb, which occurred in andesitic liquids at Mount Mazama. Strontium abundances in plagioclase phenocrysts constrain the petrogenesis of preclimactic and climactic rhyodacites. Phenocryst cores derived from high-Sr and low-Sr magmas have different Sr contents which can be resolved by microprobe. Partition coefficients for plagioclase in andesitic to rhyolitic glasses range from 2 to 7, and increase as glass %SiO₂ increases. Evolved Pleistocene rhyodacites (∼30–25,000 BP) and rhyodacites of the Holocene Llao Rock center (7015±45 BP) contain Sr-poor plagioclase and are derivatives from low-Sr magma. Rhyodacites of the Pleistocene Sharp Peak domes, Holocene Cleetwood flow (∼6850 BP), and climactic ejecta contain discrete Sr-rich and Sr-poor plagioclase phenocryst populations and are hybrids produced by mixing low-Sr rhyodacite (containing Sr-poor plag + opx + aug) with a more mafic high-Sr magma (with Sr-rich plag [ + hb?]). The data reinforce the conclusions of crystal-liquid mixing calculations (Bacon and Druitt 1988), and suggest some important refinements to the magma chamber model.     

The Petrology and Geochemistry of the Aniakchak Caldera-forming Ignimbrite, Aleutian Arc, Alaska

Aniakchak caldera, Alaska, produced a compositionally heterogeneousignimbrite 3400 years ago, which changes from rhyodacitic atthe base to andesitic at the top of the eruptive sequence. Interpretationsof compositionally heterogeneous ignimbrites typically includeeither in situ fractional crystallization of mafic magma andgeneration of a stratified magma body or replenishment of asilicic magma chamber by mafic inputs. Another possibility,silicic replenishment of a more mafic chamber, exists. Geochemicalcharacteristics of the caldera-forming rhyodacite and severallate pre-caldera rhyodacites indicate independent origins foreach, within a maximum of 5000 years prior to caldera formation.Isotopic considerations preclude derivation of the caldera-formingrhyodacite from the caldera-forming andesite. However, the caldera-formingrhyodacite can be explained as the residual liquid of a mostlycrystallized basalt, with addition of crustal material. TheAniakchak andesite probably formed in a shallow chamber by successivemixing events involving small volumes of basalt and rhyodacite,together with contamination. The pre-caldera rhyodacites representerupted portions of intruding silicic magma, whereas anotherportion homogenized with the resident mafic magma. The caldera-formingevent reflects a large influx of rhyodacite, which erupted beforesignificant mixing occurred and also triggered draining of muchof the andesitic magma from the chamber. KEY WORDS: Aniakchak; caldera-forming eruption; geochemistry; ignimbrite; silicic replenishment     

Distinguishing lower and upper crustal processes in magmas erupted during the buildup to the 7.7 ka climactic eruption of Mount Mazama, Crater Lake, Oregon, using 238U–230Th disequilibria

Uranium-series isotope ratios determined for 35 volcanic rocks and 4 glass separates erupted from ~36 to 4.8 ka at Mt. Mazama, Crater Lake, Oregon, identify both ²³⁰Th-excess and ²³⁸U-excess components. U–Th isotope compositions cover a wide range, exceeding those previously measured for the Cascade arc. Age-corrected (²³⁰Th/²³²Th) and (²³⁸U/²³²Th) activity ratios range from 1.113 to 1.464 and from 0.878 to 1.572 (44.4 % ²³⁰Th-excess to 8.8 % ²³⁸U-excess), respectively. The most distinctive aspect of the data set is the contrast in U–Th isotope ratios between low and high Sr (LSr, HSr) components that have been previously identified in products of the 7.7 ka caldera-forming climactic eruption and preclimactic rhyodacite lavas. The LSr component exclusively contains ²³⁸U-excess, but the HSr component, as well as more primitive lavas, are marked by ²³⁰Th-excess. ²³⁰Th-excesses such as those recorded at Mt. Mazama are commonly observed in the Cascades. Melting models suggest that high ²³⁰Th-excesses observed in the more primitive lavas evolved through mixing of a mantle melt with a partial melt of a mafic lower crustal composition that contained garnet in the residuum that was produced through dehydration melting of amphibolite that was initially garnet free. Dehydration melting in the lower crust offers a solution to the “hot-slab paradox” of the Cascades, where low volatile contents are predicted due to high slab temperatures, yet higher water contents than expected have been documented in erupted lavas. The ²³⁸U-excess observed at Mt. Mazama is rare in Cascade lavas, but occurs in more than half of the samples analyzed in this study. Traditionally, ²³⁸U-excess in arc magmas is interpreted to reflect slab fluid fluxing. Indeed, ²³⁸U-excess in arcs is common and likely masks ²³⁰Th-excess resulting from lower crustal interaction. Isotopic and trace element data, however, suggest a relatively minor role for slab fluid fluxing in the Cascades. We propose that ²³⁸U-excess reflects melting and assimilation of young, hydrothermally altered upper crust. The processes related to generating ²³⁸U-excess are likely important features at Mt. Mazama that accompanied development of a large-scale silicic magma chamber that led to the caldera-forming eruption.     

Origin of Phenocrysts and Compositional Diversity in Pre-Mazama Rhyodacite Lavas, Crater Lake, Oregon

Phenocrysts in porphyritic volcanic rocks may originate in avariety of ways in addition to nucleation and growth in thematrix in which they are found. Porphyritic rhyodacite lavasthat underlie the eastern half of Mount Mazama, the High Cascadeandesite/dacite volcano that contains Crater Lake caldera, containevidence that bears on the general problem of phenocryst origin.Phenocrysts in these lavas apparently formed by crystallizationnear the margins of a magma chamber and were admixed into convectingmagma before eruption. About 20 km³ of pre-Mazama rhyodacite magma erupted during arelatively short period between400 and 500 ka; exposed pre-Mazamadacites are older and less voluminous. The rhyodacites formedas many as 40 lava domes and flows that can be assigned to threeeruptive groups on the basis of composition and phenocryst content.Phenocryst abundance decreases (from 32 to 8 vol.%) and SiO₂content increases (from 68 to 73 wt.%) in the apparent orderof eruption. Phenocrysts (plagioclase, orthopyroxene, augite,and Fe-Ti oxides) are commonly fragmental or form polycrystallineaggregates with interstitial glass. Discrete phenocrysts withcomplete euhedral outlines are rare except for small elongatedcrystals. The abundance of discrete phenocrysts increases withthat of aggregates. The grain-size of minerals in the aggregatescovers the range of discrete phenocrysts (0.2–4.2 mm).Rim compositions of phenocrysts and the range of chemical zoningare almost uniform among the three rhyodacite groups, regardlessof whether crystals are discrete or in aggregates. However,a small fraction of phenocrysts, especially small elongatedcrystals, have different compositions: plagioclase with Fe-richcores and augite with Wo-poor cores, both of which are characteristicof crystals in undercooled andesite enclaves in the rhyodacites.The majority of phenocrysts were derived by disintegration ofpolycrystalline aggregates; rare, small phenocrysts crystallizedin andesitic magma similar to that represented by the andesiteenclaves. The modal and chemical compositions of the rhyodacites can beexplained by different degrees of admixing of crystals, representedby the aggregates, into magma having 4 vol.% ‘true’phenocrysts, mainly plagioclase. The aggregates may be partsof the rind formed by in situ crystallization near the walland roof of the magma chamber. The rind was disrupted duringor just before eruption, and pieces were variably disaggregatedand incorporated into erupting magma. The amount of rind incorporateddeclined during the sequence of eruptions. Owing to vesiculationof interstitial liquid and shearing during flow, crystals inthe aggregates were separated and became phenocrysts. Pre-Mazamarhyodacite was erupted dominantly as lava, as opposed to thecompositionally similar rhyodacite pumice of the Holocene caldera-formingeruption of Mount Mazama, apparently because its source chamberwas crystallizing inward rather than actively growing.     

Petrology of Volcán Tequila,Jalisco, Mexico: disequilibrium phenocryst assemblages and evolution of the subvolcanic magma system

Volcán Tequila is an extinct stratovolcano in the western Mexican Volcanic Belt that has erupted lavas ranging from andesite to rhyolite during the last 0.9 Ma. Following an early period of rhyolitic volcanism, the main edifice of the volcano was constructed by central vent eruptions that produced 25 km³ of pyroxene-andesite. At about 0.2 Ma central activity ceased and numerous flows of hornblende-bearing andesite, dacite, and rhyodacite erupted from vents located around the flanks of the volcano. Bimodal plagioclase phenocryst rim compositions in lavas from both the main edifice and the flanks indicate that magma mixing commonly occurred shortly prior to or during eruption. Compositions of endmember magmas involved in mixing, as constrained by whole-rock major and trace element abundances, phenocryst compositions, and mineral-melt exchange equilibria, are similar to those of some lavas erupted from the central vent and on the flanks of the volcano. Estimated pre-eruptive temperatures for hornblende-bearing lavas (970°–830°C) are systematically lower than for lavas that lack hornblende (1045°–970°C), whereas magmatic H₂O contents are systematically higher for hornblende-bearing lavas. In addition to stabilizing hornblende, high magmatic water contents promoted crystallization of calcic plagioclase (An_70–82). Frequent injections of magma into the base of the subvolcanic plumbing system followed by eruption of mixed magma probably prevented formation of large volumes of silicic magma, which have caused paroxysmal, caldera-forming eruptions at other stratovolcanoes in western Mexico. The later stages of volcanic activity, represented by the flank lavas, indicate a change from a large magma storage reservoir to numerous small ones that developed along a NW-trending zone parallel to regional fault trends. Sr and Nd isotopic data for lavas from the Tequila region and other volcanoes in western Mexico demonstrate that differentiated calc-alkaline magmas are formed primarily through crystal fractionation of mantle-derived calc-alkaline basalt coupled with assimilation of crustal material. Present Address:Department of the Geophysical Sciences The University of Chicago, Chicago IL, 60637, USA     

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.     

Multiple Isotopic Components in Quaternary Volcanic Rocks of the Cascade Arc near Crater Lake, Oregon

Quaternary lavas and pyroclastic rocks of Mount Mazama, CraterLake caldera, and the surrounding area have variable Sr, Nd,and Pb isotopic compositions. High-alumina olivine tholeiites(HAOT) have ⁸⁷Sr/⁸⁶Sr ratios of 0.70346–0.70364; basalticandesite, 0–70349–0.70372; shoshonitic basalticandesite, 0.70374–0.70388; and andesite, 0.70324–0.70383.Dacites of Mount Mazama have ⁸⁷Sr/⁸⁶Sr ratios of 0.70348–0.70373.Most rhyodacites converge on 0.7037. However, rhyodacite ofthe caldera-forming, climactic eruption has ⁸⁷Sr/⁸⁶Sr=0.70354because of an admixed low-⁸⁷Sr/⁸⁶Sr component. Andesitic tomafic-cumulate scoriae of the climactic eruption, and enclavesin preclimactic rhyodacites, cluster in two groups but shownearly the entire ⁸⁷Sr/⁸⁶Sr range of the data set, confirmingpreviously suggested introduction of diverse parental magmasinto the growing climactic chamber. Pb and Nd isotope ratiosdisplay less variation (²⁰⁶Pb/²⁰⁴Pb= 18.838–18.967, ²⁰⁷Pb/²⁰⁴Pb=15.556–15.616,²⁰⁸Pb/²⁰⁴Pb=38.405–38.619; _Nd= +3.9 to +6.1) and generallycovary with ⁸⁷Sr/⁸⁶ Sr ratios. Radiogenic isotope data fromCrater Lake plot with published data for other Cascade volcanoeson isotope ratio correlation diagrams. The isotopic data for the Crater Lake area require sources ofprimitive magmas to consist of depleted mantle and a subductioncomponent, introduced in variable quantity to the depleted mantlewedge. Variable degrees of melting of this heterogeneous mantle,possibly at different depths, produced the diversity of isotopiccompositions and large-ion lithophile element (LILE) abundancesin primitive magmas. Trace element ratios do not indicate presenceof an ocean island basalt (OIB) source component that has beenreported in lavas of some other Cascade volcanoes. Crustal contamination may have affected isotope ratios and LILEconcentrations in evolved HAOT, where initial LILE concentrationswere low. Contamination is more difficult to detect in the calcalkalinelavas because of their higher LILE concentrations and the smallisotopic contrast with likely contaminants, such as mid- tolower-crustal rocks thought to be equivalents of igneous rocksof the Klamath Mountains and associated lower crust. Crustalassimilation appears to be required for calcalkaline rocks onlyby ¹⁸O values, which vary from lows of +5.6 to + 6.0% in HAOTand primitive basaltic andesites to a high of +7.0% in dacite,a range that is too high to be explained by plagioclase-dominatedclosed-system fractional crystallization. Elevated ¹⁸O valuesof differentiated lavas may be attributed to interaction withrelatively ¹⁸O-rich, ⁸⁷Sr-poor crustal rocks. Variably fused granitoid blocks ejected in the climactic eruption,and rarely in late Pleistocene eruptive units, have ¹⁸Opl of–3.4 to +6.5% and ¹⁸O_qz of –2.2 to +8.0% but haveSr, Nd, and Pb isotope ratios similar to volcanic rocks (e.g.⁸⁷Sr/⁸⁶Sr0.7037). Rb and Sr data for glass separates from granodioritessuggest that the source pluton is Miocene. Glass from granodioritehas ⁸⁷Sr/⁸⁶Sr ratios as high as 0.70617. Oxygen isotope fractionationbetween quartz, plagioclase, and glass indicates requilibrationof O isotopes at magmatic temperatures, after ¹⁸O/¹⁶O had beenlowered by exchange with meteoric hydrothermal fluids. Unmeltedgranodiorite xenoliths from pre-climactic eruptive units have¹⁸O values that are consistent with onset of hydrothermal exchangeearly during growth of the climactic magma chamber. Assimilationof such upper-crustal granodiorite apparently lowered ¹⁸O valuesof rhyodacites without significantly affecting their magmaticcompositions in other ways.     

Trace element geochemistry of high alumina basalt - Andesite - Dacite - Rhyodacite lavas of the Main Volcanic Series of Santorini Volcano,Greece

Trace element systematics throughout the cal-calkaline high alumina basalt — basaltic andesite — andesite — dacite — rhyodacite lavas and dyke rocks of the Main Volcanic Series of Santorini volcano, Greece are consistent with the crystal fractionation of observed phenocryst phases from a parental basaltic magma as the dominant mechanism involved in generating the range of magmatic compositions. Marked inflection points in several variation trends correspond to changes in phenocryst mineralogy and divide the Main Series into two distinct crystallisation intervals — an early basalt to andesite stage characterised by calcic plagioclase+augite+olivine separation and a later andesite to rhyodacite stage generated by plagioclase augite+hypersthene+magnetite+apatite crystallisation. Percent solidification values derived from ratios of highly incompatible trace elements agree with previous values derived from major element data using addition-subtraction diagrams and indicate that basaltic andesites represent 47–69%; andesites 70–76%; dacites ca. 80% and rhyodacite ca. 84% crystallisation of the initial basalt magma. Least squares major element mixing calculations also confirm that crystal fractionation of the least fractionated basalts could generate derivative Main Series lavas, though the details of the least squares solutions differ significantly from those derived from highly incompatible element and addition-subtraction techniques. Main Series basalts may result from partial melting of the mantle asthenosphere wedge followed by limited olivine+pyroxene+Cr-spinel crystallisation on ascent through the sub-Aegean mantle and may fractionate to more evolved compositions at pressures close to the base of the Aegean crust. Residual andesitic to rhyodacite magmas may stagnate within the upper regions of the sialic Aegean crust and form relatively high level magma chambers beneath the southern volcanic centres of Santorini. The eruption of large volumes of basic lavas and silicic pyroclastics from Santorini may have a volcanological rather than petrological explanation.     

Pre-eruptive magma mixing in ash-flow deposits of the Tertiary Rum Igneous Centre,Scotland

The Northern Marginal Zone of the Rum Igneous Centre is a remnant of an early caldera and its infill. It is composed of intra-caldera breccias and various small-volume pyroclastic deposits, overlain by prominent rhyodacite ash-flow sheets of up to 100 m thickness. The ash-flows were fed from a feeder system near the caldera ring-fault, and intrusive rhyodacite can locally be seen grading into extrusive deposits. A variety of features suggest that the ash-flows were erupted from a magma chamber that contemporaneously hosted felsic and mafic magmas: (i) chilled basaltic inclusions in rhyodacite; (ii) formerly glassy basaltic to andesitic enclaves with fluid-fluid relationships; (iii) feldspars with thick reaction rims enclosed in the basaltic to andesitic inclusions, yet with cores chemically resembling those of the rhyodacite: (iv) trace element compositions of the rhyodacite and the mafic enclaves form a mixing line between the end-member rhyodacite and basalt compositions. Additionally, textural and chemical features in the rhyodacite feldspar phenocrysts are consistent with magma mixing; (v) feldspars with resorption embayments cutting through internal zonation of the crystals; (vi) complexly zoned crystals with sieve-textured zones that are overgrown with euhedral zones; (vii) oscillatory zonation of feldspar phenocrysts in the rhyodacite, showing sharp increases in anorthite (An 10%) followed by gradual decrease in An-content (An 4%). This evidence points to eruption of ash-flows from a felsic magma chamber that was periodically replenished by injection of mafic magma. Diffusional mixing between the two magmas was permitted by temperature and compositional differences, but was slow due to the contrast in viscosities and densities. The Fe–Ti–P-enriched basic magma that replenished the chamber was degassing on entering the lower temperature environment and soon equilibrated thermally, followed by chemical exchange between the two end-member magmas. This process formed hybrid andesite enclaves enriched in trace elements beyond that caused by simple mixing, implying trace element diffusion in addition to bulk mixing. Eruption was caused by replenishment with, and degassing of, the basic magma and the chamber partially evacuated while the process of hybridisation was underway. The erupted products record magma mixing by chamber replenishment, blending of two magmas and elemental exchange in the magma chamber, and also physical mingling in the eruptive conduit.     

Rhyodacite magma storage conditions prior to the 3430 yBP caldera-forming eruption of Aniakchak volcano, Alaska

This study presents a pre-eruptive magma storage model for the rhyodacite and andesite magmas erupted during the 3430 yBP caldera-forming eruption of Aniakchak volcano, Alaska, derived from phase equilibria experiments and petrological data. The compositions of Fe–Ti oxide pairs from the early erupted Plinian rhyodacite pumice yield core temperatures of 871–900°C, with rims up to ∼942°C, and fO₂ from −10.6 to −11.8 log units. Melt inclusions entrapped in plagioclase phenocrysts have H₂O contents between 3 and 5 wt%, estimated by FTIR and electron microprobe volatiles by difference methods, with no detectable CO₂. Assuming water saturation, this corresponds to entrapment pressures between ∼65 and 150 MPa. Phase equilibria results reproduce the natural phase assemblages at of 95–150 MPa at 870–880°C, assuming water saturation. A mismatch in experimental versus natural glass SiO₂ and Al₂O₃, and MELTS models for H₂O-undersaturated conditions indicate that the rhyodacite may not have been H₂O saturated. MELTS models with and P _total of 125–150 MPa at 870–880°C reproduce the natural groundmass glass Al₂O₃ composition best, indicating the magma may have been slightly H₂O undersaturated. Those pressures correspond to storage at 4.5–5.4 km depth in the crust. MELTS models and VBD estimates from melt inclusions in titanomagnetite grains from the andesite indicate pre-eruptive conditions of ∼1,000°C and > 110 MPa, corresponding to a minimum residence depth of ∼4.1 km assuming water saturation or greater if the magma was H₂O undersaturated. Previous geochemical studies indicate separate histories of the two magmas, though they retain some evidence that they are ultimately related through fractional crystallization processes. Analogous to the 1912 Novarupta magmas, the rhyodacite and andesite presumably originated within the same crystal mush zone beneath the edifice, yet were separated laterally and underwent different degrees of crustal assimilation. The andesite must have resided in close proximity, with ascent occurring in response to movement of the rhyodacite, and resulting in extensive syn-eruptive mingling.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

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     

Polybaric Evolution of Calc-alkaline Magmas from Nisyros, Southeastern Hellenic Arc, Greece

The lavas of Nisyros were erupted between about 0?2 m.y B.P.and 1422 A.D., and range in composition from basaltic andesiteto rhyodacite. Most were erupted prior to caldera collapse (exactdate unknown), and the post-caldera lavas are petrographically(presence of strongly resorbed phenocrysts) and chemically (lowerTiO₂ K₂O, P₂O₅, and LIL elements) distinct from the pre-calderalavas. The pre-caldera lavas do not form a continuous seriessince lavas with SiO₂ contents between 60 and 66 wt.% are absent.Nevertheless, major element variations demonstrate that fractionalcrystalliz ation (involving removal of olivine, dinopyroxene,plagioclase, and Fe-Ti oxide from the basaltic andesites andandesites and plagioclase, clinopyroxene, hypersthene, Ti-magnetite,ilmenite, apatite, and zircon from the dacites and rhyodacites)played a major role in the evolution of the pre-caldera lavas.Several lines of evidence indicate that other processes werealso important in magma evolution: (1) Quantitative modelingof major element data shows that phenocryst phases of unlikelycomposi tion or unrealistic assemblages of phenocryst phasesare required to relate the dacites and rhyodacites to the basalticandesites and andesites; (2) The proportions of olivine andclinopyroxene required in quantitative models for the initialstages of evolution differ from those observed petrographicallyand this is not likely to reflect either differential ratesof crystal settling or the curvature of cotectics along whichliquids of basaltic andesite to andesite composition lie; (3)The concentrations of Rb, Cs, Ba, La, Sm, Eu, and Th in therhyod.acites are too high for these lavas to be related to thedacites by fractional crystallization alone; and (4) ⁸⁷Sr/⁸⁶Srratios for the andesites and rhyodacites are higher than thosefor the basaltic andesites and dacites, respectively. It isshown that fractional crystallization was accompanied by assimilation,and that magma mixing played a minor role (if any) in the evolutionof the pre-caldera lavas. Trace element and isotopic data indicatethat the andesites evolved from the basaltic andesites by AFCinvolving average crust or upper crust, whereas the rhyodacitesevolved from the dacites by AFC involving lower crust. Additionalevidence for polybaric evolution is provided by the occurrenceof distinct Ab-rich cores of plagioclase phenocrysts in thedacites and rhyodacites, which record a period of high pressurecrystallization, and by the occurrence of both normal and reverse-zonedphenocrysts in the basaltic andesites and andesites. Furthermore,calculated pressures of crystallization are {small tilde}8 kbfor the dacites and rhyodacites and 3?5–4 kb for the basalticandesites and andesites. It is concluded that the dacites andrhyodacites evolved via AFC from basaltic andesites and andesiteslargely in chambers sited near the base of the crust whereasthe basaltic andesites and andesites mostly evolved in chamberssited at mid-crustal levels. Eruption from different chambersexplains the compositional gap in the chemistry of the pre-calderalavas since eruptive products represent a more or less randomsampling of residual liquids which separate (via filter pressing)from bodies of crystallizing magma at various depths. Magmamixing was important in the evolution of the post-caldera lavas,but geochemical data require that these magmas evolved fromparental magmas which were derived from a more refractory sourcethan the parental magmas to the pre-caldera lavas. ^*Present address: Netherlands Energy Research Foundation (ECN), P.O. Box 1, 1755 ZG Petten, The Netherlands     

Origin of rhyolite and rhyodacite lavas and associated mafic inclusions of Cape Akrotiri,Santorini: the role of wet basalt in generating calcalkaline silicic magmas

Amphibole-bearing mafic inclusions (low to medium-K high-alumina basalt to basaltic andesite) comprise 4.1 vol% of calc-alkaline rhyolite and rhyodacite lavas on Akrotiri Peninsula, Santorini, Greece. Physical features indicate a magmatic origin for the inclusions, involving mingling with the host silicic magma and quenching. Water contents of the mafic magmas are estimated to have been above 4% at water pressures of 1.8 kbars or more at temperatures of approximately 950–1,000 °C. Three evolutionary stages are inferred in their petrogenesis. In the first stage infiltration of slab fluids promotes partial melting in the mantle to generate primitive wet basaltic magmas enriched in LREE, LILE, Th and U in comparison to N-type MORB. In the second stage storage and crystal differentiation of primitive magmas occurred in the lithospheric mantle or deep crust, involving olivine, spinel and clinopyroxene followed by amphibole and plagioclase. In the third stage differentiated mafic magma intrudes into porphyritic silicic magma at shallower crustal levels (estimated at 7–10 km). Mingling and quenching of the mafic magmas within the silicic host causes chemical or physical interactions between the inclusions and the host prior to and during eruption. The silicic lavas have geochemical affinities with the mafic inclusions, but are relatively depleted in MREE, HREE and Y and enriched in Rb relative to Ba and K. These observations are consistent with involvement of amphibole in magma genesis due either to crystal differentiation from wet basalt or to partial melting of mafic rocks with residual amphibole. Crystallization of wet basalt in the deep crust is preferred on the basis of physical considerations.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.Editorial responsibility: I. Parsons     

Enrichment of basalt and mixing of dacite in the rootzone of a large rhyolite chamber: inclusions and pumices from the Rattlesnake Tuff, Oregon

 A variety of cognate basalt to basaltic andesite inclusions and dacite pumices occur in the 7-Ma Rattlesnake Tuff of eastern Oregon. The tuff represents ∼280 km³ of high-silica rhyolite magma zoned from highly differentiated rhyolite near the roof to less evolved rhyolite at deeper levels. The mafic inclusions provide a window into the processes acting beneath a large silicic chamber. Quenched basaltic andesite inclusions are substantially enriched in incompatible trace elements compared to regional primitive high-alumina olivine tholeiite (HAOT) lavas, but continuous chemical and mineralogical trends indicate a genetic relationship between them. Basaltic andesite evolved from primitive basalt mainly through protracted crystal fractionation and multiple cycles (≥10) of mafic recharge, which enriched incompatible elements while maintaining a mafic bulk composition. The crystal fractionation history is partially preserved in the mineralogy of crystal-rich inclusions (olivine, plagioclase ± clinopyroxene) and the recharge history is supported by the presence of mafic inclusions containing olivines of Fo₈₀. Small amounts of assimilation (∼2%) of high-silica rhyolite magma improves the calculated fit between observed and modeled enrichments in basaltic andesite and reduces the number of fractionation and recharge cycles needed. The composition of dacite pumices is consistent with mixing of equal proportions of basaltic andesite and least-evolved, high-silica rhyolite. In support of the mixing model, most dacite pumices have a bimodal mineral assemblage with crystals of rhyolitic and basaltic parentage. Equilibrium dacite phenocrysts are rare. Dacites are mainly the product of mingling of basaltic andesite and rhyolite before or during eruption and to a lesser extent of equilibration between the two. The Rattlesnake magma column illustrates the feedback between mafic and silicic magmas that drives differentiation in both. Low-density rhyolite traps basalts and induces extensive fractionation and recharge that causes incompatible element enrichment relative to the primitive input. The basaltic root zone, in turn, thermally maintains the rhyolitic magma chamber and promotes compositional zonation. Received: 1 June 1998 / Accepted: 5 February 1999     

Petrologic and experimental evidence for the movement and heating of the pre-eruptive Minoan rhyodacite (Santorini, Greece)

Hydrothermal experiments combined with petrologic observations form the basis for a new two-stage model for the evolution of the pre-eruption Minoan magma chamber at Santorini, Greece. Ninety-nine percent of the erupted volume is two-pyroxene, rhyodacitic magma that had been stored at a temperature of ∼885 °C, based on magnetite-ilmenite and QUILF geothermometry. The rest of the volume is basaltic to andesitic magma, which occurs as <10 cm scoria clasts and as small inclusions in rhyodacite pumices. Petrologic observations show that these magmas mixed at different scales and at different times (i.e., multiple batches of mafic magma). Hydrothermal experiments were carried out on samples of rhyodacite and a mafic scoria in order to determine magma storage conditions and the mixing history of the two magmas. At 885 °C, the rhyodacite must have been stored at water-saturated pressures of ∼50 MPa, based on its phase assemblage, matrix-glass composition, and crystal content. However, glass inclusions inside rhyodacitic plagioclase phenocrysts contain more than 6 wt% H₂O, indicating they formed at pressures >200 MPa. In addition, the composition of the plagioclase hosts (An_{56 ± 6}) of the inclusions require temperatures of 825 ± 25 °C at pressures >200 MPa. This demonstrates that the Minoan rhyodacitic magma underwent a two-stage evolution, first crystallizing at ∼825^∘C and >200 MPa, and then rinsing to a shallow ∼50 MPa storage region with a concomitant rise in temperature to ∼885 °C. We suggest that the episodic intrusion of mafic magmas provided the necessary heat and perhaps contributed to the ascent of the magma to shallow crustal depths where it reequilibrated before the cataclysmic eruption. Phase equilibria suggest that much of the heating of the rhyodacite occurred in the shallow storage region. Thermal budget calculations suggest that the rhyodacite magma could have been heated by intrusions of basalt rising at reasonable upwelling rates and injected into the storage zone over several hundred years. Preservation of amphibole in the mafic scoria indicate that injection of mafic magma continued up until days before the cataclysmic eruption, perhaps triggering the event. Received: 30 September 1997 / Accepted: 5 October 1998     

Rates of Thermal and Chemical Evolution of Magmas in a Cooling Magma Chamber: a Chronological and Theoretical Study on Basaltic and Andesitic Lavas from Rishiri Volcano, Japan

Rates of magmatic processes in a cooling magma chamber wereinvestigated for alkali basalt and trachytic andesite lavaserupted sequentially from Rishiri Volcano, northern Japan, bydating of these lavas using ²³⁸U–²³⁰Th radioactive disequilibriumand ¹⁴C dating methods, in combination with theoretical analyses.We obtained the eruption age of the basaltic lavas to be 29·3± 0·6 ka by ¹⁴C dating of charcoals. The eruptionage of the andesitic lavas was estimated to be 20·2 ±3·1 ka, utilizing a whole-rock isochron formed by U–Thfractionation as a result of degassing after lava emplacement.Because these two lavas represent a series of magmas producedby assimilation and fractional crystallization in the same magmachamber, the difference of the ages (i.e. 9 kyr) is a timescaleof magmatic evolution. The thermal and chemical evolution ofthe Rishiri magma chamber was modeled using mass and energybalance constraints, as well as quantitative information obtainedfrom petrological and geochemical observations on the lavas.Using the timescale of 9 kyr, the thickness of the magma chamberis estimated to have been about 1·7 km. The model calculationsshow that, in the early stage of the evolution, the magma cooledat a relatively high rate (>0·1°C/year), and thecooling rate decreased with time. Convective heat flux fromthe main magma body exceeded 2 W/m² when the magma was basaltic,and the intensity diminished exponentially with magmatic evolution.Volume flux of crustal materials to the magma chamber and rateof convective melt exchange (compositional convection) betweenthe main magma and mush melt also decreased with time, from 0·1 m/year to 10^–3 m/year, and from 1 m/yearto 10^–2 m/year, respectively, as the magmas evolved frombasaltic to andesitic compositions. Although the mechanism ofthe cooling (i.e. thermal convection and/or compositional convection)of the main magma could not be constrained uniquely by the model,it is suggested that compositional convection was not effectivein cooling the main magma, and the magma chamber is consideredto have been cooled by thermal convection, in addition to heatconduction. KEY WORDS: convection; magma chamber; heat and mass transport; timescale; U-series disequilibria     

The Cold Bay Volcanic Center,Aleutian Volcanic Arc

The Cold Bay Volcanic Center has experienced two major stages of eruptive activity. Early (M-Series) acitivity produced bimodal Hi-Alumina basalt and calc-alkaline andesite lavas while later (FPK-Series) activity produced only calc-alkaline andesite. The spectrum of basalt compositions is believed to be due to high pressure (8 kb) fractionation at or near the base of the crust. Abundant mineralogical and geochemical evidence support a lower pressure mixing origin for all andesites. Inspection of the mineralogical data has shown that the earliest (M-Series) andesites were produced by mixing of basalt (<53 wt% SiO₂) and silicic andesite (60.5 to 62.5 wt%) while later (FPK-Series) andesites resulted from the mixing of basaltic-andesite (53 to 56 wt%) and less silicic andesite (58.5 to 60.0 wt%). The major element and trace element geochemical data are consistent with a low pressure fractionation origin for the silicic endmember magmas and support the temporal variations in both mafic and silicic endmember compositions. The complete lack of crustal inclusions in all lavas is taken as evidence for a minimal crustal melting and/or assimilation role in the origin of the silicic endmembers. Many of the features of all andesites, including the important long term convergence of endmember magma compositions, are consistent with the process of liquid fractionation, accompanied by large scale magma mixing. A deduced upper limit of 62.5 wt% SiO₂ for the silicic endmember magmas suggests that liquid fractionation, in the absence of major crustal melting, cannot produce more silicic magmas. A possible explanation is the presence of a rheological barrier, based on the concept of critical crystallinity (Marsh 1981), which prohibits more silicic liquids from being extracted from a crystal-liquid suspension.     

Geology and petrogenesis of the edgecumbe volcanic field,SE Alaska: The interaction of basalt and sialic crust

The Edgecumbe volcanic field is a Holocene volcanic province located on Kruzof Island, SE Alaska. Exposed within the 260 km² field are basalt, andesitic basalt, andesite, dacite and rhyodacite. The rhyodacites were erupted after the basalts and before the andesites. The volcanics, which are Al-rich (14–18 wt%) and lack an iron enrichment trend, range from tholeiites (47 wt% SiO₂) through rhyodacites (72%), but a compositional gap of approximately 9 wt% separates the dacites and rhyodacites. Initial ⁸⁷Sr/⁸⁶Sr ranges from 0.70297 in the basalts to 0.70440 in a pyroxene andesite. δ ¹⁸O increases across the suite: 5.8‰ to 7.9‰. Plagioclase (An_32–86) is the dominant phenocryst in all but one lava. Olivine (Fo_58–86) occurs in the basic lavas (<53 wt% SiO₂), but is replaced by orthopyroxene (En_43–73) and clinopyroxene (En₃₁Wo₄₁-En₄₈Wo₄₀) in the more siliceous volcanics. In the basalts and rhyodacites, plagioclase is weakly zoned, but extreme zoning (<30 mole% An) is characteristic of phenocrysts in the intermediate lavas. Fractionation of the observed phenocryst assemblages could not have produced the more silicous volcanics. Instead they were generated by partial melting of intrusive basement (⁸⁷Sr/⁸⁶Sr=0.70487; δ ¹⁸O: 8.7–9.3) by basaltic magma and subsequent assimilation. Mass balance calculations show the rhyodacites are almost pure partial melt (<5% basaltic component) whereas the intermediate lavas contain between 30 and 60% partial melt.     

Andesites from northeastern Kanaga Island,Aleutians

Kanaga island is located in the central Aleutian island arc. Northeastern Kanaga is a currently active late Tertiary to Recent calc-alkaline volcanic complex. Basaltic andesite to andesite lavas record three episodes (series) of volcanic activity. Series I and Series II lavas are all andesite while Series III lavas are basaltic andesite to andesite. Four Series II andesites contain abundant quenched magmatic inclusions ranging in composition from high-MgO low-alumina basalt to low-MgO highalumina basalt. The spectrum of lava compositions is due primarily to fractional crystallization of a parental low-MgO high-alumina basalt but with variable degrees of crustal contamination and magma mixing. The earliest Series I lavas represent mixing between high-alumina basalt and silicic andesite with maximum SiO₂ contents of 65–67 wt %. Later Series I and all Series II lavas are due to mixing of andesite magmas of similar composition. The maximum SiO₂ content of the pre-mixed andesites magmas is estimated at 60–63 wt %. The youngest lavas (Series III) are all non-mixed and have maximum estimated SiO₂ contents of 59 wt %. The earliest Series I lavas contain a significant crustal component while all later lavas do not. It is concluded that the maximum SiO₂ contents of silicic magmas, the contribution of crustal material to silicic magma generation, and the role of magma mixing all decrease with time. Furthermore, silicic magmas generated by fractional crystallization at this volcanic center have a maximum SiO₂ content of 63 wt %. All of these features have also been documented at the central Aleutian Cold Bay Volcanic Center (Brophy 1987). Based on data from these two centers a model of Aleutian calc-alkaline magma chamber development is proposed. The main features are: (1) a single low pressure magma chamber is continuously supplied by primitive low-alumina basalt; (2) non-primary high-alumina basalt is formed along the chamber margins by selective gravitational settling of olivine and clinopyroxene and retention of plagioclase; (3) sidewall crystallization accompanied by crustal melting produces buoyant silicic (>63 wt % SiO₂) liquids that pond at the top of the chamber, and; (4) continued sidewall crystallization, now isolated from the chamber wall, produces silicic liquids with 63 wt % SiO₂ that increase the thickness and lowers the overall SiO₂ content of the upper silicic zone. It is suggested that the maximum SiO₂ content of 63% imposed on fractionation-generated magmas is due to a rheological barrier that prohibits the extraction of more silicic liquids from a crystal-liquid mush along the chamber wall.     

A model of reverse differentiation at Dikii Greben' Volcano,Kamchatka: progressive basic magma vesiculation in a silicic magma chamber

Dikii Greben' Volcano is the largest modern volcano with silicic rocks in the Kurile-Kamchatka island arc. It consists of many domes and lava flows of rhyodacite, dacite and andesite which were erupted in a reverse differentiation sequence. Non-equilibrium phenocryst assemblages (quartz + Mg-rich olivine, An-rich + An-poor plagioclase etc.), abundance of chilled mafic pillows in the dacites and andesites, and linear variations of rock compositions in binary plots are considered as mineralogical, textural and geochemical evidence for mixing. Mafic pillows in volcanics have a lower density (because of high porosity) and contain the same non-equilibrium phenocryst assemblages as the host rocks. Their groundmass contains skeletal microlites of plagioclase and amphibole proving that the groundmass as well as the pillows themselves formed from a water-rich basaltic magma at depth. They are considered as supercooled, vesiculated floating drops of a hot hybrid layer in the magma chamber which formed after refilling. The lower density of the inclusions allows them to float in the host magma and to concentrate at the top of the chamber prior to eruption. Magma mingling was effected by mechanical disintegration of the inclusions in the host magma during eruption. The rhyodacitic and basic end-members of the mixing series cannot be linked by low-P fractionation though high-P, amphibole-rich fractionation is not excluded.