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     

The 1982-83 Eruption at Galunggung Volcano, Java (Indonesia): Oxygen Isotope Geochemistry of a Chemically Zoned Magma Chamber

Mount Galunggung is a historically active volcano in southwesternJava that has erupted four times in the last two centuries.During the most recent event, which occurred during a 9–monthinterval in 1982– 83, some 305 10⁶ m³ of medium–K,calc–alkaline magma was erupted. This eruption was unusualbecause of its duration, the diversity of eruption dynamicsand products, and the range of lava compositions produced. Thecomposition of juvenile material changed gradually during thecourse of the eruption from initial plagioclase (An_60–75)and two–pyrozene bearing andesites with 58% SiO₂ to finalplagioclase (An_85–90), diopside, and olivine (Fo_85–90)bearing primitive magnesium basalts with 47% SiO₂ Mineralogicaland compositional relationships indicate a magmatic evolutioninvolving differentitation of high–Mg parental melt. Theeruptive volumes of 35 10⁶ m³ andesite, 120 10⁶ m³ maficandesite, and 150 10⁶ m³ basalt are consistent with the ideathat the 1982– 83 eruption progressively tapped and draineda magma chamber that had become chemically stratified throughextensive crystal fractionation. Separates of plagioclase and pyroxene have ¹⁸O( SMO W) rangesof + 5. 6 to + 6.0 and + 5.3 to + 5.6, respectively, with ¹⁸O_plag–pxvalues of + 0.4 to + 0.6o, indicating internal O–isotopeequiliburium at temperature of 1100–850 C. The magenesianbasalts have magmatic ¹⁸O/ ¹⁶O ratios similar to those of mid–oceanridge basalt, and the O–isotope ratios of compositionallyevolved derivative melts show no evidence for contaminationof the galunggung magmas by ¹⁸O–rich crust during differentiation.Andesites and transitional mafic and sites have a more variableO–isotope character, with laves and phenocrysts havingboth higher and lower ¹⁸O values than observed in the parentalmagnesium basalts. These features are interpreted to reflectintramagma chamber processes affecting the upper portions ofthe differentiating Galunggung magma body before the 1982–83eruption.     

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     

Petrogenesis of the Swaziland and Northern Natal Rhyolites of the Lebombo Rifted Volcanic Margin, South East Africa

The Jozini and Mbuluzi rhyolites and Oribi Beds of the southernLebombo Monocline, southeastern Africa, have geochemical characteristicsthat indicate they were derived by partial melting of a mixtureof high-Ti/Zr and low-Ti/Zr Sabie River Basalt Formation types.Compositional variations within the different rhyolite typescan largely be explained by subsequent fractional crystallization.The Sr- and Nd-isotope composition of the rhyolites is uniqueamongst Gondwana silicic large igneous provinces, having _Ndvalues close to Bulk Earth (–0·94 to 0·35)and low, but more variable, initial ⁸⁷Sr/⁸⁶Sr ratios (0·7034–0·7080).Quartz phenocryst ¹⁸O values indicate that the rhyolite magmashad ¹⁸O values between 5·3 and 6·7, consistentwith derivation from a basaltic protolith with ¹⁸O values between4·8 and 6·2. The low-¹⁸O rhyolites (< 6·0)come from the same stratigraphic horizon and are overlain andunderlain by rhyolites with more ‘normal’ ¹⁸O magmavalues. These low-¹⁸O rhyolites cannot have been produced byfractional crystallization or partial melting of mantle-derivedbasaltic material. The rhyolites have low water contents, makingit unlikely that the low ¹⁸O values are the result of post-emplacementalteration. Modification of the source by fluid–rock interactionat elevated temperatures is the most plausible mechanism forlowering the ¹⁸O magma value. It is proposed that the low-¹⁸Orhyolites were derived by melting of earlier altered rhyolitein calderas situated to the east, which were not preserved afterGondwana break-up. KEY WORDS: rhyolite; Lebombo; stable and radiogenic isotopes; low-¹⁸O magmas; partial melting     

Sampling the Cape Verde Mantle Plume: Evolution of Melt Compositions on Santo Antao, Cape Verde Islands

