Rheology of weak lithosphere beneath active volcanoes of NE Japan: Insights from postseismic deformation of 2011 Tohoku-oki earthquake

Seismicity records during the 1990’s reveals that large inland earthquakes tend to concentrate near several active volcanoes in the central part of northeastern (NE) Japan (Hasegawa & Yamamoto, 1994 Tectonophysics). The inland seismicity around active volcanoes could be related to localized zones of high strain contraction detected by the GNSS measurements in 1997–2001 (Miura et al., 2004 JGR). Several studies speculated that the localized strain contraction is caused by inelastic deformation of weak lithosphere beneath the active volcanoes (Hasegawa et al., 2004 J. Seismol. Soc.). Such weak lithosphere (i.e., low-viscosity zone or LVZ) is inferred from high heat-flow observations (Tanaka et al., 2004 EPS), lithospheric strength simulation (Shibazaki et al., 2016 GRL; Muto, 2011 Tectonophysics) and seismic-velocity tomography (Hasegawa et al., 2005 JGR). However, because of complex interplay between elastic and inelastic processes during steady-state (i.e., interseismic) crustal deformation, the physical mechanism related to inelastic deformation is still poorly understood.

When the M_w9.0 2011 Tohoku-oki earthquake occurred, strong surface deformations were observed locally near the active volcanoes (Takada & Fukushima, 2013 Nat. Geosci.) and continued for several years after the mainshock (Muto et al., 2016 GRL). Past studies (e.g., Sun et al., 2014 Nature) advocated that the earthquake-related inelastic processes such as viscoelastic mantle relaxation dominates the crustal deformations in the postseismic period. In the present study, we identified localized strain contractions near the active volcanoes by extracting the short-wavelength strain rate (Meneses-Gutierrez & Sagiya, 2016 EPSL) from the GNSS observations during 2012–2014. We explained these localized strain contraction by building three-dimensional rheological models of small-scale LVZs beneath five active volcanoes of NE Japan. We simulated the volumetric deformation of viscoelastic LVZs using power-law Burgers rheology, which previously succeeded to explain the large-scale postseismic deformation of the 2011 Tohoku-oki earthquake (Agata et al., 2019 Nat. Commun.; Muto et al., 2019 Sci. Adv.; Dhar et al., 2022 GJI). The power-law Burgers rheology represents the bi-phasic nature of rock deformations (rapid transient with subsequent steady state) and power-law dependency of strain rate to evolving stress (proxy of dominating dislocation creep in high-stress mantle condition) (Muto et al., 2019 Sci. Adv. and references therein).

We found that the localized strain contraction near the active volcanoes can be explained by small-scale LVZs which have narrow tops of 10–20 km and wide roots of 60–100 km width. Our results conclude the minimum depths of the tops and roots of LVZs are 15 km and 40 km, respectively. The geometries of the LVZs vary (e.g., upright conic or inclined shape) from volcano to volcano. The effective viscosities of the LVZs are in the order of 10¹⁷ Pa·s immediately after the earthquake and increases to the order of 10¹⁸ Pa·s over the 3 years of postseismic deformation. Our results agree with the results of several past studies (Ohzono et al., 2012 EPS; Hu et al., 2014 EPS; Muto et al., 2016 GRL) who investigated the lithospheric rheology near Mt. Naruko using the postseismic surface displacements of the 2011 Tohoku-oki and 2008 Iwate-Nairiku earthquakes.