Lithospheric delamination as the driving mechanism of intermediate-depth seismicity in the Bucaramanga Nest, Colombia.

The Bucaramanga nest (BN) is an area of exceptionally strong intermediate-depth seismicity localized in a narrow zone at 150-170 km depth beneath the continental plate in Colombia. To explain the very unusual mantle seismicity cluster in this area, we built a seismic velocity model in the vicinity of BN with the use of local earthquake tomography and developed a numerical hydromechanical model. Our seismic model shows a strong high-velocity anomaly at 130-167 km coinciding with the BN seismicity. The relocated seismicity can be separated in two clusters. We propose that the upper BN cluster at ~ 130 km depth is attributed to dehydration embrittlement, whereas the lower BN cluster at ~ 150 km depth coinciding with the high-velocity body is caused by lithospheric delamination, creating a "drip" that falls down over the subducting oceanic plate, enhancing fluid release from the slab, potentially increasing seismicity.

Spatio-temporal velocity variations observed during the pre-eruptive episode of La Palma 2021 eruption inferred from ambient noise interferometry.

On Sept. 19th, 2021, a volcanic eruption began on the island of La Palma (Canary Islands, Spain). The pre-eruptive episode was characterized by seismicity and ground deformation that started only 9.5 days before the eruption. In this study, we applied seismic interferometry to the data recorded by six broadband seismic stations, allowing us to estimate velocity variations during the weeks preceding the eruption. About 9.5 days before the eruption, we observed a reduction in the seismic velocities is registered next to the eruptive centers that opened later. Furthermore, this zone overlaps with the epicenters of a cluster of volcano-tectonic earthquakes located at shallow depth (< 4 km) and detached from the main cluster of deeper seismicity. We interpret the decrease in seismic velocities and the occurrence of such a shallow earthquake cluster as the effect of hydrothermal fluid released by the ascending magma batch and reaching the surface faster than the magma itself.

Low-degree mantle melting controls the deep seismicity and explosive volcanism of the Gakkel Ridge.

The world's strongest known spreading-related seismicity swarm occurred in 1999 in a segment of the Gakkel Ridge located at 85°E as a consequence of an effusive-explosive submarine volcanic eruption. The data of a seismic network deployed on ice floes were used to locate hundreds of local earthquakes down to ∼25 km depth and to build a seismic tomography model under the volcanic area. Here we show the seismicity and the distribution of seismic velocities together with the 3D magmatic-thermomechanical numerical model, which demonstrate how a magma reservoir under the Gakkel Ridge may form, rise and trigger volcanic eruptions in the rift valley. The ultraslow spreading rates with low mantle potential temperatures appear to be a critical factor in the production of volatile-rich, low-degree mantle melts that are focused toward the magma reservoirs within narrow magmatic sections. The degassing of these melts is the main cause of the explosive submarine eruptions.

Structure of the magma plumbing system beneath Semisopochnoi Island (Aleutian Arc) inferred from seismic tomography.

Semisopochnoi Island is a remote and little-studied volcanic island in the western part of the Aleutian Arc. The existence of several active volcanic centers and a 5000-7000-year-old large caldera makes this island an important site for volcanic hazard assessment in the Northern Pacific. Based on local seismicity data recorded by six permanent seismic stations, we created the seismic tomography model, including the 3D distributions of Vp, Vs, and Vp/Vs ratios to a depth of 10 km. This model provides the first geophysical insight into the interior structure of Semisopochnoi Island and sheds light on the processes in the magma plumbing system beneath all volcanic centers on the island. At depths of 5-10 km, we observed a columnar-shaped high Vp/Vs-ratio anomaly below the caldera in the central part of the island, which likely represents the steady magma conduit. This conduit is headed by a prominent high Vp/Vs-ratio anomaly located 3-5 km directly below the caldera, which represents the magma reservoir feeding Cerberus and other Holocene-aged volcanic centers on Semisopochnoi Island.

Rapid magma ascent beneath La Palma revealed by seismic tomography.

For the first time, we obtained high-resolution images of Earth's interior of the La Palma volcanic eruption that occurred in 2021 derived during the eruptive process. We present evidence of a rapid magmatic rise from the base of the oceanic crust under the island to produce an eruption that was active for 85 days. This eruption is interpreted as a very accelerated and energetic process. We used data from 11,349 earthquakes to perform travel-time seismic tomography. We present high-precision earthquake relocations and 3D distributions of P and S-wave velocities highlighting the geometry of magma sources. We identified three distinct structures: (1) a shallow localised region (< 3 km) of hydrothermal alteration; (2) spatially extensive, consolidated, oceanic crust extending to 10 km depth and; (3) a large sub-crustal magma-filled rock volume intrusion extending from 7 to 25 km depth. Our results suggest that this large magma reservoir feeds the La Palma eruption continuously. Prior to eruption onset, magma ascended from 10 km depth to the surface in less than 7 days. In the upper 3 km, melt migration is along the western contact between consolidated oceanic crust and altered hydrothermal material.

