Three distinct regimes of volcanic tremor associated with the eruption of Shishaldin Volcano, Alaska 1999

Tremor signals associated with the eruption of Shishaldin Volcano on 19 and 23 April 1999 were the strongest recorded anywhere in the Aleutian Arc by the Alaska Volcano Observatory (AVO) in its 10-year history. Reduced displacements (D_R) reached 23 cm² on 19 April and 43 cm² on 23 April. During the activity, D_R and spectral data with a frequency resolution of 0.1 Hz were computed and put on the World Wide Web every 10 min. These data are analyzed here. The general temporal patterns of seismicity of these eruption events were similar, but the eruptions and their effects quite different. The 19 April event is known to have culminated in a sub-Plinian phase, which ejected ash to an altitude of 16 km. Despite higher amplitudes and the largest hotspot from satellite data, the 23 April event produced little ash reaching only 6 km altitude. For several hours prior to the sub-Plinian phase on 19 April, tremor with a peak frequency of 1.3 Hz intensified. During the sub-Plinian phase the peak frequency increased to 4-8 Hz. However, in 15 h after the eruption, three episodes of stronger tremor occurred with a lower 1.0-Hz peak, alternating with weaker tremor with a 1.3-Hz peak. These transitions correspond to D_R=~8 cm². Although these strong tremor episodes produced higher D_R levels than the sub-Plinian phase, data from a pressure sensor show that only strong Strombolian explosions occurred. The suite of observations suggests three distinct tremor regimes that may correspond to slug flow, bubbly flow, and sustained strong eruptions, or a cyclic change in source parameters (e.g., geometry, sound speed, or ascent rate). This behavior occurred at Shishaldin only during the April 1999 sequence, and we are not aware of similar behavior at other volcanoes.     

Geological and seismological evidence of increased explosivity during the 1986 eruptions of Pavlof volcano,Alaska

We present results of study of the best-documented eruptions of Pavlof volcano in historic time. The 1986 eruptions were mostly Strombolian in character; a strong initial phase may have been Vulcanian. The 1986 activity erupted at least 8×10⁶ m³ of feldspar-phyric basaltic andesite lava (SiO₂=53–54%), and a comparable volume of wind-borne tephra. During the course of the eruption, 5300 explosion earthquakes occurred, the largest of which was equivalent to an M _L =2.5 earthquake. Volcanic tremor was recorded for 2600 hours, and the strongest tremor was recorded out to a distance of 160 km and had an amplitude of at least 54 cm² reduced displacement. The 1986 eruptions modified the structure of the vent area for the first time in over two decades. A possible pyroclastic flow was observed on 19 June 1986, the first time such a phenomenon has been observed at the volcano. Overall, the 1986 eruptions were the strongest and longest duration eruptions in historic time, and changed a temporal pattern of activity that had persisted from 1973–1984.     

Vent conditions for expected eruptions at Vesuvius

Determining consistent sets of vent conditions for next expected eruptions at Vesuvius is crucial for the simulation of the sub-aerial processes originating the volcanic hazard and the eruption impact. Here we refer to the expected eruptive scales and conditions defined in the frame of the EC Exploris project, and simulate the dynamics of magma ascent along the volcanic conduit for sub-steady phases of next eruptions characterized by intensities of the Violent Strombolian (VS), Sub-Plinian 2 (SP2), and Sub-Plinian 1 (SP1) scale. Sets of conditions for the simulations are determined on the basis of the bulk of knowledge on the past history of Vesuvius [Cioni, R., Bertagnini, A., Santacroce, R., Andronico, D., Explosive activity and eruption scenarios at Somma–Vesuvius (Italy): towards a new classification scheme. Journal of Volcanology and Geothermal Research, this issue.]. Volatile contents (H₂O and CO₂) are parameterized in order to account for the uncertainty in their expected amounts for a next eruption. In all cases the flow in the conduit is found to be choked, with velocities at the conduit exit or vent corresponding to the sonic velocity in the two-phase non-equilibrium magmatic mixture. Conduit diameters and vent mixture densities are found to display minimum overlapping between the different eruptive scales, while exit gas and particle velocities, as well as vent pressures, largely overlap. Vent diameters vary from as low as about 5 m for VS eruptions, to 35–55 m for the most violent SP1 eruption scale. Vent pressures can be as low as less than 1 MPa for the lowest volatile content employed of 2 wt.% H₂O and no CO₂, to 7–8 MPa for highest volatile contents of 5 wt.% H₂O and 2 wt.% CO₂ and large eruptive scales. Gas and particle velocities at the vent range from 100–250 m/s, with a tendency to decrease, and to increase the mechanical decoupling between the phases, with increasing eruptive scale. Except for velocities, all relevant vent quantities are more sensitive to the volatile content of the discharged magma for the highest eruptive scales considered.     

