Akihiko Yokoo

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.

At active volcanoes, the surface temperature and its spatial distribution can indicate changes in the underlying magmatic and hydrothermal system. Surface temperature monitoring has been widely performed using thermal remote sensing with thermal infrared (TIR) cameras. One of the drawbacks of this method is that the unknown viewing orientation of TIR images inhibits quantitative evaluations of the spatial extent and distributions of thermal anomalies. Therefore, many studies have performed 3D temperature-field reconstructions by processing TIR images photogrammetrically. However, these studies have not included the correction of TIR wave atmospheric attenuation and viewing-angle variation in emissivity in the reconstruction procedure, which can result in significant temperature misestimation. We propose a simple method that incorporates the correction into the reconstruction process, which can improve the estimation of surface temperature, especially in rugged terrains. We demonstrate our method for the active crater of Aso volcano in Japan. We create digital elevation models by applying the Structure from Motion–Multi-view Stereo algorithm to aerial visible images taken at the same time as the TIR images and then project the TIR images onto the DEM with the direct georeferencing method. The correction is carried out using the geometric parameters acquired in the process. We create 1-m resolution orthorectified thermal images of the crater on two separate dates (18 August 2020 and 16 March 2022) and calculate the heat discharge rates from steaming grounds and a volcanic lake. We find that the heat discharge rate from the fumarolic field on the south crater wall showed a sevenfold increase after the phreatic explosions that occurred on 14 and 20 October 2021. Graphical Abstract

Monitoring the depth of the magma surface at open-vent volcanoes can be a practical tool to infer temporal variations in the magma supply during an eruption. We focus on the magmatic eruption of Aso volcano in 2014–2015 to estimate the temporal change in the depth of the magma surface, and show that this needs to be coupled with an understanding of the shallow conduit geometry if it is to be done in a representative manner. The eruption lasted 5 months from November 2014 and ending with a crater floor collapse in May 2015. During the eruption, we recorded seismo-acoustic waveforms related to frequent Strombolian explosions. The infrasound signals show several distinct peak frequencies derived from acoustic resonance inside the vent. We estimate the depth of the magma surface using the time delay of seismo-acoustic signals and the peak frequency of infrasound signals. In addition, the temporal variation in the shape of the conduit is constrained by the overtone frequency of the acoustic resonance. From the beginning of the eruption to early-January 2015, the magma surface was located at a depth of $\sim $ ∼ 200 m, and the conduit was a cylindrical pipe. Later, between January and February, the magma surface rose to $\sim $ ∼ 120 m, and the shape of the conduit changed to a conical frustum flaring inside. This finding indicates that the magma was injected into the shallow conduit and that it heated and weakened the conduit wall near the magma surface. Before the cessation of the magmatic eruption, the magma surface dropped by approximately 70 m. This magma drainage and, primarily, the instability of the conduit shape caused the crater floor to collapse. We show the possibility of tracking and assessing the depth of the magma surface and the shallow conduit geometry even with limited seismo-acoustic observations.

Open-vent volcanoes provide opportunities to perform various methods of observation that can be used to study shallow plumbing systems. The depth of the magma–air interface in the shallow portion of the conduit can be used as an indicator of the volcanic activity of open-vent volcanoes. Although there are many methods used to estimate the depth, most of them cannot constrain the depth to a narrow range due to other unknown parameters. To constrain the depth more accurately, we combine two methods commonly used for estimating the depth of the magma–air interface. They consider the acoustic resonant frequency and the time delay of arrivals between the seismic and infrasound signals of explosions. Both methods have the same unknown parameters: the depth of the magma–air interface and the sound velocity inside the vent. Therefore, these unknowns are constrained so that both the observed resonant frequency and time delay can be explained simultaneously. We use seismo-acoustic data of Strombolian explosions recorded in the vicinity of Aso volcano, Japan, in 2015. The estimated depths and the sound velocities are 40–200 m and 300–680 m/s, respectively. The depth range is narrower than that of a previous study using only the time delay of arrivals. However, only a small amount of the observed data can be used for the estimation, as the rest of the data cannot provide realistic depths or sound velocities. In particular, a wide distribution of the observed time delay data cannot be explained by our simple assumptions. By considering a more complicated environment of explosions, such as source positions of explosions distributed across the whole surface of a lava pond in the conduit, most of the observed data can be used for estimation. This suggests that the factor controlling the observed time delay is not as simple as generally expected. Graphic abstract

