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.     

Magma plumbing system of the 2000 eruption of Miyakejima Volcano, Japan

During the 2000 eruption at Miyakejima Volcano, two magmas with different compositions erupted successively from different craters. Magma erupted as spatter from the submarine craters on 27 June is aphyric basaltic andesite (<5 vol% phenocrysts, 51.4–52.2 wt% SiO₂), whereas magma issued as volcanic bombs from the summit caldera on 18 August is plagioclase-phyric basalt (20 vol% phenocrysts, 50.8–51.3 wt% SiO₂). The submarine spatter contains two types of crystal-clots, A-type and A-type (andesitic type). The phenocryst assemblages (plagioclase, pyroxenes and magnetite) and compositions of clinopyroxene in these clots are nearly the same, but only A-type clots contain Ca-poor plagioclase (An < 70). We consider that the A-type clots could have crystallized from a more differentiated andesitic magma than the A-type clots, because FeO*/MgO is not strongly influenced during shallow andesitic differentiation. The summit bombs contain only B-type (basaltic type) crystal-clots of Ca-rich plagioclase, olivine and clinopyroxene. The A-type and B-type clots have often coexisted in Miyakejima lavas of the period 1469–1983, suggesting that the magma storage system consists of independent batches of andesitic and basaltic magmas. According to the temporal variations of mineral compositions in crystal-clots, the andesitic magma became less evolved, and the basaltic magma more evolved, over the past 500 years. We conclude that gradually differentiating basaltic magma has been repeatedly injected into the shallower andesitic magma over this period, causing the andesitic magma to become less evolved with time. The mineral chemistries in crystal-clots of the submarine spatter and 18 August summit bombs of the 2000 eruption fall on the evolution trends of the A-type and B-type clots respectively, suggesting that the shallow andesitic and deeper basaltic magmas existing since 1469 had successively erupted from different craters. The 2000 summit collapse occurred due to drainage of the andesitic magma from the shallower chamber; as the collapse occurred, it may have caused disruption of crustal cumulates which then contaminated the ascending, deeper basalt. Thus, porphyritic basaltic magma could erupt alone without mixing with the andesitic magma from the summit caldera. The historical magma plumbing system of Miyakejima was probably destroyed during the 2000 eruption, and a new one may now form.Editorial responsibility: S Nakada, T Druitt     

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.     

The Sr, Nd and O isotopic studies of the 1991–1995 eruption at Unzen, Japan

The magma generation at Unzen volcano may be considered as the product of crustal material mixed with mantle magma accompanied by fractional crystallization (AFC). The magma in the Unzen volcano is estimated to consist of about 50–80% of residual magma (F) and about 30–70% assimilated crustal material (A) relative to the original magma. Concerning the 1991–1995 eruption, it is estimated that the magma formed as the result of mixing of about 50–60% crustal material and about 55–65% of residual magma. An alternative magma eruption model for the 1991–1995 eruption is proposed here. In the early stage, the isotopic characteristics of 1991 eruption are defined by AFC process in the deeper magma chamber. Later, the magma ascended through the conduit and quiescently stayed for a long time in a shallow reservoir before eruption. The minerals continuously crystallized as phenocrysts especially at the chilled top and outer margin in the shallow chamber. The crystallized phenocryst mush was reworked into the central part of the magma chamber by means of magma convection and rapid magma ascent. Therefore, the reaction between phenocrysts and melt occurs only in internal chemical disequilibrium in the magma chamber. In contrast, the isotopic compositions of the original magma shall be little influenced by the above processes throughout its eruptive history. The 1991–1995 eruptive rocks of the Unzen volcano show their characteristics in Sr and Nd isotopic values independent of their two previous eruptions. However, the isotopic values of early eruptive product could represent the original magma value. This result also supports the previous work of Chen et al. (1993) [Chen, C.H., DePaolo, D.J., Nakada, S., Shieh, Y.N., 1993. Relationship between eruption volume and neodymium isotopic composition at Unzen volcano. Nature 362, 831–834], that suggested the Nd of early or precursory eruptive products could be a qualitative indicator of the maximum size of a continuing or impending eruption.     

Water contents, temperatures and diversity of the magmas of the catastrophic eruption of Nevado del Ruiz, Colombia, November 13, 1985