Damage mechanics preceding the September–October 1989 flank eruption at Mount Etna volcano inferred by seismic scaling exponents

Seismicity preceding the September–October 1989 flank eruption at Mt. Etna volcano has been analysed by means of the fractal two-point correlation dimension D_c relative to epicentres distribution and the scaling exponent D inferred from the seismic b-value. Time evolution of D_c and D preceding the eruption onset revealed: a general decrease of D_c and increase of D (9–4 months prior the eruption onset), intermittent rates of D_c increase and D decrease (4 to 1 month), steady high values of D_c and low values of D (1 month) and sharp decreases of D_c and D marking the eruption onset. Results allowed to constrain the increment of seismic damage by recognising negative (mechanical hardening) and positive (mechanical softening) feedback prior the eruption onset.     

Constraining tephra dispersion and deposition from three subplinian explosions in 2011 at Shinmoedake volcano,Kyushu, Japan

Constraining physical parameters of tephra dispersion and deposition from explosive volcanic eruptions is a significant challenge, because of both the complexity of the relationship between tephra distribution and distance from the vent and the difficulties associated with direct and comprehensive real-time observations. Three andesitic subplinian explosions in January 2011 at Shinmoedake volcano, Japan, are used as a case study to validate selected empirical and theoretical models using observations and field data. Tephra volumes are estimated using relationships between dispersal area and tephra thickness or mass/area. A new cubic B-spline interpolation method is also examined. Magma discharge rate is estimated using theoretical plume models incorporating the effect of wind. Results are consistent with observed plume heights (6.4–7.3 km above the vent) and eruption durations. Estimated tephra volumes were 15–34?×?10⁶ m³ for explosions on the afternoon of 26 January and morning of 27 January, and 5.0–7.6?×?10⁶ m³ for the afternoon of 27 January; magma discharge rates were in the range 1–2?×?10⁶ kg/s for all three explosions. Clast dispersal models estimated plume height at 7.1?±?1 km above the vent for each explosion. The three subplinian explosions occurred with approximately 12-h reposes and had similar mass discharge rates and plume heights but decreasing erupted magma volumes and durations.     

Characterization of tephra deposits with limited exposure: the example of the two largest explosive eruptions at Nisyros volcano (Greece)

Explosive eruptions associated with tephra deposits that are only exposed in proximal areas are difficult to characterize. In fact, the determination of physical parameters such as column height, mass eruption rate, erupted volume, and eruption duration is mainly based on empirical models and is therefore very sensitive to the quality of the field data collected. We have applied and compared different modeling approaches for the characterization of the two main tephra deposits, the Lower Pumice (LP) and Upper Pumice (UP) of Nisyros volcano, Greece, which are exposed only within 5 km of the probable vent. Isopach and isopleth maps were compiled for two possible vent locations (on the north and on the south rim of the caldera), and different models were applied to calculate the column height, the erupted volume, and the mass eruption rate. We found a column height of about 15 km above sea level and a mass eruption rate of about 2 × 10⁷ kg/s for both eruptions regardless of the vent location considered. In contrast, the associated wind velocity for both UP and LP varied between 0 and 20 m/s for the north and south vent, respectively. The derived erupted volume for the south vent (considered as the best vent location) ranges between 2 and 27 × 10⁸ m³ for the LP and between 1 and 5 × 10⁸ m³ for the UP based on the application of four different methods (integration of exponential fit based on one isopach line, integration of exponential and power-law fit based on two isopach lines, and an inversion technique combined with an advection–diffusion model). The eruption that produced the UP could be classified as subplinian. Discrepancies associated with different vent locations are smaller than the discrepancies associated with the use of different models for the determination of erupted mass, plume height, and mass eruption rate. Proximal outcrops are predominantly coarse grained with ≥90 wt% of the clasts ranging between −6ϕ and 0ϕ. The associated total grainsize distribution is considered to result from a combination of turbulent fallout from both the plume margins and the umbrella region, and as a result, it is fines-depleted. Given that primary deposit thickness observed on Nisyros for both LP and UP is between 1 and 8 m, if an event of similar scale were to happen again, it would have a significant impact on the entire island with major damage to infrastructure, agriculture, and tourism. Neighboring islands and the continent could also be significantly affected.     

