A model for the origin of large silicic magma chambers: precursors of caldera-forming eruptions

The relatively low rates of magma production in island arcs and continental extensional settings require that the volume of silicic magma involved in large catastrophic caldera-forming (CCF) eruptions must accumulate over periods of 10 ⁵ to 10 ⁶ years. We address the question of why buoyant and otherwise eruptible high-silica magma should accumulate for long times in shallow chambers rather than erupt more continuously as magma is supplied from greater depths. Our hypothesis is that the viscoelastic behavior of magma chamber wall rocks may prevent an accumulation of overpressure sufficient to generate rhyolite dikes that can propagate to the surface and cause an eruption. The critical overpressure required for eruption is based on the model of Rubin (1995a). An approximate analytical model is used to evaluate the controls on magma overpressure for a continuously or episodically replenished spherical magma chamber contained in wall rocks with a Maxwell viscoelastic rheology. The governing parameters are the long-term magma supply, the magma chamber volume, and the effective viscosity of the wall rocks. The long-term magma supply, a parameter that is not typically incorporated into dike formation models, can be constrained from observations and melt generation models. For effective wall-rock viscosities in the range 10 ¹⁸ to 10 ²⁰ Pa s ^–1, dynamical regimes are identified that lead to the suppression of dikes capable of propagating to the surface. Frequent small eruptions that relieve magma chamber overpressure are favored when the chamber volume is small relative to the magma supply and when the wall rocks are cool. Magma storage, leading to conditions suitable for a CCF eruption, is favored for larger magma chambers (>10 ² km ³) with warm wall rocks that have a low effective viscosity. Magma storage is further enhanced by regional tectonic extension, high magma crystal contents, and if the effective wall-rock viscosity is lowered by microfracturing, fluid infiltration, or metamorphic reactions. The long-term magma supply rate and chamber volume are important controls on eruption frequency for all magma chamber sizes. The model can explain certain aspects of the frequency, volume, and spatial distribution of small-volume silicic eruptions in caldera systems, and helps account for the large size of granitic plutons, their association with extensional settings and high thermal gradients, and the fact that they usually post-date associated volcanic deposits.     

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.     

Caldera resurgence during magma replenishment and rejuvenation at Valles and Lake City calderas

A key question in volcanology is the driving mechanisms of resurgence at active, recently active, and ancient calderas. Valles caldera in New Mexico and Lake City caldera in Colorado are well-studied resurgent structures which provide three crucial clues for understanding the resurgence process. (1) Within the limits of ⁴⁰Ar/³⁹Ar dating techniques, resurgence and hydrothermal alteration at both calderas occurred very quickly after the caldera-forming eruptions (tens of thousands of years or less). (2) Immediately before and during resurgence, dacite magma was intruded and/or erupted into each system; this magma is chemically distinct from rhyolite magma which was resident in each system. (3) At least 1?km of structural uplift occurred along regional and subsidence faults which were closely associated with shallow intrusions or lava domes of dacite magma. These observations demonstrate that resurgence at these two volcanoes is temporally linked to caldera subsidence, with the upward migration of dacite magma as the driver of resurgence. Recharge of dacite magma occurs as a response to loss of lithostatic load during the caldera-forming eruption. Flow of dacite into the shallow magmatic system is facilitated by regional fault systems which provide pathways for magma ascent. Once the dacite enters the system, it is able to heat, remobilize, and mingle with residual crystal-rich rhyolite remaining in the shallow magma chamber. Dacite and remobilized rhyolite rise buoyantly to form laccoliths by lifting the chamber roof and producing surface resurgent uplift. The resurgent deformation caused by magma ascent fractures the chamber roof, increasing its structural permeability and allowing both rhyolite and dacite magmas to intrude and/or erupt together. This sequence of events also promotes the development of magmatic–hydrothermal systems and ore deposits. Injection of dacite magma into the shallow rhyolite magma chamber provides a source of heat and magmatic volatiles, while resurgent deformation and fracturing increase the permeability of the system. These changes allow magmatic volatiles to rise and meteoric fluids to percolate downward, favouring the development of hydrothermal convection cells which are driven by hot magma. The end result is a vigorous hydrothermal system which is driven by magma recharge.     

Density and viscosity gradients in zoned magma chambers, and their influence withdrawal dynamics