The volcanic history of Santo Antão, NW Cape Verde Islands,includes the eruption of basanite–phonolite series magmasbetween 7·5 and 0·3 Ma and (melilite) nephelinite–phonoliteseries magmas from 0·7 to 0·1 Ma. The most primitivevolcanic rocks are olivine ± clinopyroxene-phyric, whereasthe more evolved rocks have phenocrysts of clinopyroxene ±Fe–Tioxide ± kaersutite ± haüyne ± titanite± sanidine; plagioclase occurs in some intermediate rocks.The analysed samples span a range of 19–0·03% MgO;the most primitive have 37–46% SiO₂, 2·5–7%TiO₂ and are enriched 50–200 x primitive mantle in highlyincompatible elements; the basanitic series is less enrichedthan the nephelinitic series. Geochemical trends in each seriescan be modelled by fractional crystallization of phenocrystassemblages from basanitic and nephelinitic parental magmas.There is little evidence for mineral–melt disequilibrium,and thus magma mixing is not of major importance in controllingbulk-rock compositions. Mantle melting processes are modelledusing fractionation-corrected magma compositions; the modelssuggest 1–4% partial melting of a heterogeneous mantleperidotite source at depths of 90–125 km. Incompatibleelement enrichment among the most primitive magma types is typicalof HIMU OIB. The Sr, Nd and Pb isotopic compositions of theSanto Antão volcanic sequence and geochemical characterchange systematically with time. The older volcanic rocks (7·5–2Ma) vary between two main mantle source components, one of whichis a young HIMU type with ²⁰⁶Pb/²⁰⁴Pb = 19·88, 7/4 =–5, 8/4 0, ⁸⁷Sr/⁸⁶Sr = 0·7033 and ¹⁴³Nd/¹⁴⁴Nd= 0·51288, whereas the other has somewhat less radiogenicSr and Pb and more radiogenic Nd. The intermediate age volcanicrocks (2–0·3 Ma) show a change of sources to two-componentmixing between a carbonatite-related young HIMU-type source(²⁰⁶Pb/²⁰⁴Pb = 19·93, 7/4 = –5, 8/4 = –38,⁸⁷Sr/⁸⁶Sr = 0·70304) and a DM-like source. A more incompatibleelement-enriched component with 7/4 > 0 (old HIMU type) isprominent in the young volcanic rocks (0·3–0·1Ma). The EM1 component that is important in the southern CapeVerde Islands appears to have played no role in the petrogenesisof the Santo Antão magmas. The primary magmas are arguedto be derived by partial melting in the Cape Verde mantle plume;temporal changes in composition are suggested to reflect layeringin the plume conduit. KEY WORDS: radiogenic isotopes; geochemistry; mantle melting; Cape Verde     

Genesis of Ultramafic Lamprophyres and Carbonatites at Aillik Bay, Labrador: a Consequence of Incipient Lithospheric Thinning beneath the North Atlantic Craton

Numerous dykes of ultramafic lamprophyre (aillikite, mela-aillikite,damtjernite) and subordinate dolomite-bearing carbonatite withU–Pb perovskite emplacement ages of 590–555 Ma occurin the vicinity of Aillik Bay, coastal Labrador. The ultramaficlamprophyres principally consist of olivine and phlogopite phenocrystsin a carbonate- or clinopyroxene-dominated groundmass. Ti-richprimary garnet (kimzeyite and Ti-andradite) typically occursat the aillikite type locality and is considered diagnosticfor ultramafic lamprophyre–carbonatite suites. Titanianaluminous phlogopite and clinopyroxene, as well as comparativelyAl-enriched but Cr–Mg-poor spinel (Cr-number < 0.85),are compositionally distinct from analogous minerals in kimberlites,orangeites and olivine lamproites, indicating different magmageneses. The Aillik Bay ultramafic lamprophyres and carbonatiteshave variable but overlapping ⁸⁷Sr/⁸⁶Sr_i ratios (0·70369–0·70662)and show a narrow range in initial _Nd (+0·1 to +1·9)implying that they are related to a common type of parentalmagma with variable isotopic characteristics. Aillikite is closestto this primary magma composition in terms of MgO (15–20wt %) and Ni (200–574 ppm) content; the abundant groundmasscarbonate has ¹³C_PDB between –5·7 and –5,similar to primary mantle-derived carbonates, and ¹⁸O_SMOW from9·4 to 11·6. Extensive melting of a garnet peridotitesource region containing carbonate- and phlogopite-rich veinsat 4–7 GPa triggered by enhanced lithospheric extensioncan account for the volatile-bearing, potassic, incompatibleelement enriched and MgO-rich nature of the proto-aillikitemagma. It is argued that low-degree potassic silicate to carbonatiticmelts from upwelling asthenosphere infiltrated the cold baseof the stretched lithosphere and solidified as veins, therebycrystallizing calcite and phlogopite that were not in equilibriumwith peridotite. Continued Late Neoproterozoic lithosphericthinning, with progressive upwelling of the asthenosphere beneatha developing rift branch in this part of the North Atlanticcraton, caused further veining and successive remelting of veinsplus volatile-fluxed melting of the host fertile garnet peridotite,giving rise to long-lasting hybrid ultramafic lamprophyre magmaproduction in conjunction with the break-up of the Rodinia supercontinent.Proto-aillikite magma reached the surface only after coatingthe uppermost mantle conduits with glimmeritic material, whichcaused minor alkali loss. At intrusion level, carbonate separationfrom this aillikite magma resulted in fractionated dolomite-bearingcarbonatites (¹³C_PDB –3·7 to –2·7)and carbonate-poor mela-aillikite residues. Damtjernites maybe explained by liquid exsolution from alkali-rich proto-aillikitemagma batches that moved through previously reaction-lined conduitsat uppermost mantle depths. KEY WORDS: liquid immiscibility; mantle-derived magmas; metasomatism, Sr–Nd isotopes; U–Pb geochronology     

The Gronnedal-Ika Carbonatite-Syenite Complex, South Greenland: Carbonatite Formation by Liquid Immiscibility

