Vulcanian Research Articles

Explosive eruption processes inferred from high-frequency seismic waveforms of eruption tremor and explosion events

SUMMARY We investigated the relation between high-frequency seismic signals and eruption size and duration using seismic data of eruption tremor and explosion events generated during sub-Plinian and Vulcanian eruptions, respectively, at various volcanoes. We estimated source amplitude functions from seismic envelope seismograms in the 5–10 Hz band, in which S waves are assumed to radiate isotropically. Because seismic data associated with explosive eruptions can be contaminated by infrasound signals, we confirmed that contamination did not significantly affect the source amplitude functions quantified from our analysed waveforms. We approximated the source amplitude functions of eruption tremor and explosion events by quadrilateral and triangular shapes. For eruption tremor, the durations of the source amplitude functions increased with decreasing slope of the initial phase, that is between onset and maximum amplitude. For explosion events, both the maximum and cumulative amplitudes of the source amplitude functions increased with increasing slope of the initial phase, but the overall durations clustered around a typical value. Moreover, the initial phase durations of eruption tremor were longer than those of explosion events. Based on eruption models proposed by previous studies, Vulcanian and sub-Plinian eruptions have been thought to be triggered by accumulation of magma at a shallow part in a conduit and mixing of cool mushy magma with hot fresh magma in a reservoir, respectively. The above differences between the source amplitude functions of eruption tremor and explosion events can be explained by the distinct eruption triggering processes of sub-Plinian and Vulcanian eruptions. Our results suggest that source amplitude functions are useful for investigating eruption processes and estimating eruption sizes and durations for seismic eruption monitoring.

New evidence of Holocene pyroclastic density currents at Galeras volcano, Colombia

The Galeras Volcanic Complex (GVC) is a composite volcano located between the Central and Western cordilleras in southern Colombia. It is considered one of the most active volcanoes in Colombia, with typically Vulcanian eruptions. The complex has been divided into different stages based on stratigraphic and geochronological differences. The most recent one is called the Galeras stage and is considered to have started ca. 4500 years. This study presents the results of new stratigraphic and geochronological analysis of the deposits that outcrop along the El Barranco river valley, NW of the GVC, which allowed us to identify newer and older deposits within the most recent stage (i.e., the Galeras stage). The deposits were analyzed via distribution, stratigraphy, geochronology, componentry and granulometry, as well as vesicularity and microtexture of the pumice fragments. We identified three previously unreported deposits, which results in a total of eight pyroclastic density current (PDC) deposits in the valley, which discordantly overlie andesitic lava flows of the previous stage (i.e., Genoy stage; 150–40 ka) of the GVC. They were named from Unit U1 to Unit U8, with units U1, U2 and U8 representing events not previously recorded in the eruptive history of Galeras volcano. Units U1 (8303 ± 97 and 8284.5 ± 90.5 cal BP), U2 (7667 ± 78 cal BP), U3 (5082 ± 198, 5096.5 ± 226.5, 5801.5 ± 481.5 and 5966 ± 235 cal BP) and U6 (2215.5 ± 101.5 and 2033 ± 91 cal BP) are interpreted as formed by PDCs generated by pumice-rich flows linked to the collapse of Subplinian eruptive columns, while units U5 (3201.5 ± 129 cal BP), U7 (1066.5 ± 109.5 and 944.5 ± 115.4 cal BP) and U8 (390.5 ± 80.5 cal BP) are interpreted as formed by PDCs generated by block and ash flows linked to the destruction of intracrater domes and the collapse of Vulcanian eruptive columns. U4 (4625 ± 181 cal BP) is associated with a “blast” type event, caused by the overpressure exerted by the gases of a viscous magma. The above mentioned record allows us to propose that the eruptive history of Galeras volcano began at the beginning of the Holocene (∼8300 years BP), and not 4500 BP, as previously established, which has hazard implications. Unit U8 in the town of La Florida as well as the nearby recorded deposits, indicates a probable PDC hazard for this population.