The petrology of the highly phyric two-pyroxene andesitic to dacitic pyroclastic rocks of the November 13, 1985 eruption of Nevado del Ruiz, Colombia, reveals evidence of: (1) increasingly fractionated bulk compositions with time; (2) tapping of a small magma chamber marginally zoned in regard to H₂O contents (1 to 4%), temperature (960–1090°C), and amount of residual melt (35 to 65%); (3) partial melting and assimilation of degassed zones in the hotter less dense interior of the magma chamber; (4) probable heating, thermal disruption and mineralogic and compositional contamination of the magma body by basaltic magma “underplating”; and (5) crustal contamination of the magmas during ascent and within the magma chamber. Near-crater fall-back or “spill-over” emitted in the middle of the eruptive sequence produced a small pyroclastic flow that became welded in its central and basal portions because of ponding and thus heat conservation on the flat glaciated summit near the Arenas crater. The heterogeneity of Ruiz magmas may be related to the comparatively small volume (0.03 km³) of the eruption, nearly ten times less than the 0.2 km³ of the Plinian phase of Mount St. Helens, and probable steep thermal and P_H2O gradients of a small source magma chamber, estimated at 300 m long and 100 m wide for an assumed ellipsoidal shape.     

Origin of differentiated lavas at Kilauea Volcano,Hawaii: Implications from the 1955 eruption

Kilauea's 1955 eruption was the first major eruption (longer than 2 days) on its east rift zone in 115 years. It lasted 88 days during which 108 × 10⁶ m³ of lava was erupted along a discontinuous, 15-km-long system of fissures. A wide compositional range of lavas was erupted including the most differentiated lavas (5.0 wt% MgO) from a historic Kilauea eruption. Lavas from the first half of the eruption are strongly differentiated (5.0–5.7 wt% MgO); later lavas are weakly to moderately differentiated (6.2–6.7 wt% MgO). Previous studies using only major-element compositions invoked either crystal fractionation (Macdonald and Eaton 1964) or magma mixing (Wright and Fiske 1971) as models to explain the wide compositional variation in the lavas. To further evaluate these models detailed petrographic, mineralogical, and whole-rock, major, and trace element XRF analyses were made of the 1955 lavas. Plagioclase and clinopyroxene in the early and late lavas show no petrographic evidence for magma mixing. Olivines from both the early and late lavas show minor resorption, which is typical of tholeiitic lavas with low MgO contents. Core-to-rim microprobe analyses across olivine, augite, and plagioclase mineral grains give no evidence of disequilibrium features related to mixing. Instead, plots of An/Ab vs distance from the core (D) and %Fo vs (D)^4.5 generated essentially linear trends indicative of simple crystal fractionation. Least-squares, mass-balance calculations for major- and trace-element data using observed mineral compositions yield excellent results for crystal fractionation (sum of residuals squared <0.01 for major elements, and <5% for trace elements); magma mixing produced less satisfactory results especially for Cr. Furthermore, trace-element plots of Zr vs Sr, Cr, and A1₂O₃ generate curved trends indicative of crystal fractionation processes. There is no evidence that mixing occurred in the 1955 lavas. Instead, the data are best explained by crystal fractionation involving a reservoir that extends at least 15 km along Kilauea's east rift zone. A dike was intruded into the rift zone from the summit reservoir eight days after the eruption started. Instead of causing magma mixing, the dike probably acted as a hydraulic plunger forcing more of the stored magma to be erupted.     

Temporal variations in magma composition at Merapi Volcano (Central Java, Indonesia): magmatic cycles during the past 2000 years of explosive activity

Merapi Volcano (Central Java, Indonesia) has been frequently active during Middle to Late Holocene time producing basalts and basaltic andesites of medium-K composition in earlier stages of activity and high-K magmas from 1900 ¹⁴C yr BP to the present. Radiocarbon dating of pyroclastic deposits indicates an almost continuous activity with periods of high eruption rates alternating with shorter time spans of distinctly reduced eruptive frequency since the first appearance of high-K volcanic rocks. Geochemical data of 28 well-dated, prehistoric pyroclastic flows of the Merapi high-K series indicate systematic cyclic variations. These medium-term compositional variations result from a complex interplay of several magmatic processes, which ultimately control the periodicity and frequency of eruptions at Merapi. Low eruption rates and the absence of new influxes of primitive magma from depth allow the generation of basaltic andesite magma (56–57 wt% SiO₂) in a small-volume magma reservoir through fractional crystallisation from parental mafic magma (52–53 wt% SiO₂) in periods of low eruptive frequency. Magmas of intermediate composition erupted during these stages provide evidence for periodic withdrawal of magma from a steadily fractionating magma chamber. Subsequent periods are characterised by high eruption rates that coincide with shifts of whole-rock compositions from basaltic andesite to basalt. This compositional variation is interpreted to originate from influxes of primitive magma into a continuously active magma chamber, triggering the eruption of evolved magma after periods of low eruptive frequency. Batches of primitive magma eventually mix with residual magma in the magmatic reservoir to decrease whole-rock SiO₂ contents. Supply of primitive magma at Merapi appears to be sufficiently frequent that andesites or more differentiated rock types were not generated during the past 2000 years of activity. Cyclic variations also occurred during the recent eruptive period since AD 1883. The most recent eruptive episode of Merapi is characterised by essentially uniform magma compositions that may imply the existence of a continuously active magma reservoir, maintained in a quasi-steady state by magma recharge. The whole-rock compositions at the upper limit of the total SiO₂ range of the Merapi suite could also indicate the beginning of another period of high eruption rates and shifts towards more mafic compositions.     