The Grønnedal-Ika complex is dominated by layered nephelinesyenites which were intruded by a xenolithic syenite and a centralplug of calcite to calcite–siderite carbonatite. Aegirine–augite,alkali feldspar and nepheline are the major mineral phases inthe syenites, along with rare calcite. Temperatures of 680–910°Cand silica activities of 0·28–0·43 weredetermined for the crystallization of the syenites on the basisof mineral equilibria. Oxygen fugacities, estimated using titanomagnetitecompositions, were between 2 and 5 log units above the fayalite–magnetite–quartzbuffer during the magmatic stage. Chondrite-normalized REE patternsof magmatic calcite in both carbonatites and syenites are characterizedby REE enrichment (La_CN–Yb_CN = 10–70). Calcite fromthe carbonatites has higher Ba (5490 ppm) and lower HREE concentrationsthan calcite from the syenites (54–106 ppm Ba). This isconsistent with the behavior of these elements during separationof immiscible silicate–carbonate liquid pairs. _Nd(T =1·30 Ga) values of clinopyroxenes from the syenites varybetween +1·8 and +2·8, and _Nd(T) values of whole-rockcarbonatites range from +2·4 to +2·8. Calcitefrom the carbonatites has ¹⁸O values of 7·8 to 8·6and ¹³C values of –3·9 to –4·6. ¹⁸Ovalues of clinopyroxene separates from the nepheline syenitesrange between 4·2 and 4·9. The average oxygenisotopic composition of the nepheline syenitic melt was calculatedbased on known rock–water and mineral–water isotopefractionation to be 5·7 ± 0·4. Nd and C–Oisotope compositions are typical for mantle-derived rocks anddo not indicate significant crustal assimilation for eithersyenite or carbonatite magmas. The difference in ¹⁸O betweencalculated syenitic melts and carbonatites, and the overlapin _Nd values between carbonatites and syenites, are consistentwith derivation of the carbonatites from the syenites via liquidimmiscibility. KEY WORDS: alkaline magmatism; carbonatite; Gardar Province; liquid immiscibility; nepheline syenite     

A Melt Inclusion Record of Volatiles, Trace Elements and Li-B Isotope Variations in a Single Magma System from the Plat Pays Volcanic Complex, Dominica, Lesser Antilles

Glass inclusions in plagioclase and orthopyroxene from daciticpumice of the Cabrits Dome, Plat Pays Volcanic Complex in southernDominica reveal a complexity of element behavior and Li–Bisotope variations in a single volcanic center that would gounnoticed in a whole-rock study. Inclusions and matrix glassesare high-silica rhyolite with compositions consistent with about50% fractional crystallization of the observed phenocrysts.Estimated crystallization conditions are 760–880°C,200 MPa and oxygen fugacity of FMQ + 1 to +2 log units (whereFMQ is the fayalite–magnetite–quartz buffer). Manyinclusion glasses are volatile-rich (up to 6 wt % H₂O and 2900ppm Cl), but contents range down to 1 wt % H₂O and 2000 ppmCl as a result of shallow-level degassing. Sulfur contents arelow throughout, with <350 ppm S. The trace element compositionof inclusion glasses shows enrichment in light rare earth elements(LREE; (La/Sm)_n = 2·5–6·6) and elevatedBa, Th and K contents compared with whole rocks and similaror lower Nb and heavy REE (HREE; (Gd/Yb)_n = 0·5–1·0).Lithium and boron concentrations and isotope ratios in meltinclusions are highly variable (20–60 ppm Li with ⁷Li= +4 to +15 ± 2; 60–100 ppm B with ¹¹B = +6 to+13 ± 2) and imply trapping of isotopically heterogeneous,hybrid melts. Multiple sources and processes are required toexplain these features. The mid-ocean ridge basalt (MORB)-likeHREE, Nb and Y signature reflects the parental magma(s) derivedfrom the mantle wedge. Positive Ba/Nb, B/Nb and Th/Nb correlationsin inclusion glasses indicate coupled enrichment in stronglyfluid-mobile (Ba, B) and less-mobile (Th, Nb) trace elements,which can be explained by fractional crystallization of plagioclase,orthopyroxene and Fe–Ti oxides. The ⁷Li and ¹¹B valuesare at the high end of known ranges for other island arc magmas.We attribute the high values to a ¹¹B and ⁷Li-enriched slabcomponent derived from sea-floor-altered oceanic crust and possiblyfurther enriched in heavy isotopes by dehydration fractionation.The heterogeneity of isotope ratios in the evolved, trappedmelts is attributed to shallow-level assimilation of older volcanicrocks of the Plat Pays Volcanic Complex. KEY WORDS: subduction; volcanic arcs; igneous processes; melt inclusions; SIMS; trace elements; lithium and boron isotopes; diffusion     

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     

Liquidus Equilibria in the System K2O-Na2O-Al2O3-SiO2-F2O-1-H2O to 100 MPa: I. Silicate-Fluoride Liquid Immiscibility in Anhydrous Systems