The liquid being sampled from a draining reservoir of density-stratified fluid, such as an erupting zoned magma chamber, is derived from a relatively thin withdrawal layer adjacent to the level of the chamber outlet. This is a consequence of the buoyancy force associated with the density gradient inhibiting vertical motion so that the opportunity for widely separated density levels (compositions) to be tapped and mingled syneruptively is suppressed.Density gradients in zoned chambers of 0.02 – 10 kgm⁻³/m are suggested by data from caldera-forming eruptions. Viscosity gradients can be specified for a given density gradient using calculated relationships between viscosity and density. Published compositional and geothermometric data are used to show that zoned high-silica rhyolites decrease in viscosity upward because of the roofward concentration of dissolved volatiles. Other zoned calc-alkaline magmas increase in viscosity upward because of decreasing temperature and concentration of network modifying cations.A method is developed of calculating the scale of the withdrawal layer thickness, δ, for given kinematic viscosity, eruption rate, and density and viscosity gradients. The method is systematized by the identification of specific flow regimes describing the action of either viscous or inertial forces in balancing the buoyancy force. Thin withdrawal layers are favoured by small eruption rates, small viscosity, and by large density gradients. For particularly steep density gradients, however, the consequently large viscosity gradient plays a role in determining the withdrawal layer thickness. Withdrawal layer thicknesses of the order of 100 m are calculated for typical pyroclastic eruptions of zoned acid magma, and are mostly independent of the viscosity gradient.The vertical scale at which a zoned chamber is instantaneously being tapped during an eruption is equal to the scale of the withdrawal layer thickness. Thus, an eruption that causes collapse of a caldera block through a height that is less than that of the withdrawal layer scale will produce magmas from deeper levels than that to which the chamber roof sinks. In this case the eruption is said to oversample the chamber with respect to the amount of caldera collapse and will produce an essentially constant range of compositions throughout. Alternatively, if the caldera collapse distance is much greater than δ then the selective withdrawal process leads to successive levels of the chamber being “skimmed off” (on a scale δ). This allows the compositional stratigraphy of the chamber to be inverted by the eruptive process, with little opportunity for syneruptive mixing between diverse magma compositions. The geological record shows that most calderas associated with zoned magmas collapsed through vertical distances in excess of 100 m (the characteristic estimate for δ) and, in agreement with our modelling of selective withdrawal, show smooth correlations between composition, or temperature, and the order of eruption.     

Ignimbrites of the Cerro Galan caldera, NW Argentina

The 35 × 20 km Cerro Galán resurgent caldera is the largest post-Miocene caldera so far identified in the Andes. The Cerro Galán complex developed on a late pre-Cambrian to late Palaeozoic basement of gneisses, amphibolites, mica schists and deformed phyllites and quartzites. The basement was uplifted in the early Miocene along large north-south reverse faults, producing a horst-and-graben topography. Volcanism began in the area prior to 15 Ma with the formation of several andesite to dacite composite volcanoes. The Cerro Galán complex developed along two prominent north-south regional faults about 20 km apart. Dacitic to rhyodacitic magma ascended along these faults and caused at least nine ignimbrite eruptions in the period 7-4 Ma (K-Ar determinations). These ignimbrites are named the Toconquis Ignimbrite Formation. They are characterised by the presence of basal plinian deposits, many individual flow units and proximal co-ignimbrite lag breccias. The ignimbrites also have moderate to high macroscopic pumice and lithic contents and moderate to low crystal contents. Compositionally banded pumice occurs near the top of some units. Many of the Toconquis eruptions occurred from vents along a north-south line on the western rim of the young caldera. However, two of the ignimbrites erupted from vents on the eastern margin. Lava extrusions occurred contemporaneously along these north-south lines. The total D.R.E. volume of Toconquis ignimbrite exceeds 500 km³.Following a 2-Ma dormant period a single major eruption of rhyodacitic magma formed the 1000-km³ Cerro Galán ignimbrite and the caldera. The ignimbrite (age 2.1 Ma on Rb-Sr determination) forms a 30–200-m-thick outflow sheet extending up to 100 km in all directions from the caldera rim. At least 1.4 km of welded intracaldera ignimbrite also accumulated. The ignimbrite is a pumice-poor, crystal-rich deposit which contains few lithic clasts. No basal plinian deposit has been identified and proximal lag breccias are absent. The composition of pumice clasts is a very uniform rhyodacite which has a higher SiO₂ content but a lower K₂O content than the Toconquis ignimbrites. Preliminary data indicate no evidence for compositional zonation in the magma chamber. The eruption is considered to have been caused by the catastrophic foundering of a cauldron block into the magma chamber.Post-caldera extrusions occurred shortly after eruption along both the northern extension of the eastern boundary fault and the western caldera margin. Resurgence also occurred, doming up the intracaldera ignimbrite and sedimentary fill to form the central mountain range. Resurgent doming was centred along the eastern fault and resulted in radial tilting of the ignimbrite and overlying lake sediments.     