Volcanic activity in the Azores. Report for 1959–1964

During the period 1959 to 1964 recorded volcanic activity in the Azores islands (Atlantic Ocean) was manifested by fumarolic fields in S. Miguel, Terceira, Graciosa, Pico and Faial and by two submarine eruptions off Pico and S. Jorge. Fumarolic activity increased in the Faial Capelinhos vent (in eruption 1957–58) after December 1963; it is concluded that the phenomenon is related to the 1963 and 1964 eruptions in the S. Jorge Channel. During February 1964 an important seismic swarm of more than 500 earthquakes, with volcanic tremor, occurred in S. Jorge.

     

Patterns in open vent, strombolian behavior at Fuego volcano, Guatemala, 2005–2007

Fuego volcano, Guatemala is a high (3,800 m) composite volcano that erupts gas-rich, high-Al basalt, often explosively. It spends many years in an essentially open vent condition, but this activity has not been extensively observed or recorded until now. The volcano towers above a region with several tens of thousands of people, so that patterns in its activity might have hazard mitigation applications. We conducted 2 years of continuous observations at Fuego (2005–2007) during which time the activity consisted of minor explosions, persistent degassing, paroxysmal eruptions, and lava flows. Radiant heat output from MODIS correlates well with observed changes in eruptive behavior, particularly during abrupt changes from passive lava effusion to paroxysmal eruptions. A short-period seismometer and two low-frequency microphones installed during the final 6 months of the study period recorded persistent volcanic tremor (1–3 Hz) and a variety of explosive eruptions. The remarkable correlation between seismic tremor, thermal output, and daily observational data defines a pattern of repeating eruptive behavior: 1) passive lava effusion and subordinate strombolian explosions, followed by 2) paroxysmal eruptions that produced sustained eruptive columns, long, rapidly emplaced lava flows, and block and ash flows, and finally 3) periods of discrete degassing explosions with no lava effusion. This study demonstrates the utility of low-cost observations and ground-based and satellite-based remote sensing for identifying changes in volcanic activity in remote regions of underdeveloped countries.     

Seismicity of Ruapehu volcano,New Zealand, 1971–1996: a review

From 1971 until 1995, the style of seismicity at Ruapehu changed little, reflecting a period of relatively low eruptive activity and consequent long-term stability within the vent system. Volcanic earthquakes and volcanic tremor were both dominated by a frequency of about 2 Hz. Volcanic earthquakes accompanied all phreatic and phreatomagmatic eruptions, but not small hydrothermal eruptions that originated within Crater Lake. Furthermore, more than half of the M_L>3 volcanic earthquakes and changes in the reduced displacement of 2 Hz volcanic tremor by as much as a factor of 20 occurred without any accompanying eruptive activity. Three and 7 Hz volcanic tremor were also recorded, although never at lower-elevation seismometers. At times, this tremor was stronger at the summit seismometer than the 2 Hz tremor. Their source regions were independent of the 2 Hz source, and located at shallower depths. Volcano-tectonic earthquakes were generally unrelated to eruptive activity. The seismicity accompanying the 1995–1996 eruptive activity was significantly different from that of the period 1971 to 1995, and included volcanic tremor with a frequency of less than 1 Hz, simultaneous changes in the amplitude of the previously independent 2 Hz and 7 Hz volcanic tremor, and finally a change in the frequency content of volcanic earthquakes and volcanic tremor from 2 Hz to wideband. Path transmission effects play an important role in determining the characteristics of seismograms at Ruapehu. The presence of Crater Lake affects both the style of eruptions and the accompanying seismicity.     

