The effect of assimilation, fractional crystallization, and ageing on U-series disequilibria in subduction zone lavas

Although most arc lavas have experienced significant magma differentiation, the effect of the differentiation process on U-series disequilibria is still poorly understood. Here we present a numerical model for simulating the effect of time-dependent magma differentiation processes on U-series disequilibria in lavas from convergent margins. Our model shows that, in a closed system with fractional crystallization, the ageing effect can decrease U-series disequilibria via radioactive decay while in an open system, both ageing and bulk assimilation of old crustal material serve to reduce the primary U-series disequilibria. In contrast, with recharge of refresh magma, significant ²²⁶Ra excess in erupted lavas can be maintained even if the average residence time is longer than 8000 years.The positive correlations of (²²⁶Ra/²³⁰Th) between Sr/Th or Ba/Th in young lavas from convergent margins have been widely used as evidence of fluid addition generating the observed ²²⁶Ra excess in subduction zones. We assess to what extent the positive correlations of (²²⁶Ra/²³⁰Th) with Sr/Th and Ba/Th observed in the Tonga arc could reflect AFC process. Results of our model show that these positive correlations can be produced during time-dependent magma differentiation at shallow crustal levels. Specifically, fractional crystallization of plagioclase and amphibole coupled with contemporaneous decay of ²²⁶Ra can produce positive correlations between (²²⁶Ra/²³⁰Th) and Sr/Th or Ba/Th (to a lesser extent). Therefore, the correlations of (²²⁶Ra/²³⁰Th) with Sr/Th and Ba/Th cannot be used to unambiguously support the fluid addition model, and the strength of previous conclusions regarding recent fluid addition and ultra-fast ascent rates of arc magmas is significantly lessened.     

Arc dacite genesis pathways: Evidence from mafic enclaves and their hosts in Aegean lavas

Mafic enclaves are commonly found in intermediate arc magmas, and their occurrence has been linked to eruption triggering by pre-eruptive magma mixing processes. New major, trace, Sr–Nd and U–Th isotope data of rocks from Nisyros in the Aegean volcanic arc are presented here. Pre-caldera samples display major and trace element trends that are consistent with fractionation of magnetite and apatite within intermediate compositions, and zircon within felsic compositions, and preclude extensive hybridization between mafic and felsic magmas. In contrast, post-caldera dacites form a mixing trend towards their mafic enclaves. In terms of U-series isotopes, most samples show small ²³⁸U excesses of up to  10%. Mafic enclaves have significantly higher U/Th ratios than their dacitic host lavas, precluding simple models that relate the mafic and felsic magmas by fractionation or aging alone. A more complicated petrogenetic scenario is required. The post-caldera dacites are interpreted to represent material remobilized from a young igneous protolith following influx of fresh mafic magma, consistent with the U–Th data and with Sr–Nd isotope constraints that point to very limited (< 10%) assimilation of old crust at Nisyros. When these results are compared to data from Santorini in the same arc, there are many geochemical similarities between the two volcanic centers during the petrogenesis of the pre-caldera samples. However, striking differences are apparent for the post-caldera lavas: in Nisyros, dacites show geochemical and textural evidence for magma mixing and remobilization by influx of mafic melts, and they erupt as viscous lava domes; in Santorini, evidence for geochemical hybridization of dacites and mafic enclaves is weak, dacite petrogenesis does not involve protolith remobilization, and lavas erupt as less viscous flows. Despite these differences, it appears that mafic enclaves in intermediate Aegean arc magmas consistently yield timescales of at least 100 kyrs between U enrichment of the mantle wedge and eruption, on the upper end of those estimated for the eruptive products of mafic arc volcanoes. Finally, the data presented here provide constraints on the rates of differentiation from primitive arc basalts to dacites (less than  140 kyrs), and on the crustal residence time of evolved igneous protoliths prior to their remobilization by mafic arc magmas (greater than  350 kyrs).     

Lead isotopes by LA-MC-ICPMS: Tracking the emergence of mantle signatures in an evolving silicic magma system

