Effusive eruption of viscous silicic magma triggered and driven by recharge: a case study of the Cerro Chascon-Runtu Jarita Dome Complex in Southwest Bolivia

 The Cerro Chascon-Runtu Jarita Complex is a group of ten Late Pleistocene (∼85 ka) lava domes located in the Andean Central Volcanic Zone of Bolivia. These domes display considerable macroscopic and microscopic evidence of magma mixing. Two groups of domes are defined chemically and geographically. A northern group, the Chascon, consists of four lava bodies of dominantly rhyodacite composition. These bodies contain 43–48% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole in a microlite-poor, rhyolitic glass. Rare mafic enclaves and selvages are present. Mineral equilibria yield temperatures from 640 to 750  °C and log ƒO₂ of –16. Geochemical data indicate that the pre-eruption magma chamber was zoned from a dominant volume of 68% to minor amounts of 76% SiO₂. This zonation is best explained by fractional crystallization and some mixing between rhyodacite and more evolved compositions. The mafic enclaves represent magma that intruded but did not chemically interact much with the evolved magmas. A southern group, the Runtu Jarita, is a linear chain of six small domes (<1 km³ total volume) that probably is the surface expression of a dike. The five most northerly domes are composites of dacitic and rhyolitic compositions. The southernmost dome is dominantly rhyolite with rare mafic enclaves. The composite domes have lower flanks of porphyritic dacite with ∼35 vol.% phenocrysts of plagioclase, orthopyroxene, and hornblende in a microlite-rich, rhyodacitic glass. Sieve-textured plagioclase, mixed populations of disequilibrium plagioclase compositions, xenocrystic quartz, and sanidine with ternary composition reaction rims indicate that the dacite is a hybrid. The central cores of the composite domes are rhyolitic and contain up to 48 vol.% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole. This is separated from the dacitic flanks by a banded zone of mingled lava. Macroscopic, microscopic, and petrologic evidence suggest scavenging of phenocrysts from the silicic lava. Mineral equilibria yield temperatures of 625–727  °C and log ƒO₂ of –16 for the rhyolite and 926–1000  °C and log ƒO₂ of –9.5 for the dacite. The rhyolite is zoned from 73 to 76% SiO₂, and fractionation within the rhyolite composition produced this variation. Most of the 63–73% SiO₂ compositional range of the lava in this group is the result of mixing between the hybrid dacite and the rhyolite. Eruption of both groups of lavas apparently was triggered by mafic recharge. A paucity of explosive activity suggests that volatile and thermal exchanges between reservoir and recharge magmas were less important than volume increase and the lubricating effects of recharge by mafic magmas. For the Runtu Jarita group, the eruption is best explained by intrusion of a dike of dacite into a chamber of crystal-rich rhyolite close to its solidus. The rhyolite was encapsulated and transported to the surface by the less-viscous dacite magma, which also acted as a lubricant. Simultaneous effusion of the lavas produced the composite domes, and their zonation reflects the subsurface zonation. The role of recharge by hotter, more fluid mafic magma appears to be critical to the eruption of some highly viscous silicic magmas. Received: 23 August 1998 / Accepted: 10 March 1999     

Sr–Nd isotopic and chemical characteristics of the silicic magma reservoir of the Aira pyroclastic eruption, southern Kyushu, Japan

Sr and Nd isotope and geochemical investigations were performed on a remarkably homogeneous, high-silica rhyolite magma reservoir of the Aira pyroclastic eruption (22,000 years ago), southern Kyushu, Japan. The Aira caldera was formed by this eruption with four flow units (Osumi pumice fall, Tsumaya pryoclastic flow, Kamewarizaka breccia and Ito pyroclastic flow). Quite narrow chemical compositions (e.g., 74.0–76.5 wt% of SiO₂) and Sr and Nd isotopic values (⁸⁷Sr/⁸⁶Sr=0.70584–0.70599 and _Nd=−5.62 to −4.10) were detected for silicic pumices from the four units, with the exception of minor amounts of dark pumices in the units. The high Sr isotope ratios (0.7065–0.7076) for the dark pumices clearly suggest a different origin from the silicic pumices. Andesite to basalt lavas in pre-caldera (0.37–0.93 Ma) and post-caldera (historical) eruptions show lower ⁸⁷Sr/⁸⁶Sr (0.70465–0.70540) and higher _Nd (−1.03 to +0.96) values than those of the Aira silicic and dark pumices. Both andesites of pre- and post-caldera stages are very similar in major- and trace-element characteristics and isotope ratios, suggesting that the both andesites had a same source and experienced the same process of magma generation (magma mixing between basaltic and dacitic magmas). Elemental and isotopic signatures deny direct genetic relationships between the Aira pumices and pre- and post-caldera lavas. Relatively upper levels of crust (middle–upper crust) are assumed to have been involved for magma generation for the Aira silicic and dark pumices. The Aira silicic magma was derived by partial melting of a separate crust which had homogeneous chemistry and limited isotope compositions, while the magma for the Aira dark pumice was generated by AFC mixing process between the basement sedimentary rocks and basaltic parental magma, or by partial melting of crustal materials which underlay the basement sediments. The silicic magma did not occupy an upper part of a large magma body with strong compositional zonation, but formed an independent magma body within the crust. The input and mixing of the magma for dark pumices to the base of the Aira silicic magma reservoir might trigger the eruptions in the upper part of the magma body and could produce a slight Sr isotope gradient in the reservoir. An extremely high thermal structure within the crust, which was caused by the uprise and accumulation of the basaltic magma, is presumed to have formed the large volume of silicic magma of the Aira stage.     

Koolau shield basalt as xenoliths entrained during rejuvenated-stage eruptions: perspectives on magma mixing

