Eruption Intensity Research Articles

Evolution of monogenetic rift-related alkaline magmatism in south Egypt: Insight from stratigraphy, geochronology, and geochemistry of the Natash volcanics

The Mesozoic Natash volcanics are rift-related monogenetic volcanic complex that erupted in south Eastern Desert of Egypt. The architecture of their building has evolved through three phases namely: (1) tuff cone consisting of explosive -dominated pyroclastics that stand for an intensity eruption peak; (2) the effusive lava-conquered phase, characterized by prevalent basalt-trachyte-rhyolite fissure eruptions; and (3) late-stage subvolcanic plugs/domes. The first phase began with the emplacement of spatter-rich ignimbrites followed by grain-supported lapilli-tuffs and hyaloclastics in a shallow-water environment (phreatomagmatic phase) with a vent-opening phase. This initial phase is followed by deposition of clast-supported conglomerates of lahar facies that represent significant erosional period, separating the pyroclastics and effusive phases. The second phase comprises multiphases displaying compositional stratification (primitive beneath evolved) of pahoehoe Hawaiian and rubbly vulcanian lavas, overlying the pyroclastic deposits. The third phase is the youngest and comprises plugs of aphyric trachytes. The volcanics occasionally have vesicle sheets, cylinder-elongate vesicles, gas blister, and inflation clefts, reflecting degassing, rupturing, brecciation/quenching history, and pressurized flows. U-Pb zircon ages from two samples of trachyte volcanics gave concordant age of 92 ± 0.9 Ma. One trachyte plug yielded concordant age of 86.6 ± 1.3 Ma. These ages represent two events belonging to the major Senonian volcanic activities. The investigated volcanics have alkaline affinity and display a continuous composition ranging from basalt to trachyte and rhyolite. Their chemical compositions display linear, curvilinear, and scattered trends on the Harker diagrams matched with the petrographic data and reflect comagmatic genesis for the entire sequence. Crystal fractionation of olivine + pyroxene + plagioclase + Fe-Ti oxides + apatite plays a role in the volcanic evolution. The origin of the mafic rocks is derived from deep mantle source in the garnet facies via low degree of partial melting. Such magma source is affected by fluid-rich metasomatism generated by dehydration subducted crust on the EM I, EMII, and HIMU mantle reservoir during Neoproterozoic time. The basalt-trachyte-rhyolite suites and their pyroclastics from the Natash area, together with Mesozoic alkaline volcanics from the surrounding terrains in Africa, constitute a remnant of extensive magmatisms, having large extension, high eruption rate, and OIB-type chemical composition.

Combining textural and geochemical investigations to explore the dynamics of magma ascent during Plinian eruptions: a Somma–Vesuvius volcano (Italy) case study

Trigger mechanisms and syn-eruptive processes of Plinian eruptions are poorly understood especially in the case of mafic powerful events. In the last decades, the combined geochemical and textural studies on volcanic rocks have proven to be fundamental tools for exploring the dynamics of magma ascent in volcanic conduits and for improving our ability to interpret volcano-monitoring signals and assess hazard. In this case study, we quantitatively investigate 2D and 3D micro-textural, geochemical, and isotopic features of pyroclastic rocks erupted during the Pomici di Base Plinian eruption (22 ka), the generally acknowledged first and most powerful event of the Somma–Vesuvius volcano. A peculiar aspect of this eruption is its high intensity that remained stable during the entire Plinian phase despite the strong magma compositional variation towards mafic terms. We infer that the transfer of magma towards the surface was intensified by the occurrence of rapid vesiculation pulses driven by limestone assimilation (skarn recycling) during magma ascent through the carbonatic bedrock. We conclude that limestone assimilation can hence be a syn-eruptive process, able to trigger further gas nucleation with deep impact on the eruption intensity, particularly crucial in the case of mafic/intermediate magma compositions.

Total grain-size distribution of four subplinian–Plinian tephras from Hekla volcano, Iceland: Implications for sedimentation dynamics and eruption source parameters