Probabilistic, scenario-based hazard assessment for pyroclastic density currents at Tungurahua volcano, Ecuador

We assess the volcanic hazard posed by pyroclastic density currents (PDCs) at Tungurahua volcano, Ecuador, using a probabilistic approach based on the analysis of calibrated numerical simulations. We address the expected variability of explosive eruptions at Tungurahua volcano by adopting a scenario-based strategy, where we consider three cases: violent Strombolian to Vulcanian eruption (VEI 2), sub-Plinian eruption (VEI 3), and sub-Plinian to Plinian eruption (VEI 4–5). PDCs are modeled using the branching energy cone model and the branching box model, considering reproducible calibration procedures based on the geological record of Tungurahua volcano. The use of different calibration procedures and reference PDC deposits allows us to define uncertainty ranges for the inundation probability of each scenario. Numerical results indicate that PDCs at Tungurahua volcano propagate preferentially toward W and NW, where a series of catchment ravines can be recognized. Two additional valleys of channelization are observed in the N and NE flanks of the volcano, which may affect the city of Baños. The mean inundation probability calculated for Baños is small (6 ± 3%) for PDCs similar to those emplaced during recent VEI 2 eruptions (July 2006, February 2008, May 2010, July 2013, February 2014, and February 2016), and on the order of 13 ± 4% for a PDC similar to that produced during the sub-Plinian phase of the August 2006 eruption (VEI 3). The highest intensity scenario (VEI 4–5), for which we present and implement a novel calibration procedure based on a few control points, produces inundation areas that nearly always include inhabited centers such as Baños, Puela, and Cotaló, among others. This calibration method is well suited for eruptive scenarios that lack detailed field information, and could be replicated for poorly known active volcanoes around the world.

The color systematics of volcanic ashfall samples in estimating eruption sequences: a case study of the 2017–2018 eruption at Shinmoe-dake, Kirishima volcano, Southwest Japan

The color of pyroclasts is fundamental, because it reflects various magma properties and eruption processes, including particle morphology, chemistry, and petrological characteristics. However, deriving the componentry ratio (CR) of pyroclasts for ongoing eruption monitoring remains challenging due to the lack of a robust quantitative standard. The derivation of the CR, as well as other petrological analyses, is too laborious and time-consuming to introduce as a sustainable monitoring method. To address this, we employed spectroscopic colorimetry to rapidly and quantitatively describe eruptive product colors, enabling CR derivation based on clear, objective standards for ash particle classification. Through color spectroscopy of bulk and sieved ash samples, we analyzed the major size fraction for time-series samples during the waxing stage of the 2017–2018 Shinmoe-dake eruption in Kirishima volcano, Southwest Japan. Our findings reveal that the color changes in bulk ash systematically changed with the evolution of componentry. This temporal color change was due to an increase in the amount of vesicular particles with clear glass against dark angular lava particles, as well as a grain size change, which we interpret as an indication of a transition from phreatic/phreatomagmatic to magmatic eruption. Subsequently, the color of the ash changed when the amount of different lava particles increased gradually, coinciding with a shift toward a more dominant effusion of lava. As the lava effusion continued, a slight reddening of the ash, indiscernible to the naked eye, was clearly detected by the spectrometer before the onset of intermittent Vulcanian eruptions. We interpreted this reddening as oxidation resulting from decreased ascent speed and caprock formation, which accumulates overpressure for Vulcanian explosions. These results highlight the potential of rapid, objective color value and componentry derivation for sustainable quasi-real-time monitoring of ongoing eruption materials.Graphical

Newly identified small vulcanian eruptions during the caldera-forming stage of Towada Volcano, Northeast Japan