Rejuvenated-stage tuff cones (Honolulu Volcanics) on Koolau volcano, Oahu, Hawaii, contain xenoliths of Koolau shield basalt. Because Koolau subaerial shield lavas represent a Hawaiian geochemical 'end member', and submarine shield lavas have compositions with some affinities to Mauna Loa and Kilauea, we analyzed 28 xenolithic basalts from Salt Lake and Koko Head cones to determine how these seemingly random samplings of the Koolau profile compare to established Koolau geochemistry. Analyses reveal that 24 are shield tholeiitic basalt—the focus of this study—and 4 are rejuvenated-stage basaltic rocks. The tholeiitic xenoliths represent largely upper Koolau shield lavas, as these samples (8.3 to 5.8 wt% MgO) have, with one exception, overall major- and trace-element compositions that overlap those of Koolau subaerial shield lavas. Secondary processes, however, created some distinctions—namely, enrichments/depletions in K, Ba, Sr, SiO₂, and FeO, and, due to zeolitization (chabazite with attending okenite and apophyllite), elevated CaO. One xenolithic basalt with 8.2 wt% MgO has higher Ti, Zr, Nb, and Sc, and lower Zr/Nb than subaerial lavas, and appears to represent relatively early, deeper shield—thereby reinforcing that the Koolau shield source varied temporally. Olivine, orthopyroxene, and plagioclase are the phenocrysts (clinopyroxene is rare), and their core compositions range widely across the suite—Fo_87.8–72, orthopyroxene Mg#s 85–72, and An_74–60. Several xenolithic basalts have both normally and reversely zoned orthopyroxene and plagioclase with a variety of core compositions (e.g., orthopyroxene-core Mg#s 82, 77, and 72, all in one sample). These compositions and zonations record evidence for wide compositional ranges of replenishment (MgO ~13–8 wt%) and reservoir (MgO ~7 to <5 wt%) magmas mixing in varying proportions; however, extreme MgO lavas (~13 and <5 wt%) are not observed as either subaerial or xenolithic basalt, but are indicated by phenocryst cores of Fo_87.8 and orthopyroxene-Mg# 72. The Koolau magma-mixing history resembles that of Kilauea, and is unlike the 'steady-state' mixing known for Mauna Loa. Finally, these basalt samples show that any xenolithic occurrence of Koolau lava is subject to the zeolitization prevalent in the tuff-cone hosts.Editorial handling: M. Carroll     

New multibeam mapping and geochemistry of the 30°–35° S sector,and overview,of southern Kermadec arc volcanism

New multibeam mapping and whole-rock geochemistry establish the first order definition of the modern submarine Kermadec arc between 30° and 35° S. Twenty-two volcanoes with basal diameters > 5 km are newly discovered or fully-mapped for the first time; Giggenbach, Macauley, Havre, Haungaroa, Kuiwai, Ngatoroirangi, Sonne, Kibblewhite and Yokosuka. For each large volcano, edifice morphology and structure, surficial deposits, lava fields, distribution of sector collapses, and lava compositions are determined. Macauley and Havre are large silicic intra-oceanic caldera complexes. For both, concentric ridges on the outer flanks are interpreted as recording mega-bedforms associated with pyroclastic density flows and edifice foundering. Other stratovolcanoes reveal complex histories, with repeated cycles of tectonically controlled construction and sector collapse, extensive basaltic flow fields, and the development of summit craters and/or small nested calderas.Combined with existing data for the southernmost arc segment, we provide an overview of the spatial distribution and magmatic heterogeneity along ∼780 km of the Kermadec arc at 30°–36°30′ S. Coincident changes in arc elevation and lava composition define three volcano–tectonic segments. A central deeper segment at 32°20′–34°10′ S has basement elevations of > 3200 m water-depth, and relatively simple stratovolcanoes dominated by low-K series, basalt–basaltic andesite. In contrast, the adjoining arc segments have higher basement elevations (typically < 2500 m water-depth), multi-vent volcanic centres including caldera complexes, and erupt sub-equal proportions of dacite and basalt–basaltic andesite. The association of silicic magmas with higher basement elevations (and hence thicker crust), coupled with significant inter- and intra-volcano heterogeneity of the silicic lavas, but not the mafic lavas, is interpreted as evidence for dehydration melting of the sub-arc crust. Conversely, the crust beneath the deeper arc segments is thinner, initially cooler, and has not yet reached the thermal requirements for anatexis. Silicic calderas with diameters > 3 km coincide with the shallower arc segments. The dominant mode of large caldera formation is interpreted as mass-discharge pyroclastic eruption with syn-eruptive collapse. Hence, the shallower arc segments are characterized by both the generation of volatile-enriched magmas from crustal melting and a reduced hydrostatic load, allowing magma vesiculation and fragmentation to initiate and sustain pyroclastic eruptions. Proposed initiation parameters for submarine pyroclastic eruptions are water-depths < 1000 m, magmas with 5–6 wt.% water and > 70 wt.% SiO₂, and a high discharge rate.     

Changes in magma composition at Arenal volcano,Costa Rica, 1968–1985: Real-time monitoring of open-system differentiation

Arenal volcano in Costa Rica has been erupting nearly continuously, but at a diminishing rate, since 1968, producing approximately 0.35 km³ of lavas and tephras that have shown consistent variations in chemistry and mineralogy. From the beginning of the eruption in July 1968 to early 1970 (stage 1, vol.=0.12 km³) tephras and lavas became richer in Ca, Mg, Ni, Cr, Fe, Ti, V, and Sc and poorer in Al₂O₃ and SiO₂. Concentrations of incompatible trace elements (including Sr) decreased by 5%–20%. Phenocryst contents increased 20–50 vol%. During stage 2 (1970–1973, vol. = 0.13 km³) concentrations of compatible trace elements rose, and concentrations of incompatible trace elements either remained constant or also rose. Al₂O₃ contents decreased by 1 wt%. Phenocryst content increased slightly, principally due to increased orthopyroxene. During stage 3 (mid-1974 to the present, vol.= 0.10 km³) concentrations of SiO₂ increased by 1 wt%, compatible trace elements decreased slightly, and incompatible trace element concentrations increased by 5% to 10%. Although crystals increased in size during stage 3, their overall abundance stayed roughly constant.Our modeling suggests that early stage-1 magmas were produced by boundary layer fractionation under high-p H₂O conditions of an unseen basaltic andesitic magma that intruded into the Arenal system after approximately 500 B.P. Changes in composition during stage 2 resulted from mixing of this more mafic original magma with new magma that had a similar SiO₂ content, but higher compatible and incompatible element concentrations. The changes during stage 3 resulted from continued influx of the same magma plus crystal removal.We conclude that the eruption proceeded in the following way. Before 1968 zoned stage-1 magma resided in the deep crust below Arenal. A new magma intruded into this chamber in July 1968 causing ejection of the stage-1 magmas. The intruding magma mixed with mafic portions of the original chamber producing the mixed lavas of stage 2. Continued mixing plus crystal fractionation along the chamber and conduit walls produced stage-3 lavas. The time scales of crustal level magmatic processes at Arenal range 10⁰–10³ years, which are 3–6 orders of magnitude shorter than those of larger, more silicic systems.     

