Eruptive versus non-eruptive behaviour of large calderas: the example of Campi Flegrei caldera (southern Italy)

Caldera eruptions are among the most hazardous of natural phenomena. Many calderas around the world are active and are characterised by recurrent uplift and subsidence periods due to the dynamics of their magma reservoirs. These periods of unrest are, in some cases, accompanied by eruptions. At Campi Flegrei caldera (CFc), which is an area characterised by very high volcanic risk, the recurrence of this behaviour has stimulated the study of the rock rheology around the magma chamber, in order to estimate the likelihood of an eruption. This study considers different scenarios of shallow crustal behaviour, taking into account the earlier models of CFc ground deformation and caldera eruptions, and including recent geophysical investigations of the area. A semi-quantitative evaluation of the different factors that lead to magma storage or to its eruption (such as magma chamber size, wall-rock viscosity, temperature, and regional tectonic strain rate) is reported here for elastic and viscoelastic conditions. Considering the large magmatic sources of the CFc ignimbrite eruptions (400–2,000 km³) and a wall-rock viscosity between 10¹⁸ and 10²⁰ Pa s, the conditions for eruptive failure are difficult to attain. Smaller source dimensions (a few cubic kilometres) promote the condition for fracture (eruption) rather than for the flow of wall rock. We also analyse the influence of the regional extensional stress regime on magma storage and eruptions, and the thermal stress as a possible source of caldera uplift. The present study also emphasises the difficulty of distinguishing eruption and non-eruption scenarios at CFc, since an unambiguous model that accounts for the rock rheology, magma-source dimensions and locations and regional stress field influences is still lacking.     

Magma contamination by direct wall rock interaction: constraints from xenoliths from the walls of a carbonate-hosted magma chamber (Vesuvius 1944 eruption)

During the 1944 eruption of Vesuvius different types of xenoliths were ejected. They represent fragments of the walls of a low volume (<0.5 km³) shallow (3–4 km depth) magma chamber. The study of these xenoliths enables us to estimate the amount of contamination occurring at the boundary of a high-T alkaline magma chamber hosted in carbonate rocks. The process of contamination of the magma by carbonates can be modelled, using isotopic and chemical data, as a mixing between magma and marbles. Mass exchanges occur at the boundary between the crystallizing magma and marble wall rocks, where endoskarn forms. The contamination of the solidification front of the chamber is very limited. The solidification front and the skarn shell effectively isolate the interior of the magma chamber from new inputs of contaminants from the carbonate wall rocks. Therefore, the main volume of magma, hosted in the magma chamber, did not undergo any significant mass exchange with the wall rocks.     

Role of precursory less-viscous mixed magma in the eruption of phenocryst-rich magma: evidence from the Hokkaido-Komagatake 1929 eruption

During the 1929 activity of Hokkaido-Komagatake volcano, the Plinian eruption of a phenocryst-rich andesite was preceded by a small eruption of more mafic magma formed by magma mixing. A similar eruption sequence has been reported for some other eruptions (Pallister et al. 1996; Venezky and Rutherford 1997), suggesting that eruption of a mixed magma is a precursor of phenocryst-rich magmas. For the purpose of understanding the tapping processes of the phenocryst-rich magma chamber, we investigated the temporal variation in the erupted magma and estimated the viscosity and density of the end-member and mixed magmas with constraints drawn from petrography. For the precursory mixed magma we estimate 33dž vol.% phenocrysts, andesitic-dacitic melt composition, 3 wt.% H₂O content, and temperature of 1040°C. In comparison, for the climactic, silicic end-member magma we estimate 48Dž vol.% phenocryst, high-silica rhyolitic melt, 3 wt.% H₂O, and temperature of 950°C, respectively. The mafic end-member magma, which was not erupted, is thought to be an almost aphyric basaltic-andesitic magma, based on mass balance calculation of the phenocryst content. The proportion of the mafic end-member magma component in the mixed magma was calculated to be 20-40 wt.%. On the basis of these data, we estimate magma viscosities of 10^3.9, 10^6.9, and 10^2.0 Pa s for the mixed, silicic end-member, and mafic end-member magmas, respectively. The calculated density differences among these magmas are inconsequential when possible errors are considered. We calculate the minimum excess pressure required for dike propagation to be 31 MPa for the silicic end-member magma and 8 MPa for the mixed magma, using the estimated viscosity and dike propagation model of Rubin (1995). If we assume that excess pressure is limited by the wall rock strength of the magma chamber, excess pressure retainable in the magma chamber is less than ca. 20 MPa. This suggests that the mixed magma was able to ascend to the surface without freezing, whereas the viscous silicic end-member magma could not. The formation and precursory eruption of the mixed magma are, therefore, effective and necessary initiation processes for the phenocryst-rich, viscous magma eruption.     