The size distribution of the population of particles injected into the atmosphere during a volcanic explosive eruption, i.e., the total grain-size distribution (TGSD), can provide important insights into fragmentation efficiency and is a fundamental source parameter for models of tephra dispersal and sedimentation. Recent volcanic crisis (e.g. Eyjafjallajökull 2010, Iceland and Córdon Caulle 2011, Chile) and the ensuing economic losses, highlighted the need for a better constraint of eruption source parameters to be used in real-time forecasting of ash dispersal (e.g., mass eruption rate, plume height, particle features), with a special focus on the scarcity of published TGSD in the scientific literature. Here we present TGSD data associated with Hekla volcano, which has been very active in the last few thousands of years and is located on critical aviation routes. In particular, we have reconstructed the TGSD of the initial subplinian–Plinian phases of four historical eruptions, covering a range of magma composition (andesite to rhyolite), eruption intensity (VEI 4 to 5), and erupted volume (0.2 to 1km3). All four eruptions have bimodal TGSDs with mass fraction of fine ash (<63μm; m63) from 0.11 to 0.25. The two Plinian dacitic-rhyolitic Hekla deposits have higher abundances of fine ash, and hence larger m63 values, than their andesitic subplinian equivalents, probably a function of more intense and efficient primary fragmentation. Due to differences in plume height, this contrast is not seen in samples from individual sites, especially in the near field, where lapilli have a wider spatial coverage in the Plinian deposits. The distribution of pyroclast sizes in Plinian versus subplinian falls reflects competing influences of more efficient fragmentation (e.g., producing larger amounts of fine ash) versus more efficient particle transport related to higher and more vigorous plumes, displacing relatively coarse lapilli farther down the dispersal axis.

Conduit stability effects on intensity and steadiness of explosive eruptions

Conduit geometry affects magma ascent dynamics and, consequently, the style and evolution of volcanic eruptions. However, despite geological evidences support the occurrence of conduit widening during most volcanic eruptions, the factors controlling conduit enlargement are still unclear, and the effects of syn-eruptive variations of conduit geometry have not been investigated in depth yet. Based on numerical modeling and the application of appropriate stability criteria, we found out a strong relationship between magma rheology and conduit stability, with significant effects on eruptive dynamics. Indeed, in order to be stable, conduits feeding dacitic/rhyolitic eruptions need larger diameters respect to their phonolitic/trachytic counterparts, resulting in the higher eruption rates commonly observed in dacitic/rhyolitic explosive events. Thus, in addition to magma source conditions and viscosity-dependent efficiency for outgassing, we suggest that typical eruption rates for different magma types are also controlled by conduit stability. Results are consistent with a compilation of volcanological data and selected case studies. As stability conditions are not uniform along the conduit, widening is expected to vary in depth, and three axisymmetric geometries with depth-dependent radii were investigated. They are able to produce major modifications in eruptive parameters, suggesting that eruptive dynamics is influenced by syn-eruptive changes in conduit geometry.

Interventions that may prevent or mollify supervolcanic eruptions

A supervolcanic eruption of 1015 kg could block the sun for years, causing mass starvation or even extinction of some species, including humans. Despite awareness of this problem for several decades, only five interventions have been proposed. In this paper, we increase the number of total possible interventions by more than an order of magnitude. The 64 total interventions involve changing magma characteristics, venting magma, strengthening the cap (rock above the magma), putting more pressure on the magma, stopping an eruption in progress, containing the erupted material, disrupting the plume, or provoking a less intense eruption. We provide qualitative evaluations of the feasibility and risk of 38 of the more promising interventions. The two most promising interventions involve putting more pressure on the magma and delaying the eruption with water dams or soil over the magma chamber. We perform a technical analysis, accurate to within an order of magnitude, and find that water dams and soil and could statistically delay the eruption for a century with 1 and 15 years of effort, respectively. All actions require essentially untested geoengineering challenges along with economic, political and general public acceptance. Further work is required to refine the science, provide cost estimates, and compare cost effectiveness with interventions focusing on adapting to a supereruption.