New insights into the 1999 eruption of Shishaldin volcano,Alaska, based on acoustic data

Data collected by a pressure sensor provide new insights into the 1999 eruption of Shishaldin volcano, Unimak Island, Alaska. On 19 April 1999, after 3 months of unrest and an extended period of low-level Strombolian activity, Shishaldin experienced a Subplinian eruption (ash plume to >16 km), followed by several episodes of strong Strombolian explosions. Acoustic data from the pressure sensor allow us to investigate the details of an eruption which was instrumentally well recorded, but with few visual observations. In the 12 h prior to the Subplinian phase, the pressure sensor detected a series of small, repeated pulses with a constant spectral peak at 2–3 Hz. The amplitude and occurrence rate of the pulses both grew such that the signal became a nearly continuous hum just before the Subplinian eruption. This humming signal may represent gas release from rising magma. The main Subplinian phase was heralded by (1) the abrupt end of the humming signal, (2) several pulses of low-frequency sound interpreted as ash bursts, and (3) a dramatic increase in seismic tremor amplitude. The change in acoustic signature at this time allows us to precisely time the start of the Subplinian eruption, previously approximated as the time of strongest tremor increase. The 50-min Subplinian phase actually contained several bursts of signal, each of which may represent a discrete volume of magma passing through the system. Following the Subplinian event, the pressure sensor recorded four discrete episodes of Strombolian gas explosions on 19–20 April and another on 22–23 April. Four of the five episodes were accompanied by strong seismic tremor; the fifth has not been previously recognized and was not associated with anomalous tremor amplitudes. In time series these events are similar to explosions recorded at other volcanoes but in general they are much larger, with maximum amplitudes of >65 Pa at 6.5 km from the vent, and they have low (0.7–1.5 Hz) peak frequencies. These large explosions occurred at rates of 3–20 per minute for 1–5 h in each episode. The explosions were accompanied by a small (<5 km above sea level) ash plume and only minor amounts of ejecta were produced. Thus, the explosion activity was dominated by gas release.     

Long-term forecasting of the extreme eruptions of Etna

Until now, the widely used predictive long-term hazard models of the most intense natural phenomena do not yield predictions other than the mean return period. To achieve less uncertain estimates, an alternative model is proposed. The analysis of the eruptions of Etna from 1868 to 1993 shows that the time interval D_t between consecutive large eruptions is clearly random, whereas the total volume of tephra ejected during D_t is roughly constant. Taking into account the above evidence, the model yields significant long-term forecasts both of the time and of the intensity of future large eruptions. The large volume of tephra ejected by the volcano since the last intense eruption in 1980 allows the model to estimate that the next large eruption will not to take place soon.     

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.     

Reconstruction of a kimberlite eruption,using an integrated volcanological,geochemical and numerical approach: A case study of the Fox Kimberlite,NWT, Canada