Effusive eruptions from a large silicic magma chamber: the Bearhead Rhyolite,Jemez volcanic field,NM

Large continental silicic magma systems commonly produce voluminous ignimbrites and associated caldera collapse events. Less conspicuous and relatively poorly documented are cases in which silicic magma chambers of similar size to those associated with caldera-forming events produce dominantly effusive eruptions of small-volume rhyolite domes and flows. The Bearhead Rhyolite and associated Peralta Tuff Member in the Jemez volcanic field, New Mexico, represent small-volume eruptions from a large silicic magma system in which no caldera-forming event occurred, and thus may have implications for the genesis and eruption of large volumes of silicic magma and the long-term evolution of continental silicic magma systems.⁴⁰Ar/³⁹Ar dating reveals that most units mapped as Bearhead Rhyolite and Peralta Tuff (the Main Group) were erupted during an ∼540 ka interval between 7.06 and 6.52 Ma. These rocks define a chemically coherent group of high-silica rhyolites that can be related by simple fractional crystallization models. Preceding the Main Group, minor amounts of unrelated trachydacite and low silica rhyolite were erupted at ∼11–9 and ∼8 Ma, respectively, whereas subsequent to the Main Group minor amounts of unrelated rhyolites were erupted at ∼6.1 and ∼1.5 Ma.The chemical coherency, apparent fractional crystallization-derived geochemical trends, large areal distribution of rhyolite domes (∼200 km²), and presence of a major hydrothermal system support the hypothesis that Main Group magmas were derived from a single, large, shallow magma chamber. The ∼540 ka eruptive interval demands input of heat into the system by replenishment with silicic melts, or basaltic underplating to maintain the Bearhead Rhyolite magma chamber.Although the volatile content of Main Group magmas was within the range of rhyolites from major caldera-forming eruptions such as the Bandelier and Bishop Tuffs, eruptions were smaller volume and dominantly effusive. Bearhead Rhyolite domes occur at the intersection of faults, and are cut by faults, suggesting that the magma chamber was structurally vented preventing volatiles from accumulating to levels high enough to trigger a caldera-forming eruption.     

Conditions for the arrest of a vertical propagating dyke

Magma ascent towards the Earth’s surface occurs through dyke propagation in the vast majority of cases. We investigate two purely mechanical effects unrelated to cooling or solidification that lead to the arrest of propagation, so that no eruption occurs. The first is that the input of magma from the source is not maintained continuously, such that a fixed volume of magma is released. Laboratory experiments show that, in this case, the dyke stops at a finite distance from the source. This behaviour is specific to the fracturing process in 3-D. We derive a relationship for the minimum magma volume required for an eruption as a function of magma buoyancy and source depth. When large magma volumes are available, eruption may also be prevented by a thick low density layer in the upper crust. Numerical studies of dyke propagation show that the dyke continues to rise even though it is negatively buoyant. Magma accumulates in a swollen nose region at the interface between the low density layer and the dense basement. Magma overpressure is largest at this interface and increases with increasing penetration into the upper layer. It may become large enough to induce horizontal fractures in the dyke walls and lateral intrusion of a sill, which prevents eruption. This requires that the thickness of the low density layer exceeds a threshold value that depends on the density contrast between magma and host rock. If the magma volume is smaller than a threshold value, neither sill intrusion nor eruption are possible and magma gets stored in a horizontal blade-shaped dyke straddling the interface. Scaling laws for variations of ascent rate and for the minimum magma volume allow diagnosis of a failed eruption.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

Some thermal problems associated with magma migration