Understanding the detailed eruptive history and conditions during the buildup to catastrophic caldera-forming eruptions is essential for understanding the evolution of caldera volcanoes and predicting eruptive hazards. Towada Volcano is an active caldera volcano in the northern part of the Northeast Japan Arc. At least two catastrophic caldera-forming eruptions (eruptive episodes N and L) occurred during its caldera-forming stage (61–15.7 ka). A detailed geological survey identified three small vulcanian tephra layers that were erupted during the caldera-forming stage. The tephra layers are blue–gray ash fall deposits that consist mainly of fresh blocky dacite–rhyolite fragments. In ascending order, these deposits are assigned to eruptive episodes O′, N′, and M′. Each eruption had a volume of <0.11 km3, which is much smaller than that of other eruptions during the caldera-forming stage. The eruptive ages of episodes O′, N′, and M′ are estimated to be ca. 40, 23.1, and 17.2 ka, respectively, based on new 14C ages and paleosol thicknesses data. At least three cycles of a medium- to large-scale (3–20 km3) explosive eruption, preceded by a small vulcanian eruption (<0.11 km3), occurred during the latter half of the caldera-forming stage. The small vulcanian eruptions (episodes O′, N′, and M′) occurred after relatively long periods of quiescence (4000–14,000 years) and were followed by medium- to large-scale explosive eruptions (episodes N, M, and L) 1500–4000 years later. This cyclic activity probably reflects changes in overpressure in the magma reservoir. As the overpressure increased, episodes O′, N′, and M′ probably occurred when the overpressure had not increased sufficiently to trigger a medium- to large-scale explosive eruption, and episodes N, M, and L occurred at higher overpressures. These cycles characterize the fermentation phase of Towada Volcano, and terminated with the formation of the Towada Caldera during episode L at 15.7 ka. The magma composition and eruption frequency changed abruptly after episode L, and a small stratovolcano was formed through intermittent eruptions of mafic magma. This change suggests that the caldera collapse during episode L changed the entire shallow magmatic system that existed during the fermentation phase, and the system shifted to a recovery (post-caldera) phase. Based on eruptive volumes and frequencies, the present Towada Volcano has not yet reached the conditions that existed during the late caldera-forming stage and is therefore unlikely to produce a catastrophic caldera-forming eruption in the near future.

Magma intrusion process during pre-magmatic period (2010−2013) of Sinabung volcano as revealed by seismicity of volcano-tectonic and hybrid earthquakes

Sinabung volcano in Indonesia reactivated, with a series of phreatic eruptions, in August and September 2010. These eruptions were followed by various types of magmatic eruption at the summit, including lava dome growth, pyroclastic density currents, lava flow, and vulcanian eruptions, starting in December 2013. Prior to the magmatic eruptions, Sinabung exhibited two sequences of phreatic eruption and an increase in the seismicity of volcanic earthquakes. This study clarifies the progress of magma intrusion during the pre-magmatic period based on the hypocenter distribution, waveform similarity, seismic moment, rupture length, and stress drop for volcano-tectonic (VT) and hybrid earthquakes. It was found that the hypocenters of VT earthquakes are distributed north and northwest of the summit at a depth range of 1–10 km below sea level. Deep seismicity (4–10 km) alternated with shallow seismicity (2–4 km). The last shallow seismicity, from July to mid-December 2013, was different from previous seismicity in terms of higher intensity, as determined from the seismic moment, a temporary increase in the rupture length, and the migration of hypocenters towards the summit. The seismicity of VT earthquakes was altered by a swarm of hybrid earthquakes in mid-December. The hybrid earthquakes were smaller, in terms of seismic moment (mostly <3 × 1010 Nm), than the VT earthquakes. Their hypocenters were concentrated in the shallowest depth range (−0.5 to 1.5 km below sea level) directly below the summit, their source process was repeatable (they could be grouped into six earthquake families), their dominant peak was in a low-frequency range (2.5–4.5 Hz), and they had a relatively low stress drop (<0.15 MPa). This suggests that the swarm of hybrid earthquakes was induced by a frequent repeated fracture of fluid-filled cracks due to the intrusion of magma up to directly below the summit. The transition of the earthquake family and the change in its source parameters, namely an increase in the stress drop prior to the appearance of a lava dome and then a slight decrease, may reflect a gradual change in internal pressure in the hypocentral zone through the magma intrusion process.