Liquidus relations in the four-component system Na₂O–Al₂O₃–SiO₂–F₂O_–1were studied at 0· 1 and 100 MPa to define the locationof fluoride–silicate liquid immiscibility and outlinedifferentiation paths of fluorine-bearing silicic magmas. Thefluoride–silicate liquid immiscibility spans the silica–albite–cryoliteand silica–topaz–cryolite ternaries and the haplogranite-cryolitebinary at greater than 960°C and 0· 1–100 MPa.With increasing Al₂O₃ in the system and increasing aluminum/alkalication ratio, the two-liquid gap contracts and migrates fromthe silica liquidus to the cryolite liquidus. The gap does notextend to subaluminous and peraluminous melt compositions. Forall alkali feldspar–quartz-bearing systems, the miscibilitygap remains located on the cryolite liquidus and is thus inaccessibleto differentiating granitic and rhyolitic melts. In peralkalinesystems, the magmatic differentiation is terminated at the albite–quartz–cryoliteeutectic at 770°C, 100 MPa, 5 wt % F and cation Al/Na =0· 75. The addition of topaz, however, significantlylowers melting temperatures and allows strong fluorine enrichmentin subaluminous compositions. At 100 MPa, the binary topaz–cryoliteeutectic is located at 770°C, 39 wt % F, cation Al/Na 0·95, and the ternary quartz–topaz–cryolite eutecticis found at 740°C, 32 wt % F, 30 wt % SiO₂ and cation Al/Na 0· 95. Such location of both eutectics enables fractionationpaths of subaluminous quartz-saturated systems to produce fluorine-rich,SiO₂-depleted and nepheline-normative residual liquids. KEY WORDS: silicate melt; granite; rhyolite; fluorine; liquid immiscibility     

Volcan Popocatepetl, Mexico. Petrology, Magma Mixing, and Immediate Sources of Volatiles for the 1994-Present Eruption

Volcán Popocatépetl has been the site of voluminousdegassing accompanied by minor eruptive activity from late 1994until the time of writing (August 2002). This contribution presentspetrological investigations of magma erupted in 1997 and 1998,including major-element and volatile (S, Cl, F, and H₂O) datafrom glass inclusions and matrix glasses. Magma erupted fromPopocatépetl is a mixture of dacite (65 wt % SiO₂, two-pyroxenes+ plagioclase + Fe–Ti oxides + apatite, 3 wt % H₂O, P= 1·5 kbar, f_O2 = NNO + 0·5 log units) and basalticandesite (53 wt % SiO₂, olivine + two-pyroxenes, 3 wt % H₂O,P = 1–4 kbar). Magma mixed at 4–6 km depth in proportionsbetween 45:55 and 85:15 wt % silicic:mafic magma. The pre-eruptivevolatile content of the basaltic andesite is 1980 ppm S, 1060ppm Cl, 950 ppm F, and 3·3 wt % H₂O. The pre-eruptivevolatile content of the dacite is 130 ± 50 ppm S, 880± 70 ppm Cl, 570 ± 100 ppm F, and 2·9 ±0·2 wt % H₂O. Degassing from 0·031 km³ of eruptedmagma accounts for only 0·7 wt % of the observed SO₂emission. Circulation of magma in the volcanic conduit in thepresence of a modest bubble phase is a possible mechanism toexplain the high rates of degassing and limited magma productionat Popocatépetl. KEY WORDS: glass inclusions; igneous petrology; Mexico; Popocatépetl; volatiles     

Oxygen Isotope Evidence for Chemical Interaction of Ki lauea Historical Magmas with Basement Rocks

Klauea historical summit lavas have a wide range in matrix ¹⁸O_VSMOWvalues (4·9–5·6) with lower values in rockserupted following a major summit collapse or eruptive hiatus.In contrast, ¹⁸O values for olivines in most of these lavasare nearly constant (5·1 ± 0·1). The disequilibriumbetween matrix and olivine ¹⁸O values in many samples indicatesthat the lower matrix values were acquired by the magma afterolivine growth, probably just before or during eruption. BothMauna Loa and Klauea basement rocks are the likely sources ofthe contamination, based on O, Pb and Sr isotope data. However,the extent of crustal contamination of Klauea historical magmasis probably minor (< 12%, depending on the assumed contaminant)and it is superimposed on a longer-term, cyclic geochemicalvariation that reflects source heterogeneity. Klauea's heterogeneoussource, which is well represented by the historical summit lavas,probably has magma ¹⁸O values within the normal mid-ocean ridgebasalt mantle range (5·4–5·8) based on thenew olivine ¹⁸O values. KEY WORDS: Hawaii; Klauea; basalt; oxygen isotopes; crustal contamination     

Petrographic, Chemical and B-Isotopic Insights into the Origin of Tourmaline-Rich Rocks and Boron Recycling in the Martinamor Antiform (Central Iberian Zone, Salamanca, Spain)

