Geochemistry of silicic magmas in the Macolod Corridor, SW Luzon, Philippines: evidence of distinct, mantle-derived, crustal sources for silicic magmas

Silicic volcanic deposits (>65 wt% SiO₂), which occur as domes, lavas and pyroclastic deposits, are relatively abundant in the Macolod Corridor, SW Luzon, Philippines. At Makiling stratovolcano, silicic domes occur along the margins of the volcano and are chemically similar to the silicic lavas that comprise part of the volcano. Pyroclastic flows are associated with the Laguna de Bay Caldera and these are chemically distinct from the domes and lavas at Makiling stratovolcano. As a whole, samples from the Laguna de Bay Caldera contain lower concentrations of MgO and higher concentrations of Fe₂O_3(t) than the samples from domes and lavas. The Laguna de Bay samples are more enriched in incompatible trace elements. The silicic rocks from the domes, Makiling Volcano and Laguna de Bay Caldera all contain high alkalis and high K₂O/Na₂O ratios. Melting experiments of primitive basalts and andesites demonstrate that it is difficult to produce high K₂O/Na₂O silicic magmas by fractional crystallization or partial melting of a low K₂O/Na₂O source. However, recent melting experiments (Sisson et al., Contrib Mineral Petrol 148:635–661, 2005) demonstrate that extreme fractional crystallization or partial melting of K-rich basalts can produce these silicic magmas. Our model for the generation of the silicic magmas in the Macolod Corridor requires partial melting of mantle-derived, evolved, moderate to K-rich, crystallized calc-alkaline magmas that ponded and crystallized in the mid-crust. Major and trace element variations, along with oxygen isotopes and ages of the deposits, are consistent with this model. Electronic Supplementary Material Supplementary material is available for this article at     

地壳硅质岩浆演化的动力学机制SCIEI北大核心CSCD

硅质岩浆的成因及演化机制研究是认识大陆地壳成分结构演化机制的关键途径。地壳硅质岩浆储库的基本性质及其演化的动力学过程是制约硅质岩浆演化及其多样性的根本原因,也是受到广泛关注的前沿和热点问题。本文总结了硅质岩浆储库的基本性质及其演化机制研究的相关进展,重点探讨了岩浆体系演化的动力学机制,即其物质成分及温度压力条件的动态演化过程,亦即岩浆储库在补给驱动下的存留、活化、分异、喷发过程。经由累积生长的方式形于地壳浅部的硅质岩浆储库以晶粥为主要赋存形式,具有较低的温度压力范围和高硅高结晶度高粘度的特征。岩浆储库的基本性质导致其自身不具有持续演化的能力,只有在基性岩浆补给的驱动下才能长时间存留分异或快速活化喷发。基性岩浆的补给及其携带的热及流体/挥发份深刻地改变了岩浆储库的物质成分及物理条件,是导致硅质岩浆分异演化或喷发的根本原因。本文还结合研究进展探讨了东南沿海地区中生代火山-侵入岩的成因演化机制及相关问题。     

Zircon from historic eruptions in Iceland: reconstructing storage and evolution of silicic magmas

Zoning patterns, U-Th disequilibria ages, and elemental compositions of zircon from eruptions of Askja (1875 AD), Hekla (1158 AD), ?r?faj?kull (1362 AD) and Torfaj?kull (1477 AD, 871 AD, 3100 BP, 7500 BP) provide insights into the complex, extended, histories of silicic magmatic systems in Iceland. Zircon compositions, which are correlated with proximity to the main axial rift, are distinct from those of mid-ocean ridge environments and fall at the low-Hf edge of the range of continental zircon. Morphology, zoning patterns, compositions, and U-Th ages all indicate growth and storage in subvolcanic silicic mushes or recently solidified rock at temperatures above the solidus but lower than that of the erupting magma. The eruptive products were likely ascending magmas that entrained a zircon ??cargo?? that formed thousands to tens of thousands of years prior to the eruptions.     

浙东晚白垩世酸性岩浆的自混合作用及其意义

