Influence of a low viscosity zone on the evolution of post-eruption deformation: A case study of the crustal deformation of Aira Caldera after the 1914 eruption of Sakurajima Volcano

Spatial viscosity variations in the crust can be inferred from geophysical imaging and may be essential for interpreting volcano deformations that show more complex behaviours than what simple uni-viscous models predict. In southern Kyushu, Japan, recent seismic imaging implies the presence of a low viscosity zone (LVZ) beneath Aira Caldera, and geodetic data in and around the caldera require uni-viscous model to have lower and higher viscosities earliest and later, respectively, in the period following the 1914 eruption of Sakurajima Volcano. Here, we use a 3D finite element model comprising an elastic layer underlain by a linear Maxwell viscoelastic layer to examine the influence of a LVZ on ground surface displacement in response to two different deformation source modes, i.e., instantaneous source deflation during a major eruption and subsequent continuous source inflation due to magma recharge. A LVZ is introduced into the viscoelastic layer by gradually reducing its viscosity towards the centre. The behaviour of the LVZ model quantified by comparison with the behaviour of a uniform viscosity (UNV) model reveals that, for a given LVZ structure, an apparent UNV model that best represents the LVZ model displacement in response to the instantaneous deflation has lower viscosity than that in response to gradual inflation, i.e., the rate-controlling viscosities of the LVZ model are those in the inner and outer parts of the LVZ for the instantaneous source deflation and subsequent continuous inflation, respectively. Such LVZ model behaviour, for a LVZ spatial extent comparable with the imaged low velocity anomaly, explains well the geodetic data in and around Aira Caldera at any stage after the 1914 eruption, as the predominant viscoelastic response changes with time from that to the instantaneous deflation to that to the subsequent continuous inflation. This study highlights the need for interdisciplinary investigations that integrate geodetic and geophysical datasets to better understand volcanic unrest.