Deep crustal earthquakes in North Tanzania, East Africa: Interplay between tectonic and magmatic processes in an incipient rift

International audienceIn this study, we explore the origin of lower crustal seismicity and the factors controlling rift propagation using seismological data recorded within the youngest part of the East African Rift System, the North Tanzanian Divergence (NTD). Most earthquakes below Lake Manyara occur at depth ranging between 20 and 40 km and have a swarm-like distribution. Focal mechanisms of 26 events indicate a combination of strike-slip and normal faulting involving Archaean basement structures and forming a relay zone. The derived local stress regime is transtensive and the minimum principal stress is oriented N110°E. Crustal seismic tomography reveals low-velocity anomalies below the rifted basins in the NTD, interpreted as localized thermomechanical perturbations promoting fluid release and subsequent seismicity in the lower crust. SKS splitting analysis in the NTD indicates seismic anisotropy beneath 17 stations most likely due to aligned magma lenses and/or dikes beneath the rift and to the lithospheric fabrics. Our results favor a strain pattern intermediate between purely mechanical and purely magmatic. We suggest that melt products arising from a large asthenospheric thermal anomaly enhance lithospheric weakening and facilitate faulting and creeping on critically oriented inherited structures of the Precambrian lower crust. Although the crust is unlikely weakened at a point comparable to other parts of the East African Rift System, this deep-seated thermomechanical process is efficient enough to allow slow rift propagation within the eastern Tanzanian cratonic edge

The crustal and upper mantle isotropic and anisotropic velocity structure beneath the Rwenzori Mountains from travel-time tomography