岩浆混合作用是造成火成岩多样性的主要原因之一,也是诱发火山喷发的重要机制。以往的研究多集中于基性和酸性岩浆之间的混合作用,但近年来酸性岩浆之间的混合作用受到越来越多的关注和研究。本文报道了浙东小雄破火山一个次级火山口内粗面质和流纹质两种酸性岩浆之间的混合现象。野外调查及岩相学研究显示,粗面质岩浆多呈大小不一的条带状以及透镜体状分布于流纹质岩浆内,局部发生扩散,粗面岩中斑晶大多为粗大的正长石斑晶,强烈熔蚀且聚斑结构普遍;在副矿物聚晶(由钛磁铁矿+磷灰石+锆石组成)的周围常可见反应边结构。流纹岩的斑晶主要由正长石、透长石及石英组成,晶体粒径较小,且熔蚀现象不发育。全岩主、微量元素特征及其他地质证据均显示,两种酸性岩浆之间以机械混合为主,其地球化学成分变化趋势主要受结晶分异过程控制。粗面质及流纹质岩浆在矿物组成、结构等方面的差异表明两者来源于同一层状岩浆房内的不同部位,其中粗面质岩浆应代表岩浆房底部及边部富晶体、贫熔体的粥状层部分(正长石+磁钛铁矿+锆石+磷灰石);而分异程度较高的流纹质岩浆则聚集于岩浆房上部形成富熔体、贫晶体的部分。两种酸性岩浆的混合现象是它们在地壳浅部层状岩浆房内自混合的结果,这一过程可能受岩浆房底部基性岩浆的聚集作用所控制,当更热、更基性的岩浆聚集时,岩浆房下部晶粥区内的粗面质岩浆迅速升温、活化,从而向上运移并与上部富熔体贫晶体的流纹质岩浆发生自混合作用。这一发现为我们理解中国东南沿海地区晚中生代大规模酸性火山喷发及岩浆演化机制、岩浆房结构提供了重要的参考,同时也为认识地壳浅部岩浆房内岩浆之间的自混合作用提供了可靠的例证。     

Experimental study of liquid evolution in an Fe-rich,layered mafic intrusion: constraints of Fe-Ti oxide precipitation on the T-fO 2 and T-ϱ paths of tholeiitic magmas

The Newark Island layered intrusion, a composite intrusion displaying a similar fractionation sequence to the Skaergaard, has both dikes which preserved liquids fed into the intrusion and chilled pillows of liquids resident in the chamber. This study reports experimentally determined one atmosphere liquid lines of descent of these compositions as a function of oxygen fugacity which varies from QFM (quartz-fayalite-magnetite) to 0.5 log₁₀ units above IW (iron-wustite). These experiments reveal a strong oxygen fugacity dependence on the order of appearance and relative abundances of the Fe–Ti oxide minerals. Titanomagnetite saturates prior to ilmenite at QFM, but the order is reversed at lower oxygen fugacities. In the layered series of the Newark Island intrusion, ilmenite arrives shortly before titanomagnetite and the titanomagnetite/ilmenite ratio decreases monotonically after the cumulus appearance of titanomagnetite. Comparison of the crystallization sequence in the intrusion with that of the experiments requires that the oxygen fugacity in the intrusion increased relative to QFM before titanomagnetite saturation and decreased afterward, but always remained between the QFM and IW buffers. Similar trends in the modes of the Fe–Ti oxides (ilmenite and titanomagnetite) in the Skaergaard, Kiglapait, and Somerset Dam intrusions along with Fe₂O₃/FeO ratios in MORBs suggest that such a temperature-oxygen fugacity path may be typical of tholeiitic magma differentiation. Calculations of the temperature-density paths of the experimental liquids indicate that, at all possible oxygen fugacities, the density must have decreased abruptly after Fe–Ti oxide saturation. Accordingly, liquids replenishing the intrusion after Fe–Ti oxide saturation should pond at the bottom of the chamber, quenching against older cumulates. Field observation at the Newark Island intrusion confirm this prediction. The similarities in the fractionation paths of several other layered intrusions to that of the Newark Island intrusion suggest that the density of the liquids in these intrusions also decreased after Fe–Ti oxide saturation. Experiments on a suggested initial Skaergaard liquid are consistent with this model.     

