Effects of relative motion between gas and liquid on 1-dimensional steady flow in silicic volcanic conduits: 1. An analytical method

We investigate the effects of vertical relative motion between gas and liquid on eruption styles by formulating a model for 1-dimensional steady flow in volcanic conduits. As magma ascends and decompresses, volatiles exsolve and volume fraction of gas increases. As a result, magma fragmentation occurs and the flow changes from bubbly flow to gas-pyroclast flow. In our model, a transitional region (‘permeable flow region’) is introduced between the bubbly flow region and the gas-pyroclast flow region. In this region, both the gas and the liquid are continuous phases, allowing the efficient vertical escape of gas through the permeable structure. We describe the features of conduit flow with relative motion of gas and liquid using non-dimensional numbers α, γ and ε. The parameter α represents the ratio of effects of wall friction to gravitational load, and is proportional to magma flow rate. The parameter γ represents the degree of decompression for the gas-pyroclast flow to reach the sound velocity at α = 1, and is proportional to r_c²/μ for given magma temperature and initial volatile content, where r_c is conduit radius and μ is liquid viscosity. The parameter ε is defined as the ratio of liquid–wall friction force to liquid–gas interaction force in the permeable flow region, and represents the efficiency of gas escape from magma. The values of γ and ε are determined only by magma properties and geological conditions such as liquid viscosity, magma permeability and conduit radius. We formulate a 1-dimensional steady-state conduit flow model to find non-dimensional magma flow rate α as a function of magma properties and geological conditions (e.g., γ and ε) under given boundary conditions. When the relative motion is taken into account with the assumption that magma fragmentation occurs when the gas volume fraction reaches some critical values, the pressure at the fragmentation level (P_f) decreases as the magma flow rate (α) decreases or the efficiency of gas escape (ε) increases, because gas escape suppresses the increase in the gas volume fraction accompanying magma ascent. When ε is so large that P_f is below the atmospheric pressure (P_a), the flow reaches the vent before fragmentation at low α. On the other hand, when ε is so small that P_f is greater than P_a, the flow reaches the vent after fragmentation at high α. These steady-state solutions of the flow at low and high α correspond to effusive and explosive eruptions, respectively. We present a graphical method to systematically find α. On the basis of the graphical method, a simple regime map showing the relationship between the assemblage of the solutions of conduit flow and the magma properties or the geological conditions is obtained.     

Shifting styles of basaltic explosive activity during the 2002–03 eruption of Mt. Etna, Italy

The 2002–03 flank eruption of Etna was characterized by two months of explosive activity that produced copious ash fallout, constituting a major source of hazard and damage over all eastern Sicily. Most of the tephra were erupted from vents at 2750 and 2800 m elevation on the S flank of the volcano, where different eruptive styles alternated. The dominant style of explosive activity consisted of discrete to pulsing magma jets mounted by wide ash plumes, which we refer to as ash-rich jets and plumes. Similarly, ash-rich explosive activity was also briefly observed during the 2001 flank eruption of Etna, but is otherwise fairly uncommon in the recent history of Etna. Here, we describe the features of the 2002–03 explosive activity and compare it with the 2001 eruption in order to characterize ash-rich jets and plumes and their transition with other eruptive styles, including Strombolian and ash explosions, mainly through chemical, componentry and morphology investigations of erupted ash. Past models explain the transition between different styles of basaltic explosive activity only in terms of flow conditions of gas and liquid. Our findings suggest that the abundant presence of a solid phase (microlites) may also control vent degassing and consequent magma fragmentation and eruptive style. In fact, in contrast with the Strombolian or Hawaiian microlite-poor, fluidal, sideromelane clasts, ash-rich jets and plumes produce crystal-rich tachylite clasts with evidence of brittle fragmentation, suggesting that high groundmass crystallinity of the very top part of the magma column may reduce bubble movement while increasing fragmentation efficiency.     

Three types of pyroclastic currents and their deposits: examples from the Vulsini Volcanoes, Italy