In this study, I investigate the crustal and upper mantle velocity structure beneath the Rwenzori Mountains in western Uganda. This mountain range is situated within the western branch of the East African Rift and reaches altitudes of more than 5000 m. I use four different approaches that belong to the travel-time tomography method. The first approach is based on the isotropic tomographic inversion of local data, which contain information about 2053 earthquakes recorded by a network of up to 35 stations covering an area of 140×90 km2. The LOTOS-09 algorithm described here is used to realize this approach. The second approach is based on the anisotropic tomographic inversion of the same local dataset. This method employs the tomographic code ANITA, developed with my participation, which provides 3D anisotropic P and isotropic S velocity distributions based on P and S travel-times from local seismicity. For the P anisotropic model, four parameters for each parameterization cell are determined. This represents an orthorhombic anisotropy with one vertically-oriented predefined direction. Three of the parameters describe slowness variations along three horizontal orientations with azimuths of 0°, 60° and 120°, and one is a perturbation along the vertical axis. The third approach is based on tomographic inversion of the teleseismic data, which contain information about the traveltimes of P-waves coming from 284 teleseismic events recorded by the seismic network stations. The TELELOTOS code, which is my own modification of the LOTOS-09 algorithm, is used in this approach. The TELELOTOS code is designed to iteratively invert the local and/or teleseismic datasets. Finally, I present the results of the new tomographic approach, which is based on the simultaneous inversion of the joint local and teleseismic data. The simultaneous use of these datasets for the tomographic inversion has several advantages. In this case, the velocity structure in the study area can be resolved as deep as in the teleseismic approach. At the same time, in the upper part of the study volume, the resolution of the obtained models is as good as in the local tomography. The TELELOTOS algorithm is used to perform the joint tomographic inversion. Special attention is paid in this work to synthetic testing. A number of different synthetic and real data tests are performed to estimate the resolution ability and robustness of the obtained models. In particular, synthetic tests have shown that the results of the anisotropic tomographic inversion of the local data have to be considered as unsatisfactory. For all approaches used in this study, I present synthetic models that reproduce the same pattern of anomalies as that obtained by inverting the real data. These models are used to interpret the results and estimate the real amplitudes of the obtained anomalies. The obtained models exhibit a relatively strong negative P anomaly (up to -10%) beneath the Rwenzori Mountains. Low velocities are found in the northeastern part of the array at shallower depths and are most likely related to sedimentary deposits, while higher velocities are found beneath the eastern rift shoulder and are thought to be related to old cratonic crust. The presence of low velocities in the northwestern part of the array may be caused by a magmatic intrusion beneath the Buranga hot springs. Relatively low velocities were observed within the lower crust and upper mantle in the western and southern parts of the study area (beneath the rift valley and the entire length of the Rwenzori range). The higher amplitude of the low-velocity anomaly in the south can be related to the thinner lithosphere in the southern part of the Albertine rift. In the center of the study area, a small negative anomaly is observed, with the intensity increasing with depth. This anomaly is presumably related to a fluids rising up from a plume branch in the deeper part of the mantle. According to the interpretation of the local earthquake distribution, the Rwenzori Mountains are located between two rift valleys with flanks marked by normal faults. The Rwenzori block is bounded by thrust faults that are probably due to compression.In der vorliegenden Arbeit untersuche ich die seismischen Geschwindigkeiten der Erdkruste und des oberen Mantels unter dem Rwenzori Gebirge im westlichen Uganda. Das Untersuchungsgebiet befindet sich im Albert Rift, einem Teil des westlichen Astes des Ostafrikanischen Grabensystems (EARS, East African Rift System). Im Gegensatz zu anderen Gebirgen im Bereich des EARS, liegen die Rwenzori Berge direkt im Riftgraben. Sie erreichen eine Höhe von über 5000 m und werden aus präkambrischen Gesteinen gebildet. Die Ursachen für die Hebung eines solchen Gebirgsmassivs innerhalb eines Riftgrabens mit überwiegend Extensionstektonik werden vielfach diskutiert. Ein Vorschlag sieht die Ursache der Gebirgsbildung in der Kollision zweier rotierender Krustenblöcke. Eine weitere, meiner Meinung nach wahrscheinlichere Theorie besagt, dass die Hebung des Rwenzori durch eine Delamination der Lithosphäre zwischen zwei benachbarten Rift Segmenten hervorgerufen wird. Die Rwenzori Berge liegen in einem der seismisch aktivsten Gebiete des Ostafrikanischen Grabensystems. In meiner Untersuchung verwende ich die Daten eines temporären seismischen Netzwerkes, das im Rahmen des RiftLink-Projektes und des BGR-GEOTHERM-Projektes im Gebiet der Rwenzoris installiert worden war von 2006 bis 2007. Es bestand aus insgesamt 35 Stationen, die über ein Gebiet von ungefähr 140×90 km2 verteilt waren. Im Nordteil des Netzwerkes variieren die Stationsabstände von 3 bis 12 km und im Süden von 20 bis 50 km. Für die Analyse der seismischen Registrierungen und die vorläufige Lokalisierung der lokalen Beben wurde das SEISAN Programmpaket verwendet. Ich präsentiere die Ergebnisse von vier tomographischen Verfahren zu Untersuchung der Geschwindigkeitsstruktur in Kruste und oberem Erdmantel unterhalb des Rwenzori Gebirges. Alle Verfahren basieren auf der Auswertung von natürlichen Quellen, wobei besondere Aufmerksamkeit auf synthetische Testmethoden gelegt wurde. Zunächst werden Daten von 2053 lokalen Erdbeben für eine tomographische Inversion mit dem LOTOS-09 Algorithmus verwendet (Koulakov 2009). Hierbei gehe ich zunächst von isotropen Geschwindigkeitsvariationen aus. Der Algorithmus basiert im Wesentlichen auf drei Schritten: (i) simultane Optimierung des 1D Geschwindigkeitsmodells und vorläufige Lokalisierung der Beben, (ii) Lokalisierung der Beben mit einem 3D Geschwindigkeitsmodell und (iii) simultane Inversion der Herdparameter und des Geschwindigkeitsmodells mit verschiedenen Modellgitterorientierungen. Die Schritte (ii) und (iii) werden iterativ durchlaufen. Um den Einfluss der Startwerte von Bebenverteilung und 1D Geschwindigkeitsmodell auf die Ergebnisse der tomographischen Inversion zu überprüfen, führe ich mehrere Checkerboard-Tests mit synthetischen 3D Anomalien durch. Diese Tests zeigen, dass die Geschwindigkeitsstruktur im nördlichen Teil der Rwenzories, wo die Stationsdichte am höchsten ist, zuverlässig bis in eine Tiefe von etwa 15 km aufgelöst werden kann. Im Südosten des Stationsnetzes ist das räumliche Auflösungsvermögen aufgrund horizontaler und vertikaler Schmiereffekte am schlechtesten. Die Auflösung des abgeleiteten Geschwindigkeitsmodells ist in horizontaler Richtung höher als in vertikaler Richtung. Wie die synthetischen Tests auch zeigen, beträgt die durchschnittliche Genauigkeit der Bebenlokalisierung etwa 1 bis 1.5 km. Bei Beben, deren Epizentren außerhalb des Stationsnetzes liegen, sind die Lokalisierungsfehler kleiner als 3 km, lediglich bei 49 Beben steigt dieser Wert auf über 5km. Die Inversion von verschiedenen Beben-Untergruppen sowie des kompletten Datensatzes führt zu guten Übereinstimmung der berechneten P- und S-Modelle. Dies bedeutet, dass die Ergebnisse stabil sind und nicht wesentlich von Zufallsfaktoren (wie etwa Fehler beim Ablesen von P- und SEinsätzen) beeinflusst werden. Im zweiten Teil wird der gleiche Datensatz unter Berücksichtigung anisotroper Geschwindigkeitsverteilung invertiert. Hierfür wird der von mir mitentwickelte ANITA Algorithmus verwendet (Koulakov et. al. 2009). Das Programm berechnet die anisotrope P- und die isotrope S-Geschwindigkeitsverteilung aus den beobachteten lokalen P- und S-Laufzeiten. Für die P-Geschwindigkeit wird orthorombische Anisotropie mit einer vorgegebenen Achse in vertikaler Richtung zugrunde gelegt. Das Programm berechnet für jede Modellzelle vier Parameter. Drei davon beschreiben die Slowness Variationen entlang horizontaler Achsen mit Azimuten von 0°, 60° und 120°. Der vierte Parameter entspricht der Perturbation der Slowness in vertikaler Richtung. Der Algorithmus basiert wieder im Wesentlichen auf drei Schritten: (i) vorläufige Lokalisierung der Beben mit einem 1D Geschwindigkeitsmodell, (ii) Lokalisierung der Beben mit einem 3D Geschwindigkeitsmodell und (iii) simultane Inversion der Herdparameter des anisotropen P-Geschwindigkeitsmodells und des isotropen SGeschwindigkeitsmodells mit verschiedenen Modellgitterorientierungen. Die Schritte (ii) und (iii) werden iterativ durchlaufen. Mehrere Tests mit synthetischen und realen Daten wurden durchgeführt, um die Zuverlässigkeit und Stabilität der berechneten Modelle zu untersuchen. Es zeigt sich, dass das anisotrope Modell nur im Nordwesten, zwischen 0.4° und 0.9° Breite sowie zwischen 30.0° und 30.3° Länge, eine mehr oder weniger gute Auflösung besitzt. Das Gebiet, in dem die isotrope Komponente des anisotropen Modells eine gute Auflösung besitzt, ist viel größer und entspricht den Ergebnissen der isotropen Inversion. Die synthetischen Tests zeigen allerdings einen kräftigen trade-off Effekt zwischen den isotropen und den anisotropen Parametern. Modelle, die aus der Inversion von isotropen Modellen berechnet wurden, besitzen im gesamten Messgebiet eine kräftige, zufällig orientierte, anisotrope Komponente. Andererseits hat ein ursprünglich anisotropes Modell bei der Inversion offensichtlich keinen Einfluss auf die isotropen Komponenten. Die Ergebnisse von gerade/ungerade Tests belegen, dass die anisotrope Komponente des berechneten Modells sehr stark vom Rauschen und der Strahlgeometrie beeinflusst wird. Die Ergebnisse variieren deshalb in Abhängigkeit von den Datenuntergruppen, die bei der Inversion benutzt wurden. Die isotrope Inversion des gleichen Datensatzes mit dem ANITA Algorithmus ergibt kleinere P-Residuen als bei der Inversion mit anisotropem P- und isotropem S-Modell. Dies kann allerdings damit erklärt werden, dass im zweiten Fall eine größere Anzahl von Parametern in die Inversion eingeht. Diese Überlegung wird dadurch gestützt, dass die berechneten Geschwindigkeitsverteilungen bei einer rein isotropen Inversion einen glätteren Verlauf aufweisen. Aufgrund dieser Tatsachen komme ich zu dem Schluss, dass eine anisotrope Inversion des bestehenden Datensatzes nicht zu zufriedenstellenden Ergebnissen führt. Es sollte allerdings betont werden, dass die isotropen Anteile des berechneten anisotropen 3D Modells stabil sind und gut mit den Resultaten der isotropen Tomographie übereinstimmen. Im dritten Teil werden die Laufzeiten teleseismischer Beben ausgewertet, um die P-Geschwindigkeitsstruktur in der Kruste und im oberen Mantel zu untersuchen. Dabei habe ich selbst die seismischen Phasen gepickt und alle Stufen der Datenanalyse durchgeführt. Zunächst war es notwendig, die teleseismischen Ereignisse in den kontinuierlichen Datenstreams zu identifizieren und zu extrahieren. Dazu wurden alle Beben aus dem ISC-Katalog gewählt, auf die folgende Bedingungen zutrafen: (i) Magnitude größer als 5.0 und (ii) Epizentralentfernung zum Zentrum des Messgebiets (30.2°O, 0.4°N) zwischen 20° und 95°. Danach wurden die theoretischen Ankunftszeiten der P-Wellen berechnet. Schließlich wurden nur solche Beben ausgesucht, die in den Seismogrammen identifiziert werden. Die Korrelation der Seismogramme aller Stationen, die ein Ereignis registriert hatten, diente zur sicheren Identifizierung der Phaseneinsätze. Insgesamt gingen 3159 P-Phasen von 284 Fernbeben in die tomographische Inversion ein. Die Berechnungen erfolgten mit dem TELELOTOS Code, den ich aus dem LOTOS-09 Algorithmus weiterentwickelt habe. Mit diesem Programm ist es möglich, lokale und teleseismische Daten kombiniert zu invertieren. Solange nur Informationen von Lokalbeben vorhanden sind, läuft das Programm im Wesentlichen wie der LOTOS-09 Algorithmus. Bei einer kombinierten Auswertung von lokalen und teleseismischen Daten werden die Lokalbeben nur benutzt, um das 1D-Modell zu optimieren. Gleichzeitig werden sie relokalisiert, während die Koordinaten der Fernbeben während des Inversionsprozesses fixiert bleiben. Werden nur teleseismische Daten für die Analyse benutzt, erfolgt weder eine Optimierung des 1D-Modells noch werden die Beben noch werden mit einem 3D-Geschwindigkeitsmodell relokalisiert. Tests mit synthetischen und realen Daten haben gezeigt, dass das invertierte 3DGeschwindigkeitsmodell von Kruste und oberem Erdmantel unterhalb der Rwenzoris robust ist und von Zufallsfaktoren nicht signifikant beeinflusst wird. Im letzten Teil meiner Untersuchung präsentiere ich die Ergebnisse dieses neuen Verfahrens. Die kombinierte tomographische Inversion von lokalen und teleseismischen Daten mit dem TELELOTOS Algorithmus hat den Vorteil, dass das invertierte Geschwindigkeitsmodell einerseits die Auflösung der teleseismischen Tomographie besitzt, andererseits aber in geringen Tiefen die bessere Auflösung der lokalen Tomographie. Mehrere synthetische Tests wurden durchgeführt, um das Auflösungsvermögen des seismischen Netzwerks in unterschiedlichen Bereichen des Messgebiets zu untersuchen. Die Ergebnisse zeigen eine gute Übereinstimmung mit den Resultaten der isotropen Lokalbebentomographie und der teleseismischen Tomographie, die in dieser Arbeit vorgestellt wurden und können wie folgt charakterisiert werden: Im Westen und Süden des Messgebiets dominieren negative Geschwindigkeitsanomalien während im Osten erhöhte Geschwindigkeiten beobachtet werden. Dies gilt für alle untersuchten Tiefenbereiche. Im Nordwesten des Messgebietes wurde eine kräftige negative Anomalie mit extrem hohen Vp/Vs Verhältnissen lokalisiert, an deren nördlichem Rand die heißen Quellen von Buranga liegen. Geochemische Untersuchungen lassen hier einen magmatischen Körper vermuten, der die Wärmequelle dieser Quellen darstellen könnte (Ochmann 2007). Weitere Zonen erniedrigter Geschwindigkeit befinden sich in geringer Tiefe nahe bei einigen östlichen Stationen und sind möglicherweise Gestainen der Paleoproterozoischen Einhaiten zuzuordnen. In der Zentral- und Südenregionen der Rwenzoris beobachten wir zwei relativ kleine Zonen erhöhter Geschwindigkeit, die das Gebirge überquert und zu teilen scheint. Interessanterweise ergeben die Untersuchungen relativ niedrige Geschwindigkeiten direkt fast unterhalb des Rwenzori Massivs. Im Wesentlichen befinden sich hier Gneise, Schiefer und Amphibolite Gesteinstypen (Shlüter 2008), denen üblicherweise erhöhte seismische Geschwindigkeiten zugeordnet werden. Es ist zu vermuten, dass es durch intensive Faltung und die rasche Hebung des Rwenzori Gebirges zu einer starken Scherung und Verzerrung des Gesteinsgefüges kam, was wiederum zu reduzierten seismischen Geschwindigkeiten führte. Es ist bemerkenswert, dass die niedrigen Geschwindigkeitsgebiete an der geringen Tiefe unter den Gebirgen gut Korrelat mit dem Gesteinsausbiss der Paleoproterozoischen Einheiten sind. Es ist allerdings zu beachten, dass sich das Gebirge außerhalb des Stationsnetzes befindet, in einem Gebiet, in dem das schlechte Auflösungsvermögen eine detaillierte Aussage über die Geschwindigkeitsstruktur verhindert. Im Westen und Süden des Messgebietes werden in Unterkruste und oberem Erdmantel relativ niedrige Geschwindigkeiten beobachtet, die sich unterhalb des Riftgrabens über die gesamte Länge der Rwenzori Region verteilen. Die stärkeren negativen Anomalien im Süden könnten mit einer ausgedünnten Lithosphäre in diesem Bereich des Albert Rifts erklärt werden. Erhöhte Geschwindigkeiten befinden sich unterhalb der östlichen Riftschulter im Bereich alter kratonischer Kruste. Im Zentrum des Messgebietes konnte eine negative Anomalie lokalisiert werden, deren Intensität mit zunehmender Tiefe anwächst. Diese Anomalie liegt möglicherweise im Bereich von aufsteigendem heißem Material, das von einem großen Mantel-Plume unterhalb des EARS herrührt (e.g., Nyblade et. al. 2000, Park and Nyblade 2006, Koulakov 2007). Ein weiteres Indiz dafür sind die Vulkanfelder der Ankole Provinz, die sich in der Nähe der beobachteten Anomalie befinden (Boven et. al. 1998). Die Verteilung der Beben in der Rwenzori Region stellt eine wichtige Information über tektonische Prozesse in der Erdkruste dar. Die Position der relokalisierten Bebenherde zeigt, dass sich die Rwenzoris zwischen zwei Riftgräben befinden, deren Ränder durch Abschiebungsstörungen markiert werden. Der Gebirgsblock selbst wird begrenzt von Aufschiebungen, die wahrscheinlich die Folge von Kompressionskräften sind

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

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

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