Proterozoic VanDieland in Central Victoria: ages,compositions and source depths for late devonian silicic magmas

The Proterozoic to Cambrian VanDieland microcontinent was accreted to mainland Australia at ca 400?Ma, and its northern tip, the Selwyn Block, forms the basement in central Victoria. Here, mainly Late Devonian, silicic magmas were derived from the Selwyn Block and intruded into the shallow crust. We use the phase petrology of Late Devonian, S-type rhyolitic ignimbrites and a xenolith of pelitic migmatite, together with Nd-model ages for the silicic magmatic rocks to constrain the lithological characteristics of the metasedimentary component of the Selwyn Block, to infer minimum depths and temperature conditions here in the Late Devonian, and the likely ages of the source rocks for the S-type magmas. The most abundant source rocks are inferred to be volcaniclastic metagreywackes, with minor metadacites, meta-andesites and metapelites. The metapelitic xenolith cannot have been the source for any of the silicic magmas but constrains the upper amphibolite-facies part of the Selwyn Block to depths around 17?km, where temperatures reached ～775?°C. The older ignimbrite magma was formed by partial melting at perhaps 770?°C and a depth of at least 33?km, while the younger ignimbrite magma formed at ～23?km and 900?°C. These depths suggest source rocks in the Paleoproterozoic to Mesoproterozoic lower parts of the Selwyn Block. Nd-model ages of the silicic magmatic rocks confirm a dominance of Mesoproterozoic to Paleoproterozoic sources. If the inferred rock types in the Mesoproterozoic formations were as current correlations suggest, the sources for the Late Devonian silicic magmas would have to lie in the undocumented Paleoproterozoic basement of the Selwyn Block. Rock types here must include andesitic to dacitic volcanic components as well as volcaniclastic greywackes and minor pelites, which suggests a continental arc setting. The Late Devonian magmatism in the region may record the progression from amphibolite- to granulite-facies conditions during post-orogenic extension, with heat advected to the crust by mantle-derived mafic magmas. These processes would have resulted in mafitisation of the deep Selwyn Block.     

18O-Enrichment of silicic magmas caused by crystal fractionation at the Galapagos Spreading Center

An extremely differentiated suite of unaltered volcanic rocks dredged from the Galapagos Spreading Center ranges in ¹⁸O from 5.7 to 7.1 At 95°W, low K-tholeiites, FeTi-basalts, andesites and rhyodacites were recovered. Their lithologic and major element geochemical variation can be accounted for by crystal fractionation of plagioclase, pyroxenes, olivine and titanomagnetite in the same proportions and amounts needed to model the ¹⁸O variation by simple Rayleigh fractionation. More complicated behaviour was observed in a FeTi-basalt suite from 85°W. This study shows that 90% fractionation only enriches the residual melt by about 1.2 in ¹⁸O. It also implies that the magma chambers along parts of the Galapagos Spreading Center were static and isolated such that extreme differentiation could occur.     

Evolution of silicic magmas in the Kos-Nisyros volcanic center,Greece: a petrological cycle associated with caldera collapse

Multiple eruptions of silicic magma (dacite and rhyolites) occurred over the last ~3 My in the Kos-Nisyros volcanic center (eastern Aegean sea). During this period, magmas have changed from hornblende-biotite-rich units with low eruption temperatures (≤750–800°C; Kefalos and Kos dacites and rhyolites) to hotter, pyroxene-bearing units (>800–850°C; Nisyros rhyodacites) and are transitioning back to cooler magmas (Yali rhyolites). New whole-rock compositions, mineral chemistry, and zircon Hf isotopes show that these three types of silicic magmas followed the same differentiation trend: they all evolved by crystal fractionation and minor crustal assimilation (AFC) from parents with intermediate compositions characterized by high Sr/Y and low Nb content, following a wet, high oxygen fugacity liquid line of descent typical of subduction zones. As the transition between the Kos-Kefalos and Nisyros-type magmas occurred immediately and abruptly after the major caldera collapse in the area (the 161 ka Kos Plateau Tuff; KPT), we suggest that the efficient emptying of the magma chamber during the KPT drew out most of the eruptible, volatile-charged magma and partly solidified the unerupted mush zone in the upper crust due to rapid unloading, decompression, and coincident crystallization. Subsequently, the system reestablished a shallow silicic production zone from more mafic parents, recharged from the mid to lower crust. The first silicic eruptions evolving from these parents after the caldera collapse (Nisyros units) were hotter (up to >100°C) than the caldera-forming event and erupted from reservoirs characterized by different mineral proportions (more plagioclase and less amphibole). We interpret such a change as a reflection of slightly drier conditions in the magmatic column after the caldera collapse due to the decompression event. With time, the upper crustal intermediate mush progressively transitioned into the cold-wet state that prevailed during the Kefalos-Kos stage. The recent eruptions of the high-SiO₂ rhyolite on Yali Island, which are low temperature and hydrous phases (sanidine, quartz, biotite), suggest that another large, potentially explosive magma chamber is presently building under the Kos-Nisyros volcanic center.     