Seismic tremor reveals active trans-crustal magmatic system beneath Kamchatka volcanoes.

The occurrence and the style of volcanic eruptions are largely controlled by the ways in which magma is stored and transported from the mantle to the surface through the crust. Nevertheless, our understanding of the deep roots of volcano-magmatic systems remains very limited. Here, we use the sources of seismovolcanic tremor to delineate the active part of the magmatic system beneath the Klyuchevskoy Volcanic Group in Kamchatka, Russia. The tremor sources are distributed in a wide spatial region over the whole range of crustal depths connecting different volcanoes of the group. The tremor activity is characterized by rapid vertical and lateral migrations explained by fast pressure transients and dynamic permeability. Our results support the conceptual model of extended and highly dynamic trans-crustal magmatic systems.

Transition from continental rifting to oceanic spreading in the northern Red Sea area.

Lithosphere extension, which plays an essential role in plate tectonics, occurs both in continents (as rift systems) and oceans (spreading along mid-oceanic ridges). The northern Red Sea area is a unique natural geodynamic laboratory, where the ongoing transition from continental rifting to oceanic spreading can be observed. Here, we analyze travel time data from a merged catalogue provided by the Egyptian and Saudi Arabian seismic networks to build a three-dimensional model of seismic velocities in the crust and uppermost mantle beneath the northern Red Sea and surroundings. The derived structures clearly reveal a high-velocity anomaly coinciding with the Red Sea basin and a narrow low-velocity anomaly centered along the rift axis. We interpret these structures as a transition of lithospheric extension from continental rifting to oceanic spreading. The transitional lithosphere is manifested by a dominantly positive seismic anomaly indicating the presence of a 50-70-km-thick and 200-300-km-wide cold lithosphere. Along the forming oceanic ridge axis, an elongated low-velocity anomaly marks a narrow localized nascent spreading zone that disrupts the transitional lithosphere. Along the eastern margins of the Red Sea, several low-velocity anomalies may represent crustal zone of massive Cenozoic basaltic magmatism.

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

Strong explosive eruptions of volcanoes throw out mixtures of gases and ash from high-pressure underground reservoirs. Investigating these subsurface reservoirs may help to forecast and characterize an eruption. In this study, we compare seismic tomography results with remote sensing and petrology data to identify deep and subaerial manifestations of pre-eruptive processes at Bezymianny volcano in Kamchatka shortly before its violent explosion on December 20, 2017. Based on camera networks we identify precursory rockfalls, and based on satellite radar data we find pre-eruptive summit inflation. Our seismic network recorded the P and S wave data from over 500 local earthquakes used to invert for a 3D seismic velocity distribution beneath Bezymianny illuminating its eruptive state days before the eruption. The derived tomography model, in conjunction with the presence of the high-temperature-stable SiO2 polymorph Tridymite in juvenile rock samples , allowed us to infer the coexistence of magma and gas reservoirs revealed as anomalies of low (1.5) and high (2.0) Vp/Vs ratios, respectively, located at depths of 2-3 km and only 2 km apart. The reservoirs both control the current eruptive activity: while the magma reservoir is responsible for episodic dome growth and lava flow emplacements, the spatially separated gas reservoir may control short but powerful explosive eruptions of Bezymianny.

Ongoing formation of felsic lower crustal channel by relamination in Zagros collision zone revealed from regional tomography.

Complex interaction of rheologically contrasting layers within the lithosphere during the collision of continental plates leads to active faulting, which represents a serious hazard to the population and infrastructure. One of the collision scenarios presumes the existence of a middle-lower crustal channel composed of subducted silicic upper crustal rocks, which is thought to exist in the Tibetan-Himalayan system. Based on the results of seismic tomography, we argue that a similar mechanism of crustal channeling takes place beneath the Zagros mountain system in southwestern Iran. The 3D seismic velocity model reveals an inverted crustal architecture of the collision zone, in which the low-velocity felsic (granitic and sedimentary) upper crustal rocks of the Arabian plate form a seismically inactive lower crustal channel below the higher-velocity mafic (basaltic) middle-upper crustal layer of the Iranian crust. Based on existing numerical models, we suggest that the formation of the felsic channel is likely governed by separation (delamination) of the weak felsic upper crust of the subducting Arabian lithosphere and its ductile underplating under rheologically stronger upper-middle crust of the Iranian plate.

Magma plumbing system and seismicity of an active mid-ocean ridge volcano.