Eruptive history of Mount Mazama and Crater Lake Caldera, Cascade Range, U.S.A.

New investigations of the geology of Crater Lake National Park necessitate a reinterpretation of the eruptive history of Mount Mazama and of the formation of Crater Lake caldera. Mount Mazama consisted of a glaciated complex of overlapping shields and stratovolcanoes, each of which was probably active for a comparatively short interval. All the Mazama magmas apparently evolved within thermally and compositionally zoned crustal magma reservoirs, which reached their maximum volume and degree of differentiation in the climactic magma chamber 7000 yr B.P.The history displayed in the caldera walls begins with construction of the andesitic Phantom Cone 400,000 yr B.P. Subsequently, at least 6 major centers erupted combinations of mafic andesite, andesite, or dacite before initiation of the Wisconsin Glaciation 75,000 yr B.P. Eruption of andesitic and dacitic lavas from 5 or more discrete centers, as well as an episode of dacitic pyroclastic activity, occurred until 50,000 yr B.P.; by that time, intermediate lava had been erupted at several short-lived vents. Concurrently, and probably during much of the Pleistocene, basaltic to mafic andesitic monogenetic vents built cinder cones and erupted local lava flows low on the flanks of Mount Mazama. Basaltic magma from one of these vents, Forgotten Crater, intercepted the margin of the zoned intermediate to silicic magmatic system and caused eruption of commingled andesitic and dacitic lava along a radial trend sometime between 22,000 and 30,000 yr B.P. Dacitic deposits between 22,000 and 50,000 yr old appear to record emplacement of domes high on the south slope. A line of silicic domes that may be between 22,000 and 30,000 yr old, northeast of and radial to the caldera, and a single dome on the north wall were probably fed by the same developing magma chamber as the dacitic lavas of the Forgotten Crater complex. The dacitic Palisade flow on the northeast wall is 25,000 yr old. These relatively silicic lavas commonly contain traces of hornblende and record early stages in the development of the climatic magma chamber.Some 15,000 to 40,000 yr were apparently needed for development of the climactic magma chamber, which had begun to leak rhyodacitic magma by 7015 ± 45 yr B.P. Four rhyodacitic lava flows and associated tephras were emplaced from an arcuate array of vents north of the summit of Mount Mazama, during a period of 200 yr before the climactic eruption. The climactic eruption began 6845 ± 50 yr B.P. with voluminous airfall deposition from a high column, perhaps because ejection of 4−12 km³ of magma to form the lava flows and tephras depressurized the top of the system to the point where vesiculation at depth could sustain a Plinian column. Ejecta of this phase issued from a single vent north of the main Mazama edifice but within the area in which the caldera later formed. The Wineglass Welded Tuff of Williams (1942) is the proximal featheredge of thicker ash-flow deposits downslope to the north, northeast, and east of Mount Mazama and was deposited during the single-vent phase, after collapse of the high column, by ash flows that followed topographic depressions. Approximately 30 km³ of rhyodacitic magma were expelled before collapse of the roof of the magma chamber and inception of caldera formation ended the single-vent phase. Ash flows of the ensuing ring-vent phase erupted from multiple vents as the caldera collapsed. These ash flows surmounted virtually all topographic barriers, caused significant erosion, and produced voluminous deposits zoned from rhyodacite to mafic andesite. The entire climactic eruption and caldera formation were over before the youngest rhyodacitic lava flow had cooled completely, because all the climactic deposits are cut by fumaroles that originated within the underlying lava, and part of the flow oozed down the caldera wall.A total of 51−59 km³ of magma was ejected in the precursory and climactic eruptions, and 40−52 km³ of Mount Mazama was lost by caldera formation. The spectacular compositional zonation shown by the climactic ejecta — rhyodacite followed by subordinate andesite and mafic andesite — reflects partial emptying of a zoned system, halted when the crystal-rich magma became too viscous for explosive fragmentation. This zonation was probably brought about by convective separation of low-density, evolved magma from underlying mafic magma. Confinement of postclimactic eruptive activity to the caldera attests to continuing existence of the Mazama magmatic system.     

Generation of Icelandic rhyolites: silicic lavas from the Torfajökull central volcano

The Torfajökull central volcano in south-central Iceland contains the largest volume of exposed silicic extrusives in Iceland (225 km³). Within SW-Torfajökull, postglacial mildly alkalic to peralkalic silicic lavas and lava domes (67–74 wt.% SiO₂) have erupted from a family of fissures 1–2.5 km apart within or just outside a large caldera (12×18 km). The silicic lavas show a fissure-dependent variation in composition, and form five chemically distinct units. The lavas are of low crystallinity (0–7 vol.%) and contain phenocrysts in the following order of decreasing abundance: plagioclase (An_10-40), Na-rich anorthoclase (<Or₂₃), clinopyroxene (Fs_37-20), FeTi oxides (Usp_32-60; Ilm_93-88), hornblende (edenitic–ferroedenitic) and olivine (Fo_22-37), with apatite, pyrrhotite and zircon as accessory phases. The phenocryst assemblage (0.2–4.0 mm) consistently exhibits pervasive disequilibrium with the host melt (glass). Xenoliths include sparse, disaggregated, and partially fused leucocratic fragments as well as amphibole-bearing rocks of broadly intermediate composition. The values of the silicic lavas are in the range 3.6–4.4, and these are lower than the values of comagmatic, contemporaneous basaltic extrusives within SW-Torfajökull, implying that the former can not be derived from the latter by simple fractional crystallization. FeTi-oxide geothermometry reveals temperatures as low as 750–800°C. To explain the fissure-dependent chemical variations, depletions, low FeTi-oxide temperatures and pervasive crystal-melt disequilibrium, we propose the extraction and collection of small parcels of silicic melt from originally heterogeneous basaltic crustal rock through heterogeneous melting and wall rock collapse (solidification front instability, SFI). The original compositional heterogeneity of the source rock is due to (1) silicic segregations, in the form of pods and lenses characteristically formed in the upper parts of gabbroic intrusives, and (2) extreme isostatic subsidence of the earlier, less differentiated lavas of the Torfajökull central volcano. Ridge migration into older crustal terranes, coupled with establishment of concentrated volcanism at central volcanoes like Torfajökull due to propagating regional fissure swarms, supplies the heat source for this overall process. Continued magmatism in these fissures promotes extensive prograde heating of older crust and the progressive vitality and rise of the central volcano magmatic system that leads to, respectively, SFI and subsidence melting. The ensuing silicic melts (with relict crystals) are extracted, collected and extruded before reaching complete internal equilibrium. Chemically, this appears as a two-stage process of crystal fractionation. In general, the accumulation of high-temperature basaltic magmas at shallow depths beneath the Icelandic rift zones and major central volcanoes, coupled with unique tectonic conditions, allows large-scale reprocessing and recycling of the low- , hydrothermally altered Icelandic crust. The end result is a compositionally bimodal proto-continental crust.     