Variation of volatile concentration in a magma system of Satsuma-Iwojima volcano deduced from melt inclusion analyses

Chemical analyses of 30 melt inclusions from Satsuma-Iwojima volcano, Japan, were carried out to investigate volatile evolution in a magma chamber beneath the volcano from about 6300 yr BP to the present. Large variations in volatile concentrations of melts were observed. (1) Water concentration of rhyolitic melts decreases with time; 3–4.6 wt.% at the time of latest caldera-forming eruption of Takeshima pyroclastic flow deposit (ca. 6300 yr BP), 3 wt.% for small pyroclastic flow (ca. 1300 yr BP) of Iwodake, post-caldera rhyolitic dome, and 0.7–1.4 wt.% for submarine lava eruption (Showa-Iwojima) in 1934. (2) Rhyolitic melts of the Takeshima and Iwodake eruptions contained CO₂ of less than 40 ppm, while the Showa-Iwojima melt has higher CO₂ concentration of up to 140 ppm. (3) Water and CO₂ concentrations of basaltic to andesitic melt of Inamuradake, a post-caldera basaltic scoria cone, are 1.2–2.8 wt.% and ≤290 ppm, respectively.Volatile evolution in the magma chamber is interpreted as follows: (1) the rhyolitic magma at the time of the latest caldera-forming eruption (ca. 6300 yr BP) was gas-saturated due to pressure variation in the magma chamber because the large variation in water concentration of the melt was attributed to exsolution of volatile in the magma prior to the eruption. Iwodake eruption (ca. 1300 yr BP) was caused by a remnant of the caldera-forming rhyolitic magma, suggested from the similarity of major element composition between these magmas. (2) Volatile composition of the Showa-Iwojima rhyolitic melt agrees with that of magmatic gases presently discharging from a summit of Iwodake, indicating the low pressure degassing condition. (3) The degassing of the magma chamber by magma convection in a conduit of Iwodake during non-eruptive but active degassing period for longer than 800 years decreased water concentration of the rhyolitic magma. (4) Geological and petrological observations indicate that a stratified magma chamber, which consists of a lower basaltic layer and an upper rhyolitic layer, might have existed during the post-caldera stage. Addition of CO₂ from the underlying basaltic magma to the upper gas-undersaturated (degassed) rhyolitic magma increased CO₂ concentration of the rhyolitic 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 dacitic magmatism at batur volcano, Bali, Indonesia: Implications for the origins of stratovolcano calderas

Batur is an active stratovolcano on the island of Bali, Indonesia, with a large, well-formed caldera whose formation is correlated with the eruption about 23,700 years ago of a thick ignimbrite sheet. Our study of the volcanic stratigraphy and geochemistry of Batur shows the formation of the caldera was signalled by a change in the composition of the erupting material from basaltic and andesitic to dacitic. The dacitic rocks are glassy, possess equilibrium phenocryst assemblages, and display compositional characteristics consistent with an origin by crystal-liquid fractionation from more mafic parent magmas in a shallow chamber, possibly at 1.5 km depth and 1000–1070°C.However, although separated by a gap of 6 wt.% SiO₂, the dacitic rocks are clearly related in their minor- and trace-element geochemistry to those basalts and basaltic andesites erupted after the caldera was formed rather than to the andesites erupted immediately before the dacites first appeared. We infer from this and published experimental modelling of the possible crystallization behaviour of basaltic magma chambers that a magmatic cycle involving caldera formation began independently of the previous activity of Batur by formation of a new, closed-system magma chamber beneath the volcano. Fractional crystallization, possibly at the walls of the chamber, led to the early production of derivative siliceous magmas and, consequently, to caldera formation, while most of the magma retained its original composition. The postcaldera Batur basalts represent the largely undifferentiated core liquids of this chamber.This model contrasts with the traditional evolutionary model for stratovolcano calderas but may be applicable to the origins of calderas similar to that of Batur, particularly those in volcanic island arcs.     

