Eruption Column Height Research Articles

Unravelling eruption dynamics at El Reventador, Ecuador: linking the physiochemical properties of volcanic ash with geophysical, geochemical and satellite remote sensing data

The physiochemical properties of volcanic ash are determined by magma ascent and eruption dynamics and provide important insights into controls on the timing and style of volcanic eruptions. However, linking petrological observations to monitoring parameters remains challenging. Here we investigate the relationships between geophysical, geochemical and satellite remote sensing data with the physiochemical properties of volcanic ash to better understand eruptive dynamics at El Reventador volcano, Ecuador. Between 2016 and 2019, eruptive activity at El Reventador was characterised by frequent explosions interspersed with effusive activity from two summit vents. We found that: (1) periods of predominantly effusive behaviour were defined by ash samples with the lowest proportions of juvenile grains, frequent, more intense thermal anomalies and frequent (ca. every 25 min) but low-energy explosions that produced lower average ash column heights (< 1200 m above the summit) and; (2) periods of predominantly explosive behaviour were defined by ash samples with the highest proportions of juvenile grains and also vesicular juvenile grains, low thermal anomalies and less frequent (ca. every 70 min), higher-energy explosions that produced higher average eruption column heights (> 1500 m above the summit). Our study shows that there are strong correlations between the physiochemical properties of erupted material and the multiparametric monitoring data which in turn link to the type of surface activity observed. As a result, the development of physical conduit models and interpretations of eruptive dynamics are made more robust combining both monitoring data and the physical properties of eruptive products.

Evolution of a large-scale phreatoplinian eruption: Constraints from the 40 ka caldera-forming eruption of Kutcharo volcano, eastern Hokkaido, Japan

“Phreatoplinian” is an explosive phreatomagmatic eruption style that is defined by the fragmentation of magma and widespread dispersal of the resulting fine ash and accretionary lapilli. These eruptions pose significant future risks at caldera volcanoes that host lakes and abundant groundwater. There have been no direct observations of a phreatoplinian eruption, therefore, constraining the detailed mechanisms and sequences of such events relies on studying the deposits of previous eruptions. In order to advance our understanding of these hazardous phenomena we conducted a case study of the 40 ka caldera-forming eruption (Kp I) from Kutcharo volcano in eastern Hokkaido, Japan. We subdivided Kp I eruption deposits into 7 units (Units 1 to 7 in ascending order). Units 1 to 6 are air fall deposits consisting of alternating thin pumice and thick silty ash layers with abundant spherical accretionary lapilli. Stratigraphically higher ash fall units are thicker, finer in grain-size, and more widely distributed. The maximum eruption column height and mass-discharge rate were calculated to be 40 km and 1.4 × 109 kg/s, respectively. Unit 7 is a climactic ignimbrite (76 km3), which is distributed widely over the area north of Kutcharo caldera.Unit 6 is the largest air fall unit and can be considered to have been deposited by a phreatoplinian eruption, given its abundant accretionary lapilli, wide dispersion, and high degree of fragmentation. Unit 6 had the highest mass discharge rate (1.4 × 109 kg/s), suggesting the interaction between magma and external water was most intense, and it is thought that a large eruption column covered eastern Hokkaido. In addition, Kp I eruption deposits commonly contain glass shards derived from fragmentation via both magma degassing and Molten Fuel Coolant Interaction (MFCI). To account for this observation, we infer that the conduit penetrated a large aquifer, and the margin of the ascending magma came into contact with this external water source. Due to repeated caldera-forming eruptions, intra-caldera filled deposits (hosting a large aquifer) likely played a key role in supplying external caldera lake water to a level near the fragmentation depth of H2O-saturated felsic magma. The occurrence of these intra-caldera conduit and caldera-lake systems may provide the required conditions for phreatoplinian eruptions at continental arc caldera volcanoes in Japan and globally.

The 23–24 March 2021 lava fountain at Mt Etna, Italy