This paper deals with ground-hugging, gas–pyroclast currents from explosive volcanic eruptions and their deposits. Key field observations and laboratory determinations are proposed to relate specific deposit types with flow regimes and particle concentration in the transport and depositional systems. Three relevant flow scenarios and corresponding deposit types have been recognized from a survey of pyroclastic successions of the Vulsini Volcanic District (central Italy): (1) dilute, turbulent, pyroclastic currents producing normally or multiply graded beds by direct suspension sedimentation; (2) concentrated bedload regions beneath suspension currents, depositing inversely graded beds by traction carpet sedimentation; (3) self-sustained, high particle concentration, laminar, mass flows developing massive, poorly sorted bodies, with opposite grading of coarse lithic and pumice clasts, overlying fine-grained, inversely graded, basal layers. Main distinguishing criteria include the occurrence and pattern of clast grading, clast–thickness relationships, grain size, ash matrix componentry and pyroclast size–density relationships. Downcurrent and temporal transitions among identified flow scenarios are likely to occur for changing energy conditions and gas–pyroclast ratio both on regional and local scales. The nature and efficiency of magma fragmentation, volatile content, conduit geometry (which determine the characteristics of the erupted mixture and possible lateral blast component at the vent), and the angle of incidence of the column collapse, are suggested as the main factors controlling the generation of one type over the other at flow inception. Dilute, fine-grained, overpressured eruption clouds are thought to favor the formation of low particle concentration turbulent currents. Column collapse over slightly inclined volcano slopes, causing a high degree of compression of the collapsing mixture and of gas expulsion, would favor the generation of high particle concentration pyroclastic currents.     

Thermal, seismic and infrasonic evidences of variable degassing rates at Stromboli volcano

During June 1999, we measured the amplitude and rate (number of events per second) at which gas exited the vent at Stromboli volcano as discrete gas bursts or puffs. This allowed us to identify two styles of gas burst (puffing) activity. The first is characterized by frequent, rapidly rising puffs, the second by less frequent, slowly rising puffs. Each style persisted over 5–40-min-long durations and was associated with a high and low number of strombolian explosions per hour, respectively. Each period was also associated with characteristic delay times between the arrival of the infrasonic and thermal signals during strombolian explosions; the delays were longer during vigorous puffing periods. To explain our observations, we propose a model in which the degassing process cycles between vigorous and weak degassing phases. During vigorous degassing phases, bubble layers ascend the conduit at a frequency of 0.5–1.0 s⁻¹. This high degassing level reflects a gas-rich magma column and leads to an increased rate in the formation of shallow foams and, hence, an increase in puffing and explosive activity, as well as a higher free surface level and/or gas jet velocity. During weak phases, bubble layers ascend the conduit at a reduced frequency of 0.2–0.3 s⁻¹. During such times the magma column is poor in gas. This leads to a decreased rate of foam layer formation and hence a reduction in puffing and explosive activity, as well as a lower free surface level and/or gas jet velocity. Variations in puffing activity can thus be used to track changes in the rate at which the shallow system is supplied by fresh, gas-rich magma. Our observations indicate that the two degassing styles last from 5 to 40 min and that the switch from one to the other occurs over a matter of minutes.     

Source mechanisms of volcanic explosion earthquakes: single force and implosive sources

Volcanic explosion earthquakes accompanying explosive eruptions are viewed as representing an abrupt release of pressurized magma beneath a volcano through a vent. I examine the source mechanisms of explosion earthquakes assuming two extreme cases of magma property. When magma is assumed to be a perfect gas which migrates upward isentropically through a vent acting as a nozzle, the seismic source is expressed as a single force. On the other hand, when magma flow is assumed to be incompressible, an implosive source is dominant. Both source models predict that the seismic magnitude is proportional to the cross-sectional area of the vent. I investigate observed seismic magnitudes for different volcanoes based on data in published papers and reports. The results show that the seismic magnitude of the largest event for each volcano is essentially proportional to the cross-sectional area of the vent for vent radii ranging from 10 to 600 m. The consistency found between the theoretical prediction and the observed relation suggests that the vent area plays a substantial role in the magnitude of explosion earthquakes. I further estimate that initial pressures in the reservoir are of the order of a few MPa, with a one order of magnitude uncertainty based on the single force model. I also apply the implosive source model to the observed seismic magnitudes and estimate the product of bulk modulus of magma and flow velocity to be 10¹⁰–10¹¹ J/m² s.     