The eruptive history of the Tequila volcanic field, western Mexico: ages, volumes, and relative proportions of lava types

The eruptive history of the Tequila volcanic field (1600 km²) in the western Trans-Mexican Volcanic Belt is based on ⁴⁰Ar/³⁹Ar chronology and volume estimates for eruptive units younger than 1 Ma. Ages are reported for 49 volcanic units, including Volcán Tequila (an andesitic stratovolcano) and peripheral domes, flows, and scoria cones. Volumes of volcanic units 1 Ma were obtained with the aid of field mapping, ortho aerial photographs, digital elevation models (DEMs), and ArcGIS software. Between 1120 and 200 kyrs ago, a bimodal distribution of rhyolite (~35 km³) and high-Ti basalt (~39 km³) dominated the volcanic field. Between 685 and 225 kyrs ago, less than 3 km³ of andesite and dacite erupted from more than 15 isolated vents; these lavas are crystal-poor and show little evidence of storage in an upper crustal chamber. Approximately 200 kyr ago, ~31 km³ of andesite erupted to form the stratocone of Volcán Tequila. The phenocryst assemblage of these lavas suggests storage within a chamber at ~2–3 km depth. After a hiatus of ~110 kyrs, ~15 km³ of andesite erupted along the W and SE flanks of Volcán Tequila at ~90 ka, most likely from a second, discrete magma chamber located at ~5–6 km depth. The youngest volcanic feature (~60 ka) is the small andesitic volcano Cerro Tomasillo (~2 km³). Over the last 1 Myr, a total of 128±22 km³ of lava erupted in the Tequila volcanic field, leading to an average eruption rate of ~0.13 km³/kyr. This volume erupted over ~1600 km², leading to an average lava accumulation rate of ~8 cm/kyr. The relative proportions of lava types are ~22–43% basalt, ~0.4–1% basaltic andesite, ~29–54% andesite, ~2–3% dacite, and ~18–40% rhyolite. On the basis of eruptive sequence, proportions of lava types, phenocryst assemblages, textures, and chemical composition, the lavas do not reflect the differentiation of a single (or only a few) parental liquids in a long-lived magma chamber. The rhyolites are geochemically diverse and were likely formed by episodic partial melting of upper crustal rocks in response to emplacement of basalts. There are no examples of mingled rhyolitic and basaltic magmas. Whatever mechanism is invoked to explain the generation of andesite at the Tequila volcanic field, it must be consistent with a dominantly bimodal distribution of high-Ti basalt and rhyolite for an 800 kyr interval beginning ~1 Ma, which abruptly switched to punctuated bursts of predominantly andesitic volcanism over the last 200 kyrs.Electronic Supplementary Material Supplementary material is available in the online version of this article at Editorial responsility: J. Donnelly-NolanThis revised version was published online in January 2005 with corrections to Tables 1 and 3.An erratum to this article can be found at     

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.     

Assessment of a shallow magmatic system: the 1888–90 eruption,Vulcano Island,Italy

