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     

Compositional Layering and Syn-eruptive Mixing of a Periodically Refilled Shallow Magma Chamber: the AD 79 Plinian Eruption of Vesuvius

A detailed study of the pyroclastic deposits of the AD 79 ‘Pompei’Plinian eruption of Vesuvius has allowed: (1) reconstructionof the thermal, compositional and isotopic (⁸⁷Sr/⁸⁶Sr) pre-eruptivelayering of the shallow magma chamber; (2) quantitative definitionof the syn-eruptive mixing between the different magmas occupyingthe chamber, and its relationships with eruption dynamics; (3)recognition of the variability of mafic magma batches supplyingthe chamber. During the different phases of the eruption 25–30%of the magma was ejected as white K-phonolitic pumice, and 70–75%as grey K-tephri-phonolitic pumice. The white pumice resultsfrom the tapping of progressively deeper magma from a body (T= 850–900%C) consisting of two distinct layers mainlyformed by crystal fractionation. The grey pumice results fromsyn-eruptive mixing involving three main end-members: the phonolitic‘white’ magmas (salic end-member, SEM), mafic cumulates(cumulate end-member, CEM) and a crystal-poor ‘grey’phono-tephritic magma (mafic end-member, MEM), which was nevererupted without first being mixed with ‘white’ magma.Evidence is provided that mixing occurred within the chamberand was characterized by a transition with time from physicalmixing at a microscopic scale to chemical hybridization. TheMEM magma had a homogeneous composition and constant ⁸⁷Sr⁸⁶Srisotopic ratio, possibly as a result of sustained convection.No unambiguous liquidus phases were found, suggesting that theMEM magma was superheated (T = 1000–1100C); its verylow viscosity was a main cause in the establishment of a physicaldiscontinuity separating the white and the grey magmas. Thewhite-grey boundary layer possibly consisted of a multiply diffusiveinterface, periodically broken and recreated, supplying thephonolitic body through mixing of moderate amounts of fractionatedgrey melts with the overlying white magma. The presence of alarge overheated mass indicates the young, growing stage ofthe AD 79 chamber, whose main engine was the periodic arrivalof hot mafic magma batches. These were characterized by K-tephriticto K-basanitic compositions, high temperatures (>1150C),high volatile contents (20–25% H₂O +Cl+F+S), low viscosities[(1+2 10² poises)] and relatively low densities (2500–2600kg/m³). The birth of the Pompei chamber followed the repeatedarrival of these batches (on average characterized by ⁸⁷Sr/⁸⁶Sr070729)into a reservoir containing a tephriticphonolitic, crystal-enriched,magma, a residue from the preceding ‘Avellino’ Plinianeruption (3400 BP).In fact, about half of magma ejected duringthe AD 79 eruption could have been inherited from pre-Avellinotimes. KEY WORDS: Vesuvius; magma chamber; magma mixing; compositional layering phonolites; magma supply; potassic magmas ^*Correponding author     

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

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

Evidence of Magma Mixing in the 'Daly Gap' of Alkaline Suites: a Case Study from the Enclaves of Pantelleria (Italy)

The island of Pantelleria consists of trachytes, pantelleritesand minor mildly alkaline basalts. Rocks of intermediate composition(falling in the so-called ‘Daly Gap’) such as mugearites,benmoreites and mafic trachytes occur only in the form of enclavesin trachytes and pantellerites inside the main caldera of theisland (Caldera ‘Cinque Denti’), which collapsedduring the ‘Green Tuff’ ignimbrite eruption at 50ka. The enclaves include volcanic, subvolcanic and intrusiverock types. The enclaves in host trachyte contain traces ofglass; devitrified glass occurs within enclaves in host pantellerites.Minerals in the enclaves show regular compositional variationswith whole-rock silica content. Glass present in the medium-grainedsamples is interpreted to be the result of incipient melting.The major and trace element compositions of the enclaves showregular and linear variations between an evolved mafic magma(hawaiite) and a felsic end-member similar to the ‘GreenTuff’ trachyte. Fractional crystallization modelling ofcompatible and incompatible trace elements (V, Ni, Zr, La, Sm,Lu, Nb, Y, Th) does not reproduce the observed trends. Rocksof intermediate composition within the ‘Daly Gap’can be explained only by magma mixing between an already differentiatedmafic magma (hawaiite) and an anorthoclase-rich trachytic meltin the lower and higher parts, respectively, of a stratifiedmagmatic chamber. Medium-grained enclaves are interpreted asthe result of fragmentation of solidified mixing layers in theroof of the magma chamber during the eruption of the ‘GreenTuff’, when the collapse of the caldera took place. Diffusioncalculations suggest a residence time of <5 days for theenclaves in their host magmas. KEY WORDS: Daly Gap; enclaves; magma mixing; Pantelleria     

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

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

