Calcareous inclusions in lavas and agglomerates of Santorini volcano

Thermally metamorphosed and metasomatised fragments of basement actinolite-chlorite-calcite-quartz schists and quartz-bearing marbles are found as inclusions in Quaternary agglomerates and historic (197 B. C.—1950) dacitic lavas of Santorini volcano, Greece.Inclusions in agglomerates preserve the structure of parent schists in the alternation of bands rich in diopside or salite with bands rich in plagioclase. By contrast, inclusions in historic dacites are not banded. Most develop a thin zone of hybrid material at the contact with enclosing lava. The assemblage calcic clinopyroxene-wollastonite-plagioclase is commonly developed. The clinopyroxene is a Fe³⁺-rich salite or ferrosalite. Andradite-rich garnet and sphene are accessory minerals. Most examples carry interstitial siliceous glass of distinctive chemical composition, and several show minor olivine, augite, hypersthene and calcic plagioclase of magmatic origin.Other inclusions exhibit the assemblage anhydrite-calcic clinopyroxene, the latter mineral ranging widely in Al content. A single example has been observed to develop two distinct assemblages, the first coarsely crystalline melilite-wollastonite-magnetite, the second finely intergrown melilite-wollastonite-andraditic garnet (-xonotlite).Stability data for hedenbergite and andradite as constituents of skarn assemblages suggest that the clinopyroxene-rich assemblages of inclusions in historic dacites formed at temperatures near to or above 800° C and oxygen fugacity (fO₂) considerably greater than that which could be imposed upon the inclusions by dacite magma (T 900° C, fO₂10^–13 atm.). Thermal breakdown of original carbonates of the inclusions probably supplied the necessary oxygen. T-fO₂ data for the reaction 4 Magnetite+18 Wollastonite 6 Andradite indicate that the assemblage melilite-wollastonite-magnetite of the last inclusion described formed at higher T and/or lower fO₂ than the assemblage melilite-wollastonite-garnet. The latter assemblage undoubtedly formed during inclusion of the fragment by dacite magma, while metamorphism by a more basic, high temperature magma may have produced the former. Temperature data for reactions limiting the stability of melilite in the system CaO-Al₂O₃-SiO₂-H₂O indicate a minimum temperature of around 800° C for formation of both assemblages.     

Formation of shoshonites from calcalkaline basalt magmas: geochemical and experimental constraints from the type locality

Independence volcano, Montana is a major center of the Absaroka volcanic field, from which absarokite, shoshonite, and banakite were originally defined. One magmatic trend at Independence volcano, from high-alumina tholeiitic basalt through shoshonite to high-K dacite, may be modeled by fractional crystallization of observed phenocryst phases (plagioclase, hypersthene, augite, and magnetite). Trace-element and Sr and Nd isotopic compositions of rocks are consistent with this model.Compositions of partial melts from experiments on four rocks at 1 atm and at 10 kbar demonstrate that rock compositions represent a nearly-anhydrous liquid line of descent at a pressure much closer to 10 kbar than to 1 atm. The line of descent involves crystallization of orthopyroxene, not olivine, resulting in strong enrichment in K₂O with little increase in SiO₂. Crystallization at either lower pressures or with water present, involving olivine, results in enrichment in both SiO₂ and K₂O.High-pressure (10 kbar) fractional crystallization of basaltic magma, resulting in formation of shoshonites, may occur at the base of thick crust (e.g., in continental interiors or in very mature arcs). At least a portion of the relationship between K₂O content of arc-related magmas and depth to the Benioff Zone may be attributed to thickening of crust towards the back-arc, resulting in higher pressures of fractionation in Moho-level chambers.     

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.     

Mid-Pleistocene lavas from the Seguam volcanic center,central Aleutian arc: closed-system fractional crystallization of a basalt to rhyodacite eruptive suite