The mechanisms by which magma migrates from the point in the earth's interior where melting occurs to the earth's surface are poorly understood. In this paper several aspects of this problem are examined. Magma can migrate upward due to its differential buoyancy on the scale of crystalline grains or as large dispairs. Magma transport is an effective means of heat transport. Magma transport at a rate of 0.15 cm/yr is equivalent to a heat flow of 10^-6 cal/cm² s. If magma encounters country rock witha lower melting point the original magma is likely to solidify while melting the country rock. This would be an effective mechanism of purging silicic rocks and incomparable elements from the lower crust. Under some circumstances magma must penetrate up to 100 km or more of cold lithospheric rock. In order for this magma to reach the surface without solidification a heated path must be provided. The heating of this path requires the solidification of some magma. It is estimated that magma penetrates the lithosphere in about 5000 years and that the crack is lined by several hundred meters of frozen basalt.     

火山活动状态及其机理的初步探讨

根据火山喷发实例总结了火山喷发在不同阶段的活动状态,并探讨了可能的物理机理。火山活动从岩浆补给到岩浆喷发的物理过程可分为3个阶段:1)岩浆补给阶段,岩浆囊压力差或过剩压力的大小决定了火山活动是否休眠或扰动,岩浆补给速率对压力差起了决定性的作用;2)通道形成阶段,当过剩压力超过围岩破裂强度时,围岩开始破裂,之后水热活动起了重要的作用;3)岩浆运移与失稳喷发阶段,主要是岩浆运移与地壳盖层的相互作用与失稳的过程。文中还讨论了火山活动状态与火山喷发危险性等级之间的关系,7个危险性等级分别对应于火山活动的7种状态,即休眠、平静、扰动、动荡、临界、活动、灾变     

Magma mixing and magma plumbing systems in island arcs

Petrographic features of mixed rocks in island arcs, especially those originating by the mixing of magmas with a large compositional and temperature difference, such as basalt and dacite, suggest that the whole mixing process from their first contact to the final cooling (= eruption) has occurred continuously and in a relatively short time period. This period is probably less than several months, considerably shorter than the whole volcanic history. There may also be a prolonged quiescent interval, lasting longer than several days, between the magmas' contact and the mechanical mixing. This interval will allow the basic magma to cool and produce a semi-solidified boundary which is later disrupted by flow movements to produce basic inclusions.Mixing of magmas of contrasting chemical composition need not be the inevitable consequence of the contact of the magmas. It is, however, made more probable by forced convection caused by motive force such as the injection of a basic magma into an acidic magma chamber. A short interval between their initial contact and the final eruption requires that the acid magma chamber has a small volume, of the same order or less than that the introduced basic magma.The volcanic activity of Myoko volcano, central Japan, of the last 100,000 years shows alternate eruptions of hybrid andesite by mixing of basaltic and dacitic magmas, and non-mixed basalt to basaltic andesite. There was a repose period of 20,000 to 30,000 years between eruptions. The acidic chamber, eventually producing the mixed andesite activity, is formed during the repose period by the « in situ » solidification of the original basic magma against its wall. The volume of the chamber is very small, probably about 10^–2 km³. Basaltic magma with constant chemical composition is supplied to the shallow chamber from another deep seated basaltic chamber. The volume of the shallow magma chamber may be critical to the characteristics of volcanic activity and its products.     

Voluminous lava-like precursor to a major ash-flow tuff: low-column pyroclastic eruption of the Pagosa Peak Dacite, San Juan volcanic field, Colorado