An integrated approach involving volcanology, geochemistry and numerical modelling has enabled the reconstruction of the volcanic history of the Fox kimberlite pipe. The observed deposits within the vent include a basal massive, poorly sorted, matrix supported, lithic fragment rich, eruption column collapse lapilli tuff. Extensive vent widening during the climactic magmatic phase of the eruption led to overloading of the eruption column with cold dense country rock lithic fragments, dense juvenile pyroclasts and olivine crystals, triggering column collapse. > 40% dilution of the kimberlite by granodiorite country rock lithic fragments is observed both in the physical componentry of the rocks and in the geochemical signature, where enrichment in Al₂O₃ and Na₂O compared to average values for coherent kimberlite is seen. The wide, deep, open vent provided a trap for a significant proportion of the collapsing column material, preventing large scale run-away in the form of pyroclastic flow onto the ground surface, although minor flows probably also occurred. A massive to diffusely bedded, poorly sorted, matrix supported, accretionary-lapilli bearing, lithic fragment rich, lapilli tuff overlies the column collapse deposit providing evidence for a late phreatomagmatic eruption stage, caused by the explosive interaction of external water with residual magma. Correlation of pipe morphology and internal stratigraphy indicate that widening of the pipe occurred during this latter stage and a thick granodiorite cobble-boulder breccia was deposited. Ash- and accretionary lapilli-rich tephra, deposited on the crater rim during the late phreatomagmatic stage, was subsequently resedimented into the vent. Incompatible elements such as Nb are used as indicators of the proportion of the melt fraction, or kimberlite ash, retained or removed by eruptive processes. When compared to average coherent kimberlite the ash-rich deposits exhibit ~ 30% loss of fines whereas the column collapse deposit exhibits ~ 50% loss. This shows that despite the poorly sorted nature of the column collapse deposit significant elutriation has occurred during the eruption, indicating the existence of a high sustained eruption column. The deposits within Fox record a complex eruption sequence showing a transition from a probable violent sub-plinian style eruption, driven by instantaneous exsolution of magmatic volatiles, to a late phreatomagmatic eruption phase. Mass eruption rate and duration of the sub-plinian phase of the eruption have been determined based on the dimensions of milled country-rock boulders found within the intra-vent deposits. Calculations show a short lived eruption of one to eleven days for the sub-plinian magmatic phase, which is similar in duration to small volume basaltic eruptions. This is in general agreement with durations of kimberlite eruptions calculated using entirely different approaches and parameters, such as predictions of magma ascent rates in kimberlite dykes.     

Reventador Volcano 2005: Eruptive activity inferred from seismo-acoustic observation

Reventador Volcano entered an eruptive phase in 2005 which included a wide variety of seismic and infrasonic activity. These are described and illustrated: volcano-tectonic, harmonic tremor, drumbeats, chugging and spasmodic tremor, long period and very long period events. The recording of this simultaneous activity on an array of three broadband, seismo-acoustic instruments provides detailed information of the state of the conduit and vent during this phase of volcanic eruption. Quasi-periodic tremor at Reventador is similar to that observed at other volcanoes and may be used as an indicator of vent aperture. Variations in the vibration modes of the volcano, frequency fluctuations and rapid temporal fluctuations suggest the influx of new material, choking of the vent and possible modification of the conduit geometry during explosions and effusion over a period of six weeks.     

A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions

During volcanic eruptions, volcanic ash transport and dispersion models (VATDs) are used to forecast the location and movement of ash clouds over hours to days in order to define hazards to aircraft and to communities downwind. Those models use input parameters, called “eruption source parameters”, such as plume height H, mass eruption rate Ṁ, duration D, and the mass fraction m₆₃ of erupted debris finer than about 4 or 63 μm, which can remain in the cloud for many hours or days. Observational constraints on the value of such parameters are frequently unavailable in the first minutes or hours after an eruption is detected. Moreover, observed plume height may change during an eruption, requiring rapid assignment of new parameters. This paper reports on a group effort to improve the accuracy of source parameters used by VATDs in the early hours of an eruption. We do so by first compiling a list of eruptions for which these parameters are well constrained, and then using these data to review and update previously studied parameter relationships. We find that the existing scatter in plots of H versus Ṁ yields an uncertainty within the 50% confidence interval of plus or minus a factor of four in eruption rate for a given plume height. This scatter is not clearly attributable to biases in measurement techniques or to well-recognized processes such as elutriation from pyroclastic flows. Sparse data on total grain-size distribution suggest that the mass fraction of fine debris m₆₃ could vary by nearly two orders of magnitude between small basaltic eruptions ( 0.01) and large silicic ones (> 0.5). We classify eleven eruption types; four types each for different sizes of silicic and mafic eruptions; submarine eruptions; “brief” or Vulcanian eruptions; and eruptions that generate co-ignimbrite or co-pyroclastic flow plumes. For each eruption type we assign source parameters. We then assign a characteristic eruption type to each of the world's  1500 Holocene volcanoes. These eruption types and associated parameters can be used for ash-cloud modeling in the event of an eruption, when no observational constraints on these parameters are available.     