At Long Valley (LV) model Sr isotope phenocryst ages and absolute U-Pb zircon ages from precaldera Glass Mountain (GM) and caldera-related Bishop Tuff (BT) rhyolites show that these crystals track >1 Myr of evolution of a voluminous rhyolite magmatic system. In detail, strong disparities between the different age populations complicate ideas for a unified model for rhyolite generation, differentiation, and storage. To better elucidate the age discrepancies a new in situ Pb isotope technique has been developed to measure the compositions of 113 individual LV feldspars (mainly sanidine) and their host glasses by UV laser ablation MC-ICPMS. Given sufficient signal the accuracy and precision of this technique approaches that of double-spike thermal ionization mass spectrometry. The utility of our technique for many geologic materials is, however, limited to determining Pb isotope ratios that include ²⁰⁶Pb, ²⁰⁷Pb, and ²⁰⁸Pb, but exclude ²⁰⁴Pb. New zircon ²³⁸U-²⁰⁶Pb crystallization ages were also obtained for two older Glass Mountain domes.A >1.5‰ difference between the Pb isotope compositions of feldspars from older (1.7-2.2 Ma) precaldera Glass Mountain (GM) rhyolites and younger LV rhyolites, including the BT, is found. The Pb isotope data for feldspars and their host glasses lie along a regional trend line between young basalts and evolved crust compositions, spanning ∼15% of that isotopic difference, and show a secular change towards increasing mantle contribution. Most feldspars have Pb isotope compositions that are similar to their host glasses and, as such, there persists an apparent >100 k.y. difference between Sr model feldspar ages and zircon ages for some GM rhyolites. Collectively, the feldspars define a Sr-Pb isotope mixing curve. Evidence for mixing complicates the interpretation that the Sr isotope data solely reflect radiogenic ingrowth. Where isotopically heterogeneous feldspar populations occur, there is greater uncertainty about the veracity of the Sr model ages. Specifically, we find no Pb isotope evidence that BT feldspars grew from older GM-like magmas.The distinct Pb isotope signatures for individual rhyolites and their feldspars support evidence based on zircon dating that LV volcanism did not erupt from a single long-lived magma chamber but rather tapped a number of different magmas. Moreover, contrary to the conventional model of gradual build-up prior to cataclysmic eruption, secular changes in the U-Pb age constraints on magma residence times and the magmas’ distinct Pb isotopic compositions suggest that, at Long Valley, eruptive volumes increase with shorter magma residence time and correlate with greater mantle input. Evidently, the plumbing and therefore activity at Long Valley was influenced by the evolving interaction between source and crustal magma system.     

Open System Magmatic Evolution of the Taos Plateau Volcanic Field, Northern New Mexico: II. The Genesis of Cryptic Hybrids

Tholeiitic lavas of the Servilleta Basalt exhibit only subtletextural and mineralogical evidence for a hybrid origin, butelemental and isotopic analyses of these basalts are best modeledin terms of mixing Servilleta parent magma with a range of contemporaneousandesite and dacite magmas. Cryptic compositional heterogeneitiesin some flows interpreted as hybrids apparently reflect incompletehomogenization following pre-emptive magma mingling. The generalscarcity of mixing-related textural disequilibrium is ascribedin part to mixing of mineralogically similar end-members. Eradicationof some phenocrysts during post-mixing residence and evolutionin a convecting magma body may be an even more important factor. Eruptions of hybrid magmas may frequently be triggered by magmamixing events (i.e. injection or replenishment), and minglingof compositionally diverse magmas may ensue as a consequenceof tapping a compositionally graded or layered magma chamber.These hybrids are instantly recognizable by the preservationof disequilibrium textures and mineral assemblages, and by discontinuouscompositional heterogeneities. Cryptic hybrids, which have notpreserved this record, will be recognizable as mixed magmasprimarily by geochemical evidence for open system evolution.     

Chemical versus Temporal Controls on the Evolution of Tholeiitic and Calc-alkaline Magmas at Two Volcanoes in the Alaska-Aleutian Arc

The Alaska–Aleutian island arc is well known for eruptingboth tholeiitic and calc-alkaline magmas. To investigate therelative roles of chemical and temporal controls in generatingthese contrasting liquid lines of descent we have undertakena detailed study of tholeiitic lavas from Akutan volcano inthe oceanic Aleutian arc and calc-alkaline products from Aniakchakvolcano on the continental Alaskan Peninsula. The differencesdo not appear to be linked to parental magma composition. TheAkutan lavas can be explained by closed-system magmatic evolution,whereas curvilinear trace element trends and a large range in⁸⁷Sr/⁸⁶Sr isotope ratios in the Aniakchak data appear to requirethe combined effects of fractional crystallization, assimilationand magma mixing. Both magmatic suites preserve a similar rangein ²²⁶Ra–²³⁰Th disequilibria, which suggests that thetime scale of crustal residence of magmas beneath both thesevolcanoes was similar, and of the order of several thousandyears. This is consistent with numerical estimates of the timescales for crystallization caused by cooling in convecting crustalmagma chambers. During that time interval the tholeiitic Akutanmagmas underwent restricted, closed-system, compositional evolution.In contrast, the calc-alkaline magmas beneath Aniakchak volcanounderwent significant open-system compositional evolution. Combiningthese results with data from other studies we suggest that differentiationis faster in calc-alkaline and potassic magma series than intholeiitic series, owing to a combination of greater extentsof assimilation, magma mixing and cooling. KEY WORDS: uranium-series; Aleutian arc; magma differentiation; time scales     

Crystallization history of rhyolites at Long Valley, California, inferred from combined U-series and Rb-Sr isotope systematics