In contrast to adjacent volcanic centers of the modern central Aleutian arc, Seguam Island developed on strongly extended arc crust. K-Ar dates indicate that mid-Pleistocene, late-Pleistocene, and Holocene eruptive phases constitute Seguam. This study focuses on the petrology of the mid-Pleistocene, 1.07–07 Ma, Turf Point Formation (TPF) which is dominated by an unusual suite of porphyritic basalt and basaltic andesite lavas with subordinate phenocryst-poor andesite to rhyodacite lavas. Increasing whole-rock FeO^*/MgO from basalt to dacite, the anhydrous Plag+Ol+Cpx±Opx±Mt phenocryst assemblage, groundmass pigeonite, and the reaction Ol+Liq=Opx preserved in the mafic lavas indicate a tholeiitic affinity. Thermometry and comparison to published phase equilibria suggests that most TPF basalts crystallized Plag+Ol+Cpx±Mt at 1160°C between about 3–5 kb (±1–2% H₂O), andesites crystallized Plag+Cpx+Opx±Mt at 1000°C between 3–4 kb with 3–5% H₂O, and dacites crystallized Plag +Cpx±Opx±Mt at 1000°C between 1–2 kb with 2–3% H₂O. All lavas crystallized at f _{o 2} close to the NNO buffer. Mineral compositions and textures indicate equilibrium crystallization of the evolved lavas; petrographic evidence of open-system mixing or assimilation is rare. MgO, CaO, Al₂O₃, Cr, Ni, and Sr abundances decrease and K₂O, Na₂O, Rb, Ba, Zr, and Pb increase with increasing SiO₂ (50–71%). LREE enrichment [(Ce/Yb)_n=1.7±0.2] characterizes most TPF lavas; total REE contents increase and Eu anomalies become more negative with increasing SiO₂. Relative to other Aleutian volcanic centers, TPF basalts and basaltic andesites have lower K₂O, Na₂O, TiO₂, Rb, Ba, Sr, Zr, Y, and LREE abundances. ⁸⁷Sr/⁸⁶Sr ratios (0.70361–0.70375) and ratios of ²⁰⁶Pb/²⁰⁴Pb (18.88–18.97), ²⁰⁷Pb/²⁰⁴Pb (15.58–15.62), ²⁰⁸Pb/²⁰⁴Pb (38.46–38.55) are the highest measured for any suite of lavas in the oceanic portion of the Aleutian arc. Conversely, Nd values (+5.8 to+6.7) are among the lowest from the Aleutians. Sr, Nd, and Pb ratios are virtually constant from basalt through rhyodacite, whereas detectable isotopic heterogenity is observed at most other Aleutian volcanic centers. Major and trace element, REE, and Sr, Nd, and Pb isotopic compositions are consistent with the basaltic andesitic, andesitic, dacitic, and rhyodacitic liquids evolving from TPF basaltic magma via closed-system fractional crystallization alone. Fractionation models suggest that removal of 80 wt% cumulate (61% Plag, 17% Cpx, 12% Opx, 7% Ol, and 3% Mt) can produce 20 wt% rhyodacitic residual liquid per unit mass of parental basaltic liquid. Petrologic and physical constraints favor segregation of small batches of basalt from a larger mid-crustal reservoir trapped below a low-density upper crustal lid. In these small magma batches, the degree of cooling, crystallization, and fractionation are functions of the initial mass of basaltic magma segregated, the thermal state of the upper crust, and the magnitude of extension. Tholeiitic magmas erupted at Seguam evolved by substantially different mechanisms than did calc-alkaline lavas erupted at the adjacent volcanic centers of Kanaga and Adak on unextended arc crust. These variable differentiation mechanisms and liquid lines of descent reflect contrasting thermal and mechanical conditions imposed by the different tectonic environments in which these centers grew. At Seguam, intra-arc extension promoted eruption of voluminous basalt and its differentiates, unmodified by interaction with lower crustal or upper mantle wallrocks.     

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     