Reconstruction of a major caldera-forming eruption from pyroclastic deposit characteristics: Kos Plateau Tuff, eastern Aegean Sea

The 161 ka explosive eruption of the Kos Plateau Tuff (KPT) ejected a minimum of 60 km³ of rhyolitic magma, a minor amount of andesitic magma and incorporated more than 3 km³ of vent- and conduit-derived lithic debris. The source formed a caldera south of Kos, in the Aegean Sea, Greece. Textural and lithofacies characteristics of the KPT units are used to infer eruption dynamics and magma chamber processes, including the timing for the onset of catastrophic caldera collapse.The KPT consists of six units: (A) phreatoplinian fallout at the base; (B, C) stratified pyroclastic-density-current deposits; (D, E) volumetrically dominant, massive, non-welded ignimbrites; and (F) stratified pyroclastic-density-current deposits and ash fallout at the top. The ignimbrite units show increases in mass, grain size, abundance of vent- and conduit-derived lithic clasts, and runout of the pyroclastic density currents from source. Ignimbrite formation also corresponds to a change from phreatomagmatic to dry explosive activity. Textural and lithofacies characteristics of the KPT imply that the mass flux (i.e. eruption intensity) increased to the climax when major caldera collapse was initiated and the most voluminous, widespread, lithic-rich and coarsest ignimbrite was produced, followed by a waning period. During the eruption climax, deep basement lithic clasts were ejected, along with andesitic pumice and variably melted and vesiculated co-magmatic granitoid clasts from the magma chamber. Stratigraphic variations in pumice vesicularity and crystal content, provide evidence for variations in the distribution of crystal components and a subsidiary andesitic magma within the KPT magma chamber. The eruption climax culminated in tapping more coarsely crystal-rich magma. Increases in mass flux during the waxing phase is consistent with theoretical models for moderate-volume explosive eruptions that lead to caldera collapse.     

The intensity of plinian eruptions

Peak intensities (magma discharge rate) of 45 Pleistocene and Holocene plinian eruptions have been inferred from lithic dispersal patterns by using a theoretical model of pyroclast fallout from eruption columns. Values range over three orders of magnitude from 1.6 × 10⁶ to 1.1 × 10⁹ kg/s. Magnitudes (total erupted mass) also vary over about three orders of magnitude from 2.0 × 10¹¹ to 6.8 × 10¹⁴ kg and include several large ignimbrite-forming events with associated caldera formation. Intensity is found to be positively correlated with the magnitude when total erupted mass (tephra fall, surges and pyroclastic flows) is considered. Initial plinian fall phases with intensities in excess of 2.0 × 10⁸ kg/s typically herald the onset of major pyroclastic flow generation and subsequent caldera collapse. During eruptions of large magnitude, the transition to pyroclastic flows is likely to be the result of high intensity, whereas the generation of pyroclastic flows in small magnitude eruptions may occur more often by reduction of magmatic volatile content or some transient change in magma properties. The correlation between plinian fall intensity and total magnitude suggests that the rate of magma discharge is related to the size of the chamber being tapped. A simple model is presented to account for the variation in intensity by progressive enlargement of conduits and vents and excess pressure at the chamber roof caused by buoyant forces acting on the chamber as it resides in the crust. Both processes are fundamentally linked to the absolute size of the pre-eruption reservoir. The data suggest that sustained eruption column heights (i.e. magma discharge rates) are indicators of eventual eruption magnitude, and perhaps eruptive style, and thus are key parameters to monitor in order to assess the temporal evolution of plinian eruptions.     

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.     

Formation and development of normal-fault calderas and the initiation of large explosive eruptions