The Pagosa Peak Dacite is an unusual pyroclastic deposit that immediately predated eruption of the enormous Fish Canyon Tuff (5000 km³) from the La Garita caldera at 28 Ma. The Pagosa Peak Dacite is thick (to 1 km), voluminous (>200 km³), and has a high aspect ratio (1:50) similar to those of silicic lava flows. It contains a high proportion (40–60%) of juvenile clasts (to 3–4 m) emplaced as viscous magma that was less vesiculated than typical pumice. Accidental lithic fragments are absent above the basal 5–10% of the unit. Thick densely welded proximal deposits flowed rheomorphically due to gravitational spreading, despite the very high viscosity of the crystal-rich magma, resulting in a macroscopic appearance similar to flow-layered silicic lava. Although it is a separate depositional unit, the Pagosa Peak Dacite is indistinguishable from the overlying Fish Canyon Tuff in bulk-rock chemistry, phenocryst compositions, and ⁴⁰Ar/³⁹Ar age.The unusual characteristics of this deposit are interpreted as consequences of eruption by low-column pyroclastic fountaining and lateral transport as dense, poorly inflated pyroclastic flows. The inferred eruptive style may be in part related to synchronous disruption of the southern margin of the Fish Canyon magma chamber by block faulting. The Pagosa Peak eruptive sources are apparently buried in the southern La Garita caldera, where northerly extensions of observed syneruptive faults served as fissure vents. Cumulative vent cross-sections were large, leading to relatively low emission velocities for a given discharge rate. Many successive pyroclastic flows accumulated sufficiently rapidly to weld densely as a cooling unit up to 1000 m thick and to retain heat adequately to permit rheomorphic flow. Explosive potential of the magma may have been reduced by degassing during ascent through fissure conduits, leading to fracture-dominated magma fragmentation at low vesicularity. Subsequent collapse of the 75×35 km² La Garita caldera and eruption of the Fish Canyon Tuff were probably triggered by destabilization of the chamber roof as magma was withdrawn during the Pagosa Peak eruption.     

Incipient magma chamber formation as a result of repetitive intrusions

An analytical solution for periodic magma intrusions in conduits was developed to study the onset of shallow magma chamber formation. The solution is based on determining when a repetitive series of intrusions can cause the wall rock of a conduit to reach its melt temperature. The results show that magma chamber formation in conduits is a strong function of the volume rate of intrusion and that magma chamber formation is likely when the intrusion rate exceeds 10^?3 km³/ yr. which agrees with observations by other investigators. Once this critical value of intrusion rate is reached, magma chambers are likely to begin forming after only a few intrusive pulses (less than ten). Results for both cylindrical conduits and dikes show cylindrical conduits are more favourable for the formation of shallow magma chambers.     

Magma mixing: petrological process and volcanological tool

Magma mixing is a widespread, if not universal igneous phenomenon of variable importance. The evidence for magma mixing is found primarily in glassy tephra; the consolidation of lava obscures the evidence. Inclusions of glass in big crystals in tephra, because of their greater range in composition compared to the whole rock and the residual glass, indicate that the big crystals were derived from separate systems which mixed together prior to and during eruption. The observed or reconstructed concentration of K₂O in inclusions of glass in large crystals represent the composition of the contaminant and host systems. Selective enrichment in K₂O during entrapment of melt by growing crystals is shown to be negligible. The weight percents of K₂O in host, contaminant and residual glass and bulk rock determine the proportions of contaminant and host required to yield either the residual glass or bulk rock. In several cases the proportion of contaminant required is substantially larger than the proportion of crystals in the hybrid magma; therefore, by heat budget argument, the contaminant was partly liquid when contamination began. In some tephra individual phenocrysts contain glasses which are more silicic toward the center of the crystal indicating that the crystal grew from a melt whose composition changed in the opposite sense to that expected for progressive solidification of a closed system. Space time associations of compositionally distinct glassy tephra with contaminated magmas suggest coexistence of basaltic and silicic melts within magma systems. Evidence of contamination is present in most tephra studied so far. Magma mixing appears to be the prevalent process whereby contamination occurs. Magma mixing seems to be particularly evident in systems where there is independent evidence for a vapor-saturated magma reservoir. Probably vapor saturation promotes mixing in magma systems. Magma mixing probably is an important mechanism of compositional diversification (differentiation) of volcanic rocks from continental margin and possibly other environments.Textural evidence of the onset of magma mixing can be related to disturbance of a complex reservoir immediately before ascent and eruption. Thus, conditions before mixing can be ascribed to the reservoir. In this way it is possible to learn about the reservoir: its composition, its diversity, its depth, its walls. It is also possible to learn about the causes of eruption: whether by increase in gas pressure due to either progressive consolidation, or heating from below by an injection of hot magma, or by encounter with ground water; whether by buoyant rise. Evaluation of these problems requires also a thorough knowledge of the chronology of particular eruptions. Thus, magma mixing is a useful volcanological tool.     