A fragmentation criterion for highly viscous bubbly magmas estimated from shock tube experiments

When a highly viscous bubbly magma is sufficiently decompressed, layer-by-layer fracturing propagates through the magma at a certain speed (fragmentation speed). On the basis of a recent shock tube theory by Koyaguchi and Mitani [Koyaguchi, T., Mitani, N. K., 2005. A theoretical model for fragmentation of viscous bubbly magmas in shock tubes. Journal of Geophysical Research 110 (B10), B10202. doi:10.1029/2004JB003513.], gas overpressures at the fragmentation surface are estimated from experimental data on fragmentation speed in shock tube experiments for natural volcanic rocks with various porosities. The results show that gas overpressure at the fragmentation surface increases as initial sample pressure increases and sample porosity decreases. We propose a new fragmentation criterion to explain the relationship between the gas overpressure at the fragmentation surface, the initial pressure and the porosity. Our criterion is based on the idea that total fragmentation of highly viscous bubbly magmas occurs when the tensile stress at the midpoint between bubbles exceeds a critical value. We obtain satisfactory agreement between our simulation and experiment when we assume that the critical value is inversely proportional to the square root of bubble wall thickness. This fragmentation criterion suggests that long micro-cracks or equivalent flaws (e.g., irregular-shaped bubbles) that reach the midpoints between bubbles are a dominant factor to determine the bulk strength of the bubbly magma.     

On magma fragmentation by conduit shear stress: Evidence from the Kos Plateau Tuff,Aegean Volcanic Arc

Large silicic explosive eruptions are the most catastrophic volcanic events. Yet, the intratelluric mechanisms underlying are not fully understood. Here we report a field and laboratory study of the Kos Plateau Tuff (KPT, 161 ka, Aegean Volcanic Arc), which provides an excellent geological example of conduit processes that control magma vesiculation and fragmentation during intermediate- to large-scale caldera-forming eruptions. A prominent feature of the KPT is the occurrence of quite unusual platy-shaped tube pumice clasts in pyroclastic fall and current deposits from the early eruption phases preceding caldera collapse. On macroscopic and SEM observations, flat clast faces are elongated parallel to tube vesicles, while transverse surfaces often occur at ~ 45° to vesicle elongation. This peculiar pumice texture provides evidence of high shear stresses related to strong velocity gradients normal to conduit walls, which induced vesiculation and fragmentation of the ascending magma. Either an increasing mass discharge rate without adequate enlargement of a narrow central feeder conduit or a developing fissure-like feeder system related to incipient caldera collapse provided suitable conditions for the generation of plate tube pumice within magma volumes under high shear during the pre-climactic KPT eruption phases. This mechanism implies that the closer to the conduit walls (where the stronger are the velocity gradients) the larger was the proportion of plate vs. conventional (lensoid) juvenile fragments in the ascending gas–pyroclast mixture. Consequently, plate pumice clasts were mainly entrained in the outer portions of the jet and convecting regions of a sustained, Plinian-type, eruption column, as well as in occasional lateral blast currents generated at the vent. As a whole, plate pumice clasts in the peripheral portions of the column were transported at lower altitudes and deposited by fallout or partial collapse closer to the vent relative to lensoid ones that dominated in the inner column portions. The plate tube pumice proportion decreased abruptly up to disappearance during the emplacement of the main pyroclastic currents and lithic-rich breccias related to extensive caldera collapse at the eruption climax, as a consequence of an overall widening of the magma feeder system through the opening of multiple conduits and eruptive vents, along with fissure erosion, concomitant to the disruption of the collapsing block.     

Magma degassing and basaltic eruption styles: a case study of 2000 year BP Xitle volcano in central Mexico