The ring fractures that form most collapse calderas are steeply inward-dipping shear fractures, i.e., normal faults. At the surface of the volcano within which the caldera fault forms, the tensile and shear stresses that generate the normal-fault caldera must peak at a certain radial distance from the surface point above the center of the source magma chamber of the volcano. Numerical results indicate that normal-fault calderas may initiate as a result of doming of an area containing a shallow sill-like magma chamber, provided that the area of doming is much larger than the cross-sectional area of the chamber and that the internal excess pressure in the chamber is smaller than that responsible for doming. This model is supported by the observation that many caldera collapses are preceded by a long period of doming over an area much larger than that of the subsequently formed caldera. When the caldera fault does not slip, eruptions from calderas are normally small. Nearly all large explosive eruptions, however, are associated with slip on caldera faults. During dip slip on, and doming of, a normal-fault caldera, the vertical stress on part of the underlying chamber suddenly decreases. This may lead to explosive bubble growth in this part of the magma chamber, provided its magma is gas rich. This bubble growth can generate an excess fluid pressure that is sufficiently high to drive a large fraction of the magma out of the chamber during an explosive eruption. Received: 2 January 1997 / Accepted: 22 April 1998     

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.     

Mechanical models for correlation of ring-fracture eruptions at Pantelleria,Strait of Sicily,with glacial sea-level drawdown

Recent K-Ar dating of eruptions at Pantelleria, a peralkaline volcanic island in the Strait of Sicily, shows a correlation between eruption of pantellerite lavas from caldera ring fractures and low stands of sea level as determined from ¹⁸O stratigraphy. Post-caldera pantellerite lavas associated with an 114-ky-old caldera erupted along the ring-fracture zone during a major low stand of sea level at about 67 Ka. The most recent episode of lava-flow emplacement began about 20 ky ago during the last glacial maximum. Magma vented along the ring fault of a 45-ky-old caldera, from fractures radial to the caldera, and along faults formed by intracaldera trapdoor uplift. Two mechanical models based on elasticity theory are presented to explain the correlation of post-caldera ring-fracture eruptions at Pantelleria with lowering of sea level. A simple analysis of a bending circular plate of thickness,T _r, and radius,R, representing the magma-chamber roof block, shows that tensile stress is concentrated by a factor of 0.75R ²/T _r ² at the lower perimeter of the plate when sea level drops. Stress changes may be even greater ifT _r is effectively less than the stratigraphic thickness due to layering of rocks in the roof block. Calculated stress changes due to a 100-m drawdown of sea level are similar in magnitude to stresses associated with dike propagation. More realistic model geometries, including different chamber shapes, a conical volcanic edifice, and sea-level drawdown beyond the surface projection of the magma chamber, were tested using the boundary-element method. Lowering sea level generates a horizontal tensile stress above the chamber, even when sea water is removed outboard of the magma chamber. For some chamber geometries the magnitude of the tensile stress maximum is greater than the 1 MPa pressure of the 100 m of removed water and is of the right order of magnitude for dike propagation. Dikes initiated by the change of the stress field may originate and propagate along fractures inboard of the chamber margin. The magnitudes of tensile maxima along the top of the chamber decrease as original sea level is moved outboard of the chamber margin and as the chamber thickness decreases. When the depth to the top of the magma chamber reaches a critical value, dependent on chamber geometry, the propagation of dikes to the surface is inhibited.     

Stratigraphy of the Toba Tuffs and the evolution of the Toba Caldera Complex,Sumatra, Indonesia

During the past 1.2 m.y., a magma chamber of batholithic proportions has developed under the 100 by 30 km Toba Caldera Complex. Four separate eruptions have occurred from vents within the present collapse structure, which formed from eruption of the 2800 km³ Youngest Toba Tuff (YTT) at 74 ka. Eruption of the three older Toba Tuffs alternated from calderas situated in northern and southern portions of the present caldera. The northern caldera apparently developed upon a large andesitic stratovolcano. The calderas associated with the three older tuffs are obscured by caldera collapse and resurgence resulting from eruption of the YTT. Samosir Island and the Uluan Block are two sides of a single resurgent dome that has resurged since eruption of the YTT. Samosir Island is composed of thick YTT caldera fill, whereas the Uluan Block consists mainly of the Oldest Toba Tuff (OTT). In the past 74000 years lava domes have been extruded on Samosir Island and along the caldera's western ring fracture. This part of the ring fracture is the site of the only current activity at Toba: updoming and fumarolic activity. The Toba eruptions document the growth of the laterally continuous magma body which eventually erupted the YTT. Repose periods between the four Toba Tuffs range between 0.34 and 0.43 m.y. and give insights into pluton emplacement and magmatic evolution at Toba.     