In this study, we present ⁸⁷Rb/⁸⁶Sr and ²³⁰Th/²³⁸U isotope analyses of glasses and phenocrysts from postcaldera rhyolites erupted between 150 to 100 ka from the Long Valley magmatic system. Both isotope systems indicate complex magma evolution with preeruptive mineral crystallization and magma fractionation, followed by extended storage in a silicic magma reservoir. Glass analyses yield a Rb-Sr isochron of 257 ± 39 ka, which can be explained by a feldspar-fractionation event ∼150 ky before eruption. Individual feldspar-glass pairs confirm this age result. A mineral ²³⁰Th-²³⁸U isochron in a low-silica rhyolite from the Deer Mountain Dome defines an age of 236 ± 1 ka, but the glass and whole rock do not lie on the isochron. U-Th fractionation of the rocks is controlled by the accessory minerals zircon and probably allanite, which crystallized at 250 ± 3 ka and 187 ± 9 ka, respectively. All major mineral phases contain accessory mineral phases; therefore, the mineral isochron represents a mixture of zircon and allanite populations. A precision of ±1 ka for the mixing array implies that the minor phases must have crystallized within this timescale. Longer periods of crystal growth would cause the mixing array to be less well defined. U-series data from other low- and high-silica rhyolites indicate younger accessory mineral crystallization events at ∼200 and 140 ka, probably related to the thermal evolution of the magma reservoir. These crystallization events are, however, only documented by the accessory minerals and had no further influence on bulk magma compositions. We interpret the indistinguishable age results from both isotope systems (∼250 ka) to record the fractionation of small magma batches by filter pressing from a much larger underlying magma volume, followed by physical isolation and extended storage at the top of the magma reservoir for up to 150 ky.     

Plume-lithosphere interaction beneath Mt. Cameroon volcano, West Africa: Constraints from U-Th-Ra and Sr-Nd-Pb isotope systematics

Precise measurements of ²³⁸U-²³⁰Th-²²⁶Ra disequilibria in lavas erupted within the last 100 yr on Mt. Cameroon are presented, together with major and trace elements, and Sr-Nd-Pb isotope ratios, to unravel the source and processes of basaltic magmatism at intraplate tectonic settings. All samples possess ²³⁸U-²³⁰Th-²²⁶Ra disequilibria with ²³⁰Th (18-24%) and ²²⁶Ra (9-21%) excesses, and there exists a positive correlation in a (²²⁶Ra/²³⁰Th)-(²³⁰Th/²³⁸U) diagram. The extent of ²³⁸U-²³⁰Th-²²⁶Ra disequilibria is markedly different in lavas of individual eruption ages, although the (²³⁰Th/²³²Th) ratio is constant irrespective of eruption age. When U-series results are combined with Pb isotope ratios, negative correlations are observed in the (²³⁰Th/²³⁸U)-(²⁰⁶Pb/²⁰⁴Pb) and (²²⁶Ra/²³⁰Th)-(²⁰⁶Pb/²⁰⁴Pb) diagrams. Shallow magma chamber processes like magma mixing, fractional crystallization and wall rock assimilation do not account for the correlations. Crustal contamination is not the cause of the observed isotopic variations because continental crust is considered to have extremely different Pb isotope compositions and U/Th ratios. Melting of a chemically heterogeneous mantle might explain the Mt. Cameroon data, but dynamic melting under conditions of high D_U and D_U/D_Th, long magma ascent time, or disequilibrium mineral/melt partitioning, is required. The most plausible scenario to produce the geochemical characteristics of Mt. Cameroon samples is the interaction of melt derived from the asthenospheric mantle with overlying sub-continental lithospheric mantle which has elevated U/Pb (>0.75) and Pb isotope ratios (²⁰⁶Pb/²⁰⁴Pb > 20.47) due to late Mesozoic metasomatism.     

Origin of hybrid ferrolatite lavas from Magic Reservoir eruptive center,Snake River Plain,Idaho

The mineralogy and geochemical characteristics of intermediate composition ferrolatites and related lavas from the Magic Reservoir eruptive center (central Snake River Plain) have been investigated to evaluate the origin and petrologic significance of these hybrid lavas. The ferrolatites are chemically uniform, but contain a disequilibrium phenocryst/xenocryst assemblage derived in part from mixed rhyolitic and basaltic magmas that are closely represented by extrusive units in the area. The hybrid lavas also contain xenoliths of Archean granulites and have high ⁸⁷Sr/ ⁸⁶Sr and low ¹⁴³Nd/¹⁴⁴Nd ratios, all of which suggest significant magma-crust interaction. Quantitative models including magma mixing, minor crystal fractionation, and crustal contamination very closely reproduce the observed compositions of these ferrolatites; closed system fractionation and (or) simple bulk contamination models are not as successful and can be ruled out. It appears that preexisting mafic and silicic magmas from distinct sources (e.g., mantle and crust) encounter one another in crustal-level magma chambers under conditions where intimate mixing may occur despite wide differences in the physical properties of these liquids.     

A Ce/Nd isotope study of crustal contamination processes affecting Palaeocene magmas in Skye,Northwest Scotland