A regular solution model for met-aluminous silicate liquids: Applications to geothermometry,immiscibility, and the source regions of basic magmas

Adopting a set of multioxide components and using published compositional data on olivineand plagioclase-liquid equilibria we have developed a 17 component regular solution model for met-aluminous silicate liquids. The partial molar excess free energies predicted from this model can be used together with phenocryst compositions as an effective geothermometer, with an approximate error of 20 °C (30 °C for olivine, 12 °C for plagioclase). The regular solution formulation is also successful in predicting liquid immiscibility at (1) high mole fractions of silica commonly observed in phase diagrams, and at (2) lower temperatures in lunar basalts and intermediate lavas. The model yields activities of silica which are consistent with those obtained from solid-liquid silica buffers in rocks which contain olivine and enstatite or quartz. From predicted activities of KAlSi₃O₈ in liquids coexisting with plagioclase a value is obtained for the limiting Henry's law activity coefficient of KAlSi₃O₈ in the solid. This coefficient agrees well with that inferred from plagioclase-sanidine equilibrium phenocryst assemblages in rhyolites. The activities of silica obtained from this model are used to place constraints on the pressure-temperature regions where various types of basic magmas are generated. In conjunction with plagioclase geothermometry an application is given where the pressure, temperature, and water content of an olivine andesite is predicted from the activity of silica.     

Petrogenesis of mafic and associated silicic end-member magmas for calc-alkaline mixed rocks in the Shirataka volcano,NE Japan

Eruptive products of the Shirataka volcano (0.9–0.7 Ma) in NE Japan are calc-alkaline andesite–dacite, and are divisible into six petrologic groups (G1–G6). Shirataka rocks possess mafic inclusions—basalt–basaltic andesite, except for G3 and G4. All rocks show mixing and mingling of the mafic and silicic end-members, with trends defined by hosts and inclusions divided into high-Cr and low-Cr types; both types coexist in G1, G2, and G5. Estimated mafic end-members are high-Cr (1120–1170°C, 48–51% SiO₂, olv ± cpx ± plg) and low-Cr type magmas (49–52% SiO₂, cpx ± plg) except for the Sr isotopic composition. In contrast, the silicic end-members of both types have similar petrologic features (790–840°C, 64–70% SiO₂, hbl ± qtz ± px + plg). High-Cr type mafic and corresponding silicic end-members have lower ⁸⁷Sr/⁸⁶Sr ratios than the low-Cr ones in each group. The trace element model calculations suggest that the low-Cr type mafic end-member magma is produced through ca. 20% fractional crystallization (olv ± cpx ± plg) from the high-Cr type one with assimilation of granitoids (r = 0.02–0.05). The silicic magmas are producible through <30% partial remelting of previously emplaced basaltic magma with assimilation of crustal components. The compositional difference between the low-K and medium-K basalts in the Shirataka volcano is mainly attributed to the different degrees of the effect of subduction derived fluid by dehydration of phlogopite. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Experiment in K/Ar geothermometry and evaluation of thermodynamic parameters in formation of intrusive bodies

Baking of host rocks by a hot dike leads to a decrease in the argon-potassium ratio in biotite and an increase in this ratio in plagioclase. The extent of such redistribution of radiogenic argon and the depths affected by thermal action of the dike accord indications of the temperature of the magma and the duration of its baking action, with a surmise as to the pressure at the time of its crystallization. —IGR Staff     