Non-explosive,constructional evolution of the ice-filled caldera at Volcán Sollipulli,Chile

 A radar and gravity survey of the ice-filled caldera at Volcán Sollipulli, Chile, indicates that the intra-caldera ice has a thickness of up to 650 m in its central part and that the caldera harbours a minimum of 6 km³ of ice. Reconnaissance geological observations show that the volcano has erupted compositions ranging from olivine basalt to dacite and have identified five distinct volcanic units in the caldera walls. Pre- or syn-caldera collapse deposits (the Sharkfin pyroclastic unit) comprise a sequence which evolved from subglacial to subaerial facies. Post-caldera collapse products, which crop out along 17 of the 20 km length of the caldera wall, were erupted almost exclusively along the caldera margins in the presence of a large body of intra-caldera ice. The Alpehué crater, formed by an explosive eruption between 2960 and 2780 a. BP, in the southwest part of the caldera is shown to post date formation of the caldera. Sollipulli lacks voluminous silicic pyroclastic rocks associated with caldera formation and the collapse structure does not appear to be a consequence of a large-magnitude explosive eruption. Instead, lateral magma movement at depth resulting in emptying of the magma chamber may have generated the caldera. The radar and gravity data show that the central part of the caldera floor is flat but, within a few hundred metres of the caldera walls, the floor has a stepped topography with relatively low-density rock bodies beneath the ice in this region. This, coupled with the fact that most of the post-caldera eruptions have taken place along the caldera walls, implies that the caldera has been substantially modified by subglacial marginal eruptions. Sollipulli caldera has evolved from a collapse to a constructional feature with intra-caldera ice playing a major role. The post-caldera eruptions have resulted in an increase in height of the walls and concomitant deepening of the caldera with time. Received: 12 June 1995 / Accepted: 7 December 1995     

Stratigraphy and eruptive dynamics of a pulsating Plinian eruption of Somma-Vesuvius: the Pomici di Mercato (8900 years B.P.)

New volcanological studies allow reconstruction of the eruption dynamics of the Pomici di Mercato eruption (ca 8,900 cal. yr B.P.) of Somma-Vesuvius. Three main Eruptive Phases are distinguished based on two distinct erosion surfaces that interrupt stratigraphic continuity of the deposits, indicating that time breaks occurred during the eruption. Absence of reworked volcaniclastic deposits on top of the erosion surfaces suggests that quiescent periods between eruptive phases were short perhaps lasting only days to weeks. Each of the Eruptive Phases was characterised by deposition of alternating fall and pyroclastic density current (PDC) deposits. The fallout deposits blanketed a wide area toward the east, while the more restricted PDC deposits inundated the volcano slopes. Eruptive dynamics were driven by brittle magmatic fragmentation of a phonolitic magma, which, because of its mechanical fragility, produced a significant amount of fine ash. External water did not significantly contribute either to fragmentation dynamics or to mechanical energy release during the eruption. Column heights were between 18 and 22 km, corresponding to mass discharge rates between 1.4 and 6 × 10⁷ kg s⁻¹. The estimated on land volume of fall deposits ranges from a minimum of 2.3 km³ to a maximum of 7.4 km³. Calculation of physical parameters of the dilute pyroclastic density currents indicates speeds of a few tens of m s⁻¹ and densities of a few kg m⁻³ (average of the lowermost 10 m of the currents), resulting in dynamic pressures lower than 3 kPa. These data suggest that the potential impact of pyroclastic density currents of the Pomici di Mercato eruption was smaller than those of other Plinian and sub-Plinian eruptions of Somma-Vesuvius, especially those of 1631 AD and 472 AD (4–14 kPa), which represent reference values for the Vesuvian emergency plan. The pulsating and long-lasting behaviour of the Pomici di Mercato eruption is unique in the history of large explosive eruptions of Somma-Vesuvius. We suggest an eruptive scheme in which discrete magma batches rose from the magma chamber through a network of fractures. The injection and rise of the different magma batches was controlled by the interplay between magma chamber overpressure and local stress. The intermittent discharge of magma during a large explosive eruption is unusual for Somma-Vesuvius, as well as for other volcanoes worldwide, and yields new insights for improving our knowledge of the dynamics of explosive eruptions.     

40Ar/39Ar Ages and stratigraphy of the Latera caldera,Italy