Scientists Discover Stromboli-Like Eruption on Volcanic Moon

Jupiter’s moon Io is known for its lava fountains and roiling lava lakes, but scientists had never seen such an intense eruption in their data until now.

Earthquake induced variations in extrusion rate: A numerical modeling approach to the 2006 eruption of Merapi Volcano (Indonesia)

Extrusion rates during lava dome-building eruptions are variable and eruption sequences at these volcanoes generally have multiple phases. Merapi Volcano, Java, Indonesia, exemplifies this common style of activity. Merapi is one of Indonesia's most active volcanoes and during the 20th and early 21st centuries effusive activity has been characterized by long periods of very slow (<0.1 m3 s−1) extrusion rate interrupted every few years by short episodes of elevated extrusion rates (1–4 m3 s−1) lasting weeks to months. One such event occurred in May–July 2006, and previous research has identified multiple phases with different extrusion rates and styles of activity. Using input values established in the literature, we apply a 1D, isothermal, steady-state numerical model of magma ascent in a volcanic conduit to explain the variations and gain insight into corresponding conduit processes. The peak phase of the 2006 eruption occurred in the two weeks following the May 27 Mw 6.4 earthquake 50 km to the south. Previous work has suggested that the peak extrusion rates observed in early June were triggered by the earthquake through either dynamic stress-induced overpressure or the addition of CO2 due to decarbonation and gas escape from new fractures in the bedrock. We use the numerical model to test the feasibility of these proposed hypotheses and show that, in order to explain the observed change in extrusion rate, an increase of approximately 5–7 MPa in magma storage zone overpressure is required. We also find that the addition of ∼1000 ppm CO2 to some portion of the magma in the storage zone following the earthquake reduces water solubility such that gas exsolution is sufficient to generate the required overpressure. Thus, the proposed mechanism of CO2 addition is a viable explanation for the peak phase of the Merapi 2006 eruption. A time-series of extrusion rate shows a sudden increase three days following the earthquake. We explain this three-day delay by the combined time required for the effects of the earthquake and corresponding CO2 increase to develop in the magma storage system (1–2 days), and the time we calculate for the affected magma to ascend from storage zone to surface (40 h). The increased extrusion rate was sustained for 2–7 days before dissipating and returning to pre-earthquake levels. During this phase, we estimate that 3.5 million m3 DRE of magma was erupted along with 11 ktons of CO2. The final phase of the 2006 eruption was characterized by highly variable extrusion rates. We demonstrate that those changes were likely controlled by failure of the edifice that had been confining the dome to Merapi's crater and subsequent large dome collapses. The corresponding reductions in confining pressure caused increased extrusion rates that rapidly rebuilt the dome and led to further collapses, a feedback cycle that prolonged the eruption. In a more general sense, this study demonstrates that both internal changes, such as magma volatile content and overpressure, and external forces, such as edifice collapse and regional earthquakes, can affect variations in eruption intensity. Further, we also demonstrate how these external forces can initiate internal changes and how these parameters may interact with one another in a feedback scenario.

Evolution of submarine eruptive activity during the 2011–2012 El Hierro event as documented by hydroacoustic images and remotely operated vehicle observations

AbstractSubmarine volcanic eruptions are frequent and important events, yet they are rarely observed. Here we relate bathymetric and hydroacoustic images from the 2011 to 2012 El Hierro eruption with surface observations and deposits imaged and sampled by ROV. As a result of the shallow submarine eruption, a new volcano named Tagoro grew from 375 to 89 m depth. The eruption consisted of two main phases of edifice construction intercalated with collapse events. Hydroacoustic images show that the eruptions ranged from explosive to effusive with variable plume types and resulting deposits, even over short time intervals. At the base of the edifice, ROV observations show large accumulations of lava balloons changing in size and type downslope, coinciding with the area where floating lava balloon fallout was observed. Peaks in eruption intensity during explosive phases generated vigorous bubbling at the surface, extensive ash, vesicular lapilli and formed high‐density currents, which together with periods of edifice gravitational collapse, produced extensive deep volcaniclastic aprons. Secondary cones developed in the last stages and show evidence for effusive activity with lava ponds and lava flows that cover deposits of stacked lava balloons. Chaotic masses of heterometric boulders around the summit of the principal cone are related to progressive sealing of the vent with decreasing or variable magma supply. Hornitos represent the final eruptive activity with hydrothermal alteration and bacterial mats at the summit. Our study documents the distinct evolution of a submarine volcano and highlights the range of deposit types that may form and be rapidly destroyed in such eruptions.