In 2021, more than 50 paroxysmal episodes occurred at the South-East Crater (SEC) of Mt Etna, Italy. The 23–24 March lava fountain was one of the longest episodes and began with weak Strombolian explosions, gradually transitioning to lava fountaining. The eruption intensity then dropped more slowly than in previous episodes, resulting in pulsating Strombolian explosions dominated by ash emission. Thirty-four tephra samples were used to reconstruct the fallout dispersal and estimate the total erupted mass. Grain size, textural, petrological and geochemical analyses indicate different features and were compared with the gas phase (SO2 and HCl) in the volcanic plume. By applying stochastic global optimization to simulations of the temporal evolution of the eruption column height and tephra dispersal and deposition, the total erupted mass retrieved (6.76 × 108 kg) matches well the total erupted mass estimation by the ground-based deposit (8.03 ± 2.38 × 108 kg), reducing the column height throughout the episode from 6.44 to 4.5 km above sea level and resulting in a mass eruption rate ranging from 1.96 × 105 to 8.18 × 103 kg/s. The unusual duration of the March episode and the characteristics of the erupted products point to the change in explosive style and magma fragmentation from fountaining to ash emission phases, associated with a slower magma supply inducing a change in magma rheology and a final, prolonged ash generation. Furthermore, this study showed that using observational data and the variation in eruption source parameters for numerical simulations can improve the accuracy of predicting the dispersal plume, thus mitigating the potential impact of longer paroxysmal episodes.

Radiocarbon dating, magnitude, and catastrophic impact of Popocatépetl's Terminal Preclassic “Lorenzo pumice” Plinian eruption and paleomagnetic age constraints on the Nealtican lava flow field (central Mexico)

The “Lorenzo pumice” (LP) fallout was produced by a Plinian (VEI = 4) eruption of Popocatépetl volcano ∼40 cal BCE during the Terminal Preclassic period. It is part of a pyroclastic sequence that consists in stratigraphic order of fall, pyroclastic density current, and lahar deposits, which in turn are partly covered by the 70 km2 andesitic-dacitic “Pedregal de Nealtican” lava flow field on the eastern slope of the volcano.The composition of the pumice varies insignificantly (SiO2 = 61.0 to 63.5 wt%) with phenocrysts of euhedral plagioclase > orthopyroxene > clinopyroxene ± olivine ± hornblende ± opaques.Sizes of lithics and the isopach map indicate an eruption column height of ∼28 km and a dispersal axis toward the ENE, where an area of ∼600 km2 that included pre-Hispanic dispersed villages was buried under >10 cm of pumice. The minimum volume for the LP fallout was calculated at ∼0.6 km3 corresponding to ∼0.2 km3 of dense rock equivalent (DRE).The short Plinian phase ended when the column collapsed abruptly, generating lethal pyroclastic density currents that were emplaced around the volcano. It was followed by periodic flooding by rain lahars that originated from the slopes of Popocatépetl and from neighboring Iztaccíhuatl, severely affecting ravines and large portions of the plains of Puebla in the east.After the cataclysmic phase, the El Ombligo vent alignment opened at ∼4000 m a.s.l. on Popocatépetl's NE upper slope from which the lavas that form the Nealtican lava flow field emanated, an effusive activity that lasted >30 years. Paleomagnetic data for these lavas and stratigraphic evidence indicate that the lava flows with a total volume of 4.2 km3 were emplaced shortly after the explosive Plinian phase and represent the degassed portion of the LP magma.This eruption had a considerable impact on pre-Hispanic settlements around the volcano, whose population exodus and relocation probably contributed to the rise of important cities in central Mexico, such as Teotihuacán, Cholula, and Cantona.

An automated ash dispersion forecast system: case study Popocatépetl volcano, Mexico

