Explore the vast range of titles available.

The dynamics of the May 18, 1980 lateral blast at Mount St. Helens, Washington (USA), were studied by means of a three‐dimensional multiphase flow model. Numerical simulations describe the blast flow as a high‐velocity pyroclastic density current generated by a rapid expansion (burst phase, lasting less than 20 s) of a pressurized polydisperse mixture of gas and particles and its subsequent gravitational collapse and propagation over a rugged topography. Model results show good agreement with the observed large‐scale behavior of the blast and, in particular, reproduce reasonably well the front advancement velocity and the extent of the inundated area. Detailed analysis of modeled transient and local flow properties supports the view of a blast flow led by a high‐speed front (with velocities between 100 and 170 m/s), with a turbulent head relatively depleted in fine particles, and a trailing, sedimenting body. In valleys and topographic lows, pyroclasts accumulate progressively at the base of the current body after the passage of the head, forming a dense basal flow depleted in fines (less than 5 wt.%) with total particle volume fraction exceeding 10−1 in most of the sampled locations. Blocking and diversion of this basal flow by topographic ridges provides the mechanism for progressive current unloading. On ridges, sedimentation occurs in the flow body just behind the current head, but the sedimenting, basal flow is progressively more dilute and enriched in fine particles (up to 40 wt.% in most of the sampled locations). In the regions of intense sedimentation, topographic blocking triggers the elutriation of fine particles through the rise of convective instabilities. Although the model formulation and the numerical vertical accuracy do not allow the direct simulation of the actual deposit compaction, present results provide a consistent, quantitative model able to interpret the observed stratigraphic sequence. Key Points The May 18, 1980 Mount St Helens blast has been simulated with a numerical model Numerical results correctly reproduce the flow front velocity and inundated area The model describes the non‐equilibrium sedimentation across topography

It is widely believed that andesitic magmas erupted at arc‐volcanoes are stored in shallow reservoirs prior to eruption, but high‐resolution images of focused regions of magma in the shallow crust are rare. We integrate seismic tomography with numerical models of magma chamber growth to constrain the magma chamber beneath Soufrière Hills Volcano, Montserrat. Our approach reveals the characteristics and dynamics of the magmatic system with a level of detail that no single method has yet achieved. The integrated analysis suggests that a magma chamber of 13 km3with over 30% melt fraction formed between 5.5 and at least 7.5 km depth, a significantly higher melt fraction than inferred from the seismic data alone. The magma chamber may have formed by incremental sill intrusion over a few thousand years and is likely to be a transient, geologically short‐lived feature. These volume and geometry estimates are critical parameters to model eruption dynamics, which in turn are key to hazard assessment and eruption forecasting. Key Points We image an active magma chamber with seismic tomography in unique detail Integration of tomography with thermal models gives new insight into magmatism Shallow magma chambers beneath arc volcanoes are transient and short lived

The CALIPSO collaborative volcano monitoring system on the Caribbean island of Montserrat includes observations of strain at depths ∼200 m using Sacks‐Evertson strainmeters. Strain data for the March 2004 explosion of the Soufrière Hills Volcano are characterized by large, roughly equal but opposite polarity changes at the two near sites and much smaller changes at a more distant site. The strain amplitudes eliminate a spherical pressure (Mogi‐type) source as the sole contributor. The initial changes are followed by smaller recoveries, but with differing relative recovery magnitudes. This dissimilarity requires a minimum of two pressure sources, which we model as a deep spherical pressure source and a shallow dike. The spherical source is fixed at the location derived from data for the massive dome collapse in July 2003. We solve for the best fitting dike plus sphere source combination. The dike geometry is consistent with earlier interpretations of dikes based on GPS data and other lines of evidence.

Several analytical and numerical eruption models have provided insight into volcanic eruption behaviour 1 , 2 , 3 , 4 , 5 , but most address plinian-type eruptions where vent conditions are quasi-steady. Only a few studies have explored the physics of short-duration vulcanian explosions 6 , 7 , 8 , 9 with unsteady vent conditions and blast events 10 , 11 . Here we present a technique that links unsteady vent flux of vulcanian explosions to the resulting dispersal of volcanic ejecta, using a numerical, axisymmetric model with multiple particle sizes. We use observational data from well documented explosions in 1997 at the Soufrière Hills volcano in Montserrat, West Indies, to constrain pre-eruptive subsurface initial conditions and to compare with our simulation results. The resulting simulations duplicate many features of the observed explosions, showing transitional behaviour where mass is divided between a buoyant plume and hazardous radial pyroclastic currents fed by a collapsing fountain 12 . We find that leakage of volcanic gas from the conduit through surrounding rocks over a short period (of the order of 10 hours) or retarded exsolution can dictate the style of explosion. Our simulations also reveal the internal plume dynamics and particle-size segregation mechanisms that may occur in such eruptions.