A simple thermodynamic model for estimating the solubility of H2O in magmas

A symmetrical, strictly regular solution model is used to estimate H₂O solubilities in silicate melts. The standard state chemical potential of dissolved H₂O and the adjustable parameter in the activity coefficient are determined by least squares analyses of data on H₂O solubility in silicate melts. The adjustable parameter in the expression for the activity coefficient (In) is a function only of the anhydrous melt composition and eleven values are provided for melts ranging in composition from picrite to rhyolite. At the 95% confidence level, the model should estimate H₂O contents to within 4.8% of the amount present if the amount present is less than 10 wt.%. This compares to the reproducibility of 2.25% of the amount present for experimental determinations. To apply the model to rocks and magmas estimates ofT, P, and the fugacity of H₂O are required.Variation of the H₂O content of the melt changes the activity of other components. Knowledge of this variation removes the requirement that the fugacity of H₂O be estimated. Application of the properties of exact differentials to the Gibbs function for the hydrous melt provides an expression relating the chemical potential of a feldspar component to the H₂O content of the melt. This expression contains a second adjustable parameter which depends on the anhydrous melt composition. Using this second expression, the H₂O content can be estimated ifT, P, and feldspar composition are known. Data are too meagre to evaluate the quantitative success of the second method.     

Evidence for extreme fractionation of peralkaline silicic magmas, the Boseti volcanic complex, Main Ethiopian Rift

Matrix glass and melt inclusions in phenocrysts from pantellerite lavas of the Boseti volcanic complex, Ethiopia, record extreme fractionation of peralkaline silicic magma, with Al₂O₃ contents as low as 2.3?wt.%, FeO* contents up to 17?wt.% and SiO₂ contents ~65?wt.%. The new data, and published data for natural and experimental glasses, suggest that the effective minimum composition for peralkaline silicic magmas has ~5?wt.% Al₂O₃, 13?wt.% FeO* and 66?±?2?wt.% SiO₂. The dominant fractionating assemblage is alkali feldspar?+?fayalite?+?hedenbergite?+?oxides?±?quartz. Feldspar – melt relationships indicate that the feldspar is close to the minimum on the albite-orthoclase solid solution loop through the entire crystallization history. There is petrographic, mineralogical and geochemical evidence that magma mixing may have been a common process in the Boseti rhyolites.     

东天山似斑状角闪辉长岩类与铁钛氧化物矿床的关系

似斑状角闪辉长岩类是东天山镁铁-超镁铁杂岩的一部分,以含有嵌晶角闪石为特征,其中部分岩体的似斑状角闪辉长岩赋存铁钛氧化物矿床。为探讨东天山似斑状角闪辉长岩类和铁钛氧化物矿床形成的内在成因关系,本次研究对二红洼、牛毛泉和香山西似斑状角闪辉长岩类的斜长石、单斜辉石、斜方辉石以及嵌晶角闪石进行了系统的电子探针成分分析,并对其形成物理化学条件进行了探讨。结果表明,东天山似斑状角闪辉长岩类的母岩浆氧逸度变化范围为NNO-0.37~NNO+1.15,和攀西地区最大的红格岩体的氧逸度一致,侵入于中下地壳,深度介于8.9~15.7km,结晶温度大于914℃。东天山似斑状角闪辉长岩中铁钛氧化物矿床的富集得益于母岩浆的高氧逸度,但铁钛氧化物矿床的形成和规模则受控于各个岩体的岩浆演化程度。二红洼、牛毛泉和香山西的岩浆的演化程度系统增高,演化程度最低的二红洼岩体并未形成铁钛氧化物矿床,而分异演化程度逐渐增高的牛毛泉和香山西似斑状角闪辉长岩类则分别形成了本区一个小型的和一个大型的铁钛氧化物矿床。     

FEG-EPMA mapping and Fe-Ti oxide mineral chemistry of Brahmaputra River sediments in Bangladesh: provenance and petrogenetic implications