An operational volcanic ash dispersion forecast system was developed for Popocatépetl. It runs automatically every day developing 108 possible scenarios of ash dispersion for the following 36 h. Scenarios are simulated for three eruption column heights: 3 km, 5 km, and 10 km above the volcano’s crater level, every hour for eruptions lasting 1 h. For each hypothetical eruption that starts every hour, the dispersion during the following 8 h is modelled. The system uses the Weather Research and Forecasting (WRF) model for weather data and the Fall3D model. It includes a visualization website that displays, among other products: ground accumulation, deposit load, and concentration at relevant flight levels. Popocatépetl volcano, located ~ 60 km from Mexico Megacity was selected as a case study. A comparison from ash forecast system results and satellite observations is presented. The system developed and tested here can be adapted to be operative at any volcano.

Estimation of Volume and Column Height from Pumiceous Tephra-Fall Deposits of Mt. Raung, East Java, Indonesia

Pumiceous tephra-fall deposits observed mainly in the western flank of Mt. Raung are indicating past explosive eruption event. The volume and the eruption column height of the eruptive event is estimated using isopach and isopleth maps. These maps are constructed from seven observation points. The pumiceous tephra-fall deposits can be divided into 9 layers from bottom to top, White Pumice Fall 1–7, Yellow Pumice Fall, and Lithic-Rich Scoria Fall. These layers not separated by any soil or weathered layer indicates no significant time gap between them. Isopach and isopleth maps are generated from White Pumice Fall 6, Yellow Pumice Fall, and Lithic-Rich Scoria Fall considering the widespread distribution of these deposits. White Pumice Fall, Yellow Pumice Fall, and Lithic-Rich Scoria Fall are resulted from eruptive event with ~25–33 km, ~35 km, and ~29 km column height and 1.3 km3, 0.73 km3, and 0.32 km3 volume, respectively. Total volume of 2.35 km3 and the column height of 25-35 km suggest that Plinian eruption type of VEI 5 magnitude produces this pumiceous tephra-fall deposit sequence.

Modeling Tephra Fall and Sediment-Water Flows to Assess Their Impacts on a Vulnerable Building Stock in the City of Arequipa, Peru

Arequipa, Peru’s second economic center hosting c. 1,110,000 inhabitants, is the largest South American city exposed to a large variety of natural hazards. At least 200,000 live in areas likely to be affected by hazards from El Misti volcano, located 17 km to the NE. A multidisciplinary project aims to address the impacts of tephra fall and frequent mass flows on the vulnerable building stock and roofs along two ravines that cross the city, enabling decision-makers to undertake retrofitting projects and improve urban risk planning. Two recent eruptions, that is, the 1440–1470 CE Vulcanian event and c. 2070 years BP Plinian eruption, were chosen as references for probable scenarios of potential tephra fall impacts from El Misti on the building roofs. Tephra fall impacts on the city depend on the eruptive style, column height, and patterns of wind directions and velocities over south Peru and roof mechanical resistance. Estimates of potential damage levels and cost range values rely on nine structural types and four classes of vulnerable roofs. Simulation runs of hyperconcentrated flows (HCF) and debris flows (DF), using three depth-averaged flow models (Titan2F, VolcFlow, and Flo-2D) along two drainage basins on the SW flank of El Misti and across Arequipa, examined three scenarios from a database of 39 recent events and other historical lahars. Simulation results showcase the extent toward the city, inundation depths ≤4.6 m, flow velocities ≤9 m/s, and dynamic pressure up to 100 kPa from three different magnitude HCFs and DFs. In both ravines, overbank flows occurred in key urban areas due to channel sinuosity and constrictions near bridges. Potential impacts on habitat stem from ranges of flow dynamic pressure and measurements of construction material. We estimated the monetary loss of buildings according to hyperconcentrated flows and debris flows scenarios to contribute to retrofitting procedure, implementation of defense work, and relocation policy.

Data assimilation of volcanic aerosol observations using FALL3D+PDAF