The magmatic system feeding the last eruption of the volcano La Fossa, Vulcano Island, Italy was studied. The petrogenetic mechanisms controlling the differentiation of erupted rocks were investigated through petrography, mineral chemistry, major, trace and rare earth element and Sr, Nd and Pb isotopic geochemistry. In addition, melt inclusion and fluid inclusion data were collected on both juvenile material and xenolithic partially melted metamorphic clasts to quantify the P-T conditions of the magma chamber feeding the eruption. A regular and continuous chemical zoning has been highlighted: rhyolites are the first erupted products, followed by trachytes and latites, whereas rhyolitic compositions were also found in the upper part of the sequence. The chemical and isotopic composition of the rhyolites indicates that they originated by fractional crystallization from latitic magmas plus the assimilation of crustal material; the trachytes represent hybrid magmas resulting from the mixing of latites and rhyolites, contaminated in the shallow magmatic system. The erupted products, primarily compositionally zoned from latites to rhyolites, are heterogeneous due to syn-eruptive mingling. The occurrence of magmacrust interaction processes, evidenced by isotopic variations (⁸⁷Sr/⁸⁶Sr=0.70474±3 to 0.70511±3; ¹⁴³Nd/¹⁴⁴Nd=0.512550±6 to 0.512614±8; ²⁰⁶Pb/²⁰⁴Pb=19.318–19.489; ²⁰⁷Pb/²⁰⁴Pb=15.642–15.782; ²⁰⁸Pb/²⁰⁴Pb=39.175–39.613), is confirmed by the presence of partially melted metamorphic xenoliths, with ⁸⁷Sr/⁸⁶Sr=0.71633±6 to 0.72505±2 and ¹⁴³Nd/¹⁴⁴Nd=0.51229±7, in rhyolites and trachytes. AFC calculations indicate a few percentage contribution of crustal material to the differentiating magmas. Thermometric measurements on melt inclusions indicate that the crystallization temperatures of the latites and trachytes were in the range of 1050–1100° C, whereas the temperature of the rhyolites appears to have been around 1000°C at the time of the eruption. Compositional data on melt inclusions reveal that the magmas involved in the eruption contained about 1–1.5 wt.% dissolved H₂O in pre-eruptive conditions. Secondary fluid inclusions found in metamorphic xenoliths give low equilibration pressure data (30–60 MPa), giving the location of the higher portions of the chamber at around 1500–2000 m of depth.     

Temporal variation in the lavas of Mayon volcano, Philippines

Chemical and petrographic analyses of 51 sequential lava flows from the central vent of Mayon volcano show cyclical variation. In the two most recent cycles, from 1800 to 1876 and from 1881 to the present, one to three basaltic flows are followed by six to ten andesitic flows. Modal and whole-rock chemical parameters show the most regular cyclical variation; calculated groundmass chemical parameters vary less regularly. There is also a long-term trend, over approximately 1700 years of exposed section, toward more basic compositions.The cyclical variation in modes and the chemical composition of the lavas apparently results from periodic influxes of basaltic magma from depth into a shallow magma system. Fractional crystallization of olivine, augite, hypersthene, calcic plagioclase, magnetite and pargasitic hornblende produces successively more andesitic lavas until the next influx of basaltic magma. Differentiation in a deep zone of magma generation is not excluded by the data, but is more likely responsible for the overall change toward more basic compositions than for the cyclical variation.Three points in a cycle — the beginning of basaltic lavas, the beginning of andesitic lavas and a leveling-off of SiO₂, K₂ O and K₂O/Na₂O values — correspond roughly to the beginning of frequent effusive eruptions (with or without an early Plinian eruption), frequent weak to moderately explosive (Strombolian) eruptions, and less frequent explosive (Vulcanian) eruptions, respectively. Recognition of the current stage in a cycle can give a qualitative indication of the nature of forthcoming eruptions. Changes in several specific parameters may precede basaltic lavas and allow early detection of basaltic influxes. These include minima in the glass inclusion/plagioclase phenocryst and phenocryst/groundmass ratios, vesicularity and groundmass TiO₂, a decrease in hypersthene phenocrysts, and constant values for the whole-rock K₂O/Na₂O ratio. The Mayon area is densely populated, making prediction of eruption type important for safety and land-use planning.     

Eruption and emplacement of a basaltic welded ignimbrite during caldera formation on Gran Canaria