Effects of relative motion between gas and liquid on 1-dimensional steady flow in silicic volcanic conduits: 2. Origin of diversity of eruption styles

We investigate the origin of diversity of eruption styles in silicic volcanoes on the basis of a 1-dimensional steady conduit flow model that considers vertical relative motion between gas and liquid (i.e., vertical gas escape). The relationship between the assemblage of steady-state solutions in the conduit flow model and magma properties or geological conditions is expressed by a regime map in the parameter space of the ratio of liquid-wall friction force to liquid–gas interaction force (non-dimensional number ε), and a normalized conduit length Λ. The regime map developed in the companion paper shows that when ε is smaller than a critical value ε_cr, a solution of explosive eruption exists for a wide range of Λ, whereas an effusive solution exists only when Λ ~ 1. On the other hand, when ε > ε_cr, an effusive solution exists for a wide range of Λ. Diversity of eruption styles observed in nature is explained by the change in ε accompanied by the change in magma viscosity during magma ascent. As magma ascends, the magma viscosity increases because of gas exsolution and crystallization, leading to the increase in ε. For the viscosity of hydrous silicic magma at magma chamber, ε is estimated to be smaller than ε_cr, indicating that an explosive solution exists for wide ranges of geological parameters. When magma flow rate is small, the viscosity of silicic magma drastically increases because of extensive crystallization at a shallow level in the conduit. In this case, ε can be greater than ε_cr; as a result, a stable effusive solution co-exists with an explosive solution.     

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 role of magma chamber-fault interaction in caldera forming eruptions

This paper examines the role of the position and orientation of a regional fault in the roof of a magma chamber on stress distribution, mechanical failure, and dyking using 2D finite element numerical simulations. The study pertains to the magma chamber behavior in the relatively short time intervals of several hundreds to thousand of years. The magma chamber is represented as an elliptical inclusion (eccentricity, a/b = 0.12) at a relative depth, H/a, of 0.9. The fault has a 45° dip and is represented by a frictionless fracture. The temperature field in the host rock is calculated assuming a quasisteady-state thermal regime that develops through periodic episodes of magma supply. The rheology of the surrounding rocks is treated using viscoelasticity with temperature activated strain-rate dependent viscosity. Strain weakening of the rocks in the ductile zone is described within the frame of the Dynamic Power Law model . The magma pressure is coupled with the deformation of the rock mass hosting the chamber, including the fault. The variation of magma pressure in response to magma supply and chamber deformation is calculated in the elastic and viscoelastic regimes. The latter corresponds to slow filling, while the former represents a filling time much less than the viscous relaxation time scale. The resulting “equation of state” for the magma chamber couples the magma pressure with the chamber volume in the elastic regime, and with the filling rate for the viscoelastic regime. Analysis of stresses is used to predict dyke propagation conditions, and the mechanical failure of the chamber roof for different fault positions and magma overpressures. Results show that an outward dipping fault located on the periphery of the chamber roof hinders the propagation of dykes to the surface, causing magma to accumulate under the footwall of the fault. At high to moderate overpressures (30–40 MPa), the fault causes localized shear failure and chamber roof collapse that might lead to the first stage of a caldera-forming eruption.     

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.     

Origin of silicic volcanic rocks in Central Costa Rica: a study of a chemically variable ash-flow sheet in the Tiribí Tuff

Chemical heterogeneities of pumice clasts in an ash-flow sheet can be used to determine processes that occur in the magma chamber because they represent samples of magma that were erupted at the same time. The dominant ash-flow sheet in the Tiribí Tuff contains pumice clasts that range in composition from 55.1 to 69.2 wt% SiO₂. It covers about 820 km² and has a volume of about 25 km³ dense-rock equivalent (DRE). Based on pumice clast compositions, the sheet can be divided into three distinct chemical groupings: a low-silica group (55.1-65.6 wt% SiO₂), a silicic group (66.2-69.2 wt% SiO₂), and a mingled group (58.6-67.7 wt% SiO₂; all compositions calculated 100% anhydrous). Major and trace element modeling indicates that the low-silica magma represents a mantle melt that has undergone fractional crystallization, creating a continuous range of silica content from 55.1-65.6 wt% SiO₂. Eu/Eu*, MREE, and HREE differences between the two groups are not consistent with crystal fractionation of the low-silica magma to produce the silicic magma. The low-silica group and the silicic group represent two distinct magmas, which did not evolve in the same magma chamber. We suggest that the silicic melts resulted from partial melting of relatively hot, evolved calc-alkaline rocks that were previously emplaced and ponded at the base of an over-thickened basaltic crust. The mingled group represents mingling of the two magmas shortly before eruption. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located at http://dx.doi.org/10.1007/s00445-001-0188-8.     

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.     