Abstract. Modelling atmospheric dispersal of volcanic ash and aerosols is becoming increasingly valuable for assessing the potential impacts of explosive volcanic eruptions on buildings, air quality, and aviation. Management of volcanic risk and reduction of aviation impacts can strongly benefit from quantitative forecasting of volcanic ash. However, an accurate prediction of volcanic aerosol concentrations using numerical modelling relies on proper estimations of multiple model parameters which are prone to errors. Uncertainties in key parameters such as eruption column height and physical properties of particles or meteorological fields represent a major source of error affecting the forecast quality. The availability of near-real-time geostationary satellite observations with high spatial and temporal resolutions provides the opportunity to improve forecasts in an operational context by incorporating observations into numerical models. Specifically, ensemble-based filters aim at converting a prior ensemble of system states into an analysis ensemble by assimilating a set of noisy observations. Previous studies dealing with volcanic ash transport have demonstrated that a significant improvement of forecast skill can be achieved by this approach. In this work, we present a new implementation of an ensemble-based data assimilation (DA) method coupling the FALL3D dispersal model and the Parallel Data Assimilation Framework (PDAF). The FALL3D+PDAF system runs in parallel, supports online-coupled DA, and can be efficiently integrated into operational workflows by exploiting high-performance computing (HPC) resources. Two numerical experiments are considered: (i) a twin experiment using an incomplete dataset of synthetic observations of volcanic ash and (ii) an experiment based on the 2019 Raikoke eruption using real observations of SO2 mass loading. An ensemble-based Kalman filtering technique based on the local ensemble transform Kalman filter (LETKF) is used to assimilate satellite-retrieved data of column mass loading. We show that this procedure may lead to nonphysical solutions and, consequently, conclude that LETKF is not the best approach for the assimilation of volcanic aerosols. However, we find that a truncated state constructed from the LETKF solution approaches the real solution after a few assimilation cycles, yielding a dramatic improvement of forecast quality when compared to simulations without assimilation.

Estimating Volcanic Eruption Column Height and Growth Rate Using X-Band Marine Radar at the Sakurajima Volcano

Estimating Volcanic Eruption Column Height and Growth Rate Using X-Band Marine Radar at the Sakurajima Volcano

Geometric estimation of volcanic eruption column height from GOES-R near-limb imagery – Part 2: Case studies

Abstract. In a companion paper (Horváth et al., 2021), we introduced a new technique to estimate volcanic eruption column height from extremely oblique near-limb geostationary views. The current paper demonstrates and validates the technique in a number of recent eruptions, ranging from ones with weak columnar plumes to subplinian events with massive umbrella clouds and overshooting tops that penetrate the stratosphere. Due to its purely geometric nature, the new method is shown to be unaffected by the limitations of the traditional brightness temperature method, such as height underestimation in subpixel and semitransparent plumes, ambiguous solutions near the tropopause temperature inversion, or the lack of solutions in undercooled plumes. The side view height estimates were in good agreement with plume heights derived from ground-based video and satellite stereo observations, suggesting they can be a useful complement to established techniques.

Geometric estimation of volcanic eruption column height from GOES-R near-limb imagery – Part 1: Methodology

Abstract. A geometric technique is introduced to estimate the height of volcanic eruption columns using the generally discarded near-limb portion of geostationary imagery. Such oblique observations facilitate a height-by-angle estimation method by offering close-to-orthogonal side views of eruption columns protruding from the Earth ellipsoid. Coverage is restricted to daytime point estimates in the immediate vicinity of the vent, which nevertheless can provide complementary constraints on source conditions for the modeling of near-field plume evolution. The technique is best suited to strong eruption columns with minimal tilting in the radial direction. For weak eruptions with severely bent plumes or eruptions with expanded umbrella clouds the radial tilt/expansion has to be corrected for either visually or using ancillary wind profiles. Validation on a large set of mountain peaks indicates a typical height uncertainty of ±500 m for near-vertical eruption columns, which compares favorably with the accuracy of the common temperature method.

