Segmentation of the Izu-Bonin and Mariana slabs based on the analysis of the Benioff seismicity distribution and regional tomography results

We present a new model of P and S velocity anomalies in the mantle down to a depth of 1300 km beneath the Izu-Bonin and Mariana (IBM) arcs. This model is derived based on tomographic inversion of global travel time data from the revised ISC catalogue. The results of inversion are thoroughly verified using a series of different tests. The obtained model is generally consistent with previous studies by different authors. We also present the distribution of relocated deep events projected to the vertical surface along the IBM arc system. Unexpectedly, the seismicity forms elongated vertical clusters instead of horizontal zones indicating phase transitions in the slab. We propose that these vertical seismicity zones mark zones of intense deformation and boundaries between semi-autonomous segments of the subducting plate. The P and S seismic tomography models consistently display the slab as prominent high-velocity anomalies coinciding with the distribution of deep seismicity. We can distinguish at least four segments which subduct differently. The northernmost segment of the Izu-Bonin arc has the gentlest angle of dipping which is explained by backward displacement of the trench. In the second segment, the trench stayed at the same location, and we observe the accumulation of the slab material in the transition zone and its further descending to the lower mantle. In the third segment, the trench is moving forward causing the steepening of the slab. Finally, for the Mariana segment, despite the backward displacement of the arc, the subducting slab is nearly vertical. Between the Izu-Bonin and Mariana arcs we clearly observe a gap which can be traced down to about 400 km in depth. Based on joint consideration of the tomography results and the seismicity distribution, we propose two different scenarios of the subduction evolution in the IBM zone during the recent time, depending on the reference frame of plate displacements. In the first case, we consider the movements in respect to the Philippine Plate, and explain the different styles of the subduction by the relative backward and forward migrations of the trench. In the second case, all the elements of the subduction system move westward in respect to the stable Asia. Different subduction styles are explained by the "anchoring" of selected segments of the slab, different physical properties of the subducting plate and the existence of buoyant rigid blocks related to sea mount and igneous provinces

Modeling the impact of melt on seismic properties during mountain building

Initiation of partial melting in the mid/lower crust causes a decrease in P-wave and S-wave velocities; recent studies imply that the relationship between these velocities and melt is not simple. We have developed a modelling approach to assess the combined impact of various melt and solid phase properties on seismic velocities and anisotropy. The modelling is based on crystallographic preferred orientation (CPO) data measured from migmatite samples, allowing quantification of the variation of seismic velocities with varying melt volumes, shapes, orientations, and matrix anisotropy. The results show non-linear behaviour of seismic properties as a result of the interaction of all of these physical properties, which in turn depend on lithology, stress regime, strain rate, pre-existing rock fabrics, and pressure-temperature conditions. This non-linear behaviour is evident when applied to a suite of samples from a traverse across a migmatitic shear zone in the Seiland Igneous Province, Northern Norway. Critically, changes in solid phase composition and CPO, and melt shape and orientation with respect to the wave propagation direction can result in huge variations in the same seismic property even if the melt fraction remains the same. A comparison with surface wave interpretations from tectonically active regions highlights the issues in current models used to predict melt percentages or partially molten regions. Interpretation of seismic data to infer melt percentages or extent of melting should, therefore, always be underpinned by robust modelling of the underlying geological parameters combined with examination of multiple seismic properties in order to reduce uncertainty of the interpretation

Ambient Seismic Noise Image of the Structurally Controlled Heat and Fluid Feeder Pathway at Campi Flegrei Caldera