The possible effects of large 19th and 20th century volcanic eruptions on zonal and hemispheric surface temperatures

The contribution of volcanic material to the stratosphere from major eruptions within the last two centuries has been estimated using volcanological criteria, including eruption type, eruption column height, volume and duration of eruption, and composition and degree of fragmentation of magma. The chronology of major explosive volcanic eruptions is compared with a record of mean surface-temperature deviation (ΔT) for the same interval constructed from all available temperature data. The temperature records are divided into 6 latitudinal zones, allowing analysis for individual zones where temperature changes induced by aerosol perturbation might be intensified.We focus on the explosive volcanic events which by our estimates injected the most material into the stratosphere. These include Tambora 1815, Krakatau 1883, Santa Maria 1902, Katmai 1912 and Quizapu 1932. Such eruptions appear to have produced a consistent but small temperature decrease on the order of 0.2° to 0.5°C on a hemispheric scale for periods ranging from one to five years, although these changes are similar to background temperature variations. The maximum change in ΔT after some of these explosions appears to lag by up to three years in going from equatorial to polar latitudes.Somewhat smaller eruptions, e.g. Agung 1963 and possibly Cosiguina 1835, seem to have produced about the same perturbation in ΔT as the larger eruptions. This suggests either a limiting mechanism on loading of the aerosol layer after a volcanic eruption or, that the composition of injected material (i.e., the ratio of silicate “dust” to volatiles, and composition of the volatiles) may significantly effect stratospheric optical depth perturbation. Temperatures do not remain depressed for a longer period after a series of closely timed eruptions (e.g., the 1881–1889 or the 1902–1903 sequences) than after single events.     

Magma supply, magma discharge and readjustment of the feeding system of mount St. Helens during 1980

The landslide and cataclysmic eruption of Mount St. Helens on May 18, 1980 triggered a sequence of explosive eruptions over the following five months. The volume of explosive products from each of these eruptions decreased uniformly over this period, and the character for each eruption progressed from a fairly continuous eruption lasting more than eight hours on May 18 to a series of short bursts, some of which were spaced 12 hours apart, on October 16–18. The transition in the character of these eruption sequences can be explained by a difference between the magma supply rate and the magma discharge rate from a shallow reservoir.The magma supply rate (MSR) is the rate with which magma is supplied to the level where disruption due to vesiculation occurs. It is determined by dividing the dense-rock-equivalent volume of eruptive products by the total duration of each eruption sequence. The magma discharge rate (MDR) is the rate with which the disrupted magma is discharged through the vent. It is determined by dividing the volume of erupted products by the duration of each explosive burst. The relative magnitude of these two quantities controls the temporal evolution of an explosive event. When MDRMSR the explosive phase of the eruption lasts for several hours as a single continuous event. When MDR>MSR, an eruption is characterized by a series of short explosive bursts at intervals of several minutes to several days. The MSR of the eruptions of 1980 decreased with time from 5500 m' s⁻¹ on May 18 to 7 m³ s⁻² on October 16–18 and approximately fits an exponential decay. The MDR for the same events remained approximately constant at 2000 m³ s⁻¹. Each explosive event has been followed by an aftershock-like series of earthquakes located beneath the volcano at depth mostly between 7 and 14 km. The seismic energy released during each of these series is proportional to the corresponding volume of erupted magma. Deformation data between June and November, 1980 indicate a subsidence of the volcanic structure which can be modeled by a volume collapse of 0.25 km³ located at 9 km depth.We propose a model in which magma is supplied from depths of 7–14 km through a narrow conduit during each eruption. It erupts to the surface at a uniform rate during each eruption. The deep seismic activity following each eruption is related to a readjustment and volume decrease in the deep feeding system. The decrease of the MSR over time is explained by an increase in the viscosity of a progressively water-depleted magma. The amount of water necessary to explain the observed decrease of the MSR is of the order of 4.6%.     

New physical characterization of the Fontana Lapilli basaltic Plinian eruption,Nicaragua