Magma Evolution of the Sete Cidades Volcano, Sao Miguel, Azores

The Sete Cidades volcano (São Miguel, Azores) is situatedat the eastern end of the ultraslow spreading Terceira riftaxis. The volcano comprises several dominantly basaltic pre-calderaeruptions, a trachytic caldera-forming stage and a post-calderastage consisting of alternating trachytic and basaltic eruptions.The post-caldera flank lavas are more primitive (>5 wt %MgO) than the pre-caldera lavas, implying extended fractionalcrystallization and longer crustal residence times for the pre-caldera,shield-building lavas. Thermobarometric estimates show thatthe ascending alkali basaltic magmas stagnated and crystallizedat the crust–mantle boundary (15 km depth), whereas themore evolved magmas mainly fractionated in the upper crust (3km depth). The caldera-forming eruption was triggered by a basalticinjection into a shallow trachytic magma chamber. Lavas fromall stages follow a single, continuous liquid line of descentfrom alkali basalt to trachyte, although slight differencesin incompatible element (e.g. Ba/Nb, La/Nb) and Sr isotope ratiosimply some heterogeneity of the mantle source. Major and traceelement data suggest similar partial melting processes throughoutthe evolution of the volcano. Slight geochemical differencesbetween post- and pre-caldera stage lavas from the Sete Cidadesvolcanic system indicate a variation in the mantle source compositionwith time. The oxygen fugacity increased from the pre-calderato the post-caldera stage lavas, probably as a result of theassimilation of crustal rocks; this is supported by the presenceof crustal xenoliths in the lavas of the flank vents. The lavasfrom the Sete Cidades volcano generally have low Sr isotoperatios; however, rocks from one post-caldera vent on the westernflank indicate mixing with magmas resembling the lavas fromthe neighbouring Agua de Pau volcano, having higher Sr isotoperatios. The different magma sources at Sete Cidades and theadjacent Agua de Pau volcano imply that, despite their closeproximity, there is only limited interaction between them. KEY WORDS: crystallization depth; fractionation; stratigraphy; Terceira rift; volcanic stages     

Zircon Hf Isotopic Evidence for Mixing of Crustal and Silicic Mantle-derived Magmas in a Zoned Granite Pluton, Eastern Australia

Zircon Hf isotopic data from a zoned pluton of the Moonbi supersuite,New England batholith, eastern Australia, are consistent withmagma mixing between two silicic melts, each derived from isotopicallydistinct sources. Although zircons from three zones within theWalcha Road pluton give a U–Pb crystallization age of249 ± 3 Ma, zircon populations from each zone have arange in _Hf. Zircons from the mafic hornblende–biotitemonzogranite pluton margin and intermediate zones have _Hf +5to +11, whereas those from the more felsic centre of the plutonhave _Hf +7 to +16, representing a total variation of 11 _Hfunits. The Lu–Hf depleted mantle model ages range from650 to 250 Ma, with the younger zircons present only in thefelsic pluton centre. The variation in _Hf indicates the involvementof silicic melts from at least two sources, one a crustal componentwith a Neoproterozoic model age and the other a primitive mantle-derivedcomponent with model ages similar to the U–Pb crystallizationage of the pluton. The zircons reflect the isotopic compositionsof the different proportions of crustal-derived silicic melt,relative to mantle-derived silicic melt, between melt generationand final pluton construction. The Walcha Road pluton is consideredto have formed by incremental assembly of progressively morefelsic melt batches resulting from mixing, replenishment andcrystal–melt separation, with final pluton constructioninvolving mechanical concentration as zones of crystal mush.The zoned pluton and, more broadly, the Moonbi supersuite provideexamples of magma mixing by which the more silicic units havemore juvenile isotopic compositions as a result of increasingproportions of residual melt from basalt fractionation, relativeto crustal partial melt. KEY WORDS: Australia; granite magma mixing; zircon; zoned pluton; Hf isotopes     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

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     