The present research deals with the FEG-EPMA mapping and Fe-Ti oxide mineral chemistry of Brahmaputra River sediments in Bangladesh. Major heavy minerals in the sediments consist of garnet (8.5–21.3%), kyanite (5.35–11.9%), monazite (2.3–5.3%), sillimanite (1.8–4.7%), zircon (3.6–9.1%), and a considerable amount of opaques mainly Fe-Ti oxide minerals (23.1–35.4%). The detrital Fe-Ti oxide minerals carry significant clues to the parent rocks or sources. In these contexts, detrital opaques (Fe-Ti oxides) have been analyzed with an electron probe microanalyzer (EPMA). These opaques (Fe-Ti oxide) display six types of textural patterns, dominantly seriate with granular, emulsion, and acicular sandwich structures and trellis type of textural patterns. These textural patterns belong to five intergrowths of Fe-Ti oxide minerals such as (1) ilmenite-hematite, (2) magnetite-ilmenite, (3) hematite-rutile, (4) ilmenite-hematite-rutile, and (5) ilmenite-rutile, where ilmenite-hematite intergrowth is common. Alteration is seen in both exsolved and unexsolved ilmenites. Textural patterns and mineral chemistry of the studied ilmenite minerals provide lines of evidence of low-temperature magmatic inheritance, later modified by diffusional processes. The estimated temperature and oxygen fugacity from the magnetite-ilmenite exsolution range from 547.6 to 558.2 °C and from 10^?21.4 to 10^?21.7, respectively. The data are also consistent with hematite-ilmenite temperature (between 537 and 540 °C) and oxygen fugacity (10^?21.7 to 10^?21.9) measurements in Cox’s Bazar beach placers. These temperatures and oxygen fugacities specify Fe-Ti oxide assemblages equilibrated in a T-fO₂ field very near to the FMQ buffer curve suggesting a crustal source (magmatic and/or metamorphic), which is modified significantly by metamorphic processes.     

Magmatic processes that generate chemically distinct silicic magmas in NW Costa Rica and the evolution of juvenile continental crust in oceanic arcs

Northwestern Costa Rica is built upon an oceanic plateau that has developed chemical and geophysical characteristics of the upper continental crust. A major factor in converting the oceanic plateau to continental crust is the production, evolution, and emplacement of silicic magmas. In Costa Rica, the Caribbean Large Igneous Province (CLIP) forms the overriding plate in the subduction of the Cocos Plate—a process that has occurred for at least the last 25 my. Igneous rocks in Costa Rica older than about 8 Ma have chemical compositions typical of ocean island basalts and intra-oceanic arcs. In contrast, younger igneous deposits contain abundant silicic rocks, which are significantly enriched in SiO₂, alkalis, and light rare-earth elements and are geochemically similar to the average upper continental crust. Geophysical evidence (high Vp seismic velocities) also indicates a relatively thick (~40 km), addition of evolved igneous rocks to the CLIP. The silicic deposits of NW Costa Rica occur in two major compositional groups: a high-Ti and a low-Ti group with no overlap between the two. The major and trace element characteristics of these groups are consistent with these magmas being derived from liquids that were extracted from crystal mushes—either produced by crystallization or by partial melting of plutons near their solidi. In relative terms, the high-Ti silicic liquids were extracted from a hot, dry crystal mush with low oxygen fugacity, where plagioclase and pyroxene were the dominant phases crystallizing, along with lesser amounts of hornblende. In contrast, the low-Ti silicic liquids were extracted from a cool, wet crystal mush with high oxygen fugacity, where plagioclase and amphibole were the dominant phases crystallizing. The hot-dry-reducing magmas dominate the older sequence, but the youngest sequence contains only magmas from the cold-wet-oxidized group. Silicic volcanic deposits from other oceanic arcs (e.g., Izu-Bonin, Marianas) have chemical characteristics distinctly different from continental crust, whereas the NW Costa Rican silicic deposits have chemical characteristics nearly identical to the upper continental crust. The transition in NW Costa Rica from mafic oceanic arc and intra-oceanic magma to felsic, upper continental crust-type magma is governed by a combination of several important factors that may be absent in other arc settings: (1) thermal maturation of the thick Caribbean plateau, (2) regional or local crustal extension, and (3) establishment of an upper crustal reservoir.     

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

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