Magmatic and Sedimentary Structure beneath the Klyuchevskoy Volcanic Group, Kamchatka, From Ambient Noise Tomography

The Klyuchevskoy Volcanic Group is a cluster of the world's most active subduction volcanoes, situated on the Kamchatka Peninsula, Russia. The volcanoes lie in an unusual off‐arc position within the Central Kamchatka Depression (CKD), a large sedimentary basin whose origin is not fully understood. Many gaps also remain in the knowledge of the crustal magmatic plumbing system of these volcanoes. We conducted an ambient noise surface wave tomography, to image the 3‐D shear wave velocity structure of the Klyuchevskoy Volcanic Group and CKD within the surrounding region. Vertical component cross correlations of the continuous seismic noise are used to measure interstation Rayleigh wave group and phase traveltimes. We perform a two‐step surface wave tomography to model the 3‐D Vsv velocity structure. For each inversion stage we use a transdimensional Bayesian Monte Carlo approach, with coupled uncertainty propagation. This ensures that our model provides a reliable 3‐D velocity image of the upper 15 km of the crust, as well as a robust assessment of the uncertainty in the observed structure. Beneath the active volcanoes, we image small slow velocity anomalies at depths of 2–5 km but find no evidence for magma storage regions deeper than 5 km—noting the 15 km depth limit of the model. We also map two clearly defined sedimentary layers within the CKD, revealing an extensive 8 km deep sedimentary accumulation. This volume of sediments is consistent with the possibility that the CKD was formed as an Eocene‐Pliocene fore‐arc regime, rather than by recent (<2 Ma) back‐arc extension.Plain Language Summary: The Klyuchevskoy Volcanic Group is a cluster of 13 volcanoes on the Kamchatkan corner of the Pacific ring of fire. The volcanoes regularly produce large eruptions, but good knowledge of the magma plumbing system beneath the surface is still lacking. Why the Klyuchevskoy Volcanic Group volcanoes lie in the location they do, in a large low‐lying depression, is also unexplained. We undertook a seismic experiment and used the data to produce a 3‐D velocity image of the subsurface beneath the volcanoes and the depression. We found that small regions of slow seismic velocity are located beneath the active volcanoes, at 2–5 km depth below sea level. This slower velocity is probably caused by magma lying within the porous fracture spaces in this rock. The seismic velocities are much faster beneath the dormant volcanoes, suggesting they have no magma beneath them. With our velocity image, we also find that the Central Kamchatka Depression is very deep, filled with over 8 km of sediments. This supports an idea that the sediments accumulated as a fore‐arc basin over many millions of years, since 40 Ma, when the active line of volcanoes was found 100 km to the west.Key Points: Three‐dimensional shear velocity structure of the Klyuchevskoy area was determined using coupled transdimensional Monte Carlo inversions. Slow velocity anomalies suggest magma storage beneath active volcanoes at 2–5 km depth (below sea level) but not in the midcrust. Sediments filling the Central Kamchatka Depression are 8 km deep, consistent with an origin of the depression as a fore‐arc basin.European Union Horizon 2020 Research and Innovation Programme http://dx.doi.org/10.13039/501100007601Russian Ministry of Education and Science http://dx.doi.org/10.13039/501100003443Alexander von Humboldt Foundation http://dx.doi.org/10.13039/10000515

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

Growth of mountain belts in central Asia triggers a new collision zone in central India

Several unusual strong earthquakes occurred in central India along the Narmada-Son Lineament (NSL) zone, far from active plate boundaries. To understand the role of collisional processes in the origin of this seismicity, we develop a numerical thermomechanical model of shortening between the Indian Plate and Asia. We show that at the final stage of collision, the shortening rate of the high mountain areas slows. The continuing convergence of India and Asia triggers the initiation of a new collision zone in continental part of India. Various geological and geophysical observations indicate that the NSL is a weakest zone with northward thrusting of the thinner central Indian lithosphere underneath the thicker northern part of the Indian Plate. We hypothesize that the NSL was reactivated during the final stage of the India Asia convergence and it will possibly form a new mountain belt within the Indian continent.ISSN:2045-232