Petrology and Fe–Ti oxide reequilibration of the 1991 Mount Unzen mixed magma

A dacitic magma (64.5 wt.% SiO₂), a mixture of phenocryst-rich rhyodacite and an aphyric mafic magma, was erupted during the recent 1991–1995 Mount Unzen eruptive cycle. The experimental and analytical results of this study reveal additional details about conditions in the premixing and postmixing magmas, and the nature of the mixing process. The preeruption rhyodacitic magma was at a temperature of 790±20°C according to Fe–Ti oxide phenocryst cores, and at a depth of 6 to 7 km (160 MPa) according to Al-in-hornblende geobarometry. The mafic magma that mixed with the rhyodacite is found as andesitic (54 to 62 wt.% SiO₂) enclaves in the erupted magma and was essentially aphyric when intruded. Phase equilibria indicate that an aphyric andesite at 160 MPa is >1030°C (H₂O-saturated) and possibly as high as 1130°C (2 wt.% H₂O). The composition of the rhyodacite which was mixed with the andesite is estimated to lie between 67 and 69 wt.% SiO₂. Using these compositions and temperatures, the temperature of the Unzen magma after mixing is estimated to be at least 850° to 870°C. The groundmass Fe–Ti oxide microphenocrysts and those in pargasite-bearing reaction zones around biotite phenocrysts both give 890±20°C temperatures; the oxide–oxide contacts give temperatures of 910±20°C. The 900±30°C postmixing temperatures are consistent with phase-equilibria experiments which show that the magma was not above 930°C at 160 MPa. Our Fe–Ti oxide reequilibration experiments suggest that the mixing of the two magmas began within a few weeks of the eruption, which is a shorter time than is calculated using available diffusion data. There is also evidence that some mixing took place much closer to the time of extrusion based on the presence of unrimmed biotite phenocrysts in the magma.     

Differentiation and magma mixing on Kilauea's east rift zone

The lavas of the 1955 east rift eruption of Kilauea Volcano have been the object of considerable petrologic interest for two reasons. First, the early 1955 lavas are among the most differentiated ever erupted at Kilauea, and second, the petrographic character and chemical composition of the lava being erupted changed significantly during the eruption. This shift, from more differentiated (MgO=5.0–5.7%) to more magnesian (MgO=6.2–6.8%) lava, has been variously interpreted, as either due to systematic excavation of a zoned, differentiated magma body, or to invasion of the differentiated magma by more primitive magma, followed by rapid mixing and eruption of the resulting hybrid magmas. Petrologic examination of several nearvent spatter samples of the late 1955 lavas shows abundant evidence for magma mixing, including resorbed and/or reversely zoned crystals of olivine, augite and plagioclase. In addition, the compositional ranges of olivine, plagioclase and groundmass sulfide are very large, implying that the assemblages are hybrid. Core compositions of olivine phenocrysts range from Fo85 to Fo77. The most magnesian olivines in these samples must have originally crystallized from a melt containing 8.0–8.5% MgO, which is distinctly more magnesian than the bulk composition of the late 1955 lavas. The majorelement and trace-element data are either permissive or supportive of a hybrid origin for the late 1955 lavas. In particular, the compositional trends of the 1955 lavas on plots of CaO vs MgO, and the virtual invariance of Al₂O₃ and Sr in these plagioclase-phyric lavas are more easily explained by magma mixing than by fractionation. The pattern of internal disequilibrium/re-equilibration in the late 1955 spatter samples is consistent with reintrusion and mixing having occurred at least twice, during the latter part of the 1955 eruption. Plagioclase zonation preserves possible evidence for additional, earlier reintrusion events. Least-squares modelling the mixing of early 1955 bulk compositions with various summit lavas±olivine pick the 1952 summit lava as most like the primitive component. The results also indicate the primitive component had MgO=7.5–8.0%, corresponding to liquidus temperatures of 1165–1175°C. The absence of Fe-Ti oxide phenocrysts in the late 1955 lavas implies that the cooler component of the hybrid had T>1110°C. Thus the thermal contrast between the two components may have been as much as 55–65°C, sufficient to produce the conspicuous disequilibrium effects visible in the spatter samples.     

Genesis of post-hotspot,A-type rhyolite of the Eastern Snake River Plain volcanic field by extreme fractional crystallization of olivine tholeiite