To investigate the relationship between volatile abundances and eruption style, we have analyzed major element and volatile (H₂O, CO₂, S) concentrations in olivine-hosted melt inclusions in tephra from the 2000 yr BP eruption of Xitle volcano in the central Trans-Mexican Volcanic Belt. The Xitle eruption was dominantly effusive, with fluid lava flows accounting for 95% of the total dense rock erupted material (1.1 km³). However, in addition to the initial, Strombolian, cinder cone-building phase, there was a later explosive phase that interrupted effusive activity and deposited three widespread ash fall layers. Major element compositions of olivine-hosted melt inclusions from these ash layers range from 52 to 58 wt.% SiO₂, and olivine host compositions are Fo_84–86. Water concentrations in the melt inclusions are variable (0.2–1.3 wt.% H₂O), with an average of 0.45±0.3 (1σ) wt.% H₂O. Sulfur concentrations vary from below detection (50 ppm) to 1000 ppm but are mostly ≤200 ppm and show little correlation with H₂O. Only the two inclusions with the highest H₂O have detectable CO₂ (310–340 ppm), indicating inclusion entrapment at higher pressures (700–900 bars) than for the other inclusions (≤80 bars). The low and variable H₂O and S contents of melt inclusions combined with the absence of less soluble CO₂ indicates shallow-level degassing before olivine crystallization and melt inclusion formation. Olivine morphologies are consistent with the interpretation that most crystallization occurred rapidly during near-surface H₂O loss. During cinder cone eruptions, the switch from initial explosive activity to effusive eruption probably occurs when the ascent velocity of magma becomes slow enough to allow near-complete degassing of magma at shallow depths within the cone as a result of buoyantly rising gas bubbles. This allows degassed lavas to flow laterally and exit near the base of the cone while gas escapes through bubbly magma in the uppermost part of the conduit just below the crater. The major element compositions of melt inclusions at Xitle show that the short-lived phase of renewed explosive activity was triggered by a magma recharge event, which could have increased overpressure in the storage reservoir beneath Xitle, leading to increased ascent velocities and decreased time available for degassing during ascent.     

Nonequilibrium flow in volcanic conduits and application to the eruptions of Mt. St. Helens on May 18, 1980, and Vesuvius in AD 79

A steady-state, one-dimensional, and nonhomogeneous two-phase flow model was developed for the prediction of local flow properties in volcanic conduits. The model incorporates the effects of relative velocity between the phases and for the variable magma viscosity. The resulting set of nonlinear differential equations was solved by a stiff numerical solver and the results were verified with the results of basaltic fissure eruptions obtained by a homogeneous two-phase flow model, before applying the model to the eruptions of Mt. St. Helens and Vesuvius volcanoes. This verification, and a study of the sensitivity of several modeling parameters, proved effective in establishing the confidence in the predicted nonequilibrium results of flow distribution in the conduits when the mass flow rate is critical or maximum. The application of the model to the plinian eruptions of Mt. St. Helens on May 18, 1980, and Vesuvius in AD 79, demonstrates the sensitivity of the magma discharge rate and distributions of pressure, volumetric fraction, and velocities of phases, on the hydrous magma viscosity feeding the volcanic conduits. Larger magma viscosities produce smaller mass discharge rates (or greater conduit diameters), smaller exit pressures, larger disequilibrium between the phases, and larger difference between the local lithostatic and fluid pressures in the conduit. This large pressure difference occurs when magma fragments and may cause a rupture of the conduit wall rocks, producing a closure of the conduit and cessation of the volcanic eruption, or water pouring into the conduit from underground aquifers leading to phreatomagmatic explosions. The motion of the magma fragmentation zone along a conduit during an eruption can be caused by the varying viscosity of magma feeding the volcanic conduit and may cause intermittent phreatomagmatic explosions during the plinian phases as different underground aquifers are activated at different depths. The variation of magma viscosity during the eruptions of Mt. St. Helens in 1980 and Vesuvius in AD 79 is normally associated with the tapping of magmas from different depths of the magma chambers. This variation of viscosity, which can include different crystal and dissolved water contents, can also produce conduit wall erosion, the onset and collapse of volcanic columns above the vent, and the onset and cessation of pyroclastic flows and surges.     

Experimental constraints on degassing and permeability in volcanic conduit flow