Generation of voluminous silicic magmas and formation of mid-Cenozoic crust beneath north-central Mexico: evidence from ignimbrites,associated lavas,deep crustal granulites,and mantle pyroxenites

 Isotopic and trace element data from mantle and granulite xenoliths are used to estimate the relative contributions of mantle and crustal components to a large ignimbrite, referred to as the upper ignimbrite, that is representative of the voluminous mid-Cenozoic rhyolites of northwestern Mexico. The study also uses data from the volcanic rocks to identify deep crustal xenoliths that are samples of new crust created by the Tertiary magmatism. The isotopic composition of the mantle component is defined by mantle-derived pyroxenites that are interpreted to have precipitated from mid-Cenozoic basaltic magmas. This component has ɛ_Nd≈+1.5, ⁸⁷Sr/⁸⁶Sr≈0.7043 and ²⁰⁶Pb/²⁰⁴Pb≈18.6. Within the upper ignimbrite and associated andesitic and dacitic lavas, initial ⁸⁷Sr/⁸⁶Sr is positively correlated with SiO₂, reaching 0.7164 in the ignimbrite. Initial ²⁰⁶Pb/²⁰⁴Pb ratios also show a positive correlation with silica, whereas ɛ_Nd values have a crude negative correlation, reaching values as low as −2. Of the four isotopically distinct crustal components identified from studies of granulite xenoliths, only the sedimentary protolith of the paragneiss xenoliths can be responsible for the high initial ⁸⁷Sr/⁸⁶Sr of the upper ignimbrite. The Nd, Sr, and Pb isotopic compositions of the upper ignimbrite can be modeled with relatively modest assimilation (≤20%) of the sedimentary component ± Proterozoic granulite. Gabbroic composition granulite xenoliths have distinctive Nd, Sr, and Pb isotope ratios that cluster closely within the range of compositions found in the andesitic and dacitic lavas. These mafic granulites are cumulates, and their protoliths are interpreted to have precipitated from the intermediate to silicic magmas at 32–31 Ma. These mafic cumulate rocks are probably representative of much of the deep crust that formed during mid-Cenozoic magmatism in Mexico. Worldwide xenolith studies suggest that the relatively great depth (≤20 km) at which assimilation-fractional crystallization took place in the intermediate to silicic magma systems of the La Olivina region is the rule rather than the exception. Oligocene ignimbrites of the southwestern United States (SWUS) have substantially lower ɛ_Nd values (e.g. <−6) than the upper ignimbrite and other rhyolites from Mexico. This difference appears to reflect a greater crustal contribution to ignimbrites of the SWUS, perhaps due to a higher temperature of the lower crust prior to the emplacement of the Oligocene basaltic magmas. Received: 16 December 1994 / Accepted: 13 September 1995     

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

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

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.     

Origin of silicic rocks of the Deccan Traps continental flood basalt province: Inferences from field observations,petrography, and geochemistry

The origin of silicic rocks (SiO₂ > 65 wt%) in Continental Flood Basalt (CFB) provinces could be attributed to complex petrogenetic processes. The 65.5–66 Ma old Deccan Traps CFB contains eight sporadic but significant silicic rock exposures that are studied here in a comprehensive framework using field observations, petrography, major oxides (n = 56), and trace element chemistry. Rhyolite and granophyre, as well as subordinate felsite, ignimbrite, trachyte, pitchstone, and microgranite coexist with volcanic and plutonic mafic rocks such as basalt, basaltic andesite, and gabbro. Multiple isolated and circular/semi-circular hills and linear dykes of silicic rocks are present in the form of lavas with prominent flow folding, rheomorphic ignimbrite, and tuffs. The ‘Rheological Agpaitic Index’ (RAI) indicates that most of the silicic rocks in the Deccan Traps are effusive in nature, except for Rajpipla, Alech, Bombay, and Osham silicic rocks, which are marked by explosive volcanism. Thermodynamic-based Rhyolite-MELTS modelling suggests that the major oxide composition of Pavagadh and Barda basalt is a likely candidate for the parental melt composition of the silicic rocks of the Deccan Traps. Ba, Sr, P, Zr, and Ti anomalies are consistent with the fractionation of K-feldspar, plagioclase, apatite, zircon, and Fe-Ti oxides, respectively. Two broad REE patterns are noticed in the Deccan Traps silicic rocks: a flat pattern for Barda, Alech, and Chogat-Chamardi silicic rocks, and a steep REE pattern for Osham, Rajula, Pavagadh, Rajpipla, and Bombay silicic rocks. Trace element modelling reveals that 5–10 % partial melting of a spinel peridotite source could produce an REE pattern and abundances similar to the associated basalts. Further extensive fractional crystallization (60–90 %) of the parental mafic melt at a deeper depth (where spinel is stable) could produce the REE composition and pattern observed in most silicic rocks except for those of Barda, Alech, and Chogat-Chamardi, which require fractional crystallization of the same parental melt at a shallower depth (where spinel is not stable). The geochemical variability of Deccan Traps silicic rocks reveals an origin from a mantle-derived parental mafic melt that evolved via the assimilation and fractional crystallization (AFC) process to form the silicic exposures, which is typical of silicic volcanism in other global CFBs.     

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

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

Geochemistry of ultra-potassic rhyodacite magmas from the area of the Orlovka Massif of Li-F granites in Eastern Transbaikalia: Evidence from study of melt inclusions in quartz