Phase Equilibria Constraints on the Chemical and Physical Evolution of the Campanian Ignimbrite

The Campanian Ignimbrite is a > 200 km³ trachyte–phonolitepyroclastic deposit that erupted at 39·3 ± 0·1ka within the Campi Flegrei west of Naples, Italy. Here we testthe hypothesis that Campanian Ignimbrite magma was derived byisobaric crystal fractionation of a parental basaltic trachyandesiticmelt that reacted and came into local equilibrium with smallamounts (5–10 wt%) of crustal rock (skarns and foid-syenites)during crystallization. Comparison of observed crystal and magmacompositions with results of phase equilibria assimilation–fractionationsimulations (MELTS) is generally very good. Oxygen fugacitywas approximately buffered along QFM + 1 (where QFM is the quartz–fayalite–magnetitebuffer) during isobaric fractionation at 0·15 GPa ( 6km depth). The parental melt, reconstructed from melt inclusionand host clinopyroxene compositions, is found to be basaltictrachyandesite liquid (51·1 wt% SiO₂, 9·3 wt%MgO, 3 wt% H₂O). A significant feature of phase equilibria simulationsis the existence of a pseudo-invariant temperature, 883 °C,at which the fraction of melt remaining in the system decreasesabruptly from 0·5 to < 0·1. Crystallizationat the pseudo-invariant point leads to abrupt changes in thecomposition, properties (density, dissolved water content),and physical state (viscosity, volume fraction fluid) of meltand magma. A dramatic decrease in melt viscosity (from 1700Pa s to 200 Pa s), coupled with a change in the volume fractionof water in magma (from 0·1 to 0·8) and a dramaticdecrease in melt and magma density acted as a destabilizingeruption trigger. Thermal models suggest a timescale of 200kyr from the beginning of fractionation until eruption, leadingto an apparent rate of evolved magma generation of about 10^–3km³/year. In situ crystallization and crystal settling in density-stratifiedregions, as well as in convectively mixed, less evolved subjacentmagma, operate rapidly enough to match this apparent volumetricrate of evolved magma production. KEY WORDS: assimilation; Campanian Ignimbrite; fractional crystallization; magma dynamics; phase equilibria     

Lower Crustal Magma Genesis and Preservation: a Stochastic Framework for the Evaluation of Basalt-Crust Interaction

We present a quantitative assessment of the thermal and dynamicresponse of an amphibolitic lower crust to the intrusion ofbasaltic dike swarms in an arc setting. We consider the effectof variable intrusion geometry, depth of intrusion, and basaltflux on the production, persistence, and interaction of basalticand crustal melt in a stochastic computational framework. Distinctmelting and mixing environments are predicted as a result ofthe crustal thickness and age of the arc system. Shallow crustal(30 km) environments and arc settings with low fluxes of mantle-derivedbasalt are likely repositories of isolated pods of mantle andcrustal melts in the lower crust, both converging on daciticto rhyodacitic composition. These may be preferentially rejuvenatedin subsequent intrusive episodes. Mature arc systems with thickercrust (50 km) produce higher crustal and residual basaltic meltfractions, reaching 0·4 for geologically reasonable basaltfluxes. The basaltic to basaltic andesite composition of bothcrustal and mantle melts will facilitate mixing as the networkof dikes collapses, and Reynolds numbers reach 10^–4–1·0in the interiors of dikes that have been breached by ascendingcrustal melts. This may provide one mechanism for melting, assimilation,storage and homogenization (MASH)-like processes. Residual mineralassemblages of crust thickened by repeated intrusion are predictedto be garnet pyroxenitic, which are denser than mantle peridotiteand also generate convective instabilities where some of thecrustal material is lost to the mantle. This reconciles thethinner than predicted crust in regions that have undergonea large flux of mantle basalt for a prolonged period of time,and helps explain the enrichment of incompatible elements suchas K₂O, typical of mature arc settings, without the associatedmass balance problem. KEY WORDS: crustal anatexis; delamination; lower crust; magma mixing; thermal model     

Shift and Rotation of Composition Trends by Magma Mixing: 1983 Eruption at Miyake-jima Volcano, Japan