This study assesses the effect of decompression rate on two processes that directly influence the behavior of volcanic eruptions: degassing and permeability in magmas. We studied the degassing of magma with experiments on hydrated natural rhyolitic glass at high pressure and temperature. From the data collected, we defined and characterized one degassing regime in equilibrium and two regimes in disequilibrium. Equilibrium bubble growth occurs when the decompression rate is slower than 0.1 MPa s^–1, while higher rates cause porosity to deviate rapidly from equilibrium, defining the first disequilibrium regime of degassing. If the deviation is large enough, a critical threshold of super-saturation is reached and bubble growth accelerates, defining the second disequilibrium regime. We studied permeability and bubble coalescence in magma with experiments using the same rhyolitic melt in open degassing conditions. Under these open conditions, we observed that bubbles start to coalesce at ~43 vol% porosity, regardless of decompression rate. Coalescence profoundly affects bubble texture and size distributions, and induces the melt to become permeable. We determined coalescence to occur on a time scale (~180 s) independent of decompression rate. We parameterized and incorporated our experimental results into a 1D conduit flow model to explore the implications of our findings on eruptive behavior of rhyolitic melts with low crystal contents stored in the upper crust. Compared to previous models that assume equilibrium degassing of the melt during ascent, the introduction of disequilibrium degassing reduces the deviation from lithostatic pressure by ~25%, the acceleration at high porosities (>50 vol%) by a factor 5, and the associated decompression rate by an order of magnitude. The integration of the time scale of coalescence to the model shows that the transition between explosive and effusive eruptive regimes is sensitive to small variations of the initial magma ascent speed, and that flow conditions near fragmentation may significantly be affected by bubble coalescence and gas escape.Editorial responsibility: D. Dingwell     

Fine ash content of explosive eruptions

In explosive eruptions, the mass proportion of ash that is aerodynamically fine enough to cause problems with jet aircraft or human lungs (< 30 to 60 μm in diameter) is in the range of a few percent to more than 50%. The proportions are higher for silicic explosive eruptions, probably because vesicle size in the pre-eruptive magma is smaller than those in mafic magmas. There is good evidence that pyroclastic flows produce high proportions of fine ash by communition and it is likely that this process also occurs inside volcanic conduits and would be most efficient when the magma fragmentation surface is well below the summit crater. Reconstructed total grain size distributions for several recent explosive eruptions indicate that basaltic eruptions have small proportions of very fine ash (~ 1 to 4%) while tephra generated during silicic eruptions contains large proportions (30 to > 50%).     

Mixed deposits of simultaneous strombolian and phreatomagmatic volcanism: Rothenberg volcano, east Eifel volcanic field

The Pleistocene basanite-tephrite Rothenberg cone complex in the East Eifel was constructed by alternating dominantly Strombolian (S1–3) and dominantly phreatomagmatic (P1–3) phases of volcanism along a NNE-SSW linear vent system. Strombolian eruptions, from the central vent of the S1 scoria cone, and phreatomagmatic eruptions, from a vent on the southern margin of the cone, occurred simultaneously during the second phreatomagmatic phase (P2). The P2 deposits are a complex sequence in which Strombolian fallout ejecta is intimately admixed with phreatomagmatic fallout and pyroclastic surge material. Every bed contains at least trace amounts of ejecta from both sources but, at every site, an alternation of Strombolian-dominant and phreatomagmatic-dominant units is recorded. Each bed also shows marked lateral changes with a progressive northward increase in the proportion of Strombolian material. The two eruptive styles produced morphologically distinct clast populations often with widely separated (5–7 φ) grain size modes. The phreatomagmatic component of the P2 deposits is inferred to be the result of shallow interaction of external water and cool, partially degassed magma which reached the surface at a time when the magma column was retreating from the northern Strombolian central vent.The Rothenberg deposits illustrate the complexity and sensitivity of controls on Strombolian and associated phreatomagmatic volcanism, and the shallow depth of fragmentation during such eruptions. During such shallow eruptions minor, ephemeral and localised variations in the rate of rise and discharge of magma, and vent geometry and hydrology significantly influence the magma:water ratio and hence eruptive style.     

Strombolian and phreatomagmatic deposits of Ohakune craters, Ruapehu, New Zealand: A complex interaction between external water and rising basaltic magma

