Abstract:Temporal variation of ground surface displacement rates at volcanoes may represent changes to the volumetric magma supply rate over time. If crustal rock responds elastically to a magmatic deformation source, the volume change of the source has a simple one-to-one relation to the surface deformation. For viscoelastic crustal rock, however, the correlation is not always straightforward, particularly at volcanoes where crustal viscosity is lowered by a high geothermal gradient. For a continuous inflation of the deformation source, the viscoelastic surface displacement depends on the balance between the uplift due to magma supply and subsidence due to viscoelastic relaxation. In this study, using a 3D finite element model composed of an elastic layer underlain by a linear Maxwell viscoelastic layer with a spatially uniform viscosity, we estimate the temporal variation of volumetric magma supply rate into the upper crust beneath Aira caldera in southern Kyushu, Japan, after the 1914 eruption of Sakurajima volcano. It is found that the geodetic data, including levelling and GNSS displacement fields, require a sill-like magmatic deformation source at a depth of 11 km to inflate at the following rates for six different periods: (I) ~6.9–9.4 × 106 m3/yr in 1914–1934, (II) ~9.1–16.7 × 106 m3/yr in 1934–1960, (III) ~1.6–3.8 × 106 m3/yr in 1960–1968, (IV) ~8.1–11.0 × 106 m3/yr in 1968–1976, (V) ~ −1.0–2.2 × 106 m3/yr in 1976–1997, and (VI) ~5.8–9.4 × 106 m3/yr in 1997–2007. The constrained viscosity (η c ) of the viscoelastic crust ranges from ~5 × 1018 Pa s to ~1020 Pa s or more. The lowest and highest values produce the best fit to the data in the earliest period and in period V showing surface subsidence, respectively, but for the other periods any viscosity is acceptable unless it is lower than ~1018 Pa s. The lowest viscosity providing the best fit in the earliest period indicates that viscoelastic displacement in response to the 1914 eruption plays an important role in the surface deformation field. The observed surface subsidence in period V can be explained either by inflation or deflation of the deformation source, depending on whether η c is lower or higher than ~1019 Pa s, respectively. The subsequent surface recovery in period VI requires an increase in the inflation rate for any case of η c . The volume of magma accumulated since the 1914 eruption (ΔV) is predicted to be ~0.6–0.8 km3 in 2020. The estimated inflation/deflation rate of the deformation source, in comparison with the observed volumetric eruption rate, reveals that a significant eruption occurs only when ΔV is ~0.4 km3 or more. However, even if ΔV is beyond the critical value, the eruptive activity does not depend on ΔV. The temporal inflation rate inferred in this study gives an opportunity to discuss eruptive dynamics in relation to the observed eruptive activity, such as a quantitative relationship between magma supply rate into, and discharge rate from, a deformation source.

Short-period flow oscillation during eruptions of Onikobe geyser, NE Japan: Insights from thermal infrared observation and acoustic measurements

Physics of a toy geyser

Lava fountain jet noise during the 2018 eruption of fissure 8 of Kīlauea volcano

Potential fluid flow pathways and the geothermal structure of Kuju revealed by azimuthal anisotropic ambient noise tomography

Ground deformations associated with an overpressurized hydrothermal systems at Azuma volcano (Japan) revealed by InSAR data

Spatiotemporal Variations and Postseismic Relaxation Process Around Mt. Fuji, Japan, During and After the 2011 Tohoku‐Oki Earthquake

Investigating Mt. Etna lava fountains by seismic and infrasonic signals: 14-21 February 2021 case study

Very-long-period signal reveals lava lake sloshing and its interaction with a deep reservoir in Nyiragongo volcano

Experimental investigation on geysering induced intermittent flow oscillation

Hidden earthquakes unveil the dynamic evolution of a large-scale explosive eruption

Hybrid long-period volcanic events observed in off Nicobar region, the Andaman Sea from a passive OBS experiment

Variable inflation rate of a magmatic deformation source beneath Aira caldera after the 1914 eruption of Sakurajima volcano: Inferences from a linear Maxwell viscoelastic model constrained by geodetic data

Fitting Jet Noise Similarity Spectra to Volcano Infrasound Data

Geophysical imaging of the shallow geyser and hydrothermal reservoir structures of Spouter Geyser, Yellowstone National Park: Geyser Dynamics I

Temporal Changes in Seismic Velocity and Attenuation at The Geysers Geothermal Field, California, from Double‐Difference Tomography

A numerical model for precursory time sequences of the phreatic eruptions of Mt. Ontake, central Japan

Electrotomography Studies in the Uzon Caldera at Shaman Geyser (Kamchatka Peninsula)

A water flow model of the active crater lake at Aso volcano, Japan: fluctuations of magmatic gas and groundwater fluxes from the underlying hydrothermal system

Time‐Resolved Trigger Processes Leading to the Plinian Eruptions at Sakurajima Volcano, Japan

Response of the underground environment of the KAGRA observatory against the air-pressure disturbance from the Tonga volcano eruption on January 15th, 2022

High Skewness Infrasound from the eruption of Nabro Volcano, Eritrea: Comparison with supersonic jet and rocket engine data