The last 1100 years of activity of La Fossa caldera, Vulcano Island (Italy): new insights into stratigraphy, chronology, and landscape evolution

A detailed study of past eruptive activity is crucial to understanding volcanic systems and associated hazards. We present a meticulous stratigraphic analysis, a comprehensive chronological reconstruction, thorough tephra mapping, and a detailed analysis of the interplay between primary and secondary volcanic processes of the post-900 AD activity of La Fossa caldera, including the two main systems of La Fossa volcano and Vulcanello cones (Vulcano Island, Italy). Our analyses demonstrate how the recent volcanic activity of La Fossa caldera is primarily characterized by effusive and Strombolian activity and Vulcanian eruptions, combined with sporadic sub-Plinian events and both impulsive and long-lasting phreatic explosions, all of which have the capacity to severely impact the entire northern sector of Vulcano island. We document a total of 30 eruptions, 25 from the La Fossa volcano and 5 from Vulcanello cones, consisting of ash to lapilli deposits and fields of ballistic bombs and blocks. Volcanic activity alternated with significant erosional phases and volcaniclastic re-sedimentation. Large-scale secondary erosion processes occur in response to the widespread deposition of fine-grained ash blankets, both onto the active cone of La Fossa and the watersheds conveying their waters into the La Fossa caldera. The continuous increase in ground height above sea level, particularly in the western sector of the caldera depression where key infrastructure is situated, is primarily attributed to long-term alluvial processes. We demonstrate how a specific methodological approach is key to the characterization and hazard assessment of low-to-high intensity volcanic activity, where tephra is emitted over long time periods and is intercalated with phases of erosion and re-sedimentation.

Vulcanian eruptions at Sakurajima Volcano: geophysical data, numerical modelling, and petrological evidence

We investigated the magma conduit system beneath three active craters that have repeatedly generated Vulcanian eruptions at Sakurajima Volcano, Japan, by analysing seismic, infrasound, and tilt data. The hypocentres of the earthquakes associated with the Vulcanian eruptions are distributed separately at shallow depths of approximately 0.5 km beneath the craters. Infrasound indicated that the swelling of the crater floor starts approximately 0.2 s after the occurrence of earthquakes and that the eruption starts about 0.3 to 0.5 s later. During an eruption, tilt vectors at a station approximately 2.7 km far from the active craters indicated a deflationary trend directed toward one of the three active craters. A 1-D conduit flow simulation indicated pressure build-up at a depth of 0.4–1.0 km beneath the crater, consistent with previously reported pressure sources during eruptions detected by geodesy. Volcanic lapilli emitted from the three craters have the same chemical composition; hence, Vulcanian eruptions of all three studied craters originate from a common magma source.

Tephra segregation profiles based on disdrometer observations and tephra dispersal modeling: Vulcanian eruptions of Sakurajima volcano, Japan

The profile of tephra concentration along a volcanic plume (i.e., the tephra segregation profile) is an important source parameter for the simulation of tephra transport and deposition and thus for the tephra sedimentation load. The most commonly-used approach is to treat an eruption as a single event (i.e., with a time-averaged mass eruption rate; MER). In this case, it is common to use pre-determined profiles that feature most of the tephra segregate at the top of the plume. However, case studies based on observations have revealed that large concentration maxima also appear at the lower part of the plume. To investigate this discrepancy, the impact of plume height on the temporal variations in the MER is examined. To this end, we use the tephra transport and dispersion model Tephra4D with MER estimates obtained from geophysical monitoring and maximum plume height observations to calculate the spatial distribution of the tephra deposit load for 39 eruptive events that consisted of explosions and quasi-steady particle emission from the Sakurajima volcano, Japan. A comparison of the model results with observations from a disdrometer network revealed that for both kinds of activity, maxima in tephra segregation can occur at heights below the reported plume height. The tephra segregation profiles of Vulcanian eruptions at Sakurajima volcano are consistent with most of the modeling studies giving profiles that feature most of the tephra segregating at the top of the plume if the temporal variation of the MER is taken into consideration to properly represent the total series of eruptive events in a sequence. This highlights that even though the activity at Sakurajima volcano is commonly characterized simply as Vulcanian eruptions, in addition to the primary plume developed due to the initial instantaneous release caused by the explosion, the subsequent continuous plume that can accompany the eruption plays an important role in particle emission. Calculations could not reproduce the simultaneous deposition of particles with a wide range of settling velocities in observations, suggesting the importance of volcanic ash fingers caused by gravitational instability in tephra transport simulations.Graphical