The Ohakune Craters form one of several parasitic centres surrounding Ruapehu volcano, at the southern end of the Taupo Volcanic Zone. An inner scoria cone and an outer, probably older, tuff ring are the principal structures in a nested cluster of four vents.The scoria cone consists of alternating lava flows and coarse, welded and unwelded, strombolian block and bomb beds. The strombolian beds consist of principally two discrete types of essential clast, vesicular bombs and dense angular blocks. Rare finer-grained beds are unusually block-rich. The tuff ring consists of alternating strombolian and phreatomagmatic units. Strombolian beds have similar grain size characteristics to scoria cone units, but contain more highly vesicular unoxidised bombs and few blocks. Phreatomagmatic deposits, which contain clasts with variable degrees of palagonitisation, consist of less well-sorted airfall deposits and very poorly sorted, crystal-rich pyroclastic surge deposits.Disruption by expanding magmatic gas bubbles was a major but relatively constant influence on both strombolian and phreatomagmatic eruptions at Ohakune. Instead, the nature of deposits was principally controlled by two other variables, vent geometry and the relative influence of external water during volcanism. During tuff-ring construction, magma is considered to have risen rapidly to the surface, and to have been ejected without sufficient residence time in the vent for non-explosive degassing. Availability of external water principally governed the eruption mechanism and hence the nature of the deposits. Essentials clasts of the scoria cone are, by comparison, dense, degassed and oxidised. It is suggested that a change in vent geometry, possibly the construction of the tuff ring itself, permitted lava ponding and degassing during scoria cone growth. During strombolian eruptions, magma remaining in the vent probably became depleted in gas, leading to the formation of an inert zone, or crust, above actively degassing magma. Subsequent explosions had therefore to disrupt both this passive crust and underlying, vesiculating magma “driving” the eruption. Cycles of strombolian eruption are thought to have stopped when the thickness of the inert crust precluded explosive eruption and only recommenced when some of this material was removed, either as a lava flow or during phreatomagmatic explosions when external water entered the vent. Such explosions probably formed the unusually fine-grained and block-rich beds in the strombolian sequence.The Ohakune deposits are an excellent example of the products of explosive eruption of fluid, gas-rich basic magma vesiculating under very near-surface conditions. A complex interplay of rate of magma rise, rate and depth of formation of gas bubbles, vent geometry, abundance of shallow external water, wind velocity and accumulation rate of ejecta determines the nature of deposits of such eruptions.     

Degassing during magma ascent in the Mule Creek vent (USA)

 The structures and textures of the rhyolite in the Mule Creek vent (New Mexico, USA) indicate mechanisms by which volatiles escape from silicic magma during eruption. The vent outcrop is a 300-m-high canyon wall comprising a section through the top of a feeder conduit, vent and the base of an extrusive lava dome. Field relations show that eruption began with an explosive phase and ended with lava extrusion. Analyses of glass inclusions in quartz phenocrysts from the lava indicate that the magma had a pre-eruptive dissolved water content of 2.5–3.0 wt% and, during eruption, the magma would have been water-saturated over the vertical extent of the present outcrop. However, the vesicularity of the rhyolite is substantially lower than that predicted from closed-system models of vesiculation under equilibrium conditions. At a given elevation in the vent, the volume fraction of primary vesicles in the rhyolite increases from zero close to the vent margin to values of 20–40 vol.% in the central part. In the centre the vesicularity increases upward from approximately 20 vol.% at 300 m below the canyon rim to approximately 40 vol.% at 200 m, above which it shows little increase. To account for the discrepancy between observed vesicularity and measured water content, we conclude that gas escaped during ascent, probably beginning at depths greater than exposed, by flow through the vesicular magma. Gas escape was most efficient near the vent margin, and we postulate that this is due both to the slow ascent of magma there, giving the most time for gas to escape, and to shear, favouring bubble coalescence. Such shear-related permeability in erupting magma is supported by the preserved distribution of textures and vesicularity in the rhyolite: Vesicles are flattened and overlapping near the dense margins and become progressively more isolated and less deformed toward the porous centre. Local zones have textures which suggest the coalescence of bubbles to form permeable, collapsing foams, implying the former existence of channels for gas migration. Local channelling of gas into the country rocks is suggested by the presence of sub-horizontal syn-eruptive rhyolitic tuffisite veins which depart from the vent margin and invade the adjacent country rock. In the central part of the vent, similar local channelling of gas is indicated by steep syn-eruption tuffisite veins which cut the rhyolite itself. We conclude that the suppression of explosive eruption resulted from gas separation from the ascending magma and vent structure by shear-related porous flow and channelling of gas through tuffisite veins. These mechanisms of gas loss may be responsible for the commonly observed transition from explosive to effusive behaviour during the eruption of silicic magma. Received: 24 May 1995 / Accepted: 13 March 1996     