Sound waves generated by erupting volcanoes can be used to infer important source dynamics, yet acoustic source‐time functions may be distorted during propagation, even at local recording distances ( 15 km). The resulting uncertainty in source estimates can be reduced by improving constraints on propagation effects. We aim to quantify potential distortions caused by wave steepening during nonlinear propagation, with the aim of improving the accuracy of volcano‐acoustic source predictions. We hypothesize that wave steepening causes spectral energy transfer away from the dominant source frequency. To test this, we apply a previously developed single‐point, frequency domain, quadspectral density‐based nonlinearity indicator to 30 acoustic signals from Vulcanian explosion events at Sakurajima Volcano, Japan, in an 8‐day data set collected by five infrasound stations in 2013 with 2.3‐ to 6.2‐km range. We model these results with a 2‐D axisymmetric finite‐difference method that includes rigid topography, wind, and nonlinear propagation. Simulation results with flat ground indicate that wave steepening causes up to 2 dB (1% of source level) of cumulative upward spectral energy transfer for Sakurajima amplitudes. Correction for nonlinear propagation may therefore provide a valuable second‐order improvement in accuracy for source parameter estimates. However, simulations with wind and topography introduce variations in the indicator spectra on order of a few decibels. Nonrandom phase relationships generated during propagation or at the source may be misinterpreted as nonlinear spectral energy transfer. The nonlinearity indicator is therefore best suited to small source‐receiver distances (e.g., 2 km) and volcanoes with simple sources (e.g., gas‐rich strombolian explosions) and topography.

We estimated the vent radius within Showa Crater of Sakurajima volcano from ejection velocity and flow rate of gas‐particle mixtures. The ejection velocity was calculated from video recordings, and the flow rate and volume from infrasound data. Based on the assumption that the vent shape does not change during an explosion, the vent radius was estimated from 201 impulsive Vulcanian eruptions at Showa Crater, yielding values of 6.4–42.3 m (median 23.8 m), which is comparable with the width of fresh lava capping the vent, as photographed from a helicopter. Long‐term changes in vent radius (i.e., over several months) show a relationship with magma accumulation within a reservoir 2–5 km beneath the crater. If the top of the conduit is assumed to be cylindrical, then the vertical extent of the gas‐rich zone is estimated to be 120 m, which may reflect the depth of gas accumulation and buildup of significant overpressure.

Transient volcanic plumes, having similar eruption duration and rise timescales, characterize many unsteady Strombolian to Vulcanian eruptions. Despite being more common, such plumes are less studied than their steady state counterpart from stronger eruptions. Here we investigate the initial dynamics of transient volcanic plumes using high‐speed (visible light and thermal) and high‐resolution (visible light) videos from Strombolian to Vulcanian eruptions of Stromboli (Italy), Fuego (Guatemala), and Sakurajima (Japan) volcanoes. Physical parameterization of the plumes has been performed by defining their front velocity, velocity field, volume, and apparent surface temperature. We also characterized the ejection of the gas‐pyroclast mixture at the vent, in terms of number, location, duration, and frequency of individual ejection pulses and of time‐resolved mass eruption rate of the ejecta's ash fraction. Front velocity evolves along two distinct trends related to the initial gas‐thrust phase and later buoyant phase. Plumes' velocity field, obtained via optical flow analysis, highlights different features, including initial jets and the formation and/or merging of ring vortexes at different scales. Plume volume increases over time following a power law trend common to all volcanoes and affected by discharge history at the vent. Time‐resolved ash eruption rates range between 10² and 10⁷ kg/s and may vary up to 2 orders of magnitude within the first seconds of eruption. Our results help detailing how the number, location, angle, duration, velocity, and time interval between ejection pulses at the vents crucially control the initial (first tens of second), and possibly later, evolution of transient volcanic plumes.