Dispersal and grain size characteristics of the May 14, 2018 Shinmoedake eruption deposit, Kirishima Volcano, Japan, based on post-eruption field survey and meteorological datasets

This study describes the dispersal and grain size characteristics of the May 14, 2018 Shinmoedake eruption deposits of Kirishima Volcano in southern Kyushu, southwestern Japan. We discuss the eruption sequence, including the temporal variations in the behavior of the plume, by combining field and meteorological datasets. Following a magmatic activity in 2011 characterized by a substantial change in the eruption style (from subplinian eruptions to lava effusion) and subsequent vulcanian explosions, the Shinmoedake crater experienced intermittent eruptions in 2018. The May 14, 2018 eruption began at 14:44 with a vulcanian eruption, with the eruption plume rising 4500 m above the crater rim. Thereafter, it transitioned to an ash eruption; the plume height decreased gradually until the eruption ceased at 16:10. The tephra fall deposits were distributed more than 27 km to the southeast of the source crater; the mass of the tephra fall deposit was approximately 2.1 × 107 kg, calculated based on an isomass map. The deposit incidence differed between the east and west sides of the major dispersal axis. The deposits found east of the main dispersal axis were primarily composed of coarse to medium sand-sized particles with no fine fraction (fine sand to silt in size). In contrast, the deposits west of the axis were finer-grained than those east of the axis. We analyzed photographs of the eruption plume, along with the regional meteorological data and the dispersal and grain-size characteristics of the deposits, and reached the following conclusion: during the May 14, 2018 eruption, the wind directions above the Shinmoedake crater fluctuated across altitudes. The westerly winds dispersed the eruption plume that rose to a higher altitude, containing coarser tephra associated with the initial vulcanian eruption, further to the east rather than along the main axis. In contrast, a lower-altitude ash eruption plume that was rich in fine materials was dispersed westward rather than along the main axis, which was influenced by northerly winds. The findings of this study can support the analysis of similar volcanic events.Graphical

Insights into the sticking probability of volcanic ash particles from laboratory experiments

Although the characterization of the sticking and aggregation probability is essential to the description of volcanic ash dispersal and sedimentation, there is still no general model describing the sticking probability of volcanic ash. Experiments of dry particle–plate collisions in an enclosed box were carried out to characterize quantitatively the sticking efficiency of volcanic particles and silica beads in a limit case scenario where the mass of one of the particles is much greater than the others. Silica beads and volcanic particles from a Sakurajima Vulcanian eruption were filmed impacting a glass plate with a High-Speed Camera. The sticking probability is calculated from an equation depending on the particle diameter, impact velocity, and two experimental parameters (a, q). Particle size was found to dominantly control the sticking probability, with small particles more likely sticking on the glass plate than large particles. These experiments represent a significant step forward in the quantification of the sticking efficiency of fine volcanic ash (< 63 μm).

What controls the formation of vulcanian bombs? A case study from the 1 February 2014 eruption of Tungurahua (Ecuador)