Rhyolites occur as a subordinate component of the basalt-dominated Eastern Snake River Plain volcanic field. The basalt-dominated volcanic field spatially overlaps and post-dates voluminous late Miocene to Pliocene rhyolites of the Yellowstone–Snake River Plain hotspot track. In some areas the basalt lavas are intruded, interlayered or overlain by ~15 km³ of cryptodomes, domes and flows of high-silica rhyolite. These post-hotspot rhyolites have distinctive A-type geochemical signatures including high whole-rock FeO_tot/(FeO_tot+MgO), high Rb/Sr, low Sr (0.5–10 ppm) and are either aphyric, or contain an anhydrous phenocryst assemblage of sodic sanidine ± plagioclase + quartz > fayalite + ferroaugite > magnetite > ilmenite + accessory zircon + apatite + chevkinite. Nd- and Sr-isotopic compositions overlap with coeval olivine tholeiites (ɛ_Nd = −4 to −6; ⁸⁷Sr/⁸⁶Sr_i = 0.7080–0.7102) and contrast markedly with isotopically evolved Archean country rocks. In at least two cases, the rhyolite lavas occur as cogenetic parts of compositionally zoned (~55–75% SiO₂) shield volcanoes. Both consist dominantly of intermediate composition lavas and have cumulative volumes of several 10’s of km³ each. They exhibit two distinct, systematic and continuous types of compositional trends: (1) At Cedar Butte (0.4 Ma) the volcanic rocks are characterized by prominent curvilinear patterns of whole-rock chemical covariation. Whole-rock compositions correlate systematically with changes in phenocryst compositions and assemblages. (2) At Unnamed Butte (1.4 Ma) the lavas are dominated by linear patterns of whole-rock chemical covariation, disequilibrium phenocryst assemblages, and magmatic enclaves. Intermediate compositions in this group resulted from variable amounts of mixing and hybridization of olivine tholeiite and rhyolite parent magmas. Interestingly, models of rhyolite genesis that involve large degrees of melting of Archean crust or previously consolidated mafic or silicic Tertiary intrusions do not produce observed ranges of Nd- and Sr-isotopes, extreme depletions in Sr-concentration, and cogenetic spectra of intermediate rock compositions for both groups. Instead, least-squares mass-balance, energy-constrained assimilation and fractional crystallization modeling, and mineral thermobarometry can explain rhyolite production by 77% low-pressure fractional crystallization of a basaltic trachyandesite parent magma (~55% SiO₂), accompanied by minor (0.03–7%) assimilation of Archean upper crust. We present a physical model that links the rhyolites and parental intermediate magmas to primitive olivine tholeiite by fractional crystallization. Assimilation, recharge, mixing and fractional melting occur to limited degrees, but are not essential parts of the rhyolite formation process. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. This paper constitutes part of a special issue dedicated to Bill Bonnichsen on the petrogenesis and volcanology of anorogenic rhyolites.     

Evolution of the crustal magma plumbing system during the build-up to the 22-ka caldera-forming eruption of Santorini (Greece)

The formation of shallow, caldera-sized reservoirs of crystal-poor silicic magma requires the generation of large volumes of silicic melt, followed by the segregation of that melt and its accumulation in the upper crust. The 21.8?±?0.4-ka Cape Riva eruption of Santorini discharged >10 km³ of crystal-poor dacitic magma, along with <<1 km³ of hybrid andesite, and collapsed a pre-existing lava shield. We have carried out a field, petrological, chemical, and high-resolution ⁴⁰Ar/³⁹Ar chronological study of a sequence of lavas discharged prior to the Cape Riva eruption to constrain the crustal residence time of the Cape Riva magma reservoir. The lavas were erupted between 39 and 25 ka, forming a ～2-km³ complex of dacitic flows, coulées and domes up to 200 m thick (Therasia dome complex). The Therasia dacites show little chemical variation with time, suggesting derivation from one or more thermally buffered reservoirs. Minor pyroclastic layers occur intercalated within the lava succession, particularly near the top. A prominent pumice fall deposit correlates with the 26-ka Y-4 ash layer found in deep-sea sediments SE of Santorini. One of the last Therasia lavas to be discharged was a hybrid andesite formed by the mixing of dacite and basalt. The Cape Riva eruption occurred no more than 2,800?±?1,400 years after the final Therasia activity. The Cape Riva dacite is similar in major element composition to the Therasia dacites, but is poorer in K and most incompatible trace elements (e.g. Rb, Zr and LREE). The same chemical differences are observed between the Cape Riva and Therasia hybrid andesites, and between the calculated basaltic mixing end-members of each series. The Therasia and Cape Riva dacites are distinct silicic magma batches and are not related by shallow processes of crystal fractionation or assimilation. The Therasia lavas were therefore not simply precursory leaks from the growing Cape Riva magma reservoir. The change 21.8 ky ago from a magma series richer in incompatible elements to one poorer in those elements is one step in the well documented decrease with time of incompatibles in Santorini magmas over the last 530 ky. The two dacitic magma batches are interpreted to have been emplaced sequentially into the upper crust beneath the summit of the volcano, the first (Therasia) then being partially, or wholly, flushed out by the arrival of the second (Cape Riva). This constrains the upper-crustal residence time of the Cape Riva reservoir to less than 2,800?±?1,400 years, and the associated time-averaged magma accumulation rate to >0.004 km³ year^-1. Rapid ascent and accumulation of the Cape Riva dacite may have been caused by an increased flux of mantle-derived basalt into the crust, explaining the occurrence of hybrid andesites (formed by the mixing of olivine basalt and dacite in approximately equal proportions) in the Cape Riva and late Therasia products. Pressurisation of the upper crustal plumbing system by sustained, high-flux injection of dacite and basalt may have triggered the transition from prolonged, largely effusive activity to explosive eruption and caldera collapse.     

The geology of Mount Hasan stratovolcano, central Anatolia, Turkey

Mount Hasan is a double-peaked stratovolcano, located in Central Anatolia, Turkey. The magmas erupted from this multi-caldera complex range from basalt to rhyolite, but are dominated by andesite and dacite. Two terminal cones (Big Mt. Hasan and Small Mt. Hasan) culminate at 3253 m and 3069 m respectively. There are four evolutionary stages in the history of the volcanic complex (stage 1: Kecikalesi volcano, 13 Ma, stage 2: Palaeovolcano, 7 Ma, stage 3: Mesovolcano and stage 4: Neovolcano). The eruptive products consist of lava flows, lava domes, and pyroclastic rocks. The later include ignimbrites, phreatomagmatic intrusive breccias and nuées ardentes, sometimes reworked as lahars. The total volume is estimated to be 354 km³, the area extent 760 km². Textural and mineralogical data suggest that both magma mixing and fractional crystallization were involved in the generation of the andesites and dacites. The magmas erupted from the central volcanoes show a transition with time from tholeite to calc-alkaline. Three generations of basaltic strombolian cones and lava flows were emplaced contemporaneously with the central volcanoes. The corresponding lavas are alkaline with a sodic tendency.     