The TIDES-COST Action (STSM-ES1401-34011) provided a travel grant to framework the research project. The Japan Society for the Promotion of Science - Short-Term Fellowship (JSPS/OF215/022) financed the work, undertaken at Tohoku University and concluded at the University of Aberdeen. We thank Giuseppe Vilardo and Eliana Bellucci Sessa for providing the geomorphological maps, and Simona Petrosino and Paola Cusano for the P- and S-wave pickings used to localise the seismicity. Informal revisions from Guido Ventura, Nick Rawlinson and Chris Kilburn helped us improving the analyses and interpretation, respectively. We acknowledge the help of Naveed Khan in parallelising the codes and two anonymous reviewers who improved the analysis, interpretation, and readibility with their comments. All data to reproduce the maps can be downloaded from the World Data Center PANGAEA data repository, permanent link https://doi.pangaea.de/10.1594/PANGAEA.890238.Peer reviewedPublisher PD

Post-supereruption recovery at Toba Caldera

Large calderas, or supervolcanoes, are sites of the most catastrophic and hazardous events on Earth, yet the temporal details of post-supereruption activity, or resurgence, remain largely unknown, limiting our ability to understand how supervolcanoes work and address their hazards. Toba Caldera, Indonesia, caused the greatest volcanic catastrophe of the last 100 kyr, climactically erupting ~74 ka. Since the supereruption, Toba has been in a state of resurgence but its magmatic and uplift history has remained unclear. Here we reveal that new 14 C, zircon U-Th crystallization and (U-Th)/He ages show resurgence commenced at 69.7±4.5 ka and continued until at least ~2.7 ka, progressing westward across the caldera, as reflected by post-caldera effusive lava eruptions and uplifted lake sediment. The major stratovolcano north of Toba, Sinabung, shows strong geochemical kinship with Toba, and zircons from recent eruption products suggest Toba's climactic magma reservoir extends beneath Sinabung and is being tapped during eruptions

A large magmatic sill complex beneath the Toba caldera

International audienceAn understanding of the formation of large magmatic reservoirs is a key issue for the evaluation of possible strong volcanic eruptions in the future. We estimated the size and level of maturity of one of the largest volcanic reservoirs, based on radial seismic anisotropy. We used ambient-noise seismic tomography below the Toba caldera (in northern Sumatra) to observe the anisotropy that we interpret as the expression of a fine-scale layering caused by the presence of many partially molten sills in the crust below 7 kilometers. This result demonstrates that the magmatic reservoirs of present (non-eroded) supervolcanoes can be formed as large sill complexes and supports the concept of the long-term incremental evolution of magma bodies that lead to the largest volcanic eruptions

Rapid changes in magma storage beneath the Klyuchevskoy group of volcanoes inferred from time-dependent seismic tomography

We present the results of time-dependent local earthquake tomography for the Kluchevskoy group of volcanoes in Kamchatka, Russia. We consider the time period from 1999 to 2009, which covers several stages of activity of Kluchevskoy and Bezymianny volcanoes. The results are supported by synthetic tests that recover a common 3D model based on data corresponding to different time windows. Throughout the period, we observe a robust feature below 25 km depth with anomalously high Vp/Vs values (up to 2.2). We interpret this feature as a channel bringing deep mantle materials with high fluid and melt content to the bottom of the crust. This mantle channel directly or indirectly determines the activity of all volcanoes of the Kluchevskoy group. In the crust, we model complex structure that varies over time. During the pre-eruptive period, we detected two levels of potential magma storage: one in the middle crust at 10–12 km depth and one close to the surface just below Kluchevskoy volcano. In 2005, a year of powerful eruptions of Kluchevskoy and Besymiyanny volcanoes, we observe a general increase in Vp/Vs throughout the crust. In the relaxation period following the eruption, the Vp/Vs values are generally low, and no strong anomalous zones in the crust are observed. We propose that very rapid variations in Vp/Vs are most likely due to abrupt changes in the stress and deformation states, which cause fracturing and the active transport of fluids. These fluids drive more fracturing in a positive feedback system that ultimately leads to eruption. We envision the magma reservoirs beneath the Kluchevskoy group as sponge-structured volumes that may quickly change the content of the molten phases as fluids pulse rapidly through the system