Vulcanian eruptions are very common at many volcanoes around the world that erupt intermediate to silicic magmas. This type of eruption generates a wide variability of bombs and blocks preserving information onto the conduit processes that strongly control the dynamics of these events. After 84 years of repose, a new cycle of eruptive activity of Tungurahua volcano (Ecuador) started in October 1999, consisting of recurrent low-to-moderate explosive phases, which included sporadic stronger Strombolian to sub-Plinian pulses, as well as Vulcanian outbursts. The 1 February 2014 eruption was one of the most important Vulcanian events and was characterized by highly energetic explosions resulting from a plug conduit failure that generated a ∼ 9 km-high eruptive column and fallback pyroclastic density currents. Four different types of blocks and bombs were found in the deposits of the pyroclastic density currents: dense fragments (DB), breadcrust bombs (BCB), cauliflower bombs (CFB) and foliated, banded bombs (FB). All the different types of bombs have homogenous andesitic bulk-rock compositions but different, highly evolved matrix glass compositions ranging from rhyolitic for BCB to dacitic for CFB and DB, suggesting the occurrence of contrasting shallow crystallization processes within the conduit. The wide variability of the bombs in terms of patterns of surface cracks, external morphologies, internal density and vesicularity gradients records different conditions of formation from a same magma composition.The combination of morphological measurements, compositional data and textural analysis allow us (i) to infer the pre-eruption physical and rheological state of the conduit plug, (ii) to reconstruct the short-scale lateral and vertical gradients characterizing the plug, and (iii) to clarify the specific mechanism controlling the formation conditions of CFB, still poorly defined in the literature.

Pre‐Eruptive Outgassing and Pressurization, and Post‐Fragmentation Bubble Nucleation, Recorded by Vesicles in Breadcrust Bombs From Vulcanian Activity at Guagua Pichincha Volcano, Ecuador

AbstractBreadcrust bombs formed during Vulcanian eruptions are assumed to originate from the shallow plug or dome. Their rim to core texture reflects the competition between cooling and degassing timescales, which results in a dense crust with isolated vesicles contrasting with a highly vesicular vesicle network in the interior. Due to relatively fast quenching, the crust can shed light on pre‐ and syn‐eruptive conditions prior to or during fragmentation, whereas the interior allows us to explore post‐fragmentation vesiculation. Investigation of pre‐ to post‐fragmentation processes in breadcrust bombs from the 1999 Vulcanian activity at Guagua Pichincha, Ecuador, via 2D and 3D textural analysis reveals a complex vesiculation history, with multiple, spatially localized nucleation and growth events. Large vesicles (Type 1), present in low number density in the crust, are interpreted as pre‐eruptive bubbles formed by outgassing and collapse of a permeable bubble network during ascent or stalling in the plug. Haloes of small, syn‐fragmentation vesicles (Type 2), distributed about large vesicles, are formed by pressurization and enrichment of volatiles in these haloes. The nature of the pressurization process in the plug is discussed in light of seismicity and ground deformation signals, and previous textural and chemical studies. A third population (Type 3) of post‐fragmentation small vesicles appears in the interior of the bomb, and growth and coalescence of Type 2 and 3 vesicles causes the transition from isolated to interconnected bubble network in the interior. We model the evolution of viscosity, bubble growth rate, diffusion timescales, bubble radius and porosity during fragmentation and cooling. These models reveal that thermal quenching dominates in the crust whereas the interior undergoes a viscosity quench caused by degassing, and that the transition from crust to interior corresponds to the onset of percolation and development of permeability in the bubble network.

Vulcanian eruption processes inferred from volcanic glow analysis at Sakurajima volcano, Japan

Volcanic glow, which is often observed at active craters during nighttime, contains information on high-temperature zones deep in the crater that are hidden from view. We analyzed 90 eruption videos capturing volcanic glow before the onset of Vulcanian eruptions at the Showa crater of Sakurajima volcano (Japan), and found that they show clear temporal changes. The red (R)-value and the green-to-red ratio (GR ratio) are used to track glow intensity and temperature change, respectively. We found two types of temporal change in glow intensity: a short-term change approximately 1 s before the onset of explosions, and a long-term change 3 to 330 s before an explosion. Short-term changes were observed for 28 of 90 eruptions, and coincided with an increase in the GR ratio and with a modest precursory infrasound pressure increase. Since an increasing GR ratio means a temperature increase, we consider that the change was associated with the opening of tensile cracks on the crater floor before an explosion. However, the other 62 eruptions did not show short-term changes, suggesting that the opening of cracks is sporadic. Long-term changes were observed for 73 of 90 eruptions, the GR ratio did not exhibit a significant increase. We suggest that long-term changes are caused by gradual gas leakage, which precedes initiation of Vulcanian eruptions. We hypothesize that the long-term gas leakage causes gradual decompression of the conduit, which induces supersaturation of volatile-rich magma beneath the lava plug causing explosive expansion, thereby triggering an eruption.