Dikes, stocks and/or sheet flows of felsic volcanic and subvolcanic rocks are typically observed in the vicinity of rare-metal Li-F granite massifs. Their ubiquitous spatial association to rare-metal granites and, often, geochemical affinity to them suggest their certain petrological relation. Compositionally unique ultrapotassic trachydacites enriched in many rare elements were found among these rocks within the Khangilay complex of ore deposits in Eastern Transbaikalia. Melt inclusions in rock-forming quartz were studied to reconstruct the composition and evolution of parent melt. The obtained data demonstrated the existence of a super-potassic peraluminous melt (K₂O = 6.12 wt %, Na₂O = 1.08 wt %) having elevated contents of rare lithophile elements (730 ppm Rb₂O and 900 ppm BaO). The ion-microprobe content of Li is 354.23 ppm at a relatively low F content (up to 0.5 wt %). The residual melt is characterized by the most unusual composition: extremely low contents of mafic components and basicity (< 0.5 wt % femic oxides), a high Al index (A/CNK = 1.53) at comparatively low SiO₂ (60 wt %), and high total sodic alkalinity (more than 10 wt % K₂O + Na₂O; 6.11 wt % Na₂O). Such a composition corresponds to ongonite magma. However, the melt contains no F but has a high Cl content (0.34 wt %), which corresponds to the limit Cl saturation of haplogranite melt. SHRIMP-II U-Pb zircon dating showed significant difference between rare metal granites and trachyrhyodacites of the Khangilay complex of ore deposits: 139.9 ± 1.9 Ma and 253.4 ± 2.4 Ma, respectively. The geochemical similarity of these rocks, primarily in terms of abundance of refractory elements, REE distribution patterns, and initial Sr ratio, indicates their derivation from similar protolith.     

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

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

Origin and interaction of mafic and felsic magmas in an evolving late orogenic setting: the Early Paleozoic Terra Nova Intrusive Complex, Antarctica

The Abbott Unit (∼508 Ma) and the Vegetation Unit (∼475 Ma) of the Terra Nova Intrusive Complex (northern Victoria Land, Antarctica) represent the latest magmatic events related to the Early Paleozoic Ross Orogeny. They show different emplacement styles and depths, ranging from forcible at 0.4–0.5 GPa for the Abbott Unit to passive at ∼0.2 GPa for the Vegetation Unit. Both units consist of mafic, felsic and intermediate facies which collectively define continuous chemical trends. The most mafic rocks from both units show different enrichment in trace element and Sr-Nd isotopic signatures. Once the possible effects of upper crustal assimilation-fractional crystallisation (AFC) and lower crustal coupled AFC and magma refilling processes have been taken into account the following features are recognised: (1) the modelled primary Abbott Unit magma shows a slightly enriched incompatible element distribution, similar to common continental arc basalts and (2) the modelled primary Vegetation Unit magma displays highly enriched isotope ratios and incompatible element patterns. We interpreted these major changes in magmatic affinity and emplacement style as linked to a major change in the tectonic setting affecting melt generation, rise and emplacement of the magmas. The Abbott Unit mafic melts were derived from a mantle wedge above a subduction zone, with subcontinental lithospheric mantle marginally involved in the melting column. The Vegetation Unit mafic melts are regarded as products of a different source involving an old layer of subcontinental lithospheric mantle. The crustal evolution of both types of mafic melts is marked by significant compositional contrasts in Sr and Nd isotopes between mafic and associated felsic rocks. The crustal isotope signature showed an increase with felsic character. Geochemical variations for both units can be accounted for by a similar two-stage hybridisation process. In the first stage, the most mafic magma evolved mainly by fractional crystallisation coupled with assimilation of metasedimentary rocks having crustal time-integrated Sr and Nd compositions similar to those of locally exposed metamorphic basement. The second stage involves contaminated products mixing with independently generated crustal melts. Petrographic, geochemical and isotope data also provide evidence of significant compositional differences in the felsic end-members, pointing to the involvement of metaigneous and metasedimentary source rocks for the Abbott granite and Vegetation leucogranite, respectively. Received: 31 March 1998 / Accepted: 3 May 1999     

Fe-Ti oxide geothermometry: thermodynamic formulation and the estimation of intensive variables in silicic magmas