Twelve¹³⁸Ce/¹³⁶Ce isotope determinations, 31 Nd isotope analyses, and 31 REE profiles are presented for Tertiary basic to intermediate igneous rocks from the Isle of Skye, NW Scotland. The aim of this work is to precisely identify the contamination mechanisms of basic magmas emplaced through old crust, and to test the effectiveness of Ce isotope analysis as a petrogenetic tool.Combined Ce/Nd isotope analysis enables the modelling of the light REE profiles of the mantle-derived precursors to contaminated lavas, using different crustal end-members, in order to compare these with the magmatic lineage of uncontaminated Skye lavas. The geochemical data support a contamination mechanism involving a granitic melt, produced either by large degree melting of Scourian granulitefacies acid sheets, or (possibly) by melting of intermediate gneiss out of isotopic equilibrium.Basaltic lavas showing strong isotopic contamination effects yield calculated degrees of crustal contamination by large degree granitic melts of ca. 8 or 9% based on Ce and Nd isotopic data respectively. However, for lavas with liquidus temperatures of over 1250° C, the temperature dependence of the degree of contamination is weak.The combination of this evidence with new and published Pb isotope data suggests that the bulk of crustal contamination of the Skye lavas occurred in sill complexes at distinct levels in the crust, rather than during the actual ascent of magma through the crust in dykes. It is suggested on the basis of published fluid dynamic and field evidence that the assimilation of large degree melts of acid gneiss by turbulently flowing magma is more likely than assimilation of small degree disequilibrium melts from more refractory intermediate gneisses.It is concluded that Ce isotope analysis is a viable and useful adjunct to Nd isotope data in petrogenetic studies of continental igneous rocks emplaced through old basement.     

Generation of Deccan Trap magmas

Deccan Trap magmas may have erupted through multiple centers, the most prominent of which may have been a shield volcano-like structure in the Western Ghats area. The lavas are predominantly tholeiitic; alkalic mafic lavas and carbonatites are rare. Radioisotope dating, magnetic chronology, and age constraints from paleontology indicate that although the eruption started some 68 Ma, the bulk of lavas erupted at around 65–66 Ma. Paleomagnetic constraints indicate an uncertainty of ± 500,000 years for peak volcanic activity at 65 m.y. in the type section of the Western Ghats. Maximum magma residence times were calculated in this study based on growth rates of “giant plagioclase” crystals in lavas that marked the end phase of volcanic activity of different magma chambers. These calculations suggest that the > 1.7 km thick Western Ghats section might have erupted within a much shorter time interval of ∼ 55,000 years, implying phenomenal eruption rates that are orders of magnitude larger than any present-day eruption rate from any tectonic environment. Other significant observations/conclusions are as follows: (1) Deccan lavas can be grouped into stratigraphic subdivisions based on their geochemistry; (2) While some formations are relatively uncontaminated others are strongly contaminated by the continental crust; (3) Deccan magmas were produced by 15–30% melting of a Fe-rich lherzolitic source at ∼ 3–2 GPa; (4) Parent magmas of the relatively uncontaminated Ambenali formation had a primitive composition with 16%MgO, 47%SiO₂; (5) Deccan magmas were generated much deeper and by significantly more melting than other continental flood basalt provinces; (6) The erupted Deccan tholeiitic lavas underwent fractionation and magma mixing at ∼ 0.2 GPa. The composition and origin of the crust and crust/mantle boundary beneath the Deccan are discussed with respect to the influence of Deccan magmatic episode.     

Sr isotopic evidence for a multi-source origin of the potassic magmas in the neapolitan area (S. Italy)

New Sr isotopic data on lavas and xenoliths from Somma-Vesuvius and other nearby volcanic areas (Phlegrean Fields and Ischia) are presented. Chemical and isotopic evidences show that not all the Phlegrean Fields rocks belong to the low K series, but some of them may be interpreted as low pressure differentiates of Somma magmas, i.e. as a part of the high K series. Two rock groups are defined in the Ischia low K series, which are well identified both in time and in chemical and isotopic features, and cannot be derived from the same magma source. The low K series in the studied area generally has lower Sr isotopic values than the high K series.Historical Vesuvian lavas show two distinct linear trends with negative slopes when⁸⁷Sr/⁸⁶Sr ratios are plotted against their ages of eruption. Such trends are interpreted to result from mixing of magmas in two separate reservoirs. Evidence from the Vesuvian ejecta shows that Somma-Vesuvius magmas underwent high or low pressure fractionation, in connection with different events of the Vesuvian activity. Distinct magma reservoirs developed episodically at different depths. Isotopic and geochemical evidences do not favour large scale assimilation of crustal materials by Somma-Vesuvius magmas, but instead appear to reflect mantle characteristics.A minimum of three different (inhomogeneous) source regions is necessary to account for the isotopic features of the studied rocks.     

Oxygen isotope composition of phenocrysts from Tristan da Cunha and Gough Island lavas: variation with fractional crystallization and evidence for assimilation