The Fontana Lapilli deposit was erupted in the late Pleistocene from a vent, or multiple vents, located near Masaya volcano (Nicaragua) and is the product of one of the largest basaltic Plinian eruptions studied so far. This eruption evolved from an initial sequence of fluctuating fountain-like events and moderately explosive pulses to a sustained Plinian episode depositing fall beds of highly vesicular basaltic-andesite scoria (SiO₂ > 53 wt%). Samples show unimodal grain size distribution and a moderate sorting that are uniform in time. The juvenile component predominates (> 96 wt%) and consists of vesicular clasts with both sub-angular and fluidal, elongated shapes. We obtain a maximum plume height of 32 km and an associated mass eruption rate of 1.4 × 10⁸ kg s⁻¹ for the Plinian phase. Estimates of erupted volume are strongly sensitive to the technique used for the calculation and to the distribution of field data. Our best estimate for the erupted volume of the majority of the climactic Plinian phase is between 2.9 and 3.8 km³ and was obtained by applying a power-law fitting technique with different integration limits. The estimated eruption duration varies between 4 and 6 h. Marine-core data confirm that the tephra thinning is better fitted by a power-law than by an exponential trend.     

Multifractal characterization of Vesuvio lava-flow margins and its implications

The digitized lava-flow margins of well-defined extended eruptions occurring at Vesuvio in 1760, 1794, 1861, 1906, 1929 and 1944 are found to follow fractal behaviours inside a scaling region enclosed between 50 and 400 m. Although the invariance region is well respected, the fractal dimension D varies from one lava flow to another: the more irregular the lava-flow margin, the larger the value of D. The ascertained dependence of D on the duration of premonitory activity, preceding the emission of lavas, might provide some insight into the inner volcanic processes before the eruption and into the dynamical processes operating during flow emplacement.     

The rio caliente ignimbrite: Analysis of a compound intraplinian ignimbrite from a major late quaternary Mexican eruption

The Rio Caliente ignimbrite is a multi-flow unit orcompound ignimbrite formed during a major late Quaternary explosive rhyolitic eruption of La Primavera volcano, Mexico. The eruption sequence of the ignimbrite is complex and it occurs between lower and upper plinian air-fall deposits. It is, therefore, anintraplinian ignimbrite. Air-fall layers, pyroclastic surge, mudflow and fluviatile reworked pumice deposits also occur interbedded between ignimbrite flow units. A chaotic near-vent facies of the ignimbrite includes co-ignimbrite lag breccias segregated from proximal pumice flows. The facies locates a central vent but one which could not have been associated with a well defined edifice. Many of the lithics in the exposed lag breccias and near-vent facies of the ignimbrite appear to be fragments of welded Rio Caliente ignimbrite, and indicate considerable vent widening, or migration, during the eruption. Nearer vent the ignimbrite is thickest and composed of the largest number of flow units. Here it is welded and is a simple cooling unit. Evidence suggests that it was only the larger thicker pumice flows that escaped to the outer parts of the sheet. Detailed analysis of four flow units indicates that the pumice flows were generally poorly expanded, less mobile flows which would be produced by collapse of low eruption columns. The analogy of a compound ignimbrite with a compound lava flow is, therefore, good — a compound lava flow forms instead of a simple one when the volumetric discharge rate (or intensity) is low, and in explosive eruptions this predicts lower eruption column heights. A corollary is that the ignimbrite has a high aspect ratio. The complex eruption sequence shows the reinstatement of plinian activity several times during the eruption after column collapse occurred. This, together with erosional breaks and evidence that solidified fragments of already welded ignimbrite were re-ejected, all suggest the eruption lasted a relatively significant time period. Nearly 90 km³ of tephra were erupted. The associated plinian pumice fall is one of the largest known having a volume of 50 km³ and the ignimbrite, plus a co-ignimbrite ash-fall, have a volume of nearly 40 km³. Published welding models applied to the reejected welded blocks indicate an eruption duration of 15-20d, and a maximum average magma-discharge rate of 1.4 × 10⁴ m³/s for the ignimbrite. This is low intensity when compared with available data from other ignimbrite-forming eruptions, and concurs with all the geological evidence presented. The total eruption duration was perhaps 15-31d, which is consistent with other estimates of the duration of large magnitude explosive silicic eruptions.