Climate change modulates the stratospheric volcanic sulfate aerosol lifecycle and radiative forcing from tropical eruptions

Explosive volcanic eruptions affect climate, but how climate change affects the stratospheric volcanic sulfate aerosol lifecycle and radiative forcing remains unexplored. We combine an eruptive column model with an aerosol-climate model to show that the stratospheric aerosol optical depth perturbation from frequent moderate-magnitude tropical eruptions (e.g. Nabro 2011) will be reduced by 75% in a high-end warming scenario compared to today, a consequence of future tropopause height rise and unchanged eruptive column height. In contrast, global-mean radiative forcing, stratospheric warming and surface cooling from infrequent large-magnitude tropical eruptions (e.g. Mt. Pinatubo 1991) will be exacerbated by 30%, 52 and 15% in the future, respectively. These changes are driven by an aerosol size decrease, mainly caused by the acceleration of the Brewer-Dobson circulation, and an increase in eruptive column height. Quantifying changes in both eruptive column dynamics and aerosol lifecycle is therefore key to assessing the climate response to future eruptions.

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): A new observational database to support explosive eruptive column model validation and development

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Design of parametric risk transfer solutions for volcanic eruptions: an application to Japanese volcanoes

Abstract. Volcanic eruptions are rare but potentially catastrophic phenomena, affecting societies and economies through different pathways. The 2010 Eyjafjallajökull eruption in Iceland, a medium-sized ash-fall-producing eruption, caused losses in the range of billions of dollars, mainly to the aviation and tourism industries. Financial risk transfer mechanisms such as insurance are used by individuals, companies, governments, etc., to protect themselves from losses associated with natural catastrophes. In this work, we conceptualize and design a parametric risk transfer mechanism to offset losses to building structures arising from large, ash-fall-producing volcanic eruptions. Such a transfer mechanism relies on the objective measurement of physical characteristics of volcanic eruptions that are correlated with the size of resulting losses (in this case, height of the eruptive column and predominant direction of ash dispersal) in order to pre-determine payments to the risk cedent concerned. We apply this risk transfer mechanism to the case of Mount Fuji in Japan by considering a potential risk cedent such as a regional government interested in offsetting losses to dwellings in the heavily populated prefectures of Tokyo and Kanagawa. The simplicity in determining eruptive column height and ash fall dispersal direction makes this design suitable for extrapolation to other volcanic settings worldwide where significant ash-fall-producing eruptions may occur, provided these parameters are reported by an official, reputable agency and a suitable loss model is available for the volcanoes of interest.

The Arce Tephra: Two subsequent paroxysmal Plinian eruptions from Coatepeque Caldera (El Salvador)

The Coatepeque volcanic complex in El Salvador produced at least four Plinian eruptions within the last 80 kyr. The eruption of the 72 ka old Arce Tephra formed the Coatepeque Caldera and was one of the most powerful explosive eruptions in El Salvador. Hitherto it was thought that the Arce tephra had been emplaced only by one, mostly Plinian, eruptive event that ended with the deposition of a thick ignimbrite.However, our stratigraphic, geochemical, and zircon data reveal a temporally closely- spaced double eruption separated by a gap of only a couple of hundred years, and we therefore distinguish Lower and Upper Arce Tephras. Both eruptions produced in the beginning a series of fallout units generated from fluctuating eruption columns and turning wind directions. The final phase of the Upper Arce eruption produced surge deposits by several eruption column collapses before the terminal phase of catastrophic ignimbrite eruption and caldera collapse.Mapping of the individual tephra units including the occurrences of distal marine and lacustrine ash layers in the Pacific Ocean, the Guatemalan lowlands and the Caribbean Sea, result in 25.6km3 tephra volume, areal distribution of 4×105km2 and eruption column heights between 20–33km for the Lower Arce eruption, and 40.5km3 tephra volume, including 10km3 for the ignimbrite, distributed across 6×105km2 and eruption column heights of 23–28km for the Upper Arce eruption. These values and the detailed eruptive sequence emphasize the great hazard potential of possible future highly explosive eruptions at Coatepeque Caldera, especially for this kind of double eruption.