Evidence for high fluid/melt content beneath Krakatau volcano (Indonesia) from local earthquake tomography

Within the KRAKMON project for multiparameter monitoring of Anak Krakatau volcano (Indonesia), a network of temporary stations was installed on the islands of the Krakatau complex as well as in the surrounding areas of the Sunda Strait, Sumatra and Java. The network was operated from June 2005 until January 2006. More than 700 local events were recorded during this experiment, and travel times from these events were used to perform a tomographic inversion for P and S velocities and for the Vp/Vs ratio. In this study, special attention was paid to the validation of the computed model based on different tests, such as inversion of independent data subsets and synthetic modeling. Although the network configuration and the distribution of the events are not favorable for high-quality tomographic imaging, we have obtained some important and robust features which give information about sources of volcanic activity in the Krakatau complex. The most interesting feature of this study is a zone of high Vp/Vs ratio beneath the Krakatau complex. At depths down to 4 km depth we observe anticorrelation of higher P- and lower S-velocities that leads to Vp/Vs ratio higher than 2. This is a probable indicator of the presence of partially molten and/or with high fluid content material with a composition corresponding to deeper layers. It is important that the anomaly of high Vp/Vs ratio beneath the Krakatau complex appears to be separated in two parts at a depth of 5-6 km. This fits to results of geobarometric analysis that presume the existence of several levels of magma chambers beneath Anak Krakatau.</p

Asymmetric caldera-related structures in the area of the Avacha group of volcanoes in Kamchatka as revealed by ambient noise tomography and deep seismic sounding

Avacha group includes two active and potentially dangerous volcanoes, Avachinsky and Koryaksky, located close to Petropavlovsk-Kamchatsky, the main city of Kamchatka. We present the results of two independent seismic studies of shallow crustal structures beneath the Avacha group based on passive and active source observations. The first study is based on the analysis of continuous recording by 11 seismic stations installed over the Avacha group in 2012 and 7 permanent stations in the same region. We present a series of 2D Rayleigh-wave group velocity maps based on correlation of ambient noise, that were then converted into 3D distribution of shear wave velocity. The second work was based on the reprocessing of an active source deep seismic sounding profile across the Avachinsky volcano that was shot in 1982–1984. We made the analysis of travel times of refracted waves using a 2D tomography inversion. The resulting seismic models appear to be consistent with each other and show clear low-velocity zone to the SW of the Avachinsky volcano and high velocity structures to NE. These observations also agree with the existing gravity and magnetotelluric measurements. Based on the obtained seismic models we identify two large buried calderas and large lava flows that are thought to be related to a series of large eruption episodes of Avachinsky occurred within the last 30,000 years

Asymmetric caldera-related structures in the area of the Avacha group of volcanoes in Kamchatka as revealed by ambient noise tomography and deep seismic sounding

International audienceAvacha group includes two active and potentially dangerous volcanoes, Avachinsky and Koryaksky, located close to Petropavlovsk-Kamchatsky, the main city of Kamchatka. We present the results of two independent seismic studies of shallow crustal structures beneath the Avacha group based on passive and active source observations. The first study is based on the analysis of continuous recording by 11 seismic stations installed over the Avacha group in 2012 and 7 permanent stations in the same region. We present a series of 2D Rayleigh-wave group velocity maps based on correlation of ambient noise, that were then converted into 3D distribution of shear wave velocity. The second work was based on the reprocessing of an active source deep seismic sounding profile across the Avachinsky volcano that was shot in 1982-1984. We made the analysis of travel times of refracted waves using a 2D tomography inversion. The resulting seismic models appear to be consistent with each other and show clear low-velocity zone to the SW of the Avachinsky volcano and high velocity structures to NE. These observations also agree with the existing gravity and magnetotelluric measurements. Based on the obtained seismic models we identify two large buried calderas and large lava flows that are thought to be related to a series of large eruption episodes of Avachinsky occurred within the last 30,000 years