Pre-eruption processes are investigated for magmas erupted in1983 from Miyake-jima volcano, which is one of the most activevolcanoes in Japan. The whole-rock compositional trends of theeruptive products are principally smooth and linear. Magmaserupted from some fissures have compositions that deviate fromthe main linear trend. Phenocryst contents of samples displacedfrom the linear compositional trends are significantly lowerthan those of samples on the main trends. Anorthite-rich plagioclasephenocrysts, present throughout the 1983 products, are too calcicto have crystallized from the erupted magma composition, andwere derived from a basaltic magma through magma mixing. Althoughthe linear whole-rock composition trends favor simple two-componentmagma mixing, this cannot explain the presence of samples thatdeviate from the main trend. Instead, the observed compositiontrends were formed by mixing of a homogeneous basaltic magmawith andesitic magmas exhibiting compositional diversity. Theoriginal linear composition trends of the andesitic end-membermagma were rotated and shifted to the direction of the basalticend-member magma by magma mixing. The samples out of the maintrends represent magmas with less basaltic component than thoseon the trend. The density and viscosity of the basaltic end-membermagma were comparable with those of the andesitic end-membermagmas. The basaltic magma, discharged from one magma chamberat 2 kbar pressure, was injected into a magma chamber at lowerpressure occupied by the chemically zoned andesite magma (1kbar), and possibly as a fountain. To establish the characteristicmixing trend of the 1983 magma, the basaltic component musthave been distributed systematically in the zoned andesite magma.A requirement is that the basaltic magma spread laterally andmixed with the andesite magma at various levels of ascent ofthe fountain in the host andesite magma. Analysis of compositionalzoning in titanomagnetite crystals revealed that the eruptionof the 1983 magmas was initiated soon after the replenishmentof the basaltic magma in the 1 kbar magma chamber. KEY WORDS: compositional trend; liquid–liquid blending; magma chamber; magma mixing; Miyake-jima Volcano     

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     

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

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

The Lilloise Intrusion, East Greenland: Fractionation of a Hydrous Alkali Picritic Magma

The Lilloise is an 8 km4 km layered mafic intrusion which cutsthe plateau basalts of the East Greenland Tertiary province.Lilloise was intruded at 50 Ma, 4–5 Ma after cessationof the voluminous tholeiitic magmatism which accompanied riftingof the East Greenland continental margin. Lilloise is unusualamong layered intrusions in the province because it had a hydrousalkali picrite parent magma and generated a late-stage effluxof magmatic water from the intrusion into the aureole rocks.The three major subdivisions of the layered rocks are: olivine-clinopyroxene,olivine-clinopyroxene-plagioclase and plagioclase-amphibolecumulates. Massive subsidence of the intrusion before completesolidification resulted in deformation of the internal layeringand downturn of the bedding in the surrounding basalts. A strikingfeature of the intrusion is the injection of the layered rocksby a plexus of magmatic sheets which formed at the time of subsidence.The composition of these sheets is representative of the fractionationtrend of the intrusion and ranges from hawaiite to mildly saturatedquartz trachyte. The fractionation trend is successfully explainedby extraction of cumulus minerals of the layered rocks froma parent magma represented by alkali picrite dykes of a contemporaneousregional dyke swarm. Saturated to mildly over-saturated syenitesare a major component of the East Greenland province and theLilloise intrusion is illustrative of an important magmatictrend towards such compositions at this stage in the openingof the North Atlantic. KEY WORDS: Lilloise intrusion; East Greenland; alkali picrite magma; layered intrusion; magmatic differentiation ^*Corraponding author     

Compositional Zoning of the Bishop Tuff

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

Fountains in Magma Chambers