Eruptive dynamics and caldera collapse during the Onano eruption, Vulsini, Italy

The Onano explosive eruption of the Latera Volcanic Complex (Vulsini Volcanoes, Quaternary potassic Roman Comagmatic Region, Italy) provides an interesting example of multiple changes of eruptive style that were concomitant with a late phase of collapse of the polygenetic Latera Caldera. This paper reports a reconstruction of the event based on field analysis, laboratory studies of grain size and density of juvenile clasts, and re-interpretation of available subsurface geology data. The Onano eruption took place in a structurally weak area, corresponding to a carbonate substrate high bordered by the pre-existing Latera caldera and Bolsena volcano-tectonic depression, which controlled the ascent and eruption of a shoshonitic-phonotephritic magma through intersecting rim fault systems. Temporal changes of magma vesiculation, fragmentation and discharge rate, and consequent eruptive dynamics, were strongly controlled by pressure evolution in the magma chamber and changing vent geometry. Initially, pumice-rich pyroclastic flows were emplaced, followed by spatter- and lithic-rich flows and fallout from energetic fire-fountaining. The decline of magma pressure due to the partial evacuation of the magma chamber induced trapdoor collapse of the chamber roof, which involved part of the pre-existing caldera and external volcano slopes and eventually led to the present-day caldera. The widening of the vent system and the emplacement of the main pyroclastic flow and associated co-ignimbrite lag breccia marked the eruption climax. A sudden drop of the confining pressure, which is attributed to a pseudo-rigid behaviour of the magma chamber wall rocks during a phase of rapid magma drainage, led to extensive magma vesiculation and fragmentation. The disruption of the magma chamber roof and waning magma pressure in the late eruption stage favoured the explosive interaction of residual magma with groundwater from the confined carbonate aquifer. Pulsating hydrostatic and magma pressures produced alternating hydromagmatic pyroclastic surges, strombolian fallout and spatter flows.     

Small inflation sources producing seismic and infrasonic pulses during the 2000 eruptions of Miyake-jima, Japan

During the 2000 activity of Miyake-jima volcano, Japan, we detected long period seismic signals with initial pulse widths of 1-2 s, accompanied by infrasonic pulses with almost the same pulse widths. The seismic signals were observed from 13 July 2000, a day before the second summit eruption. The occurrences of the seismic signals were intermittent with a gradual increase in their magnitudes and numbers building toward a significant explosive eruption on 18 August. After the eruption, the seismic and infrasonic events ceased. The results of a waveform inversion show that the initial motions were excited by an isotropic inflation source beneath the south edge of the caldera at a depth of 1.4 km. On the other hand, the sources of the infrasonic pulses were located in the summit caldera area. The times at which the infrasonic pulses were emitted at the surface were delayed by about 3 s from the origin times of the seismic events. It is suggested that small isotropic inflations excited seismic waves in the crust and simultaneously caused acoustic waves that traveled in the conduit and produced infrasonic pulses at the crater bottom. Considering the observed time differences and gas temperatures emitted from the vent, the conduit should have been filled with vapor mixed with SO₂ gas and volcanic ash. The change of the time differences between the seismic and infrasonic signals suggests that the seismic source became shallower within half a day before the August 18 explosive eruption. We interpret the source process as a fragmentation process of magma in which gas bubbles burst and quickly released part of the pressure that had been sustained by the tensional strength of magma.