Eruption dynamics of the 23 February 2013 event at Mt. Etna

Volcanic activity at Mt. Etna in the last decade has mostly been manifested by sequences of short paroxysmal episodes characterised by powerful lava fountains and high eruption columns. On the 23 February 2013, an exceptionally intense episode occurred at the New South-East Crater, producing a fountain >800 m high (among the highest ever recorded at Etna) and a ~9 km eruption column that dispersed ash >400 km from the vent. Textural and petrographic analyses of lapilli revealed that magma erupted during the high-intensity phase is characterised by low microlite contents (<7 area%), high vesicularity (76–83%), and high vesicle number densities (6–8.2 × 106 cm−3). The short-lived initial Strombolian explosions removed viscous magma from the conduit, enabling the rapid ascent of gas-rich, microlite-poor magma and the eruption of an 800 m high fountain and 9 km high eruption column. For the 23 February eruption, the high vesicularity and low microlite content of the pyroclasts support the hypothesis that volatile-rich magma was the driver of the high intensity lava fountain. This eruptive event, along with three other recent events at Etna over the last 15 years, can be defined as subplinian based on eruption rate and column height, but also generated incandescent 800–1000 m fountains. For these reasons, we propose to term this event, and others at Etna characterised by similar eruption features and parameters, as subplinian fountaining events.

Exploring How Eruption Source Parameters Affect Volcanic Radiative Forcing Using Statistical Emulation

AbstractThe radiative forcing caused by a volcanic eruption is dependent on several eruption source parameters such as the mass of sulfur dioxide (SO2) emitted, the eruption column height, and the eruption latitude. General circulation models with prognostic aerosol and chemistry schemes can be used to investigate how each parameter influences the volcanic forcing. However, the range of multidimensional parameter space that can be explored is restricted because such simulations are computationally expensive. Here we use statistical emulation to explore the radiative impact of eruptions over a wide covarying range of SO2 emission magnitudes, injection heights, and eruption latitudes based on only 30 simulations. We use the emulators to build response surfaces to visualize and predict the sulfate aerosol e‐folding decay time, the stratospheric aerosol optical depth, and net radiative forcing of thousands of different eruptions. We find that the volcanic stratospheric aerosol optical depth and net radiative forcing are primarily determined by the mass of SO2 emitted, but eruption latitude is the most important parameter in determining the sulfate aerosol e‐folding decay time. The response surfaces reveal joint effects of the eruption source parameters in influencing the net radiative forcing, such as a stronger influence of injection height for tropical eruptions than high‐latitude eruptions. We also demonstrate how the emulated response surfaces can be used to find all combinations of eruption source parameters that produce a particular volcanic response, often revealing multiple solutions.

Tephra From the 3 March 2015 Sustained Column Related to Explosive Lava Fountain Activity at Volcán Villarrica (Chile)