Tourmaline in the Martinamor antiform occurs in tourmalinites(rocks with >15–20% tourmaline by volume), clasticmetasedimentary rocks of the Upper Proterozoic Monterrubio formation,quartz veins, pre-Variscan orthogneisses and Variscan graniticrocks. Petrographic observations, back-scattered electron (BSE)images, and microprobe data document a multistaged developmentof tourmaline. Overall, variations in the Mg/(Mg + Fe) ratiosdecrease from tourmalinites (0·36–0·75),through veins (0·38–0·66) to granitic rocks(0·23–0·46), whereas Al increases in thesame order from 5·84–6·65 to 6·22–6·88apfu. The incorporation of Al into tourmaline is consistentwith combinations of ^xAl(NaR)_–1 and AlO(R(OH))_–1exchange vectors, where ^x represents X-site vacancy and R is(Mg + Fe²⁺ + Mn). Variations in ^x/(^x + Na) ratios are similarin all the types of tourmaline occurrences, from 0·10to 0·53, with low Ca-contents (mostly <0·10apfu). Based on field and textural criteria, two groups of tourmaline-richrocks are distinguished: (1) pre-Variscan tourmalinites (probablyCadomian), affected by both deformation and regional metamorphismduring the Variscan orogeny; (2) tourmalinites related to thesynkinematic granitic complex of Martinamor. Textural and geochemicaldata are consistent with a psammopelitic parentage for the protolithof the tourmalinites. Boron isotope analyses of tourmaline havea total range of ¹¹B values from –15·6 to 6·8;the lowest corresponding to granitic tourmalines (–15·6to –11·7) and the highest to veins (1·9to 6·8). Tourmalines from tourmalinites have intermediate¹¹B values of –8·0 to +2·0. The observedvariations in ¹¹B support an important crustal recycling ofboron in the Martinamor area, in which pre-Variscan tourmaliniteswere remobilized by a combination of mechanical and chemicalprocesses during Variscan deformation, metamorphism and anatexis,leading to the formation of multiple tourmaline-bearing veinsand a new stage of boron metasomatism. KEY WORDS: tourmalinites; metamorphic and granitic rocks; mineral chemistry; whole-rock chemistry; boron isotopes     

Magmatic Differentiation at an Island-arc Caldera: Okmok Volcano, Aleutian Islands, Alaska

Okmok volcano is situated on oceanic crust in the central Aleutianarc and experienced large (15 km³) caldera-forming eruptionsat 12 000 years BP and 2050 years BP. Each caldera-forming eruptionbegan with a small Plinian rhyodacite event followed by theemplacement of a dominantly andesitic ash-flow unit, whereaseffusive inter- and post-caldera lavas have been more basaltic.Phenocryst assemblages are composed of olivine + pyroxene +plagioclase ± Fe–Ti oxides and indicate crystallizationat 1000–1100°C at 0·1–0·2 GPain the presence of 0–4% H₂O. The erupted products followa tholeiitic evolutionary trend and calculated liquid compositionsrange from 52 to 68 wt % SiO₂ with 0·8–3·3wt % K₂O. Major and trace element models suggest that the moreevolved magmas were produced by 50–60% in situ fractionalcrystallization around the margins of the shallow magma chamber.Oxygen and strontium isotope data (¹⁸O 4·4–4·9,⁸⁷Sr/ ⁸⁶Sr 0·7032–0·7034) indicate interactionwith a hydrothermally altered crustal component, which led toelevated thorium isotope ratios in some caldera-forming magmas.This compromises the use of uranium–thorium disequilibria[(²³⁰Th/ ²³⁸U) = 0·849–0·964] to constrainthe time scales of magma differentiation but instead suggeststhat the age of the hydrothermal system is 100 ka. Modellingof the diffusion of strontium in plagioclase indicates thatmany evolved crystal rims formed less than 200 years prior toeruption. This addition of rim material probably reflects theremobilization of crystals from the chamber margins followingreplenishment. Basaltic recharge led to the expansion of themagma chamber, which was responsible for the most recent caldera-formingevent. KEY WORDS: Okmok; caldera; U-series isotopes; Sr-diffusion; time scales; Aleutian arc     