The Latera caldera is a well-exposed volcano where more than 8 km³ of mafic silica-undersaturated potassic lavas, scoria and felsic ignimbrites were emplaced between 380 and 150 ka. Isotopic ages obtained by ⁴⁰Ar/³⁹Ar analysis of single sanidine crystals indicate at least four periods of explosive eruptions from the caldera. The initial period of caldera eruptions began at 232 ka with emplacement of trachytic pumice fallout and ignimbrite. They were closely followed by eruption of evolved phonolitic magma. After roughly 25 ky, several phonolitic ignimbrites were deposited, and they were followed by phreatomagmatic eruptions that produced trachytic ignimbrites and several smaller ash-flow units at 191 ka. Compositionally zoned magma then erupted from the northern caldera rim to produce widespread phonolitic tuffs, tephriphonolitic spatter, and scoria-bearing ignimbrites. After 40 ky of mafic surge deposit and scoria cone development around the caldera rim, a compositionally zoned pumice sequence was emplaced around a vent immediately northwest of the Latera caldera. This activity marks the end of large-scale explosive eruptions from the Latera volcano at 156 ka.     

Eruptive patterns for the volcano piton de la fournaise, Reunion Island

A study of the historic record of activity of Piton de la Fournaise has revealed a cyclic pattern of eruption involving effusion of oceanite lava from major-flank centers every 20–40 years. Calculated volumes of the recent lava flows and pyroclastic ejecta have established an effusion rate of 3.9 m³ s⁻¹ since 1931 and 6.2 m³ s⁻¹ since 1951. Flank eruptions outside the present caldera define a distribution maximum which is expected to correlate with the depth range of a high-level magma reservoir.A model has been constructed which requires replenishment of a high-level magma chamber at a constant rate and regular eruption from summit and minor-flank centers, acting as “safety valves” to the magma chamber; when the magma chamber reaches its maximum expansion, a major-flank outburst of oceanitic lava occurs.The fact that calculated effusion rates are not consistent with radiometric dates implies an increase in effusion volume with time for the volcano.     

Geology of the peralkaline volcano at Pantelleria,Strait of Sicily

Situated in a submerged continental rift, Pantelleria is a volcanic island with a subaerial eruptive history longer than 300 Ka. Its eruptive behavior, edifice morphologies, and complex, multiunit geologic history are representative of strongly peralkaline centers. It is dominated by the 6-km-wide Cinque Denti caldera, which formed ca. 45 Ka ago during eruption of the Green Tuff, a strongly rheomorphic unit zoned from pantellerite to trachyte and consisting of falls, surges, and pyroclastic flows. Soon after collapse, trachyte lava flows from an intracaldera central vent built a broad cone that compensated isostatically for the volume of the caldera and nearly filled it. Progressive chemical evolution of the chamber between 45 and 18 Ka ago is recorded in the increasing peralkalinity of the youngest lava of the intracaldera trachyte cone and the few lavas erupted northwest of the caldera. Beginning about 18 Ka ago, inflation of the chamber opened old ring fractures and new radial fractures, along which recently differentiated pantellerite constructed more than 25 pumice cones and shields. Continued uplift raised the northwest half of the intracaldera trachyte cone 275 m, creating the island's present summit, Montagna Grande, by trapdoor uplift. Pantellerite erupted along the trapdoor faults and their hingeline, forming numerous pumice cones and agglutinate sheets as well as five lava domes. Degassing and drawdown of the upper pantelleritic part of a compositionally and thermally stratified magma chamber during this 18-3-Ka episode led to entrainment of subjacent, crystal-rich, pantelleritic trachyte magma as crenulate inclusions. Progressive mixing between host and inclusions resulted in a secular decrease in the degree of evolution of the 0.82 km³ of magma erupted during the episode.The 45-Ka-old caldera is nested within the La Vecchia caldera, which is thought to have formed around 114 Ka ago. This older caldera was filled by three widespread welded units erupted 106, 94, and 79 Ka ago. Reactivation of the ring fracture ca. 67 Ka ago is indicated by venting of a large pantellerite centero and a chain of small shields along the ring fault. For each of the two nested calderas, the onset of postcaldera ring-fracture volcanism coincides with a low stand of sea level.Rates of chemical regeneration within the chamber are rapid, the 3% crystallization/Ka of the post-Green Tuff period being typical. Highly evolved pantellerites are rare, however, because intervals between major eruptions (averaging 13–6 Ka during the last 190 Ka) are short. Benmoreites and mugearites are entirely lacking. Fe-Ti-rich alkalic basalts have erupted peripherally along NW-trending lineaments parallel to the enclosing rift but not within the nested calderas, suggesting that felsic magma persists beneath them. The most recent basaltic eruption (in 1891) took place 4 km northwest of Pantelleria, manifesting the long-term northwestward migration of the volcanic focus. These strongly differentiated basalts reflect low-pressure fractional crystallization of partial melts of garnet peridotite that coalesce in small magma reservoirs replenished only infrequently in this continental rift environment.     