Cyclic layering is a common feature of the ultramafic zone oflayered intrusions and is usually attributed to the entry ofnew pulses of dense magma into the chamber. Since the crystallizationof olivine and bronzite lowers the density of the magma, a newpulse of the parent magma will be denser than the fractionatedmagma in the chamber. If the new pulse enters with excess momentumit will initially rise up into the host magma to form a fountain,then fall back around the feeder when negative buoyancy forcesovercome the initial momentum of the pulse. Laboratory experimentsusing aqueous solutions with both point and line sources havebeen conducted to obtain a quantitative understanding of thefluid-dynamical processes that are important in fountains. Itis observed that convection within the fountain is highly turbulent,resulting in appreciable entrainment of the host magma. A gravity-stratifiedhybrid layer develops at the floor and this breaks up into aseries of double-diffusive convecting layers if the new pulseis hotter than the host magma. The number of layers that formdepends on a number of factors, especially R, the ratio of thecontributions of composition and heat to the total density differencebetween the host magma and the new pulse. Raising the valueof R, results in the formation of more, thinner layers. The thickness of the hybrid layer at any time t is given byH = h₀+(V₀/A)t where V₀ is the volume flux through the feederand A is the horizontal area of the chamber. h₀ is related tothe initial steady-state height of the fountain and, for a linesource, is given by h₀=CU₀^4/3 d^–1(g/)^–2/3 whereU₀ is the volume flux per unit length, g is the accelerationdue to gravity, d is the width of the feeder, is the densityof the host magma, is the density difference between the magmasand C is a constant. Calculations based on these results and the consideration ofthe flow in the feeder dykes below the chamber indicate thata fountain will rise at least 350 m in a continental magma chamberif the feeder width is greater than 10 m. This will lead toextensive mixing between the new pulse and the fractionatedmagma in the chamber, producing a zoned hybrid layer at thefloor that is commonly over 1000 m thick. If the chamber receivesmany pulses of dense magma, the resulting zoning may persistthroughout much of the life of the chamber, especially if thefirst pulse to enter becomes contaminated by light magma releasedby melting at the margins. The highest Mg/Fe ratio for olivineand pyroxenes from cyclic units from the ultramafic zones oflayered intrusions is often well below the value expected forminerals crystallizing from a melt derived directly from themantle, supporting the hypothesis that new pulses of dense magmacan mix extensively with the fractionated magma in the chamber. The feeder dykes to some oceanic magma chambers, such as theBay of Islands Ophiolite, are believed to be narrower, so thatfountains do not rise more than a few metres above the floorof the chamber. This restricts mixing between the input magmaand the host magma and can result in the formation of a hybridzone that is only a few metres thick.     

The 1996 Eruption of Karymsky Volcano, Kamchatka: Historical Record of Basaltic Replenishment of an Andesite Reservoir

The simultaneous eruption in 1996 of andesite from Karymskyvolcano and of basalt from the Academy Nauk vent 6 km away appearsto provide a case of mafic recharge of an andesite reservoirfor which the time of recharge is exactly known and direct samplesof the recharging magma are available. The explosive phreato-magmaticeruption of basalt was terminated in less than 24 h, whereasandesite erupted continuously during the following 4 years.Detailed petrological study of volcanic ash, bombs and lavasof Karymsky andesite erupted during the period 1996–1999provides evidence for basaltic replenishment at the beginningof the eruptive cycle, as well as a record of compositionalvariations within the Karymsky magma reservoir induced by basalticrecharge. Shortly after the beginning of the eruption the compositionof the matrix glass of the Karymsky tephra became more maficand then, within 2 months, gradually returned to its originalstate and remained almost constant for the following 3 years.Further evidence for basaltic replenishment is provided by thepresence of xenocrysts of basaltic origin in the andesite ofKarymsky. A conspicuous portion of the plagioclase phenocrystsin the Karymsky andesite has calcic cores, with compositionsand textures resembling those of plagioclases in the AcademyNauk basalt. Similarly, the earlier portion of the andesiteof the eruption sequence contains rare olivines, which occuras resorbed cores in pyroxenes. The composition of the olivinematches that of olivines in the Academy Nauk basalt. The sequenceof events appears to be: (1) injection of basaltic magma intothe Karymsky chamber with immediate, compensating expulsionof pre-existing chamber magma from the Karymsky central vent;(2) direct mixing of basaltic and andesitic magmas with dispersalof phenocrysts associated with the basalt throughout the andesiteso that newly mixed magma appeared at the vent within 2 months;(3) re-establishment of thermal and chemical equilibrium withinthe reservoir involving crystallization in the new hybrid liquid,which returned the melt composition to ‘normal’,formed rims on inherited calcic plagioclase, and caused theresorption of dispersed olivine xenocrysts. Taken together,these findings indicate that the Karymsky magma reservoir wasrecharged by basalt at the onset of the 1996 eruptive cycle.The rapidity and thoroughness of mixing of the basalt with thepre-existing andesite probably reflects the modest contrastin temperature, viscosity, and density between the two magmas. KEY WORDS: Karymsky; Kamchatka; magma mixing; andesite; volcanic glass; plagioclase     