At mid-ocean ridges volcanism generally decreases with spreading rate but surprisingly massive volcanic centres occur at the slowest spreading ridges. These volcanoes can host unexpectedly strong earthquakes and vigorous, explosive submarine eruptions. Our understanding of the geodynamic processes forming these volcanic centres is still incomplete due to a lack of geophysical data and the difficulty to capture their rare phases of magmatic activity. We present a local earthquake tomographic image of the magma plumbing system beneath the Segment 8 volcano at the ultraslow-spreading Southwest Indian Ridge. The tomography shows a confined domain of partial melt under the volcano. We infer that from there melt is horizontally transported to a neighbouring ridge segment at 35 km distance where microearthquake swarms and intrusion tremor occur that suggest ongoing magmatic activity. Teleseismic earthquakes around the Segment 8 volcano, prior to our study, indicate that the current magmatic spreading episode may already have lasted over a decade and hence its temporal extent greatly exceeds the frequent short-lived spreading episodes at faster opening mid-ocean ridges.

Breathing of the Nevado del Ruiz volcano reservoir, Colombia, inferred from repeated seismic tomography.

Nevado del Ruiz volcano (NRV), Columbia, is one of the most dangerous volcanoes in the world and caused the death of 25,000 people in 1985. Using a new algorithm for repeated tomography, we have found a prominent seismic anomaly with high values of the Vp/Vs ratio at depths of 2-5 km below the surface, which is associated with a shallow magma reservoir. The amplitude and shape of this anomaly changed during the current phase of unrest which began in 2010. We interpret these changes as due to the ascent of gas bubbles through magma and to degassing of the reservoir. In 2011-2014, most of this gas escaped through permeable roof rocks, feeding surface fumarole activity and leading to a gradual decrease of the Vp/Vs ratio in the reservoir. This trend was reversed in 2015-2016 due to replenishment of the reservoir by a new batch of volatile-rich magma likely to sustain further volcanic activity. It is argued that the recurring "breathing" of the shallow reservoir is the main cause of current eruptions at NRV.

The feeder system of the Toba supervolcano from the slab to the shallow reservoir.

The Toba Caldera has been the site of several large explosive eruptions in the recent geological past, including the world's largest Pleistocene eruption 74,000 years ago. The major cause of this particular behaviour may be the subduction of the fluid-rich Investigator Fracture Zone directly beneath the continental crust of Sumatra and possible tear of the slab. Here we show a new seismic tomography model, which clearly reveals a complex multilevel plumbing system beneath Toba. Large amounts of volatiles originate in the subducting slab at a depth of ∼150 km, migrate upward and cause active melting in the mantle wedge. The volatile-rich basic magmas accumulate at the base of the crust in a ∼50,000 km(3) reservoir. The overheated volatiles continue ascending through the crust and cause melting of the upper crust rocks. This leads to the formation of a shallow crustal reservoir that is directly responsible for the supereruptions.

Seismic tomography of the area of the 2010 Beni-Ilmane earthquake sequence, north-central Algeria.

The region of Beni-Ilmane (District of M'sila, north-central Algeria) was the site of an earthquake sequence that started on 14 May 2010. This sequence, which lasted several months, was triggered by conjugate E-W reverse and N-S dextral faulting. To image the crustal structure of these active faults, we used a set of 1406 well located aftershocks events and applied the local tomography software (LOTOS) algorithm, which includes absolute source location, optimization of the initial 1D velocity model, and iterative tomographic inversion for 3D seismic P- and S-wave velocities (and the Vp/Vs ratio), and source parameters. The patterns of P-wave low-velocity anomalies correspond to the alignments of faults determined from geological evidence, and the P-wave high-velocity anomalies may represent rigid blocks of the upper crust that are not deformed by regional stresses. The S-wave low-velocity anomalies coincide with the aftershock area, where relatively high values of Vp/Vs ratio (1.78) are observed compared with values in the surrounding areas (1.62-1.66). These high values may indicate high fluid contents in the aftershock area. These fluids could have been released from deeper levels by fault movements during earthquakes and migrated rapidly upwards. This hypothesis is supported by vertical sections across the study area show that the major Vp/Vs anomalies are located above the seismicity clusters.

4D Arctic: A Glimpse into the Structure and Evolution of the Arctic in the Light of New Geophysical Maps, Plate Tectonics and Tomographic Models.

Knowledge about the Arctic tectonic structure has changed in the last decade as a large number of new datasets have been collected and systematized. Here, we review the most updated, publicly available Circum-Arctic digital compilations of magnetic and gravity data together with new models of the Arctic's crust. Available tomographic models have also been scrutinized and evaluated for their potential to reveal the deeper structure of the Arctic region. Although the age and opening mechanisms of the Amerasia Basin are still difficult to establish in detail, interpreted subducted slabs that reside in the High Arctic's lower mantle point to one or two episodes of subduction that consumed crust of possibly Late Cretaceous-Jurassic age. The origin of major igneous activity during the Cretaceous in the central Arctic (the Alpha-Mendeleev Ridge) and in the proximity of rifted margins (the so-called High Arctic Large Igneous Province-HALIP) is still debated. Models of global plate circuits and the connection with the deep mantle are used here to re-evaluate a possible link between Arctic volcanism and mantle plumes.