1998 Eruption at Volcán Cerro Azul, Galápagos Islands: I. Syn-Eruptive Petrogenesis

The 1998 eruption of Volcán Cerro Azul in the Galápagos Islands produced two intra-caldera vents and a flank vent that erupted more than 1.0×10⁸ m³ of lava. Lava compositions changed notably during the 5-week eruption, and contemporaneous eruptions in the caldera and on the flank produced different compositions. Lavas erupted from the flank vent range from 6.3 to 14.1% MgO, nearly the entire range of MgO contents previously reported from the volcano. On-site monitoring of eruptive activity is linked with petrogenetic processes such that geochemical variations are evaluated in a temporal context. Lavas from the 1998 eruption record two petrogenetic stages characterized by progressively more mafic lavas as the eruption proceeded. Crystal compositions, whole rock major and trace element compositions, and isotope ratios indicate that early lavas are the product of mixing between 1998 magma and remnant magma of the 1979 eruption. Intra-caldera lavas and later lavas have no 1979 signature, but were produced by the 1998 magma incorporating olivine and clinopyroxene xenocrysts. Thus, early magma petrogenesis is characterized by mixing with the 1979 magma, followed by the magma progressively entraining wehrlite cumulate mush.Editorial Responsibility: M.R. Carroll     

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.     

The Negra Muerta Volcanic Complex, southern Central Andes: geochemical characteristics and magmatic evolution of an episodically active volcanic centre

This petrologic analysis of the Negra Muerta Volcanic Complex (NMVC) contributes to understanding the magmatic evolution of eruptive centres associated with prominent NW-striking fault zones in the southern Central Andes. Specifically, the geochemical characteristics and magmatic evolution of the two eruptive episodes of this Complex are analysed. The first one occurred as an explosive eruption at 9 Ma and is represented by a strongly welded, fiamme-rich, andesitic to dacitic ignimbrite deposit. The second commenced with an eruption of a rhyolitic ignimbrite at 7.6 Ma followed by effusive discharge of hybrid lavas at 7.3 Ma and by emplacement of andesitic to rhyodacitic dykes and domes. Both explosive and effusive eruptions of the second episode occurred within a short time span, but geochemical interpretations permit consideration of the existence of different magmas interacting in the same magma chamber. Our model involves an andesitic recharge into a partially cooled rhyolitic magma chamber, pressurising the magmatic system and triggering explosive eruption of rhyolitic magma. Chemical or mechanical evidence for interaction between the rhyolitic and andesitic magma in the initial stages are not obvious because of their difference in composition, which could have been strong enough to inhibit the interaction between the two magmas. After the initial explosive stages of the eruption at 7.6 Ma, the magma chamber become more depressurised and the most mafic magma settled in compositional layers by fractional crystallisation. Restricted hybridisation occurred and was effective between adjacent and thermally equivalent layers close to the top of the magma chamber. At 7.3 Ma, increments of caldera formation were accompanied by effusive discharge of hybrid lavas through radially disposed dykes whereby andesitic magma gained in importance toward the end of this effusive episode in the central portion of the caldera. Assimilation during turbulent ascent (ATA) is invoked to explain a conspicuous reversed isotopic signature (⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd) in the entire volcanic series. Therefore, the 7.6 to 7.3 Ma volcanic rocks of the NMVC resulted from synchronous and mutually interacting petrological processes such as recharge, fractional crystallization, hybridisation, and Assimilation during Turbulent Ascent (ATA).Geochemical characteristics of both volcanic episodes show diverse type and/or depth in the sources and variable influence of upper crustal processes, and indicate a recurrence in the magma-forming conditions. Similarly, other minor volcanic centres in the transversal volcanic belts of the Central Andes repeated their geochemical signatures throughout the Miocene.