Deformation Source Revealed From Leveling Survey in Jigokudani Valley, Tateyama Volcano, Japan

We modeled vertical deformation detected from leveling survey in Jigokudani valley, Tateyama volcano, central Japan. In Jigokudani valley, uplift of 4 cm/year was previously detected during the period from 2007 to 2010 by Interferometric Synthetic Aperture Radar (InSAR). To confirm whether this inflation has continued to the present, we conducted leveling survey in Jigokudani valley since 2015. Most bench marks showed subsidence up to 5.6 cm during the four-year period from October 2016 to September 2020, while a bench mark locates at the center of the leveling route uniquely showed uplift of 1.6 cm. We applied a dislocation source model to the deformation using a grid search method. A crack with a length of 350 m, a width of 100 m, a strike of N117°E and a dip of 61° is located at a depth of 50 m near the center of Jigokudani valley (Koya jigoku and the new fumarolic area) where highly activating recently. Closing of the crack of 344 cm yields volume decreases of 120,400 m3. Striking direction of the crack is parallel to the line of which are old explosion craters (Mikurigaike and Midorigaike ponds) and corresponds to current maximum compressive stress field in the region of Hida Mountains including Tateyama volcano. The deformation source of the previous period from 2007 to 2010 detected from InSAR was estimated to be at a depth of 50 m and a gas chamber was correspondingly found from the audio-frequency magnetotelluric (AMT) survey. The estimated crack in this study is also located at a similar position of the gas chamber which was also identified from AMT survey. During the period from 2015 to 2016, the crack opened (i.e., inflated) and the inflation stopped during the next one-year period from 2016 to 2017. During the period from 2017 to 2020, the crack turned to closing (i.e., deflation), probably because of the increase in emission of volcanic fluid or gas with a formation of a new crater at the western side of Jigokudani valley (Yahata jigoku) during the period from 2017 to 2018.

Modeling long-term volcanic deformation at Kusatsu-Shirane and Asama volcanoes, Japan, using the GNSS coordinate time series

Long-term deformation of Kusatsu-Shirane and Asama volcanoes in central Japan were investigated using Global Navigation Satellite System (GNSS) measurements. Large postseismic deformation caused by the 2011 Tohoku earthquake—which obscures the long-term volcanic deformation—was effectively removed by approximating the postseismic and other recent tectonic deformation in terms of quadrature of the geographical eastings/northings. Subsequently, deformation source parameters were estimated by the Markov Chain Monte Carlo (MCMC) method and linear inversion, employing an analytical model that calculates the deformation from an arbitrary oriented prolate/oblate spheroid. The deformation source of Kusatsu-Shirane volcano was found to be a sill-like oblate spheroid located a few kilometers northwest of the Yugama crater at a depth of approximately 4 $$\text {km}$$ km , while that of Asama was also estimated to be a sill-like oblate spheroid beneath the western flank of the edifice at a depth of approximately 12 $$\text {km}$$ km , along with the previously reported shallow east–west striking dike at a depth of approximately 1 $$\text {km}$$ km . It was revealed that (1) volume changes of the Kusatsu-Shirane deformation source and the shallow deformation source of Asama were correlated with the volcanic activities of the corresponding volcanoes, and (2) the Asama deep source has been steadily losing volume, which may indicate that the volcano will experience fewer eruptions in the near future.

A New Analysis of Caldera Unrest through the Integration of Geophysical Data and FEM Modeling: The Long Valley Caldera Case Study

The Long Valley Caldera, located at the eastern edge of the Sierra Nevada range in California, has been in a state of unrest since the late 1970s. Seismic, gravity and geodetic data strongly suggest that the source of unrest is an intrusion beneath the caldera resurgent dome. However, it is not clear yet if the main contribution to the deformation comes from pulses of ascending high-pressure hydrothermal fluids or low viscosity magmatic melts. To characterize the nature of the intrusion, we developed a 3D finite element model which includes topography and crust heterogeneities. We first performed joint numerical inversions of uplift and Electronic Distance Measurement baseline length change data, collected during the period 1985–1999, to infer the deformation-source size, position, and overpressure. Successively, we used this information to refine the source overpressure estimation, compute the gravity potential and infer the intrusion density from the inversion of deformation and gravity data collected in 1982–1998. The deformation source is located beneath the resurgent dome, at a depth of 7.5 ± 0.5 km and a volume change of 0.21 ± 0.04 km3. We assumed a rhyolite compressibility of 0.026 ± 0.0011 GPa−1 (volume fraction of water between 0% and 30%) and estimated a reservoir compressibility of 0.147 ± 0.037 GPa−1. We obtained a density of 1856 ± 72 kg/m3. This density is consistent with a rhyolite melt, with 20% to 30% of dissolved hydrothermal fluids.

Combined Geodetic and Seismological Study of the December 2020 Mw = 4.6 Thiva (Central Greece) Shallow Earthquake