Translations of volcanological terms: cross-cultural standards for teaching, communication, and reporting

When teaching at a non-English language university, we often argue that because English is the international language, students need to become familiar with English terms, even if the bulk of the class is in the native language. However, to make the meaning of the terms clear, a translation into the native language is always useful. Correct translation of terminology is even more crucial for emergency managers and decision makers who can be confronted with a confusing and inconsistently applied mix of terminology. Thus, it is imperative to have a translation that appropriately converts the meaning of a term, while being grammatically and lexicologically correct, before the need for use. If terms are not consistently defined across all languages following industry standards and norms, what one person believes to be a dog, to another is a cat. However, definitions and translations of English scientific and technical terms are not always available, and language is constantly evolving. We live and work in an international world where English is the common language of multi-cultural exchange. As a result, while finding the correct translation can be difficult because we are too used to the English language terms, translated equivalents that are available may not have been through the peer review process. We have explored this issue by discussing grammatically and lexicologically correct French, German, Icelandic, Indonesian, Italian, Portuguese, Russian, Spanish, and Japanese versions for terms involved in communicating effusive eruption intensity.

The relationship between glaciation and volcanism: Numerical simulation and the Holocene volcanism in Kamchatka

This study is concerned with numerical simulation of the strain due to glaciation and glacial melting, when a magma zone (a layer containing inclusions of magma and magma cumulates) is present at the crust–mantle boundary. According to analytical solutions of this problem that involves viscous relaxation of an uncompensated depression at the place of the molten glacier, the depth to the zone of increased shear stresses beneath the depression is proportional to its width, while the relaxation duration is proportional to viscosity of the lithosphere and is a few thousand years. These fundamental estimates are corroborated by our numerical simulation. According to it, the magma zone at the Moho boundary shields the zone of increased shear stresses, limiting it from below. The maximum values (12–25 MPa) with glacial thickness 500–1000 m are reached at the top of this layer of low viscosity. The directions of maximum compression (s1) as calculated for the time after the melting indicate that the magma that rises along dikes is displaced from the center of the magma lens toward its periphery. It is found that glacial unloading makes the dipping faults in the crust above the low-viscosity layer attractors for the rising magma. Glacial unloading accelerates, by factors of a few times, the magma generation in the mantle that occurs following the mechanism of adiabatic decompression, as well as facilitating the accumulation of mantle fluids in the zone of increased shear stresses at the boundary of the low viscosity layer. The magma traverses this deep fluid collector and increases the intensity and explosivity of eruptions at the beginning of an interglacial period. Our numerical simulation results are in general agreement with published data on Early Holocene volcanic eruptions that occurred after the second phase of the Late Pleistocene glaciation in Kamchatka.

Rapid crystallization during recycling of basaltic andesite tephra: timescales determined by reheating experiments