The 14.1 Ma old composite ignimbrite cooling unit P1 (45 km³) on Gran Canaria comprises a lower mixed rhyolite-trachyte tuff, a central rhyolite-basalt mixed tuff, and a slightly rhyolite-contaminated basaltic tuff at the top. The basaltic tuff is compositionally zoned with (a) an upward change in basalt composition to higher MgO content (4.3–5.2 wt.%), (b) variably admixed rhyolite or trachyte (commonly <5 wt.%), and (c) an upward increasing abundance of basaltic and plutonic lithic fragments and cognate cumulate fragments. The basaltic tuff is divided into three structural units: (I) the welded basaltic ignimbrite, which forms the thickest part (c. 95 vol.%) and is the main subject of the present paper; (II) poorly consolidated massive, bomb- and block-rich beds interpreted as phreatomagmatic pyroclastic flow deposits; and (III) various facies of reworked basaltic tuff. Tuff unit I is a basaltic ignimbrite rather than a lava flow because of the absence of top and bottom breccias, radial sheet-like distribution around the central Tejeda caldera, thickening in valleys but also covering higher ground, and local erosion of the underlying P1 ash. A gradual transition from dense rock in the interior to ash at the top of the basaltic ignimbrite reflects a decrease in welding; the shape of the welding profile is typical for emplacement temperatures well above the minimum welding temperature. A similar transition occurs at the base where the ignimbrite was emplaced on cold ground in distal sections. In proximal sections the base is dense where it was emplaced on hot felsic P1 tuff. The intensity of welding, especially at the base, and the presence of spherical particles and of mantled and composite particles formed by accretion and coalescence in a viscous state imply that the flow was a suspension of hot magma droplets. The flow most likely had to be density stratified and highly turbulent to prevent massive coalescence and collapse. Model calculations suggest eruption through low pyroclastic fountains (<1000 m high) with limited cooling during eruption and turbulent flow from an initial temperature of 1160°C. The large volume of 26 km³ of erupted basalt compared with only 16 km³ of the evolved P1 magmas, and the extremely high discharge rates inferred from model calculations are unusual for a basaltic eruption. It is suggested that the basaltic magma was erupted and emplaced in a fashion commonly only attributed to felsic magmas because it utilized the felsic P1 magma chamber and its ring-fissure conduits. Evolution of the entire P1 eruption was controlled by withdrawal dynamics involving magmas differing in viscosity by more than four orders of magnitude. The basaltic eruption phase was initially driven by buoyancy of the basaltic magma at chamber depth and continued degassing of felsic magma, but most of the large volume of basalt magma was driven out of the reservoir by subsidence of a c. 10 km diameter roof block, which followed a decrease in magma chamber pressure during low viscosity basaltic outflow.     

Volcano-stratigraphy and geochemistry of collision-related volcanism on the Erzurum–Kars Plateau, northeastern Turkey

The Eastern Anatolia Region exhibits one of the world's best exposed and most complete transects across a volcanic province related to a continental collision zone. Within this region, the Erzurum–Kars Plateau is of special importance since it contains the full record of collision-related volcanism from Middle Miocene to Pliocene. This paper presents a detailed study of the volcanic stratigraphy of the plateau, together with new K–Ar ages and several hundred new major- and trace-element analyses in order to evaluate the magmatic evolution of the plateau and its links to collision-related tectonic processes. The data show that the volcanic units of the Erzurum–Kars Plateau cover a broad compositional range from basalts to rhyolites. Correlations between six logged, volcano-stratigraphic sections suggest that the volcanic activity may be divided into three consecutive Stages, and that activity begins slightly earlier in the west of the plateau than in the east. The Early Stage (mostly from 11 to 6 Ma) is characterised by bimodal volcanism, made up of mafic-intermediate lavas and acid pyroclastic rocks. Their petrography and high-Y fractionation trend suggest that they result from crystallization of anhydrous assemblages at relatively shallow crustal levels. Their stratigraphy and geochemistry suggest that the basic rocks erupted from small transient chambers while the acid rocks erupted from large, zoned magma chambers. The Middle Stage (mostly from 6–5 Ma) is characterised by unimodal volcanism made up predominantly of andesitic–dacitic lavas. Their petrography and low-Y fractionation trend indicate that they resulted from crystallization of hydrous (amphibole-bearing) assemblages in deeper magma chambers. The Late Stage (mostly 5–2.7 Ma) is again characterised by bimodal volcanism, made up mainly of plateau basalts and basaltic andesite lavas and felsic domes. Their petrography and high-Y fractionation trend indicate that they resulted from crystallization of anhydrous assemblages at relatively shallow crustal levels. AFC modelling shows that crustal assimilation was most important in the deeper magma chambers of the Middle Stage. The geochemical data indicate that the parental magma changed little throughout the evolution of the plateau. This parental magma exhibits a distinctive subduction signature represented by selective enrichment in LILE and LREE thought to have been inherited from a lithosphere modified by pre-collision subduction events. The relationships between magmatism and tectonics support models in which delamination of thickened subcontinental lithosphere cause uplift accompanied by melting of this enriched lithosphere. Magma ascent, and possibly magma generation, is then strongly controlled by strike-slip faulting and associated pull-apart extensional tectonics.     

Chromite in the Paricutin lava flows (1943–1952)