We have measured the δ¹⁸O values of the major phenocrysts (olivine, clinopyroxene and plagioclase) present in lavas from Tristan da Cunha and Gough Island. These islands, which result from the same mantle plume, have enriched radiogenic isotope ratios and are, therefore, prime candidates for an oxygen isotope signature that is distinct from that of MORB. Consistent differences between the δ¹⁸O values of olivine, pyroxene and feldspar in the Gough lavas show that the phenocrysts in the mafic Gough Island lavas are in oxygen isotope equilibrium. The olivines in lavas with SiO₂ <50 wt% have a mean δ¹⁸O value of 5.19‰, consistent with crystallization from a magma having the same oxygen isotope composition as MORB. Phenocrysts in all the Gough lavas show a systematic increase in δ¹⁸O value as silica content increases, which is consistent with closed-system fractional crystallization. The lack of enrichment in δ¹⁸O of the Gough magmas suggests that the mantle source contained <2% recycled sediment. In contrast, the Tristan lavas with SiO₂ >48 wt% contain phenocrysts which have δ¹⁸O values that are systematically ∼0.3‰ lower than their counterparts from Gough. We suggest that the parental mafic Tristan magmas were contaminated by material from the volcanic edifice that acquired low δ¹⁸O values by interaction with water at high temperatures. The highly porphyritic SiO₂-poor lavas show a negative correlation between olivine δ¹⁸O value and whole-rock silica content rather than the expected positive correlation. The minimum δ¹⁸O value occurs at an SiO₂ content of about 45 wt%. Below 45 wt% SiO₂, magmas evolved via a combination of assimilation, fractionational crystallization and crystal accumulation; above 45 wt% SiO₂, magmas appeared to have evolved via closed-system fractional crystallization. Received: 23 November 1998 / Accepted: 27 September 1999     

Petrogenesis of isotopically unusual Pliocene olivine leucitites from Deep Springs Valley, California

High-K mafic alkalic lavas (5.4 to 3.2 wt% K₂O) from Deep Springs Valley, California define good correlations of increasing incompatible element (e.g., Sr, Zr, Ba, LREE) and compatible element contents (e.g., Ni, Cr) with increasing MgO. Strontium and Nd isotope compositions are also correlated with MgO; ⁸⁷Sr/⁸⁶Sr ratios decrease and ɛ_Nd values increase with decreasing MgO. The Sr and Nd isotope compositions of these lavas are extreme compared to most other continental and oceanic rocks; ⁸⁷Sr/⁸⁶Sr ratios range from 0.7121 to 0.7105 and ɛ_Nd values range from −16.9 to −15.4. Lead isotope ratios are relatively constant, ²⁰⁶Pb/²⁰⁴Pb ∼17.2, ²⁰⁷Pb/²⁰⁴Pb ∼15.5, and ²⁰⁸Pb/²⁰⁴Pb ∼38.6. Depleted mantle model ages calculated using Sr and Nd isotopes imply that the reservoir these lavas were derived from has been distinct from the depleted mantle reservoir since the early Proterozoic. The Sr-Nd-Pb isotope variations of the Deep Springs Valley lavas are unique because they do not plot along either the EM I or EM II arrays. For example, most basalts that have low ɛ_Nd values and unradiogenic ²⁰⁶Pb/²⁰⁴Pb ratios have relatively low ⁸⁷Sr/⁸⁶Sr ratios (the EM I array), whereas basalts with low ɛ_Nd values and high ⁸⁷Sr/⁸⁶Sr ratios have radiogenic ²⁰⁶Pb/²⁰⁴Pb ratios (the EM II array). High-K lavas from Deep Springs Valley have EM II-like Sr and Nd isotope compositions, but EM I-like Pb isotope compositions. A simple method for producing the range of isotopic and major- and trace-element variations in the Deep Springs Valley lavas is by two-component mixing between this unusual K-rich mantle source and a more typical depleted mantle basalt. We favor passage of MORB-like magmas that partially fused and were contaminated by potassic magmas derived from melting high-K mantle veins that were stored in the lithospheric mantle. The origin of the anomalously high ⁸⁷Sr/⁸⁶Sr and ²⁰⁸Pb/²⁰⁴Pb ratios and low ɛ_Nd values and ²⁰⁶Pb/²⁰⁴Pb ratios requires addition of an old component with high Rb/Sr and Th/Pb ratios but low Sm/Nd and U/Pb ratios into the mantle source region from which these basalts were derived. This old component may be sediments that were introduced into the mantle, either during Proterozoic subduction, or by foundering of Proterozoic age crust into the mantle at some time prior to eruption of the lavas. Received: 28 February 1997 / Accepted: 9 July 1998     

Mixing of rhyolite,trachyte and basalt magma erupted from a vertically and laterally zoned reservoir,composite flow P1, Gran Canaria