Short time scales of magma-mixing processes prior to the 2011 eruption of Shinmoedake volcano, Kirishima volcanic group, Japan

We estimated time scales of magma-mixing processes just prior to the 2011 sub-Plinian eruptions of Shinmoedake volcano to investigate the mechanisms of the triggering processes of these eruptions. The sequence of these eruptions serves as an ideal example to investigate eruption mechanisms because the available geophysical and petrological observations can be combined for interpretation of magmatic processes. The eruptive products were mainly phenocryst-rich (28 vol%) andesitic pumice (SiO₂ 57 wt%) with a small amount of more silicic pumice (SiO₂ 62–63 wt%) and banded pumice. These pumices were formed by mixing of low-temperature mushy silicic magma (dacite) and high-temperature mafic magma (basalt or basaltic andesite). We calculated the time scales on the basis of zoning analysis of magnetite phenocrysts and diffusion calculations, and we compared the derived time scales with those of volcanic inflation/deflation observations. The magnetite data revealed that a significant mixing process (mixing I) occurred 0.4 to 3 days before the eruptions (pre-eruptive mixing) and likely triggered the eruptions. This mixing process was not accompanied by significant crustal deformation, indicating that the process was not accompanied by a significant change in volume of the magma chamber. We propose magmatic overturn or melt accumulation within the magma chamber as a possible process. A subordinate mixing process (mixing II) also occurred only several hours before the eruptions, likely during magma ascent (syn-eruptive mixing). However, we interpret mafic injection to have begun more than several tens of days prior to mixing I, likely occurring with the beginning of the inflation (December 2009). The injection did not instantaneously cause an eruption but could have resulted in stable stratified magma layers to form a hybrid andesitic magma (mobile layer). This hybrid andesite then formed the main eruptive component of the 2011 eruptions of Shinmoedake.     

Pillow talk

Three distinct types of pillows and pillow lava sequences with different modes of origin have been recognized in the extrusive sequences comprising the upper parts of ophiolite complexes that represent the mafic portion of the floor of an Early Cretaceous back-arc basin in southern Chile. One type of pillow formed by non-explosive submarine effusion. A second type formed by magmatic intrusion into pre-existing aquagene tuff formed by explosive eruption. The third type of pillow occurs within dikes, forming pillowed dikes, possibly as a result of vapor streaming within a cooling dike. Where studied in southern Chile, aquagene tuffs and intrusive pillows decrease and water-lain pillows increase in relative abundance from north to south. This variation corresponds with a north-to-south decrease in both the relative volume of extrusives to extensional dikes and the range and volume of differentiated rocks, suggesting a southward increase in rate of extension relative to rate of magma supply within the spreading ridges at which the ophiolites formed. In the northern part of the original basin where the rate of extension was small relative to the rate of magma supply, magma remained in magma chambers longer, resulting in a greater range and volume of differentiated rocks. The larger volume of more differentiated, cooler and more viscous magmas, in conjunction with the likelihood of more violent eruption of volatile-rich differentiates, may have been responsible for the large volume of aquagene tuff in the northern part of the original basin. These observations in southern Chile suggest that ophiolites which contain a great abundance of aquagene tuffs and intrusive pillow lavas formed in tectonic environments in which the rate of extension was small relative to the rate of magma supply (island arcs, embryonic marginal basins). Ophiolites with predominantly water-lain pillowed and massive lavas formed in tectonic environments in which the rate of extension was large relative to the rate of magma supply (mid-ocean ridges, mature back-arc basins). Thus geologic field data may supplement geochemical data as a tool in distinguishing the original igneous-tectonic environments in which ophiolites originate.