A new thermodynamic formulation of the Fe–Ti oxide geothermometer/oxygen barometer is developed. The method is based upon recently calibrated models for spinel solid solutions in the quinary system (Fe²⁺, Mg)(Al,Fe³⁺,Cr)₂O₄–(Fe²⁺, Mg)₂TiO₄ by Sack and Ghiorso, and rhombohedral oxides in the quaternary system (Fe²⁺,Mg,Mn)TiO₃–Fe₂O₃ (this paper). The formulation is internally consistent with thermodynamic models for (Fe²⁺,Mg)-olivine and -orthopyroxene solid solutions and end-member thermodynamic properties tabulated by Berman. The constituent expressions account for compositional and temperature dependent cation ordering and reproduce miscibility gap features in all of the component binaries. The calibration does not account for the excess Gibbs energy resulting from compositional and temperature dependent magnetic ordering in either phase. This limits application of the method to assemblages that equilibrated at temperatures above 600° C. Practical implementation of the proposed geothermometer/oxygen barometer requires minimal use of projection algorthms in accommodating compositions of naturally occurring phases. The new formulation is applied to the estimation of temperature and oxygen fugacity in a wide variety of intermediate to silicic volcanic rocks. In combination with previous work on olivine and orthopyroxene thermodynamics, equilibration pressures are computed for a subset of these volcanics that contain the assemblage quartz, oxides and either ferromagnesian silicate. The calculated log₁₀ f _{O 2}-T relations are reflected in coexisting ferromagnesian mineral assemblages. Volcanics with the lowest relative oxygen fugacity (log₁₀ f _{O 2}) are characterized by the assemblage olivine-quartz, those with slightly higher log₁₀ f _{O 2} s, by the assemblage orthopyroxene-quartz. The sequence proceeds with the necessary phases biotite-feldspar, then hornblende-quartz-clinopyroxene, and finally at the highest log₁₀ f _{O 2} s, sphene-quartz-clinopyroxene. Quantitative analysis of these trends, utilizing thermodynamic data for the constituent phases, establishes that, in most cases, the T-log₁₀ f _{O 2}value computed from the oxides is consistent with the compositions of coexisting silicate phases, indicating that phenocryst equilibrium was achieved prior to eruption. There is, however, considerable evidence of oxide-silicate disequilibrium in samples collected from more slowly cooled domes and obsidians. In addition, T-log₁₀ f _{O 2}trends from volcanic rocks that contain biotite and orthopyroxene are interpreted to imply a condition of Fe²⁺–Mg exchange disequilibrium between orthopyroxene and coexisting ferromagnesian silicates and melt. It is suspected that many biotite-feldspar-quartz-orthopyroxene bearing low temperature volcanic rocks inherit orthopyroxene xenocrysts which crystallized earlier in the cooling history of the magma body.The problem is probably at least as complex as that of the feldspars... A.F. Buddington (1956)     

The evolution of young silicic lavas at Medicine Lake Volcano,California: Implications for the origin of compositional gaps in calc-alkaline series lavas

At Medicine Lake Volcano, California, the compositional gap between andesite (57–62 wt.% SiO₂) and rhyolite (73–74 wt.% SiO₂) has been generated by fractional crystallization. Assimilation of silicic crust has also occurred along with fractionation. Two varieties of inclusions found in Holocene rhyolite flows, hornblende gabbros and aphyric andesites, provide information on the crystallization path followed by lavas parental to the rhyolite. The hornblende gabbros are magmatic cumulate residues and their mineral assemblages are preserved evidence of the phases that crystallized from an andesitic precursor lava to generate the rhyolite lavas. The andesitic inclusions represent samples of a parental andesite and record the early part of the differentiation history. Olivine, plagioclase and augite crystallization begins the differentiation history, followed by the disappearance of olivine and augite through reaction with the liquid to form orthopyroxene and amphibole. Further crystallization of the assemblage plagioclase, amphibole, orthopyroxene, magnetite, and apatite from a high-SiO₂ andesite leads to rhyolite. This final crystallization process occurs on a cotectic that is nearly horizontal in temperature-composition space. Since a large amount of crystallization occurs over a limited temperature interval, a compositional gap develops between rhyolite and high SiO₂ andesite.Liquidus surfaces with shallow slopes in temperature-composition space are characteristic of several late-stage crystallization assemblages in the andesite to rhyolite compositional range. Experimentally produced plagioclase+ amphibole+orthopyroxene+magnetite and plagioclase+ augite+low-Ca pyroxene+magnetite cotectics have liquidus slopes that are nearly flat. At other calc-alkaline volcanic centers crystallization processes involving large compositional changes over small temperature intervals may also be important in the development of bimodal volcanism (i.e. the existence of a composition gap). At Mt. Mazama and Mt. St. Helens, USA and Aso Caldera and Shikotsu, Japan the amphibole-bearing assemblage was important. At Krakatau, Indonesia and Katmai, USA, an augite+orthopyroxene-bearing assemblage was important. In addition to its role in the production of a compositional gap between intermediate and rhyolitic lavas, the crystallization process increases the H₂O content of the residual liquid. This rapid increase in residual liquid volatile content which results from the precipitation of a large proportion of crystalline solids may be an important factor among several that lead to explosive silicic eruptions.     