Assessments of pyroclastic flow (PF) hazards are commonly based on mapping of PF and surge deposits and estimations of inundation limits, and/or computer models of varying degrees of sophistication. In volcanic crises a PF hazard map may be sorely needed, but limited time, exposures, or safety aspects may preclude fieldwork, and insufficient time or baseline data may be available for reliable dynamic simulations. We have developed a statistically constrained simulation model for block-and-ash type PFs to estimate potential areas of inundation by adapting methodology from Iverson et al. (Geol Soc America Bull 110:972–984, 1998 ) for lahars. The predictive equations for block-and-ash PFs are calibrated with data from several volcanoes and given by A  = (0.05 to 0.1) V 2/3 , B  = (35 to 40) V 2/3 , where A is cross-sectional area of inundation, B is planimetric area and V is deposit volume. The proportionality coefficients were obtained from regression analyses and comparison of simulations to mapped deposits. The method embeds the predictive equations in a GIS program coupled with DEM topography, using the LAHARZ program of Schilling ( 1998 ). Although the method is objective and reproducible, any PF hazard zone so computed should be considered as an approximate guide only, due to uncertainties on the coefficients applicable to individual PFs, the authenticity of DEM details, and the volume of future collapses. The statistical uncertainty of the predictive equations, which imply a factor of two or more in predicting A or B for a specified V , is superposed on the uncertainty of forecasting V for the next PF to descend a particular valley. Multiple inundation zones, produced by simulations using a selected range of volumes, partly accommodate these uncertainties. The resulting maps show graphically that PF inundation potentials are highest nearest volcano sources and along valley thalwegs, and diminish with distance from source and lateral distance from thalweg. The model does not explicitly consider dynamic behavior, which can be important. Ash-cloud surge impact limits must be extended beyond PF hazard zones and we provide several approaches to do this. The method has been used to supply PF and surge hazard maps in two crises: Merapi 2006; and Montserrat 2006–2007.

The SEA‐CALIPSO experiment was designed to image structure related to active volcanism beneath the island of Montserrat in the Caribbean. As part of that experiment, over 200 Texan recorders with 5 Hz geophones were deployed in 3 linear arrays at a nominal spacing of 100m, primarily to record an airgun source towed offshore around the island. Because the recorders were operating in continuous mode for three days, a number of shallow microearthquakes under the active summit of Soufriere Hills Volcano (SHV) were also recorded. 20 events were sufficiently well recorded and located that they could be used to identify and map reflections from deep subsurface structure. Here we report on the processing of those recordings as multichannel CMP reflection sources, with emphasis on careful statics correction and coherency enhancement. The resulting reflection images indicate subhorizontal layering at depths between 6 and 19km which we interpret as sills.

Noritic anorthosite, gabbroic anorthosite and hornblende‐gabbro xenoliths are ubiquitous in the host andesite at Montserrat. Other xenoliths include quartz diorite, metamorphosed biotite‐gabbro, plagioclase‐hornblendite and plagioclase‐clinopyroxenite. Mineral compositions suggest a majority of the xenoliths are cognate. Cumulate, hypabyssal and crescumulate textures are present. A majority of the xenoliths are estimated to have seismic velocities of 6.7–7.0 km/s for pore‐free assemblages. These estimates are used in conjunction with petrological models to constrain the SEA CALIPSO seismic data and the structure of the crust beneath Montserrat. Andesitic upper crust is interpreted to overlie a lower crust dominated by amphibole and plagioclase. Xenolith textures and seismic data indicate the presence of hypabyssal intrusions in the shallow crust. The structure of the crust is consistent with petrological models indicating that fractionation is the dominant process producing andesite at Montserrat.

Dome growth at the Soufriere Hills volcano (1996 to 1998) was frequently accompanied by repetitive cycles of earthquakes, ground deformation, degassing, and explosive eruptions. The cycles reflected unsteady conduit flow of volatile-charged magma resulting from gas exsolution, rheological stiffening, and pressurization. The cycles, over hours to days, initiated when degassed stiff magma retarded flow in the upper conduit. Conduit pressure built with gas exsolution, causing shallow seismicity and edifice inflation. Magma and gas were then expelled and the edifice deflated. The repeat time-scale is controlled by magma ascent rates, degassing, and microlite crystallization kinetics. Cyclic behavior allows short-term forecasting of timing, and of eruption style related to explosivity potential.

Five Vulcanian explosions were triggered by collapse of the Soufrière Hills Volcano lava dome in 2003. We report strainmeter data for three explosions, characterized by four stages: a short transition between the onset of disturbance and a pronounced change in strain; a quasi‐linear ramp accounting for the majority of strain change; a more gradual continued decline of strain to a minimum value; and a strain recovery phase lasting hours. Remarkable ∼800 s barometric gravity waves propagated at ∼30 m s−1. Eruption volumes estimated from plume height and strain data are 0.32–0.42 × 106, 0.26–0.49 × 106, and 0.81–0.84 × 106 m3, for Explosions 3–5 respectively, consistent with quasi‐cylindrical conduit drawdown <2 km. The duration of vigorous explosion is given by the strain signature, indicating mass fluxes of order 107 kg s−1. Conduit pressures released reflect static weight of porous gas‐charged magma, and exsolution‐generated overpressures of order 10 MPa.

We explore volcanotectonic earthquake activity during the early stages of volcanic activity (1995–96) at Soufrière Hills volcano. We focus on several zones of temporally‐confined seismic activity that comprise assemblages of small scale structural elements, >2–4 km distant from the volcano, at depths 2–4 km below sea level, and consider seismicity in relation to regional tectonics and heterogeneity of stiffness and strength as revealed by the SEA‐CALIPSO tomography experiment. The clustered seismicity and relatively aseismic zones are interpreted to reflect a broad weakened tectonic zone of ESE trend that crosses Montserrat, and the ascent of a magmatic dike of NNE trend, which altered the stress distribution to promote localized fault movements and caused localized dilatation with changes in pore‐fluid pressures, to either weaken or strengthen the rock mass depending on the local polarity of strain.