A data driven approach to mineral chemistry unveils magmatic processes associated with long-lasting, low-intensity volcanic activity

The most frequent volcanic eruptions are of low-intensity and small magnitude. They produce abundant ash-sized (< 2 mm) clasts, which are too small to establish quantitative links between magmatic processes and eruptive dynamics using classic approaches. This inhibits our ability to study the past behaviour of frequently erupting volcanoes, essential to predict their future activity and mitigate their impact. The Palizzi unit (10–13th century, Vulcano, Italy) includes a prototype sequence of ash deposits resulting from prolonged Vulcanian eruptions punctuated by those of two larger sub-Plinian events. We apply Hierarchical Clustering to chemical analyses of clinopyroxene collected along the stratigraphy to decipher magma dynamics during this eruptive period. We identify periods of magma input and we link deep magmatic processes to eruptive dynamics, also showing that our approach can be used to connect magma and eruptive dynamics in any volcanic sequence. This is essential to track the processes occurring during frequent eruptions and to identify the build-up to larger explosive events.

Seismic tremor location of 10 large paroxysmal eruptions of Tungurahua volcano, Ecuador

SUMMARYThe most recent eruptive period of Tungurahua volcano lasted 17 yr (1999–2016), generating strong eruptive phases with the release of large amounts of seismic and acoustic energies. We have selected 10 large eruptions and located their seismic sources. The location method involves applying a new method for correcting for site effects in the seismic records and identifying, in the frequency domain, the signals from the ground-coupled airwaves and pyroclastic flows. The locations are computed by minimizing the differences of the source energy rates in the [0.4, 2.5] Hz range. The results suggest that the first three eruptions have locations mainly concentrated at the [−1, 2] km depth range (asl), and the locations of the fourth eruption are mainly coincident with depths of a possible reservoir estimated from petrological studies. This fact strongly suggests that this eruption was responsible of changes in the dynamics of Tungurahua, making the volcano prone to produce Vulcanian eruptions.

Thermobarometry, vesicularity and microlite populations in pyroclastic density current deposits from the last 4500 years of activity at Galeras, Colombia

Vulcanian eruptions during the 19th and 20th Centuries at the Galeras volcano have produced small eruption columns and bomb fields around the active vent. Four pyroclastic density current (PDC) deposits from the early history of the currently active Galeras cone which have ages between 4500 and 2300 years before present have a much wider distribution. The first two (ML 1 and ML 2) were deposited from hot, dilute, pumice-rich pyroclastic density currents. Amphibole thermobarometry indicates that the magma for these eruptions tapped reservoirs between 5 and 28 km (ML 1) and 33–41 km (ML 2) below the crater. Debse, degassed magma formed a conduit plug which persisted for up to several weeks before the eruptions based on the thickness of breakdown rims around the amphibole. Amphibole-bearing pumice came from a reservoir 5–26 km below the crater and since it lacks breakdown rims it must have reached surface in <8 days. The vesicularity of the pumice indicates two phases of bubble nucleation and growth which coincides with growth of plagioclase and orthopyroxene microlites. The upper two layers (ML 3 and ML 4) are block and ash flows. ML3 shows significant hydrothermal alteration of the dense lithic clasts and has only minor vesicular scoria. The magma for the ML 3 eruption came from a reservoir 17–29 km below the crater and the presence of 20 μm thick breakdown rims on the amphibole implies that the degassing plug was stable for ∼20 days, though the degree of hydrothermal alteration of the clasts appears to require a significantly longer time. The uppermost layer (ML4) is relatively unaltered with a large proportion of pumice. Dense lithics indicate plug formation from magma that had its source at ∼17 km below the crater and remained stable during degassing for ∼20 days. The absence of amphibole from the pumice in ML 4 implies a distinct magma composition was emplaced after the plug.