Replenishment of a mafic magma in a zoned felsic magma chamber beneath Rishiri Volcano, Japan

The trachytic Tanetomi lava from Rishiri Volcano, northern Japan, provides useful information concerning how a replenished mafic magma mixes with a compositionally zoned felsic magma in a magma chamber. The Tanetomi lava was erupted in the order of Lower lava 1 (LL1, 59.2-59.8 wt.% in SiO₂), Lower lava 2 (LL2, 58.4-59.1 wt.%), and Upper lava (UL, 59.9-65.1 wt.%). Evidence for mixing with a mafic magma is observed only in the LL2, in which a greater amount of crystals derived from the mafic magma occurs in rocks with higher SiO₂ content. The whole-rock compositional trend of the Tanetomi lavas is fairly smooth except for the LL2 lava composition, which scatter along the main composition trend. There is no reasonable composition of basaltic magma on the extrapolation of the LL2 composition trend, and the trend cannot be explained by a simple two-component magma mixing. Before the replenishment, the felsic magma was zoned in composition (58-65 wt.% in SiO₂) and temperature (1030-920°C) in the magma chamber located at the pressure of ~2 kbar. The compositional variation of the main felsic magma was produced by extraction of a fractionated interstitial melt from mush zones along the chamber walls and its subsequent mixing with the main magma (boundary layer fractionation). The LL1 magma tapped the magma chamber soon after the replenishment, before the mafic magma mixed with the overall felsic magma. Then the basalt magma mixed heterogeneously with the upper part of the felsic magma by forced convection as a fountain during injection. The mixing of the basalt magma with compositionally zoned felsic magma resulted in the characteristic composition trend of the LL2. The fraction of basaltic magma in the LL2 magma is estimated to be at most 10%. Despite such a small proportion, the basalt magma was mixed completely with the felsic magma, probably because the crystallinity of undercooled basalt magma was low enough to behave as a liquid.     

Intracrustal origin of Arenal basaltic andesite in the light of solid–melt interactions and related compositional buffering

The origin of Arenal basaltic andesite can be explained in terms of fractional crystallization of a parental high-alumina basalt (HAB), which assimilates crustal rocks during its storage, ascent and evolution. Contamination of this melt by Tertiary calc-alkalic intrusives (quartz–diorite and granite, with ⁸⁷Sr/⁸⁶Sr ratios ranging 0.70381–0.70397, nearly identical with those of the Arenal lavas) occurs at upper crustal levels, following the interaction of ascending basaltic magma masses with gabbroic–anorthositic layers. Fragments of these layers are found as inclusions within Arenal lavas and tephra and may show reaction rims (1–5 mm thick, consisting of augite, hypersthene, bytownitic–anorthitic plagioclase, and granular titanomagnetite) at the gabbro–lava interface. These reaction rims indicate that complete `assimilation' was prevented since the temperature of the host basaltic magma was not high enough to melt the gabbroic materials (whose mineral phases are nearly identical to the early formed liquidus phases in the differentiating HAB). Olivine gabbros crystallized at pressure of about 5–6 kbar and equilibrated with the parental HAB at pressures of 3–6 kbar (both under anhydrous and hydrous conditions), and temperatures ranging 1000–1100°C. In particular, `deeper' interactions between the mafic inclusions and the hydrous basaltic melt (i.e., with about 3.5 wt.% H₂O) are likely to occur at 5.4 (±0.4) kbar and temperatures approaching 1100°C. The olivine gabbros are thus interpreted as cumulates which represent crystallized portions of earlier Arenal-type basalts. Some of the gabbros have been `mildly' tectonized and recrystallized to give mafic granulites that may exhibit a distinct foliation. Below Arenal volcano a zoned magma chamber evolved prior the last eruptive cycle: three distinct andesitic magma layers were produced by simple AFC of a high-alumina basalt (HAB) with assimilation of Tertiary quartz–dioritic and granitic rocks. Early erupted 1968 tephra and 1969 lavas (which represent the first two layers of the upper part of a zoned magma chamber) were produced by simple AFC, with fractionation of plagioclase, pyroxene and magnetite and concomitant assimilation of quartz–dioritic rocks. Assimilation rates were constant (r<sub>1</sub>=0.33) for a relative mass of magma remaining of 0.77–0.80, respectively. Lavas erupted around 1974 are less differentiated and represent the `primitive andesitic magma type' residing within the middle–lower part of the chamber. These lavas were also produced by simple AFC: assimilation rates and the relative mass of magma remaining increased of about 10%, respectively (r₁=0.36, and F=0.89). Ba enrichment of the above lavas is related to selective assimilation of Ba from Tertiary granitic rocks. Lava eruption occurred as a dynamic response to the intrusion of a new magma into the old reservoir. This process caused the instability of the zoned magma column inducing syneruptive mixing between portions of two contiguous magma layers (both within the column itself and at lower levels where the new basalt was intruded into the reservoir). Syneruptive mixing (mingling) within the middle–upper part of the chamber involved fractions of earlier gabbroic cumulitic materials (lavas erupted around 1970). On the contrary, within the lower part of the chamber, mixing between the intruded HAB and the residing andesitic melt was followed by simple fractional crystallization (FC) of the hybrid magma layer (lavas erupted in 1978–1980). By that time the original magma chamber was completely evacuated. Lavas erupted in 1982/1984 were thus modelled by means of `open system' AFCRE (i.e., AFC with continuous recharge of a fractionating magma batch during eruption): in this case assimilation rates were r₁=0.33 and F=0.86. Recharge rates are slightly higher than extrusion rates and may reflect differences in density (between extruded and injected magmas), together with dynamic fluctuations of these parameters during eruption. Ba and LREE (La, Ce) enrichments of these lavas can be related to selective assimilation of Tertiary granitic and quartz–dioritic rocks. Calculated contents for Zr, Y and other REE are in acceptable agreement with the observed values. It is concluded that simple AFC occurs between two distinct eruption cycles and is typical of a period of repose or mild and decreasing volcanic activity. On the contrary, magma mixing, eventually followed by fractional crystallization (FC) of the hybrid magma layer, occurs during an ongoing eruption. Open-system AFCRE is only operative when the original magma chamber has been totally replenished by the new basaltic magma, and seems a prelude to the progressive ceasing of a major eruptive cycle.     