Microcrystalline inclusions within microlite-poor matrix are surprisingly common in low intensity eruptions around the world, yet their origin is poorly understood. Inclusions are commonly interpreted as evidence of crystallization along conduit margins. Alternatively, these clasts may be recycled from low level eruptions where they recrystallize by heating within the vent. We conducted a series of experiments heating basaltic andesite lapilli from temperatures below the glass transition (~690 °C) to above inferred eruption temperatures (>1150 °C) for durations of 2 to >60 minutes. At 690 °C < T < 800 °C, crystallization is evident after heating for ~20 minutes; at T > 800 °C, crystallization occurs in <5 minutes. At T ≥ 900 °C, all samples recrystallize extensively in 2–10 minutes, with pyroxenes, Fe-oxides, and plagioclase. Experimental crystallization textures closely resemble those observed in natural microcrystalline inclusions. Comparison of inclusion textures in lapilli from the active submarine volcano NW Rota-1, Mariana arc and subaerial volcano Stromboli suggest that characteristic signatures of clast recycling are different in the two environments. Specifically, chlorine assimilation provides key evidence of recycling in submarine samples, while bands of oxides bordering microcrystalline inclusions are unique to subaerial environments. Correct identification of recycling at basaltic vents will improve (lower) estimates of mass eruption rate and help to refine interpretations of eruption dynamics.

Climatic response of extratropics to intense tropical eruptions

Climatic response of extratropics to intense tropical eruptions

Retrieval and intercomparison of volcanic SO2 injection height and eruption time from satellite maps and ground-based observations

Syneruptive gas flux time series can, in principle, be retrieved from satellite maps of SO2 collected during and immediately after volcanic eruptions, and used to gain insights into the volcanic processes which drive the volcanic activity. Determination of the age and height of volcanic plumes are key prerequisites for such calculations. However, these parameters are challenging to constrain using satellite-based techniques. Here, we use imagery from OMI and GOME-2 satellite sensors and a novel numerical procedure based on back-trajectory analysis to calculate plume height as a function of position at the satellite measurement time together with plume injection height and time at a volcanic vent location. We applied this new procedure to three Etna eruptions (12 August 2011, 18 March 2012 and 12 April 2013) and compared our results with independent satellite and ground-based estimations. We also compare our injection height time-series with measurements of volcanic tremor, which reflects the eruption intensity, showing a good match between these two datasets. Our results are a milestone in progressing towards reliable determination of gas flux data from satellite-derived SO2 maps during volcanic eruptions, which would be of great value for operational management of explosive eruptions.

Formation of obsidian pyroclasts by sintering of ash particles in the volcanic conduit

The ranges in intensity and style of volcanic eruptions, from highly explosive Plinian eruptions to quiescent lava extrusions, depend on the style and efficiency of gas loss from ascending magma. Obsidian pyroclasts – small, glassy pieces of quenched magma found in some volcanic tephra beds – may preserve valuable information about magma degassing in their vesicle textures and volatile contents. Accurate interpretation of their textures and volatiles, however, requires understanding the mechanism of formation of the pyroclasts. Obsidian pyroclasts from the ca. 1325–1350 C.E. North Mono eruption of Mono Craters (CA, USA) were analyzed and found to have H2O and CO2 contents indicating that they were formed at pressures in the approximate range of 3–40 MPa. Many also contain domains with differing vesicle textures, separated by boundaries containing xenocrystic material, indicating that they are composed of smaller fragments that have sutured together. More than half of the pyroclasts analyzed contained small (∼10 μm), highly distorted vesicles, with multi-cuspate morphology, interpreted as the remnants of interstitial gas trapped amongst sintered fragments of melt/glass. Rounded vesicles are also common and are interpreted to result from surface tension-driven relaxation of the distorted vesicles. Calculated timescales of sintering and relaxation are consistent with timescales for pyroclast formation indicated by H2O re-equilibration within the heterogeneous pyroclasts. This sintering model for the origin of obsidian pyroclasts is further supported by the observation that spherical vesicles are found mainly in H2O-rich pyroclasts, and distorted vesicles mainly in H2O-poor pyroclasts. We conclude that obsidian pyroclasts generated during the North Mono eruption were formed by cycles of fragmentation, sintering/suturing, and relaxation, over a very wide range of depths within the conduit; we find no evidence to support pumice (foam) collapse as the formation mechanism. Similar textures, and the occurrence of xenolithic material, in obsidian pyroclasts in other eruptions suggest that sintering may be generally responsible for the origin of obsidian pyroclasts. Our conceptual model indicates that volatile contents in obsidian pyroclasts reflect both degassing of bubbly magma and the composition of gas trapped between sintering particles.