Hazardous Changes in Soil CO2 Emissions at Vulcano, Italy, in 2021

AbstractThe La Fossa volcano on the Island of Vulcano, Italy, showed signs of more energetic fumarolic–solfataric activity during 2021. Several increases in volcanic gas emissions and seismicity, namely, “crisis,” punctuated the passive degassing at Vulcano that had ensued after the last 1888–1890 vulcanian eruption. Most of the gases (i.e., up to 90%) were emitted at the crater cone while the diffuse degassing of CO2 at Vulcano Porto accounted for more than 10% of the volcanic emissions. Two anomalous degassing zones at the base of the volcanic cone (i.e., Palizzi and Faraglione) showed notable changes in the gas output during the volcanic crisis. In these zones, increases of soil CO2 flux (φCO2) had several practical implications other than of volcanological interest, owing to the risk related to people's exposure to volcanic gas emissions. The results of this study reveal variations of the average φCO2 from 74 g m−2 d−1 during September 2021 to 370 g m−2 d−1 in November 2021, which were 27% and 538% higher than the statistical background since 1988 (φCO2 ≈ 58 g m−2 d−1), respectively. These observations helped in volcanic surveillance at Vulcano. The soil CO2 partitioning determined using both φCO2 and carbon isotope measurements, helped track changes in the volcanic CO2 output from 9.97 · 104 to 101.15 · 104 kg d−1. Estimates for volcanic CO2 suggest that the instability of a magmatic body caused a transition from background fumarolic–solfataric activity toward an unrest event after September 2021.

Behavior of Sn-3.0Ag-0.5Cu solder/Cu fluxless soldering via Sn steaming under formic acid atmosphere

Formic acid (FA) atmosphere is promising to achieve chemical reduction at the oxidized surfaces of the solder and matrix material during soldering, replacing chemical flux being challenging as the continuous decrease in the pitch of solder bump in a three-dimensional integrated circuit. Although a previous study observed the effect of Sn steaming on the wettability of the solder during FA soldering, it has not been understood yet by adequate studies. This study demonstrates the fluxless soldering behavior of Sn-3.0Ag-0.5Cu (SAC) solder/Cu under FA atmosphere. Resembling vulcanian eruption Sn steaming occurred with a solid-state Sn and promoted Sn–Cu intermetallic compound (IMC) formations obviously at temperature (210 ° C) lower than the melting temperature of Sn (231 ° C). As indicated by in-situ observations, the IMC formations led to the final spreading area of solder under FA soldering was larger than that under the soldering with rosin mildly activated (RMA) flux. The wettability of SAC solder under FA atmosphere could be improved with continuous heating.

Global shockwaves of the Hunga Tonga-Hunga Ha’apai volcano eruption measured at ground stations

The eruption of the Tonga volcano created globally propagating spherical shockwaves in the atmosphere. Analyses are done to data from two southern U.S. stations of the author sampling at 3-21s intervals and 189 weather stations at 1-5min intervals. The shockwaves arrived from two routes in the atmosphere: the shortest spherical arc and the longer spherical arc through the antipole. In most stations, signals up to the 6th path of shockwaves were recorded as the waves traveled around the globe multiple times. The speed of shockwaves is estimated to be 309.5± 2.9m/s, consistent with the speed of sound at the top of the troposphere where a waveguide exists. Discussion is made on the post-shockwave ringing of 4-8min as higher amplitude oscillations above the level of pre-shockwaves background noise. A theoretical wave dispersion is derived which verifies that the spherical shockwave's phase speed is the same as the speed of sound.