Open-System Magma Chamber Evolution: an Energy-constrained Geochemical Model Incorporating the Effects of Concurrent Eruption, Recharge, Variable Assimilation and Fractional Crystallization (EC-E'RA{chi}FC)

Significant petrogenetic processes governing the geochemicalevolution of magma bodies include magma Recharge (includingformation of ‘quenched inclusions’ or enclaves),heating and concomitant partial melting of country rock withpossible ‘contamination’ of the evolving magma body(Assimilation), and formation and separation of cumulates byFractional Crystallization (RAFC). Although the importance ofmodeling such open-system magma chambers subject to energy conservationhas been demonstrated, the effects of concurrent removal ofmagma by eruption and/or variable assimilation (involving imperfectextraction of anatectic melt from wall rock) have not been considered.In this study, we extend the EC-RAFC model to include the effectsof Eruption and variable amounts of assimilation, A. This model,called EC-E'RAFC, tracks the compositions (trace elements andisotopes), temperatures, and masses of magma body liquid (melt),eruptive magma, cumulates and enclaves within a composite magmaticsystem undergoing simultaneous eruption, recharge, assimilationand fractional crystallization. The model is formulated as aset of 4 + t + i + s coupled nonlinear differential equations,where the number of trace elements, radiogenic and stable isotoperatios modeled are t, i and s, respectively. Solution of theEC-E'RAFC equations provides values for the average temperatureof wall rock (T_a), mass of melt within the magma body (M_m),masses of cumulates (M_ct), enclaves (M_en) and wall rock () and the masses of anatectic melt generated () and assimilated (). In addition, t trace element concentrations and i + s isotopic ratios inmelt and eruptive magma (C_m, _m, _m), cumulates (C_ct, _m, _m), enclaves(C_en, , ) and anatectic melt (C_a, , ) as a function of magma temperature (T_m) are also computed. Input parametersinclude the (user-defined) equilibration temperature (T_eq),a factor describing the efficiency of addition of anatecticmelt () from country rock to host magma, the initial temperatureand composition of pristine host melt (, , , ), recharge melt (, , , ) and wall rock (, , , ), distribution coefficients (D_m, D_r, D_a) and their temperaturedependences (H_m, H_r, H_a), latent heats of transition (meltingor crystallization) for wall rock (h_a), pristine magma (h_m)and recharge magma (h_r) as well as the isobaric specific heatcapacity of assimilant (C_p,a), pristine (C_p,m) and recharge(C_p,r) melts. The magma recharge mass and eruptive magma massfunctions, M_r(T_m) and M_e(T_m), respectively, are specified apriori. M_r(T_m) and M_e(T_m) are modeled as either continuous orepisodic (step-like) processes. Melt productivity functions,which prescribe the relationship between melt mass fractionand temperature, are defined for end-member bulk compositionscharacterizing the local geologic site. EC-E'RAFC has potentialfor addressing fundamental questions in igneous petrology suchas: What are intrusive to extrusive ratios (I/E) for particularmagmatic systems, and how does this factor relate to rates ofcrustal growth? How does I/E vary temporally at single, long-livedmagmatic centers? What system characteristics are most profoundlyinfluenced by eruption? What is the quantitative relationshipbetween recharge and assimilation? In cases where the extractionefficiency can be shown to be less than unity, what geologiccriteria are important and can these criteria be linked to fieldobservations? A critical aspect of the energy-constrained approachis that it requires integration of field, geochronological,petrologic, and geochemical data, and, thus, the EC-ERAFC ‘systems’approach provides a means for answering broad questions whileunifying observations from a number of disciplines relevantto the study of igneous rocks. KEY WORDS: assimilation; energy conservation; eruption; open system; recharge     

An Experimental Study of Fe-Al Solubility in the System Corundum-Hematite up to 40 kbar and 1300{degrees}C

The mutual solubility in the system corundum–hematite[-(Al, Fe³⁺)₂O₃] was investigated experimentally using bothsynthetic and natural materials. Mixtures of -Al₂O₃ and -Fe₂O₃(weight ratios of 8:2 and 10:1) were used as starting materialsfor synthesis experiments in air at 800–1300°C withrun times of 7–34 days. Experiments at 8–40 kbarand 490–1100°C were performed in a piston-cylinderapparatus (run times of 0·8–7·4 days) usinga natural diasporite consisting of 60–70 vol. % diasporeand 20–30 vol. % Ti-hematite. During the diasporite–corunditetransformation, the FeTiO₃ component (12–18 mol %) ofTi-hematite only slightly increased, implying that oxygen fugacitywas maintained at high values. Run products were studied byelectron microprobe and X-ray diffraction (Rietveld) techniques.An essentially linear volume of mixing exists in the solid solutionwith a slight positive deviation at the hematite side. Up to1000°C, corundum contains <4 mol % Fe₂O₃ and hematite<10 mol % Al₂O₃; at 1200°C these amounts increase to9·3 and 17·0 mol %, respectively. At 1300°Chematite was no longer stable and coexists with the orthorhombic phase . The present results agree with corundum (solvus) compositions obtained inprevious studies but indicate a larger solubility of Al in hematite.The miscibility gap in the solution can be modelled with anasymmetric Margules equation with interaction parameters (2uncertainties): ; ; ; . Application of the corundum–hematite solution as a solvus geothermometer is limited because of thescarcity of suitable rock compositions. KEY WORDS: corundum; hematite; corundum–hematite miscibility gap; experimental study; Margules model; metabauxite     

Source of the Quaternary Alkalic Basalts, Picrites and Basanites of the Potrillo Volcanic Field, New Mexico, USA: Lithosphere or Convecting Mantle?

The <80 ka basalts–basanites of the Potrillo VolcanicField (PVF) form scattered scoria cones, lava flows and maarsadjacent to the New Mexico–Mexico border. MgO ranges upto 12·5%; lavas with MgO < 10·7% have fractionatedboth olivine and clinopyroxene. Cumulate fragments are commonin the lavas, as are subhedral megacrysts of aluminous clinopyroxene(with pleonaste inclusions) and kaersutitic amphibole. REE modellingindicates that these megacrysts could be in equilibrium withthe PVF melts at 1·6–1·7 GPa pressure. Thelavas fall into two geochemical groups: the Main Series (85%of lavas) have major- and trace-element abundances and ratiosclosely resembling those of worldwide ocean-island alkali basaltsand basanites (OIB); the Low-K Series (15%) differ principallyby having relatively low K₂O and Rb contents. Otherwise, theyare chemically indistinguishable from the Main Series lavas.Sr- and Nd-isotopic ratios in the two series are identical andvary by scarcely more than analytical error, averaging ⁸⁷Sr/⁸⁶Sr= 0·70308 (SD = 0·00004) and ¹⁴³Nd/¹⁴⁴Nd = 0·512952(SD=0·000025). Such compositions would be expected ifboth series originated from the same mantle source, with Low-Kmelts generated when amphibole remained in the residuum. ThreePVF lavas have very low Os contents (<14 ppt) and appearto have become contaminated by crustal Os. One Main Series picritehas 209 ppt Os and has a _Os value of +13·6, typical forOIB. This contrasts with published ¹⁸⁷Os/¹⁸⁸Os ratios for KilbourneHole peridotite mantle xenoliths, which give mostly negative_Os values and show that Proterozoic lithospheric mantle formsa thick Mechanical Boundary Layer (MBL) that extends to 70 kmdepth beneath the PVF area. The calculated mean primary magma,in equilibrium with Fo₈₉, has Na₂O and FeO contents that givea lherzolite decompression melting trajectory from 2·8GPa (95 km depth) to 2·2 GPa (70 km depth). Inverse modellingof REE abundances in Main Series Mg-rich lavas is successfulfor a model invoking decompression melting of convecting sub-lithosphericlherzolite mantle (Nd = 6·4; T_p 1400°C) between90 and 70 km. Nevertheless, such a one-stage model cannot accountfor the genesis of the Low-K Series because amphibole wouldnot be stable within convecting mantle at T_f 1400°C. Thesemagmas can only be accommodated by a three-stage model thatenvisages a Thermal Boundary Layer (TBL) freezing conductivelyonto the 70 km base of the Proterozoic MBL during the 20 Myrtectonomagmatic quiescence before PVF eruptions. As it grew,this was veined by hydrous small-fraction melts from below.The geologically recent arrival of hotter-than-ambient (T_p 1400°C) convecting mantle beneath the Potrillo area re-meltedthe TBL and caused the magmatism. KEY WORDS: western USA; picrites; Sr–Nd–Os isotopes; petrogenetic modelling; thermal boundary layer     

Basalt and Sediment Geochemistry and Magma Petrogenesis in a Transect from Oceanic Island Arc to Rifted Continental Margin Arc: the Kermadec--Hikurangi Margin, SW Pacific

Sediment mixing and recycling through a subduction zone canbe detected in lead isotopes and trace elements from basaltsand sediments from the Kermadec-Hikurangi Margin volcanic arcsystem and their coupled back-arc basins. Sr, Nd and Pb isotopesfrom the basalts delineate relatively simple, almost overlapping,arrays between back-arc basin basalts of the Havre Trough-NgatoroBasin (⁸⁷Sr/⁸⁶Sr = 0.70255; _Nd=+9.3; ²⁰⁶Pb/²⁰⁴Pb = 18.52; ²⁰⁸Pb/²⁰⁴Pb= 38.18), island arc basalts from the Kermadec Arc togetherwith basalts from Taupo Volcanic Zone (⁸⁷Sr/⁸⁶Sr 0.7042; _Nd= +5; ²⁰⁶Pb/²⁰⁴Pb= 18.81; ²⁰⁸Pb/²⁰⁴Pb = 38.61), and sedimentsderived from New Zealand's Mesozoic (Torlesse) basement (⁸⁷Sr/⁸⁶Sr 0.715; _Nd —4; ²⁰⁶Pb/²⁰⁴Pb 18.86; ²⁰⁸Pb/²⁰⁴Pb 38.8).Basalts from the arc front volcanoes have high Cs, Rb, Ba, Th,U and K, and generally high but variable Ba/La, Ba/Nb ratios,characteristic of subduction-related magmas, relative to typicaloceanic basalts. These signatures are diluted in the back-arcbasins, which are more like mid-ocean ridge basalts. Strongchemical correlations in plots of SiO₂ vs CaO and loss on ignitionfor the sediments (finegrained muds) are consistent with mixingbetween detrital and biogenic (carbonate-rich) components. Otherdata, such as Zr vs CaO, are consistent with the detrital componentcomprising a mixture of arc- and continent-derived fractions.In chondrite-normalized diagrams, most of the sediments havelight rare earth element enriched patterns, and all have negativeEu anomalies. The multielement diagrams have negative spikesat Nb, P and Ti and distinctive enrichments in the large ionlithophile elements and Pb relative to mantle. Isotopic measurementsof Pb, Sr and Nd reveal restricted fields of Pb isotopes butwide variation in Nd and Sr relative to other sediments fromthe Pacific Basin. Rare K-rich basalts from Clark Volcano towardthe southern end of the oceanic Kermadec Island Arc show unusualand primitive characteristics ( 2% K₂O at 50% SiO₂, Ba 600p.p.m., 9–10% MgO and Ni > 100 p.p.m.) but have highlyradiogenic Sr, Nd and Pb isotopes, similar to those of basaltsfrom the continental Taupo Volcanic Zone. These oceanic islandarc basalts cannot have inherited their isotope signatures throughcrustal contamination or assimilation—fractional crystallizationtype processes, and this leads us to conclude that source processesvia bulk sediment mixing, fluid and/or melt transfer or somecombination of these are responsible. Although our results showclear chemical gradients from oceanic island arc to continentalmargin arc settings (Kermadec Arc to Taupo Volcanic Zone), overlapbetween the data from the oceanic and continental sectors suggeststhat the lithospheric (crustal contamination) effect may beminimal relative to that of sediment subduction. Indeed, itis possible to account for the chemical changes by a decreasenorthward in the sediment flux into the zone of magma genesis.This model receives support from recent sediment dispersal studiesin the Southern Ocean which indicate that a strong bottom current(Deep Western Boundary Current) flows northward along the easterncontinental margin of New Zealand and sweeps continental derivedsediment into the sediment-starved oceanic trench system. Thetrace element and isotopic signatures of the continental derivedcomponent of this sediment are readily distinguished, but alsodiluted in a south to north direction along the plate boundary. KEY WORDS: subduction zone basalts; sediments; Sr-, Nd-, Pb-isotopes; trace elements ^*Present address: School of Earth Sciences, University of Melbourne, Parkville, Vic. 3052, Australia.     

Fluid Flow During Contact Metamorphism at Mary Kathleen, Queensland, Australia

Corella marbles in the Mary Kathleen Fold Belt were infiltratedby fluids during low-pressure (200-MPa) contact metamorphismassociated with the intrusion of the Burstall granite at 1730–1740Ma. Fluids emanating from the granite [whole-rock (WR) ¹⁸O=8.1–8.6%]produced Fe-rich massive and banded garnet—clinopyroxeneskarns [¹⁸O(WR)=9.1–11.9%]. Outside the skarn zones, marblemineralogies define an increase in temperature (500 to >575C) and XCO₂ (0.05 to >0.12) towards the granite, andmost marbles contain isobarically univariant or invariant assemblagesin the end-member CaO–MgO–Al₂O₃–SiO₂–H₂O–CO₂system. Marbles have calcite (Cc) ¹⁸O and ¹³C values of 12.3–24.6%and –1.0 to –3.9%, respectively. A lack of down-temperaturemineral reactions in the marbles suggests that pervasive fluidinfiltration did not continue after the thermal peak of contactmetamorphism. The timing of fluid flow probably correspondsto a period of high fluid production and high intrinsic permeabilitiesduring prograde contact metamorphism. The petrology and stableisotope geochemistry of the marbles suggest that these rockswere infiltrated by water-rich fluids. If fluid flow occurredup to the peak of contact metamorphism, the mineralogical andisotopic resetting is best explained by fluids flowing up-temperaturetoward the Burstall granite. However, if fluid flow ceased beforthe peak of regional metamorphism, the fluid flow directioncannot be unambiguously determined. At individual outcrops,marble ¹⁸O(Cc) values vary by several permil over a few squaremetres, suggesting that fluid fluxes varied by at least an orderof magnitude on the metre to tens-of-metre scale. Fluids werefocused across lithological layering; however, mesoscopic fracturesare not recognized. The focusing of fluids was possibly viamicrofractures, and the variation in the degree of resettingmay reflect variations in microcrack density and fracture permeability.The marble—skarn contacts represent a sharp discontinuityin both major element geochemistry and ¹⁸O values, suggestingthat, at least locally, little fluid flow occurred across thesecontacts.     

Geochemistry of South African On- and Off-craton, Group I and Group II Kimberlites: Petrogenesis and Source Region Evolution

Bulk-rock geochemical compositions of hypabyssal kimberlites,emplaced through the Archaean Kaapvaal craton and ProterozoicNamaqua–Natal belt, are used to estimate close-to-primarymagma compositions of Group I kimberlites (Mg-number = 0·82–0·87;22–28 wt % MgO; 21–30 wt % SiO₂; 10–17 wt% CaO; 0·2–1·7 wt % K₂O) and Group II kimberlites(Mg-number = 0·86–0·89; 23–29 wt %MgO; 28–36 wt % SiO₂; 8–13 wt % CaO; 1·6–4·6wt % K₂O). Group I kimberlites are distinguished from GroupII by their lower Ba/Nb (<12), Th/Nb (<1·1) andLa/Nb (<1·1) but higher Ce/Pb (>22) ratios. Thedistinct rare earth element patterns of the two types of kimberlitesindicate a more highly metasomatized source for Group II kimberlites,with more residual clinopyroxene and less residual garnet. Thesimilarity of Sr and Nd isotope ratios and diagnostic traceelement ratios (Ce/Pb, Nb/U, La/Nb, Ba/Nb, Th/Nb) of Group Ikimberlites to ocean island basalts (OIB), but more refractoryMg-numbers and Ni contents, are consistent with derivation ofGroup I kimberlites from subcontinental lithospheric mantle(SCLM) that has been enriched by OIB-like melts or fluids. Sourceenrichment ages and plate reconstructions support a direct associationof these melts or fluids with Mesozoic upwelling beneath southernAfrica of a mantle plume(s), at present located beneath thesouthern South Atlantic Ocean. In contrast, the geochemicalcharacteristics of both on- and off-craton Group II kimberlitesshow strong similarity to calc-alkaline magmas, particularlyin their Nb and Ta depletion and Pb enrichment. It is suggestedthat Group II kimberlites are derived from both Archaean andProterozoic lithospheric mantle source regions metasomatizedby melts or fluids associated with ancient subduction events,unrelated to mantle plume upwelling. The upwelling of mantleplumes beneath southern Africa during the Mesozoic, at the timeof Gondwana break-up, may have acted as a heat source for partialmelting of the SCLM and the generation of both Group I and GroupII kimberlite magmas. KEY WORDS: kimberlite; geochemistry; petrogenesis; mantle plumes; South Africa