Stress Field Control during Large Caldera-Forming Eruptions

Crustal stress field can have a significant influence on the way magma is channelled through the crust and erupted explosively at the surface. Large Caldera Forming Eruptions (LCFEs) can erupt hundreds to thousands of cubic kilometres of magma in a relatively short time along fissures under the control of a far-field extensional stress. The associated eruption intensities are estimated in the range 109 - 1011 kg/s. We analyse syn-eruptive dynamics of LCFEs, by simulating numerically explosive flow of magma through a shallow dyke conduit connected to a magma chamber that in turn is fed by a deeper magma reservoir, both under the action of an extensional far-field stress. Results indicate that huge amounts of high viscosity silicic magma can be erupted over timescales of a few to several hours. Our study provides answers to outstanding questions relating to the intensity and duration of catastrophic volcanic eruptions in the past. In addition, it presents far-reaching implications for the understanding of dynamics and intensity of large-magnitude volcanic eruptions on Earth and to highlight the necessity of a future research to advance our knowledge of these rare catastrophic events.

Ionospheric detection and localization of volcano eruptions on the example of the April 2015 Calbuco events

AbstractUsing data from ground‐based Global Navigation Satellite Systems (GNSS) receivers located in southern Chile, we study the ionospheric total electron content (TEC) response to two eruptions of the Calbuco volcano that occurred on 22–23 April 2015. In both cases, the TEC response showed quasi‐periodic signals with several consecutive wave trains. The averaged amplitude of the observed covolcanic TEC perturbations amounted 0.45 total electron content unit, 1 TECU = 1016 el m−2 (TECU) for the first eruption and 0.16 TECU for the second one. We compare amplitudes of the TEC response to volcano eruptions of different intensity from our and previously published data, and we show that both the intensity and the background ionospheric conditions define the amplitude of ionospheric covolcanic disturbances. The relative contribution, however, scales with the eruption intensity. The traveltime diagrams allowed to estimate the propagation speed of the observed covolcanic TEC perturbations as ~900–1200 m/s, which is close to the acoustic (or shock acoustic) waves speed at the ionospheric height. The spectrograms are consistent with the conclusion on the acoustic nature of the observed TEC perturbations. Finally, we use the approximation of a spherical wave propagating at a constant velocity from a point source, and for the first time, we calculate the location of the volcanic source and the onset time of the volcano eruption from ionospheric measurements. We show that even from 30 s ionospheric GPS data it is possible to “localize” the eruptive source within several degrees of latitude/longitude.

Magma decompression rates during explosive eruptions of Kīlauea volcano, Hawaii, recorded by melt embayments

The decompression rate of magma as it ascends during volcanic eruptions is an important but poorly constrained parameter that controls many of the processes that influence eruptive behavior. In this study, we quantify decompression rates for basaltic magmas using volatile diffusion in olivine-hosted melt tubes (embayments) for three contrasting eruptions of Kīlauea volcano, Hawaii. Incomplete exsolution of H2O, CO2, and S from the embayment melts during eruptive ascent creates diffusion profiles that can be measured using microanalytical techniques, and then modeled to infer the average decompression rate. We obtain average rates of ~0.05–0.45 MPa s−1 for eruptions ranging from Hawaiian style fountains to basaltic subplinian, with the more intense eruptions having higher rates. The ascent timescales for these magmas vary from around ~5 to ~36 min from depths of ~2 to ~4 km, respectively. Decompression-exsolution models based on the embayment data also allow for an estimate of the mass fraction of pre-existing exsolved volatiles within the magma body. In the eruptions studied, this varies from 0.1 to 3.2 wt% but does not appear to be the key control on eruptive intensity. Our results do not support a direct link between the concentration of pre-eruptive volatiles and eruptive intensity; rather, they suggest that for these eruptions, decompression rates are proportional to independent estimates of mass discharge rate. Although the intensity of eruptions is defined by the discharge rate, based on the currently available dataset of embayment analyses, it does not appear to scale linearly with average decompression rate. This study demonstrates the utility of the embayment method for providing quantitative constraints on magma ascent during explosive basaltic eruptions.

