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.

Model of Surface Subsidence of Pyroclastic Flow: the August 29, 2019 Shiveluch Volcano Eruption, Kamchatka

The 2019–2021 Eruptions of Klyuchevskoi Volcano, Kamchatka

Material Properties and Triggering Mechanisms of an Andesitic Lava Dome Collapse at Shiveluch Volcano, Kamchatka, Russia, Revealed Using the Finite Element Method

Modeling a Potential Collapse on Koryaksky Volcano and Assessing the Volcanic Hazard for the Elizovo–Petropavlovsk Agglomeration, Kamchatka

The April 11, 2023 Catastrophic Explosive Eruption of Sheveluch Volcano, Kamchatka

On the Geomagnetic and Ionospheric Variations after the 2023 Strong Eruption of the Shiveluch Volcano

Seismic structure beneath Ebeko Volcano and surrounding areas of Paramushir Island (Kuril Arc) inferred from local earthquake tomography

Declining magma supply to a poroelastic magma mush explains long-term deformation at Soufrière Hills Volcano

Postseismic Processes in the Region of the July 29, 2021 Chignik Earthquake, Alaska: Part I. Modeling Results

Anatomy of the Bezymianny volcano merely before an explosive eruption on 20.12.2017

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

Response of the Lower and Upper Ionosphere after the Eruption of Shiveluch Volcano on April 10, 2023

Upward Magma Migration Within the Multi-Level Plumbing System of the Changbaishan Volcano (China/North Korea) Revealed by the Modeling of 2018–2020 SAR Data

Magmatic and sedimentary structure beneath the Klyuchevskoy Volcanic Group, Kamchatka, from Ambient Noise Tomography

Seismic Velocity Changes beneath Ebeko Volcano (Kuril Islands) Based on Observations in July–August 2021

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

The Ash Discharged by Sheveluch Volcano (Kamchatka, Russia) in April 2023 As a Source of Water-Soluble Salts

Dynamics of the 2021 unrest at Changbaishan Tianchi volcano from ALOS-2/PALSAR-2 and seismic data

A Characteristic Electrostatic Structure of Eruptive Plumes Emitted by Large Explosive Eruptions of Shiveluch and Bezymianny Volcanoes, Kamchatka

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

The rebirth and evolution of Bezymianny volcano, Kamchatka after the 1956 sector collapse