The 14.1 Ma composite welded ignimbrite P1 (45 km³ DRE) on Gran Canaria is compositionally zoned from a felsic lower part to a basaltic top. It is composed of four component magmas mixed in vertically varying proportions: (1) Na-rhyolite (10 km³) zoned from crystal-poor to highly phyric; (2) a continuously zoned, evolved trachyte to sodic trachyandesite magma group (6 km³); (3) a minor fraction of Na-poor trachyandesite (<1 km³); and (4) nearly aphyric basalt (26 km³) zoned from 4.3 to 5.2 wt% MgO. We distinguish three sites and phases of mixing: (a) Mutual mineral inclusions show that mixing between trachytic and rhyolitic magmas occurred during early stages of their intratelluric crystallization, providing evidence for long-term residence in a common reservoir prior to eruption. This first phase of mixing was retarded by increasing viscosity of the rhyolite magma upon massive anorthoclase precipitation and accumulation. (b) All component magmas probably erupted through a ring-fissure from a common upper-crustal reservoir into which the basalt intruded during eruption. The second phase of mixing occurred during simultaneous withdrawal of magmas from the chamber and ascent through the conduit. The overall withdrawal and mixing pattern evolved in response to pre-eruptive chamber zonation and density and viscosity relationships among the magmas. Minor sectorial variations around the caldera reflect both varying configurations at the conduit entrance and unsteady discharge. (c) During each eruptive pulse, fragmentation and particulate transport in the vent and as pyroclastic flows caused additional mixing by reducing the length scale of heterogeneities. Based on considerations of magma density changes during crystallization, magma temperature constraints, and the pattern of withdrawal during eruption, we propose that eruption tapped the P1 magma chamber during a transient state of concentric zonation, which had resulted from destruction of a formerly layered zonation in order to maintain gravitational equilibrium. Our model of magma chamber zonation at the time of eruption envisages a basal high-density Na-poor trachyandesite layer that was overlain by a central mass of highly phyric rhyolite magma mantled by a sheath of vertically zoned trachyte-trachyandesite magma along the chamber walls. A conventional model of vertically stacked horizontal layers cannot account for the deduced density relationships nor for the withdrawal pattern.     

Petrogenesis and evolution of Mt. Vulture alkaline volcanism (Southern Italy)

Summary The Late Pleistocene Mt. Vulture strato-volcano developed at the intersection of NE-SW and NW-SE lithospheric fault systems, on the easternmost border of the Apennine compressional front overthrust onto the Apulian foreland. The initial phase of the volcanic activity is represented by pyroclastic deposits, including lava blocks, and subordinate eccentric domes, mostly phonolitic in composition. The later stages of activity formed the bulk of the strato-volcano (pyroclastic products and subordinate lavas), mostly tephritic in composition, with minor intercalations of basanite, mela-foidite and melilitite lavas and dikes. Variations in rock and mineral composition suggest that the volumetrically predominant basanite-tephrite (foidite)-phonotephrite-phonolite series can be accounted for by fractional crystallization processes starting from basanitic parental magmas, in agreement with the remarkably constant ⁸⁷Sr/⁸⁶Sr isotopes (0.70586–0.70581). Mass-balance calculations indicate that the variably differentiated magmas may have been produced by removal of wehrlite, clinopyroxenite and syenite cumulates, some of which are occasionally found as cognate xenoliths in the volcanics. Fractionation processes probably developed in multiple-zoned magma chambers, at depths of 3–5 km, corresponding to the tectonic discontinuity between the allochthonous Apennine formations and the underlying Apulian platform. Highly differentiated phonolitic magmas capping the magma chambers and their conduits thus appear to have fed the initial volcanic activity, whereas dominantly tephritic products were erupted in later stages. The least evolved mafic magmas, namely basanites, mela-foidites and melilitites, are characterized by diverse Na/K ratios and critical SiO₂-undersaturation, which indicate their derivation as independent melts generated from distinct, heterogeneously enriched mantle sources and by variable partial melting degrees. Primitive mantle-normalized incompatible element patterns of Vulture mafic lavas invariably share analogies with both orogenic subduction-related magmas (high Low Field Strength Elements/High Field Strength Elements ratios, K, Rb and Th contents and marked Ti and Nb negative anomalies) and alkaline lavas from within-plate and rift settings (high Light Rare Earth Elements, P, Zr, Nb and Na). These geochemical features may be accounted for by magma generation from deep lithospheric mantle sources, enriched in Na-alkali silicate/carbonatite anorogenic components, subsequently affected by orogenic subduction-related K-metasomatism, analogous to that which modified magma sources of the Roman Magmatic Province along the internal Apennine Chain. Received April 12, 2000; revised version accepted June 7, 2001     

The origin of high-Mg magmas in Mt Shasta and Medicine Lake volcanoes,Cascade Arc (California): higher and lower than mantle oxygen isotope signatures attributed to current and past subduction