Small euhedral chromite crystals are found in olivine macrophenocrysts (Fo_80–84) from the basaltic andesites (150 ppm Cr) erupted in 1943–1947, and in orthopyroxene macrophenocrysts of the andesites (75 ppm Cr) erupted in 1947–1952. The majority of the chromite octahedra are 5–20 μm in diameter, and some are found in clusters and linear chains of three or more oriented chromite crystals. The composition of the majority of the chromite grains within olivine and orthopyroxene macrophenocrysts is Fe²⁺/(Fe²⁺+Mg)=0.5–0.6, Cr/(Cr+Al)=0.5–0.6 and Fe³⁺/(Fe³⁺+Al+Cr)=0.2–0.3. The chromite crystals in contact with the groundmass are larger, subhedral, and grade in composition from chromite cores to magnetite rims. Comparison of the composition of chromite with those of other volcanic rocks shows that the most primitive Paricutin chromite is richer in total iron and higher in Fe³⁺/(Fe³⁺+Al+Cr) than primary chromite in most lavas. The linear chains of oriented chromite octahedra are found in olivine and orthopyroxene macrophenocrysts, and in the groundmass. These chromite chains are thought to result from diffusion-controlled crystallization because of the very high partition coefficient (1000) of Cr between chromite and melt. We conclude that chromite was a primary phase in the lavas at the time of extrusion and that magnetite only crystallized after extrusion during cooling of the lava flows. The presence of chromite microphenocrysts in andesitic lavas containing as little as 70 ppm Cr can be explained by dissolved H₂O in the melt depressing the liquidus temperature for orthopyroxene such that chromite becomes a liquidus phase. The influence of dissolved H₂O can also explain the lack of plagioclase macrophenocrysts in most of the lavas and the relatively high partition coefficient (20) of Ni between olivine and melt and the high partition coefficient (40) of Cr between orthopyroxene and melt. The liquidus temperature of the basaltic andesite is estimated to have been less than 1140°C, assuming H₂O>1 wt.%, and the log f_O2 to have been above that of the QFM buffer. The chromite and orthopyroxene liquidus temperature of the andesites, assuming H₂O>1 wt.%, is estimated to have been 1100°C or less. The derivation of the later andesites from the earlier basaltic andesites has been explained by a combination of fractional crystallization of olivine, orthopyroxene and plagioclase, and assimilation of xenoliths. The significantly lower Cr, Ni and Mg of the andesites may have been in part due to the separation of olivine macrophenocrysts plus enclosed chromite crystals from the earlier basaltic andesites.     

Apoyo caldera, Nicaragua: A major quaternary silicic eruptive center

Apoyo caldera, near Granada, Nicaragua, was formed by two phases of collapse following explosive eruptions of dacite pumice about 23,000 yr B.P. The caldera sits atop an older volcanic center consisting of lava flows, domes, and ignimbrite (ash-flow tuff). The earliest lavas erupted were compositionally homogeneous basalt flows, which were later intruded by small andesite and dacite flows along a well defined set of N—S-trending regional faults. Collapse of the roof of the magma chamber occurred along near-vertical ring faults during two widely separated eruptions. Field evidence suggests that the climactic eruption sequence opened with a powerful plinian blast, followed by eruption column collapse, which generated a complex sequence of pyroclastic surge and ignimbrite deposits and initiated caldera collapse. A period of quiescence was marked by the eruption of scoria-bearing tuff from the nearby Masaya caldera and the development of a soil horizon. Violent plinian eruptions then resumed from a vent located within the caldera. A second phase of caldera collapse followed, accompanied by the effusion of late-stage andesitic lavas, indicating the presence of an underlying zoned magma chamber. Detailed isopach and isopleth maps of the plinian deposits indicate moderate to great column heights and muzzle velocities compared to other eruptions of similar volume. Mapping of the Apoyo airfall and ignimbrite deposits gives a volume of 17.2 km³ within the 1-mm isopach. Crystal concentration studies show that the true erupted volume was 30.5 km³ (10.7 km³ Dense Rock Equivalent), approximately the volume necessary to fill the caldera. A vent area located in the northeast quadrant of the present caldera lake is deduced for all the silicic pyroclastic eruptions. This vent area is controlled by N—S-trending precaldera faults related to left-lateral motion along the adjacent volcanic segment break. Fractional crystallization of calc-alkaline basaltic magma was the primary differentiation process which led to the intermediate to silicic products erupted at Apoyo. Prior to caldera collapse, highly atypical tholeiitic magmas resembling low-K, high-Ca oceanic ridge basalts were erupted along tension faults peripheral to the magma chamber. The injection of tholeiitic magmas may have contributed to the paroxysmal caldera-forming eruptions.     