Villarrica Volcano in the Southern Andes of Chile has sustained persistent open conduit activity, associated with a dynamic lava lake since at least 1984-85. Increased seismic activity and degassing in August 2014 was followed by increased Strombolian activity starting in February 2015 and culminated with a 1.5 km-high lava fountain on 3 March 2015. This eruption produced tephra fallout, spatter agglutination, clastogenic lavas and mixed avalanche deposits. Field work carried out in 2017 and 2018 focused upon the tephra fallout deposits from the 3 March 2015 paroxysmal eruption. Tephra fall deposits produced by lava fountain and a rapid transition to sustained eruption column at Villarrica are here described for first time. Tephra dispersion was observed mostly southeast of the crater in a narrow, elongated, nearly-elliptic area. Minimum observed tephra loading was estimated to be 80-120 g/m2 at crosswind locations, while maximum load estimation was about 11,500 g/m2, at 7-9 km distance from the vent. At 6-8 km downwind, the tephra fall deposit consisted of a single black scoria layer, fall varied from a lapilli blanket at medial distances (9-19 km) to a few individual particles at distal sites (20-43 km). All samples are vitrophyric in texture, with a mineral assemblage of plagioclase (3-30 %), clinopyroxene (2-5 %), olivine (1-4 %) and scarce orthopyroxene (<1 %) phenocrysts. We estimate that 1.4 x 109 kg, equivalent to ~2.4 x 106 m3 of basaltic andesite tephra (52.8 to 54.7 SiO2 wt. %), were erupted within a period of about 16 minutes at a mass eruption rate of ~1.5 x 106 kg/s. The sustained eruption column height was estimated at ~10.8 km, with magnitude and intensity of 2.11 and 9.13 respectively. We propose that a rejuvenated, volatile-rich magmatic intrusion is the probable triggering factor for the paroxysm. The Villarrica eruption is classified as a short-lived, large-scale lava fountain, similar to intense historical lava fountains at Mt. Etna and Kilauea. As evident by devastating 20th century Villarrica eruptions which also developed similar lava fountain behavior, this eruption style should be considered as frequent and a potential source of volcanic hazard in the future.

Eruption dynamics and explosive-effusive transitions during the 1400 cal BP eruption of Opala volcano, Kamchatka, Russia

Deposits and pumice from the 1400 cal BP eruption of Opala volcano record activity that occurred at the explosive-effusive transition, resulting in intermittent, or stop-start, behavior, where explosive activity resumed following a pause. The eruption deposited distinctive, biotite-bearing rhyolite tephra across much of Kamchatka, and its stratigraphy consists of a lithic-rich pumice fall, overlain by pumice falls and pyroclastic density deposits, with the proportion of the latter increasing with height. This sequence repeats such that the middle of the total deposit is marked by a lithic-rich fall with abundant obsidian clasts. Notably, the eruptive pumice are poorly vesiculated, with vesicle textures that record fragmentation of a partially collapsed magmatic foam. The eruption vent, Baranii Amphitheater is filled with obsidian lavas of the same composition as the rhyolite tephra. Based upon the stratigraphic and compositional relations, we divide the eruption into four phases. Phase I initiated with eruption of a lithic-rich pumice fall, followed by eruption of Plinian falls and pyroclastic density currents. During Phase II, the eruption paused for at least 5–6h; in this time, microlites nucleated and began to grow in the magma. Phase III essentially repeated the Phase I sequence. Obsidian lavas were emplaced during Phase IV. The pumice textures suggest that the magma ascended very near the threshold decompression rate for the transition between explosive (fast) and effusive (slow) behavior. The pause during Phase II likely occurred as decompression slowed enough for the magma to develop sufficient permeability for gas to escape resulting in collapse of the magmatic foam, stopping the eruption and temporarily sealing the conduit. After about 5–6h, eruption resumed with, once again, magma decompressing very near the explosive-effusive transition. Phase III ended when the decompression rate slowed and lava dome emplacement began. Distributions of pumice and lithic clasts, and inclusion of data from previous workers, indicate minimum deposit volumes of ~0.75 and ~0.75–1.15km3 (DRE) and eruption column heights of ~18 and ~20km for Phases I and III, respectively. Phases I-III had a likely total duration of ~60–80h, including a pause in activity of 5–6h during Phase II. This study demonstrates that analysis of vesicle textures from numerous pumice combined with stratigraphic data can reveal syn-eruptive changes in and links between magma permeability, decompression rate, and eruption style.

A08 古記録と景観再現CGにもとづく1707年富士山宝永噴火の噴煙柱高度推定(噴火と噴煙のダイナミクス,口頭発表)

A08 古記録と景観再現CGにもとづく1707年富士山宝永噴火の噴煙柱高度推定(噴火と噴煙のダイナミクス,口頭発表)