Ground-coupled airwaves at Pavlof Volcano, Alaska, and their potential for eruption monitoring

An abnormally high number of explosion quakes were noted during the monitoring effort for the 2007 eruption of Pavlof Volcano on the Alaska Peninsula. In this study, we manually cataloged the explosion quakes from their characteristic ground-coupled airwaves. This study investigates how the ground-coupled airwaves might be used in a monitoring or analysis effort by estimating energy release and gas mass release. Over 3 × 104 quakes were recorded. The energy release from the explosions is approximated to be 3 × 1011 J, and the total gas mass (assuming 100 % water) released was 450 t. The tracking of explosion quakes has the potential to estimate relative eruption intensity as a function of time and is thus a useful component of a seismic monitoring program.

Magma reservoirs from the upper crust to the Moho inferred from high-resolution Vp and Vs models beneath Mount St. Helens, Washington State, USA

The size, frequency, and intensity of volcanic eruptions are strongly controlled by the volume and connectivity of magma within the crust. Several geophysical and geochemical studies have produced a comprehensive model of the magmatic system to depths near 7 km beneath Mount St. Helens (Washington State, USA), currently the most active volcano in the Cascade Range. Data limitations have precluded imaging below this depth to observe the entire primary shallow magma reservoir, as well as its connection to deeper zones of magma accumulation in the crust. The inversion of P and S wave traveltime data collected during the active-source component of the iMUSH (Imaging Magma Under St. Helens) project reveals a high P-wave (Vp)/S-wave (Vs) velocity anomaly beneath Mount St. Helens between depths of 4 and 13 km, which we interpret as the primary upper–middle crustal magma reservoir. Beneath and southeast of this shallow reservoir, a low Vp velocity column extends from 15 km depth to the Moho. Deep long-period events near the boundary of this column indicate that this anomaly is associated with the injection of magmatic fluids. Southeast of Mount St. Helens, an upper–middle crustal high Vp/Vs body beneath the Indian Heaven Volcanic Field may also have a magmatic origin. Both of these high Vp/Vs bodies are at the boundaries of the low Vp middle–lower crustal column and both are directly above high Vp middle–lower crustal regions that may represent cumulates associated with recent Quaternary or Paleogene–Neogene Cascade magmatism. Seismicity immediately following the 18 May 1980 eruption terminates near the top of the inferred middle–lower crustal cumulates and directly adjacent to the inferred middle–lower crustal magma reservoir. These spatial relationships suggest that the boundaries of these high-density cumulates play an important role in both vertical and lateral transport of magma through the crust.

On the relationship between eruption intensity and volcanic plume height: Insights from three-dimensional numerical simulations

Height of plumes generated during explosive volcanic eruptions is commonly used to estimate the associated eruption intensity (i.e., mass eruption rate; MER). In order to quantify the relationship between plume height and MER, we performed a parametric study using a three-dimensional (3D) numerical model of volcanic plumes for different vent sizes. The results of five simulations indicate that the flow pattern in the lower region of the plume systematically changes with vent size, and hence, with MER. For MERs <4×107kgs−1, the flow in the lower region has a jet-like structure (the jet-like regime). For MERs >108kgs−1, the flow shows a fountain-like structure (the fountain-like regime). The flow pattern of plumes with 4×107kgs−1<MERs <108kgs−1 shows transitional features between the two flow regimes. Within each of the two flow regimes, the plume height increases as the MER increases, whereas plume heights remain almost constant or even decrease as MER increases in the transitional regime; as a result, the jet-like and fountain-like regimes show distinct relationships of plume height and MER. The different relationships between the two regimes reflect the fact that the efficiency of entrainment of ambient air in the jet-like regime is substantially lower than that in the fountain-like regime. It is suggested that, in order to estimate eruption intensity from the observed plume heights, it is necessary to take the different flow regimes depending on MER into account.