Interrelations between coeval mafic and A-type silicic magmas from composite dykes in a bimodal suite of southern Israel, northernmost Arabian–Nubian Shield: Geochemical and isotope constraints

Late Neoproterozoic bimodal dyke suites are abundant in the Arabian–Nubian Shield. In southern Israel this suite includes dominant alkaline quartz porphyry dykes, rare mafic dykes, and numerous composite dykes with felsic interiors and mafic margins. The quartz porphyry chemically corresponds to A-type granite. Composite dykes with either abrupt or gradational contacts between the felsic and mafic rocks bear field, petrographic and chemical evidence for coexistence and mixing of basaltic and rhyolitic magmas. Mixing and formation of hybrid intermediate magmas commenced at depth and continued during emplacement of the dykes. Oxygen isotope ratios of alkali feldspar in quartz porphyry (13 to 15‰) and of plagioclase in trachydolerite (10–11‰) are much higher than their initial magmatic ratios predicted by equilibrium with unaltered quartz (8 to 9‰) and clinopyroxene (5.8‰). The elevation of δ¹⁸O in alkali feldspar and plagioclase, and extensive turbidization and sericitization call for post-magmatic low-temperature (≤ 100 °C) water–rock interaction. Hydrous alteration of alkali feldspar, the major carrier of Rb and Sr in the quartz–porphyry, also accounts for the highly variable and unusually high I(Sr) of 0.71253 to 0.73648.

The initial ¹⁴³Nd/¹⁴⁴Nd ratios, expressed by ε_Nd(T) values, are probably unaltered and show small variation in mafic and felsic rocks within a narrow range from + 1.4 to + 3.3. The Nd isotope signature suggests either a common mantle source for the mafic and silicic magmas or a juvenile crustal source for the felsic rocks (metamorphic rocks from the Elat area). However, oxygen isotope ratios of zircon in quartz porphyry [δ¹⁸O(Zrn) = 6.5 to 7.2‰] reveal significant crustal contribution to the rhyolite magma, suggesting that mafic and A-type silicic magmas are not co-genetic, although coeval. Comparison of ¹⁸O/¹⁶O ratios in zircon allows to distinguish two groups of A-type granites in the region: those with mantle-derived source, δ¹⁸O(Zrn) ranging from 5.5 to 5.8‰ (Timna and Katharina granitoids) and those with major contribution of the modified juvenile crustal component, δ¹⁸O(Zrn) varying from 6.5 to 7.2‰ (Elat quartz porphyry dykes and the Yehoshafat alkaline granite). This suggests that A-type silicic magmas in the northern ANS originated by alternative processes almost coevally.     

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.     

Origin of basaltic magmas in the Mojave Desert area,California

The basaltic lavas erupted throughout the Mojave Desert are basanites (SiO₂<46%, normative nepheline>5%, and K₂O>1.5%), alkali-olivine basalts (SiO₂=46–48%; ne=0–5%; and K₂O=1.0–1.5%), and low-alumina, sub-alkaline basalts (SiO₂=48–51%; ne=0; K₂O<1.0%). One volcano, Pisgah Crater, erupted five times, with lava from each successive phase containing more silica and less potash than the one proceeding it. This compositional trend is the reverse of that expected from differentiation of a single alkalic magma, and therefore, may represent a succession of magmas tapped from a zone of continuing partial melting in the mantle.These lava compositions suggest that first melting was under high water pressure and was followed by relatively dry partial melting of gamet-orthopyroxene-clinopyroxene-olivine assemblages. The successive increase in silica and alkali decrease also requires that the partial melting zone move to shallower levels.All lavas sampled in the Mojave Desert area have compositions that can best be explained by the extraction of magma from such a rising melting zone, analogous to the mantle diapirs suggested by Green and Ringwood.     

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.