Petrogenesis of an 800 m lava sequence in eastern Uruguay: insights into magma chamber processes beneath the Paraná flood basalt province

Major and trace element along with representative Sr, Nd and Pb isotope data are presented for drill core samples which intersect an 800 m lava pile in eastern Uruguay. The lavas form part of the Paraná flood basalt province, are low-Ti in composition but distinct from the low-Ti Gramado magma type, and have been termed the Treinte Y Trés magma type. The lava pile overlies a large positive gravity anomaly inferred to reflect an east–west trending, mid-crustal mafic intrusive body with a calculated volume of 35,000 km³. Smooth up-section compositional variations in the basalts are interpreted to record magma evolution within this mid-crustal magma chamber. ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb increase throughout the sequence yet Mg remains relatively constant in the lower 200 m of the sequence, suggesting a role for magma chamber recharge. Above this the lavas show a regular, up-section decrease in Mg coupled with increasing ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb and this is interpreted to reflect crystal fractionation combined with crustal contamination. The data provide further evidence that contamination of flood basalt magmas in crustal magma chambers is a common phenomenon and calculations suggest that the amount of crustal addition may be as high as 60–70%. Nevertheless, the effects of this crustal contamination do not appear able to account for the discrepancy between key incompatible trace element ratios and isotope ratios of the lavas and those of any putative mantle plume. In fact, La/Ta decreases with decreasing Mg and increasing ⁸⁷Sr/⁸⁶Sr indicating that the effects of crustal contamination were actually to reduce La/Ta and implying that the parental magmas had very high La/Ta (90). These constraints are clearly inconsistent with an asthenospheric origin for the parental magmas and so, consistent with mass balance calculations, it is inferred that they were derived from the lithospheric mantle.     

Petrology and geochronology of basalt breccia from the 1996 earthquake swarm of Loihi seamount, Hawaii: magmatic history of its 1996 eruption

 Samples of basalt were collected during the Rapid Response cruise to Loihi seamount from a breccia that was probably created by the July to August 1996 Loihi earthquake swarm, the largest swarm ever recorded from a Hawaiian volcano. ²¹⁰Po–²¹⁰Pb dating of two fresh lava blocks from this breccia indicates that they were erupted during the first half of 1996, making this the first documented historical eruption of Loihi. Sonobuoys deployed during the August 1996 cruise recorded popping noises north of the breccia site, indicating that the eruption may have been continuing during the swarm. All of the breccia lava fragments are tholeiitic, like the vast majority of Loihi's most recent lavas. Reverse zoning at the rim of clinopyroxene phenocrysts, and the presence of two chemically distinct olivine phenocryst populations, indicate that the magma for the lavas was mixed just prior to eruption. The trace element geochemistry of these lavas indicates there has been a reversal in Loihi's temporal geochemical trend. Although the new Loihi lavas are similar isotopically and geochemically to recent Kilauea lavas and the mantle conduits for these two volcanoes appear to converge at depth, distinct trace element ratios for their recent lavas preclude common parental magmas for these two active volcanoes. The mineralogy of Loihi's recent tholeiitic lavas signify that they crystallized at moderate depths (∼8–9 km) within the volcano, which is approximately 1 km below the hypocenters for earthquakes from the 1996 swarm. Taken together, the petrological and seismic evidence indicates that Loihi's current magma chamber is considerably deeper than the shallow magma chamber (∼3–4 km) in the adjoining active shield volcanoes. Received: 21 August 1997 / Accepted: 15 February 1998     

Large-scale silicic alkalic magmatism associated with the Buckhorn Caldera, Trans-Pecos Texas, USA: comparison with Pantelleria, Italy

Three major rhyolite systems in the northeastern Davis and adjacent Barrilla Mountains include lava units that bracketed a large pantelleritic ignimbrite (Gomez Tuff) in rapid eruptions spanning 300,000 years. Extensive silicic lavas formed the shields of the Star Mountain Formation (37.2 Ma-K/Ar; 36.84 Ma ³⁹Ar/⁴⁰Ar), and the Adobe Canyon Formation (37.1 Ma-K/Ar; 36.51-³⁹Ar/⁴⁰Ar). The Gomez Tuff (36.6 Ma-K/Ar; 36.74-³⁹Ar/⁴⁰Ar) blanketed a large region around the 18×24 km diameter Buckhorn caldera, within which it ponded, forming sections up to 500 m thick. Gomez eruption was preceded by pantelleritic rhyolite domes (36.87, 36.91 Ma-³⁹Ar/⁴⁰Ar), some of which blocked movement of Star Mountain lava flows. Following collapse, the Buckhorn caldera was filled by trachyte lava. Adobe Canyon rhyolite lavas then covered much of the region. Star Mountain Formation (~220 km³) is composed of multiple flows ranging from quartz trachyte to mildly peralkalic rhyolite; three major types form a total of at least six major flows in the northeastern Davis Mountains. Adobe Canyon Formation (~125 km³) contains fewer flows, some up to 180 m thick, of chemically homogenous, mildly peralkalic comendite, extending up to 40 km. Gomez Tuff (~220 km³) may represent the largest known pantellerite. It is typically less than 100 m thick in extra-caldera sections, where it shows a pyroclastic base and top, although interiors are commonly rheomorphic, containing flow banding and ramp structures. Most sections contain one cooling unit; two sections contain a smaller, upper cooling unit. Chemically, the tuff is fairly homogeneous, but is more evolved than early pantelleritic domes. Overall, although Davis Mountains silicic units were generated through open system processes, the pantellerites appear to have evolved by processes dominated by extensive fractional crystallization from parental trachytes similar to that erupted in pre- and post-caldera lavas. Comparison with the Pantelleria volcano suggests that the most likely parental magma for the Buckhorn series is transitional basalt, similar to that erupted in minor, younger Basin and Range volcanism after about 24 Ma. Roughly contemporaneous mafic lavas associated with the Buckhorn caldera appear to have assimilated or mixed with crustal melts, and, generally, may not be regarded as mafic precursors of the Buckhorn silicic rocks, They thus form a false Daly Gap as opposed to the true basalt/trachyte Daly gap of Pantelleria. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. This paper constitutes part of a special issue dedicated to Bill Bonnichsen on the petrogenesis and volcanology of anorogenic rhyolites.     

