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 km3 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.

We present a three‐dimensional density model of the lithosphere and upper mantle for the Middle East and surroundings based on seismic, gravity, and seismic tomography data and analyze the main factors responsible for the density variations. The gravity effect of the crust is calculated and removed from the observed field using the most recent crustal model. The residual gravity anomalies are jointly inverted with the residual topography to image the density distribution within the upper mantle. The inversion is constrained by an initial density model based on seismic tomography. The obtained density variations span in a large range (±60 kg/m3), revealing strong asymmetry in the density structure of the Arabian plate. The uppermost mantle layer in the Arabian Shield is relatively dense. However, below a depth of ~100 km we observe a strong low‐density anomaly. In contrast, the mantle density in the Arabian platform increases at the same depths. The most pronounced decrease of the mantle density occurs in the Gulf of Aden, Red Sea, and East African Rift. Underneath the northern Red Sea the low‐density anomaly is limited to the depth ~150 km, while in the southern part it is likely linked to a mantle plume. The densest mantle material is found under the South Caspian basin, which is likely associated with an eclogite body in the uppermost mantle. In the collision zones (the Zagros Belt and the Hellenic Arc), the high‐density lithosphere shows the location of the subducting plates.

We construct a new-generation 3D density model of the upper mantle of Asia and its surrounding areas based on a joint interpretation of several data sets. A recent model of the crust combining nearly all available seismic data is employed to calculate the impact of the crust on the gravity anomalies and observed topography and to estimate the residual mantle anomalies and residual topography. These fields are jointly inverted to calculate the density variations in the lithosphere and upper mantle down to 325 km.As an initial approximation, we estimate density variations using a seismic tomography model. Seismic velocity variations are converted into temperatures and then to density variations based on mineral physics constraints. In the Occam-type inversion, we fit both the residual mantle gravity anomalies and residual topography by finding deviations to the initial model. The obtained corrections improve the resolution of the initial model and reflect important features of the mantle structure that are not well resolved by the seismic tomography. The most significant negative corrections of the upper mantle density, found in the Siberian and East European cratons, can be associated with depleted mantle material. The most pronounced positive density anomalies are found beneath the Tarim and South Caspian basins, Barents Sea, and Bay of Bengal. We attribute these anomalies to eclogites in the uppermost mantle, which have substantially affected the evolution of the basins. Furthermore, the obtained results provide evidence for the presence of eclogites in the oceanic subducting mantle lithosphere.

The Klyuchevskoy group of volcanoes (KGV) in Kamchatka includes three presently active volcanoes (Klyuchevskoy, Bezymianny, and Tolbachik) located close together in an area of approximately 50 × 80 km. These three volcanoes have completely different compositions and eruption styles from each other. We have analyzed new data recorded by a temporary seismic network consisting of 22 seismic stations operated within the area of Tolbachik in 2014–2015 in conjunction with the data from the permanent network and the temporary PIRE network deployed at the Bezymianny volcano in 2009. The arrival times of the P and S waves were inverted using a local earthquake tomography algorithm to derive 3‐D seismic models of the crust beneath the KGV as well as accurate seismicity locations. High‐resolution structures beneath the Tolbachik volcanic complex were identified for the first time in this study. The tomography results reveal three different types of feeding system for the main KGV volcanoes. The basaltic lavas of the Klyuchevskoy volcano are supplied directly from a reservoir at a depth of 25–30 km through a nearly vertical pipe‐shaped conduit. The explosive Bezymianny volcano is fed through a dispersed system of crustal reservoirs where a lighter felsic material separates from the mafic component and ascends to the upper crust to form andesitic magma sources. For Tolbachik, low‐viscosity volatile‐saturated basalts ascend from two deep reservoirs following a system of fractures in the crust associated with the intersections of regional faults.

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.