On 2 December 2020, a moderate and shallow Mw = 4.6 earthquake occurred in Boeotia (Central Greece) near the city of Thiva. Despite its magnitude, the co-seismic ground deformation field was detectable and measurable by Sentinel-1, ascending and descending, synthetic aperture interferometry radar (InSAR) acquisitions. The closest available GNSS station to the epicenter, located 11 km west, measured no deformation, as expected. We proceeded to the inversion of the deformation source. Moreover, we reassessed seismological data to identify the activated zone, associated with the mainshock and the aftershock sequence. Additionally, we used the rupture plane information from InSAR to better determine the focal mechanism and the centroid location of the mainshock. We observed that the mainshock occurred at a shallower depth and the rupture then expanded downdip, as revealed by the aftershock distribution. Our geodetic inversion modelling indicated the activation of a normal fault with a small left-lateral component, length of 2.0 km, width of 1.7 km, average slip of 0.2 m, a low dip angle of 33°, and a SW dip-direction. The inferred fault top was buried at a depth of ~0.5 km, rooted at a depth of ~1.4 km, with its geodetic centroid buried at 1.0 km. It was aligned with the Kallithea fault. In addition, the dip-up projection of the modeled fault to the surface was located very close (~0.4 km SW) to the mapped (by existing geological observations) trace of the Kallithea fault. The ruptured area was settled in a transition zone. We suggest the installation of at least one GNSS and seismological station near Kallithea; as the activated zone (inferred by the aftershock sequence and InSAR results) could yield events with M≥5.0, according to empirical laws relating to rupture zone dimensions and earthquake magnitude.

Inflating Source Imaging and Stress/Strain Field Analysis at Campi Flegrei Caldera: The 2009–2013 Unrest Episode

In this study, we analyze the 2009–2013 uplift phenomenon at Campi Flegrei (CF) caldera in terms of temporal and spatial variations in the stress/strain field due to the effect of an inflating source. We start by performing a 3D stationary finite element (FE) modeling of X-band COSMO-SkyMed DInSAR and GPS mean velocities to retrieve the geometry and location of the deformation source. The modeling results suggest that the best-fit source is a three-axis oblate spheroid ~3 km deep, which is mostly elongated in the NE–SW direction. Furthermore, we verify the reliability of model results by calculating the total horizontal derivative (THD) of the modeled vertical velocity component; the findings emphasize that the THD maxima overlap with the projection of source boundaries at the surface. Then, we generate a 3D time-dependent FE model, comparing the spatial and temporal distribution of the shear stress and volumetric strain with the seismic swarms beneath the caldera. We found that low values of shear stress are observed corresponding with the shallow hydrothermal system where low-magnitude earthquakes occur, whereas high values of shear stress are found at depths of about 3 km, where high-magnitude earthquakes nucleate. Finally, the volumetric strain analysis highlights that the seismicity occurs mainly at the border between compression and dilatation modeled regions, and some seismic events occur within compression regions.

Episodic transport of discrete magma batches beneath Aso volcano

Magma ascent, storage, and discharge in the trans-crustal magma plumbing system are key to long-term volcanic output and short-term eruption dynamics. Petrological analytics, geodetic deformations and mechanical modeling have shaped the current understanding of magma transport. However, due to the lack of observations, how a distinct magma batch transports from a crystal-rich mush region to a crystal-poor pool with eruptible magma remains enigmatic. Through stacking of tilt and seismic waveform data, we find that episodic long-period tremors (LPTs) located near sea level beneath the Aso volcano are accompanied by a synchronous deformation event, which initiates ~50 seconds before individual LPT event and concludes seconds after. The episodic deformation source corresponds to either an inflation or a deflation event located ~3 km below sea level, with a major volumetric component (50-440 cubic meters per event) and a minor high-angle normal-fault component. We suggest that these deformation events likely represent short-lived, episodic upward transport of discrete magma batches accompanied by high-angle shear failure near the roof of the inferred magma chamber at relatively high temperature, whereas their recurrences, potentially composition dependent, are regulated by the brittle-to-ductile transition rheology under low differential stress and high strain rate due to the surge of magma from below, regulating long-term volcanic output rate. The magma ascent velocity, decompression rates, and cumulative magma output deduced from the episodic deformation events before recent eruptions in Aso volcano are compatible with retrospective observations of the eruption style, tephra fallouts, and plume heights, promising real-time evaluation of upcoming eruptions.

Modeling long-term volcanic deformations at the Kusatsu-Shirane and Asama volcanoes, Japan using the GNSS coordinate time series

Long-term deformations of the Kusatsu-Shirane and Asama volcanoes in central Japan were investigated using Global Navigation Satellite System (GNSS) measurements. Large postseismic deformations caused by the 2011 Tohoku earthquake — which obscure the long-term volcanic deformations — were effectively removed by approximating the postseismic and other recent tectonic deformations in terms of quadrature of the geographical eastings/northings. Subsequently, deformation source parameters were estimated by the Markov Chain Monte-Carlo (MCMC) method and linear inversion. The deformation source of the Kusatsu-Shirane volcano was found to be a sill-like oblate spheroid located a few kilometers northwest of the Yugama crater at a depth of approximately five km, while that of Asama was also estimated to be a sill-like oblate spheroid located at the western flank of the edifice at a depth of approximately 13 km, along with the previously reported shallow east-west striking dike at a depth of approximately 1 km. It was revealed that 1) volume changes of the Kusatsu-Shirane deformation source and the shallow deformation source of Asama were correlated with the volcanic activities of the corresponding volcanoes, and 2) the Asama deep source has been steadily losing volume, which may indicate that the volcano will experience less eruptions in the near future.