Geological evolution of a pleistocene rhyolitic center — Sierra La Primavera, Jalisco, México

The Sierra La Primavera volcanic complex consists of late Pleistocene comenditic lava flows and domes. ash-flow tuff, air-fall pumice, and cold caldera-lake sediments. The earliest lavas were erupted about 120,000 years ago, and were followed approximately 95,000 years ago by the eruption of about 20 km³ of magma as ash flows that form the compositionally-zoned Tala Tuff. Collapse of the roof zone of the magma chamber led to the formation of a shallow 11-km-diameter caldera. It soon filled with water, forming a caldera lake in which sediment began to collect. At about the same time, two central domes erupted through the middle of the lake and a “giant pumice horizon”, an important stratigraphic marker, was deposited. Shortly thereafter ring domes erupted along two parallel arcs: one along the northeast portion of the ring fracture, and the other crossing the middle of the lake. All these events occurred during a period of approximately 5,000–10,000 years. Sedimentation continued and a period of volcanic quiescence was marked by the deposition of some 30 m of fine-grained ashy sediments virtually free from pumice lapilli. Approximately 75,000 years ago, a new group of ring domes erupted at the southern margin of the lake. These domes are lapped by only 10–20 m of sediments, as uplift resulting from renewed insurgence of magma brought an end to the lake. This uplift culminated in the eruption, beginning approximately 60,000 years ago, of aphyric lavas along a southern arc. The youngest of these lavas erupted approximately 20,000–30,000 years ago.The four major fault systems in the Sierra La Primavera are related to caldera collapse or to uplift caused by the insurgence of the southern are magma. Steam vents and larga-discharge 65°C hot springs are associated with the faulting. Calculated equilibrium temperatures of the geothermal fluids are 170°C, but temperatures in excess of 240°C have been encountered in an exploratory drill hole.A seismic survey showed attenuation of both S and P waves within the caldera, P waves attenuated more severely than S waves. The greatest attenuation is associated with an area of steam vents, and the rapid lateral variations in attenuation suggest that they are produced by a shallow geothermal system rather than by underlying magma.     

The 1999 and 2000 eruptions of Mount Cameroon: eruption behaviour and petrochemistry of lava

Mount Cameroon (4,095 m high and with a volume of ~1,200 km³) is one of the most active volcanoes in Africa, having erupted seven times in the last 100 years. This stratovolcano of basanite and hawaiite lavas has an elliptical shape, with over a hundred cones around its flanks and summit region aligned parallel to its NE--SW-trending long axis. The 1999 (28 March–22 April) eruption was restricted to two sites: ~2,650 m (site 1) and ~1,500 m (site 2). Similarly, in the eruption in 2000 (28 May–19 June), activity occurred at two sites: ~4,095 m (site 1) and ~3,300 m (site 2). During both eruptions, the higher vents were more explosive, with strombolian activity, while the lower vents were more effusive. Accordingly, most of the lava (~8×10⁷ m³ in 1999 and ~6×10⁶ m³ in 2000) was emitted from the lower sites. The 1999–2000 lavas are predominantly basanites with low Ni (5–79 ppm), Cr (40–161 ppm) and mg numbers (34–40). Olivine (Fo_77–85, phenocrysts and Fo_68–72, microlites), clinopyroxene (Wo₄₇En₄₁Fs₁₀ to Wo₅₁En₃₄Fs₁₅), plagioclase (An_49–67) and titanomagnetite are the principal phenocryst and groundmass phases. The lavas contain xenocrysts of olivine and clinopyroxene, which are interpreted as fragments of intrusive rocks disrupted by magma ascent. Major and trace element characteristics point to early fractionation of olivine. The clinopyroxenes (Al₂O₃ 1.36–7.83 wt%) have high Al^iv/Al^vi ratios (1.3–1.8) and are rich in TiO₂, characteristics typical of low pressure clinopyroxenes. Geochemical differences between the 1999–2000 lavas and those from previous eruptions, such as higher Nb/Zr of the former, suggest that different eruptions discharged magmas that evolved differently in space and time. Geophysical and petrological data indicate that these fractionated magmas originated just below the geophysical Moho (at 20–22 km) in the lithospheric mantle. During ascent, the magmas disrupted intrusions and earlier magma pockets. The main ascent path is below the summit, where newly arrived magma degasses. Degassed magma simultaneously intrudes the flank rift zones where most lava is extruded.An erratum to this article can be found at