We report the oxygen isotope composition of olivine and orthopyroxene phenocrysts in lavas from the main magma types at Mt Shasta and Medicine Lake Volcanoes: primitive high-alumina olivine tholeiite (HAOT), basaltic andesites (BA), primitive magnesian andesites (PMA), and dacites. The most primitive HAOT (MgO > 9 wt%) from Mt. Shasta has olivine δ¹⁸O (δ¹⁸O_Ol) values of 5.9–6.1‰, which are about 1‰ higher than those observed in olivine from normal mantle-derived magmas. In contrast, HAOT lavas from Medicine Lake have δ¹⁸O_Ol values ranging from 4.7 to 5.5‰, which are similar to or lower than values for olivine in equilibrium with mantle-derived magmas. Other magma types from both volcanoes show intermediate δ¹⁸O_Ol values. The oxygen isotope composition of the most magnesian lavas cannot be explained by crustal contamination and the trace element composition of olivine phenocrysts precludes a pyroxenitic mantle source. Therefore, the high and variable δ¹⁸O_Ol signature of the most magnesian samples studied (HAOT and BA) comes from the peridotitic mantle wedge itself. As HAOT magma is generated by anhydrous adiabatic partial melting of the shallow mantle, its 1.4‰ range in δ¹⁸O_Ol reflects a heterogeneous composition of the shallow mantle source that has been influenced by subduction fluids and/or melts sometime in the past. Magmas generated in the mantle wedge by flux melting due to modern subduction fluids, as exemplified by BA and probably PMA, display more homogeneous composition with only 0.5‰ variation. The high-δ¹⁸O values observed in magnesian lavas, and principally in the HAOT, are difficult to explain by a single-stage flux-melting process in the mantle wedge above the modern subduction zone and require a mantle source enriched in ¹⁸O. It is here explained by flow of older, pre-enriched portions of the mantle through the slab window beneath the South Cascades.     

U-series evidence for crustal involvement and magma residence times in the petrogenesis of Parinacota volcano, Chile

In this study, we present Th–U disequilibria as well as radiogenic and trace element data for recent volcanic rocks from the Nevados de Payachata volcano which erupted through ∼70 km of continental crust in the Central Volcanic Zone of the Andes (18°S, 69°W). Both lavas and mineral separates were analyzed by mass spectrometry for ²³⁸U–²³⁰Th disequilibria. The lavas are characterized either by ²³⁰Th enrichment or depletion relative to its parent nuclide ²³⁸U. Mineral separates are used to derive U–Th isochron ages and these ages compare favorably with inferred stratigraphic ages or K–Ar ages, although in one case the U–Th age is significantly older than the stratigraphic age. Despite relatively constant Sr, Nd, and Pb isotope ratios, the lavas display inverse trends in ²³⁰Th/²³⁸U versus Ce/Yb or Ba/Hf diagrams. These trends cannot be interpreted by simple two-component mixing. Rather, there must be three (and perhaps four components) involved in the genesis of the Parinacota lavas. A mantle wedge, a slab fluid, and a lower crustal component can be identified. A sediment component is more difficult to detect as it is difficult to decipher its signature because of the strong crustal influence. The existence of binary arrays can be explained by variable amounts of crustal material. The process of crust–mantle interaction must have been short enough to preserve U–Th disequilibrium (<300 ka). Received: 21 April 1999 / Accepted: 11 March 2000     

Mineral isotope evidence for the contemporaneous process of Mesozoic granite emplacement and gneiss metamorphism in the Dabie orogen

Kinetics of isotopic equilibrium in the mineral radiometric systems of igneous and metamorphic rocks is an important issue in geochronology. It turns out that temperature is the most important factor in dictating isotopic equilibrium or disequilibrium with respect to diffusion mechanism. Contemporaneous occurrence of Mesozoic granites and gneisses in the Dabie orogen of China allows us to evaluate the thermal effect of magma emplacement and associated metamorphism on mineral radiometric systems. Zircon U-Pb, mineral Rb-Sr and O isotope analyses were carried out for a Cretaceous granite and its host gneiss (foliated granite) from North Dabie. Zircon U-Pb dating gave consistently concordant ages of 127 ± 3 Ma and 128 ± 2 Ma for the granite and the gneiss, respectively. A direct correspondence in equilibrium state is observed between the O and Rb-Sr isotope systems of both granitic and gneissic minerals. Mineral O isotope temperatures correlate with O diffusion closure temperatures under conditions of slow cooling, indicating attainment and preservation of O isotope equilibrium in these minerals. The mineral Rb-Sr isochron of granite, constructed by biotite, feldspar, apatite and whole-rock with the O isotope equilibrium, yields a meaningful age of 118 ± 3 Ma, which is in accordance with the mineral Rb-Sr isochron age of 122 ± 1 Ma for the host gneiss. The consistency in both U-Pb and Rb-Sr ages between the granite and the gneiss suggests a contemporaneous process of crystallizing the zircons and resetting the Rb-Sr radiometric systems during magma emplacement and granite foliation. Whereas the zircon U-Pb ages for both granite and gneiss are interpreted as the timing of magma crystallization, the young Rb-Sr isochron ages record the timing of Sr diffusion closure during the slow cooling. Protolith of the gneiss crystallized shortly before intrusion of the granite, so that it was able to be foliated by voluminous emplacement of coeval mafic to felsic magmas derived by anatexis of orogenic lithospheric keel. Therefore, extensional collapse of collision-thickened crust at Early Cretaceous is suggested to trigger the post-collisional magmatism, which in turn serves as an essential driving force for the contemporaneous high-T deformation/metamorphism.     