Erratum to Compositional variations and magma mixing in the 1991 eruptions of Hudson volcano,Chile

The August 1991 eruptions of Hudson volcano produced ~2.7 km³ (dense rock equivalent, DRE) of basaltic to trachyandesitic pyroclastic deposits, making it one of the largest historical eruptions in South America. Phase 1 of the eruption (P1, April 8) involved both lava flows and a phreatomagmatic eruption from a fissure located in the NW corner of the caldera. The paroxysmal phase (P2) began several days later (April 12) with a Plinian-style eruption from a different vent 4 km to the south-southeast. Tephra from the 1991 eruption ranges in composition from basalt (phase 1) to trachyandesite (phase 2), with a distinct gap between the two erupted phases from 54–60 wt% SiO₂. A trend of decreasing SiO₂ is evident from the earliest part of the phase 2 eruption (unit A, 63–65 wt% SiO₂) to the end (unit D, 60–63 wt% SiO₂). Melt inclusion data and textures suggest that mixing occurred in magmas from both eruptive phases. The basaltic and trachyandesitic magmas can be genetically related through both magma mixing and fractional crystallization processes. A combination of observed phase assemblages, inferred water content, crystallinity, and geothermometry estimates suggest pre-eruptive storage of the phase 2 trachyandesite at pressures between ~50–100 megapascal (MPa) at 972 ± 26°C under water-saturated conditions (log fO₂ –10.33 (±0.2)). It is proposed that rising P1 basaltic magma intersected the lower part of the P2 magma storage region between 2 and 3 km depth. Subsequent mixing between the two magmas preferentially hybridized the lower part of the chamber. Basaltic magma continued advancing towards the surface as a dyke to eventually be erupted in the northwestern part of the Hudson caldera. The presence of tachylite in the P1 products suggests that some of the magma was stalled close to the surface (<0.5 km) prior to eruption. Seismicity related to magma movement and the P1 eruption, combined with chamber overpressure associated with basalt injection, may have created a pathway to the surface for the trachyandesite magma and subsequent P2 eruption at a different vent 4 km to the south-southeast. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

O, S, Sr, and Pb isotope variations in volcanic rocks from the Northern Mariana Islands: implications for crustal recycling in intra-oceanic arcs

Quaternary lavas from the Northern Mariana Islands have respective O- and S-isotope ranges ofδ¹⁸O = +5.7 to +6.6 (‰ SMOW) andδ³⁴S = +2.0to+20.7 (‰ CDT). Chemically evolved andesites and dacites with meanδ¹⁸O = +6.3 ± 0.2 are slightly enriched in¹⁸O with respect to unfractionated basalts of<53%SiO₂ with meanδ¹⁸O = +6.0 ± 0.1. This¹⁸O enrichment can be explained in terms of differentiation of parental mafic magmas havingδ¹⁸O values between +5.7 to +6.2‰ through closed-system crystal fractionation because the lavas from all nine islands of the arc define a coherentδ¹⁸OSiO₂ trend. The S-isotope composition of oxidized magmas is not modified extensively through the degassing of SO₂; therefore, the meanδ³⁴S value of ca. +11‰ for the Mariana lavas is considered to be representative of their source region.The enrichment of¹⁸O and³⁴S in Mariana Arc parental magmas relative to ocean floor basalts withδ¹⁸O ca. + 5.7‰ andδ³⁴S = ca.0.3‰ is attributed to the recycling of¹⁸O- and³⁴S-rich crustal components (sediment withδ¹⁸O = ca. + 25‰ and seawater sulfate withδ³⁴S = ca. +20‰ into the upper mantle source region for these arc magmas. This interpretation is consistent with enrichments of radiogenic Sr and Pb in the same lavas relative to ocean-floor basalts erupted either side of the arc, which are presumed to share a common upper mantle source. This enrichment is considered to reflect the mixing of two components, one having a typical upper mantle composition and the other having a more radiogenic character similar to that of western Pacific pelagic sediments.     

Petrological, geochemical and chronological constraints for the tectonic setting of the Dongco ophiolite in Tibet

The Dongco ophiolite occurred in the middle-western segment of the Bangong-Nujiang suture zone. The thickness of the ophiolite suite is more than 5 km, which is composed, from bottom to top, of the mantle peridotite, mafic-ultramafic cumulates, basic sills (dykes) and basic lava and tectoni- cally emplaced in Jurassic strata (Mugagongru Group). The Dongco cumulates consist of dunite- troctolite-olivine-gabbro, being a part of DTG series of mafic-ultramafic cumulates. The basic lavas are characterized by being rich in alkali (Na2O K2O), TiO2, P2O5 and a LREE-rich type pattern dip- ping right with [La/Yb]=6.94―16.6 as well as a trace elements spider-diagram with normal anomaly of Th, Nb, Ta, Hf. Therefore, the Dongco basic lavas belong to ocean-island basalt (OIB) and dis- tinctly differ from mid-ocean ridge basalt (MORB) and island-arc basalt (IAB) formed in the plate convergence margin. The basic lavas have higher 87Sr/86Sr (0.704363―0.705007), lower 143Nd/144Nd (0.512708―0.512887) and εNd(t ) from 2.7― 5.8, indicating that they derive from a two-components mixing mantle source of depleted mantle (DM) and enriched mantle (EMI). From above it is ready to see that the Dongco ophiolite forms in oceanic island (OIB) where the mantle source is replaced by a large amount of enriched material, therefore it distinctly differs from these ophiolites formed in island-arc and mid-oecan ridge. Newly obtained SHRIMP U-Pb dating for zircon of the cumulate troctolite is 132 ± 3 Ma and whole-rock dating of ~(39)Ar/~(40)Ar for the basalt is 173.4 ± 2.7 Ma and 140.9 ± 2.8 Ma, indicating that the Dongco ophiolite formed at Early Cretaceous and the middle-western segment of the Bangong-Nujiang oceanic basin was still in the developing and evolving period at Early Cretaceous.