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.

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

Quantitative Assessment of Temporal Changes in Subaqueous Hydrothermal Activity in Active Crater Lakes During Unrest Based on a Time-Series of Lake Water Chemistry

From sink to volcanic source: Unravelling missing terrestrial eruption records by characterization and high‐resolution chronology of lacustrine volcanic density flow deposits, Lake Inawashiro‐ko, Fukushima, Japan

High spatio-temporal resolution observations of crater lake temperatures at Kawah Ijen volcano, East Java, Indonesia

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

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

Branched Conduit Structure Beneath the Active Craters of Sakurajima Volcano Inferred From Muography

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

Precursory crater contraction associated with the 2017 eruption of Shinmoe-dake volcano (Japan) detected by PALSAR-2 and Sentinel-1 InSAR

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

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

An easier approach for helium isotope flux estimation in a submerged caldera

Episodic transport of discrete magma batches beneath Aso volcano

Increase in magma supply to Sakurajima volcano’s (Japan) shallow magma chamber over the past 500 years

Asymmetric magma plumbing system beneath Axial Seamount based on full waveform inversion of seismic data

Groundwater flow in the Ischia volcanic island (Italy) and its implications for thermal water abstraction

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

Chemistry, Growth, and Fate of the Unique, Short‐Lived (2019–2020) Water Lake at the Summit of Kīlauea Volcano, Hawaii

Changes in the magmatic plumbing system associated with the Gotemba sector collapse at Mount Fuji, Japan

Altimeter-derived marine gravity variations reveal the magma mass motions within the subaqueous Nishinoshima volcano, Izu–Bonin Arc, Japan

Inter-annual variations in the chemical and thermal properties of well waters from the volcanic hydrothermal areas: Usu Volcano, Hokkaido, Japan