Geochemical variations in Andean basaltic and silicic lavas from the Villarrica-Lanin volcanic chain (39.5° S): an evaluation of source heterogeneity,fractional crystallization and crustal assimilation

At 39.5° S in the southern volcanic zone of the Andes three Pleistocene-recent stratovolcanoes, Villarrica, Quetrupillan and Lanin, form a trend perpendicular to the strike of the Andes, 275 to 325 km from the Peru-Chile trench. Basalts from Villarrica and Lanin are geochemically distinct; the latter have higher incompatible element abundances and La/Sm but lower Ba/La and alkali metal/La ratios. These differences are consistent with our previously proposed models involving: a) a west to east decrease in an alkali metal-rich, high Ba/La slab-derived component which causes an across strike decrease in degree of melting; or b) a west to east increase in the contamination of subduction-related magma by enriched subcontinental lithospheric mantle. Silicic and mafic lavas from the stratovolcanoes have overlapping Sr, Nd and O isotopic ratios. Silicic lavas also have geochemical differences that parallel those of their associated basalts, e.g., rhyolite from Villarrica has lower La/Sm and incompatible element contents than high-SiO₂ andesite from Lanin. At each volcano the most silicic lavas can be modelled by closed system fractional crystallization while andesites are best explained by magma mixing. Apparently crustal contamination was not an important process in deriving the evolved lavas. Basaltic flows from small scoria cones, 20–35 km from Villarrica volcano have high incompatible element contents and low Ba/La, like Lanin basalts, but trend to higher K/Rb (356–855) and lower ⁸⁷Sr/ ⁸⁶Sr (0.70361–0.70400) than basalts from either stratovolcano. However all basalts have similar Nd, Pb and O isotope ratios. The best explanation for the unique features of the cones is that the sources of SVZ magmas, e.g., slab-derived fluids or melts of the subcontinental lithospheric mantle, have varying alkali metal and radiogenic Sr contents. These heterogeneities are not manifested in stratovolcano basalts because of extensive subcrustal pooling and mixing. This model is preferable to one involving crustal contamination because it can account for variable Sr isotope ratios and uniform Nd and Pb isotope ratios among the basalts, and the divergence of the cones from across-strike geochemical trends defined by the stratovolcanoes.     

Petrology of mafic lavas within the Onega plateau, central Karelia: evidence for 2.0 Ga plume-related continental crustal growth in the Baltic Shield

The Onega plateau constitutes part of a vast continental flood basalt province in the SE Baltic Shield. It consists of Jatulian-Ludikovian submarine volcanic, volcaniclastic and sedimentary sequences attaining in places 4.5?km in thickness. The parental magmas of the lavas contained ～10% MgO and were derived from melts generated in the garnet stability field at depths 80–100?km. The Sm-Nd mineral and Pb-Pb whole-rock isochron ages of 1975?±?24 and 1980?±?57 Ma for the upper part of the plateau and a SHRIMP U-Pb zircon age of 1976?±?9 Ma for its lower part imply the formation of the entire sequence within a short time span. These ages coincide with those of picrites in the Pechenga-Imandra belt (the Kola Peninsula) and komatiites and basalts in the Karasjok-Kittilä belt (Norway and Finnmark). Together with lithostratigraphic, chemical and isotope evidence, these ages suggest the derivation of the three provinces from a single large (～2000?km in diameter) mantle plume. These plume-generated magmas covered ～600,000?km² of the Baltic Shield and represent a major contribution of juvenile material to the existing continental crust at 2.0 Ga. The uppermost Onega plateau lavas have high (Nb/Th)_N?=?1.4–2.4, (Nb/La)_N= 1.1–1.3, positive ?Nd(T) of +3.2 and unradiogenic Pb-isotope composition (μ₁?= 8.57), comparable with those of modern oceanic plume-derived magmas (oceanic flood basalt and ocean island basalt). These parameters are regarded as source characteristics. The lower sequences have (Nb/Th)_N= 0.58–1.2, (Nb/La)_N= 0.52–0.88 and ?Nd(T) =?2.6. They have experienced mixing with 10–30% of continental crust and resemble contaminated lavas from other continental flood basalt provinces. The estimated Nb/U ratios of 53?±?4 in the uncontaminated rocks are similar to those found in the modern mantle (～47) suggesting that by 2.0 Ga a volume of continental crust similar to the present-day value already existed.