The Petrology and Geochemistry of Oto-Zan Composite Lava Flow on Shodo-Shima Island, SW Japan: Remelting of a Solidified High-Mg Andesite Magma

The Oto-Zan lava in the Setouchi volcanic belt is composed ofphenocryst-poor, sparsely plagioclase-phyric andesites (sanukitoids)and forms a composite lava flow. The phenocryst assemblagesand element abundances change but Sr–Nd–Pb isotopiccompositions are constant throughout the lava flow. The sanukitoidat the base is a high-Mg andesite (HMA) and contains Mg- andNi-rich olivine and Cr-rich chromite, suggesting the emplacementof a mantle-derived hydrous (7 wt % H₂O) HMA magma. However,Oto-Zan sanukitoids contain little H₂O and are phenocryst-poor.The liquid lines of descent obtained for an Oto-Zan HMA at 0·3GPa in the presence of 0·7–2·1 wt % H₂Osuggest that mixing of an HMA magma with a differentiated felsicmelt can reasonably explain the petrographical and chemicalcharacteristics of Oto-Zan sanukitoids. We propose a model wherebya hydrous HMA magma crystallizes extensively within the crust,resulting in the formation of an HMA pluton and causing liberationof H₂O from the magma system. The HMA pluton, in which interstitialrhyolitic melts still remain, is then heated from the base byintrusion of a high-T basalt magma, forming an H₂O-deficientHMA magma at the base of the pluton. During ascent, this secondaryHMA magma entrains the overlying interstitial rhyolitic melt,resulting in variable self-mixing and formation of a zoned magmareservoir, comprising more felsic magmas upwards. More effectiveupwelling of more mafic, and hence less viscous, magmas througha propagated vent finally results in the emplacement of thecomposite lava flow. KEY WORDS: high-Mg andesite; sanukitoid; composite lava; solidification; remelting     

Magma Generation at a Large, Hyperactive Silicic Volcano (Taupo, New Zealand) Revealed by U-Th and U-Pb Systematics in Zircons

Young (<65 ka) explosive silicic volcanism at Taupo volcano,New Zealand, has involved the development and evacuation ofseveral crustal magmatic systems. Up to and including the 26·5ka 530 km³ Oruanui eruption, magmatic systems were contemporaneousbut geographically separated. Subsequently they have been separatedin time and have vented from geographically overlapping areas.Single-crystal (secondary ionization mass spectrometry) andmultiple-crystal (thermal ionization mass spectrometry) zirconmodel-age data are presented from nine representative eruptiondeposits from 45 to 3·5 ka. Zircon yields vary by threeorders of magnitude, correlating with the degrees of zirconsaturation in the magmas, and influencing the spectra of modelages. Two adjacent magma systems active up to 26·5 kashow wholly contrasting model-age spectra. The smaller systemshows a simple unimodal distribution. The larger system, usingdata from three eruptions, shows bimodal model-age spectra.An older 100 ka peak is interpreted to represent zircons (antecrysts)derived from older silicic mush or plutonic rocks, and a youngerpeak to represent zircons (phenocrysts) that grew in the magmabody immediately prior to eruption. Post-26·5 ka magmabatches show contrasting age spectra, consistent with a mixtureof antecrysts, phenocrysts and, in two examples, xenocrystsfrom Quaternary plutonic and Mesozoic–Palaeozoic metasedimentaryrocks. The model-age spectra, coupled with zircon-dissolutionmodelling, highlight contrasts between short-term silicic magmageneration at Taupo, by bulk remobilization of crystal mushand assimilation of metasediment and/or silicic plutonic basementrocks, and the longer-term processes of fractionation from crustallycontaminated mafic melts. Contrasts between adjacent or successivemagma systems are attributed to differences in positions ofthe source and root zones within contrasting domains in thequartzo-feldspathic (<15 km deep) crust below the volcano. KEY WORDS: zircon; U-series dating; rhyolite; Taupo Volcanic Zone; Taupo volcano