S U M M A R YWe present a new tomographic model for P-and S-velocity anomalies beneath Europe (30 • N-55 • N, 5 • W-40 • E), extending in depth up to 700 km and constrained by inversion of data from the International Seismological Center (ISC) catalogue. The algorithm uses the traveltimes from events located in the study area recorded by all available worldwide stations, as well as times from teleseismic events recorded by European stations. The events from the ISC catalogue have been relocated and combined into composite events. All the traveltimes were corrected for crustal structure using the reference model EuCRUST-07. The resulting velocity anomalies show similar large-scale patterns as observed in previous studies, but have a higher resolution, which allows detection of some features in more detail. For example, it is now possible to assess the depth extension of the small slow velocity body beneath the Eifel region and Eger graben. The P and the S model show a good consistency in the uppermost 200 km below most of the European area and in some parts even in the deeper layers (e.g. beneath the Apennines and the Hellenic arc). The new model provides clear images of some principal features, which were previously detected in a limited number of studies, while the comparison between P-and S-velocity anomalies provides novel constraints to address on their nature (e.g. the gap in the Adriatic plate subducted below the central-southern Apennines). In this paper, we pay special attention to testing the reliability of the results. The random noise effect is evaluated using a test with independent inversion of two data subsets (with odd/even events). The spatial resolution is estimated using different checkerboard tests. Furthermore, we present a synthetic model with realistic patterns, which reproduces after performing forward and inverse modelling the same shape and amplitudes of the anomalies as in the case of the real data inversion. In this case, the parameters of the model can be used to assess the amplitudes of P and S anomalies that is critical for evaluation of other petrophysical parameters (temperature, density, composition, etc.) in the upper mantle.

GF is high and melt water is present under ice cover [11][12] Greenland to explain the origin of the observed melting beneath the ice cover (Figure 1). This are controlled by a combination of GF and non-GF influences, we build our calibration 137 strategy on estimating GF required to reproduce the observed thawed basal ice conditions, 138 discounting basal ice melt rates as a proxy for GF. This has the effect that GF estimates will 139 likely be biased downwards where basal melt is rapid; nevertheless, our strategy is 140 sufficiently effective to separate out the signal of a strong and spatially extensive geothermal 141 anomaly beneath the GIS and provides a hard lower bound for GF values at the observed 142 basal melt locations. 143The anomalous GF zone lies in the area with the highest density of direct measurements. 150One potential cause of elevated GF is illustrated by seismic data that link our west-to-east GF 151anomaly with a zone of low-seismic-velocity mantle, a "negative anomaly", beneath Iceland 6- Greenland may be the expression of Iceland hotspot history. The geothermal anomaly 237 provides evidence for a more northerly hotspot track than previously proposed and will offer 238 a useful test for existing paleoreconstructions of absolute plate motion. This study advocates 239 a previously undocumented strong coupling between Greenland's present-day ice dynamics, 240 subglacial hydrology, and the remote tectonothermal history of the North Atlantic region.

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.

S U M M A R Y∼82 000 P and S arrival times from ∼3000 sources recorded by ∼250 seismic stations from the revised ISC catalogue are employed to study a circular area of 6 • radius centred on the Dead Sea. We use the linearized tomographic approach based on the rays constructed in a 1-D spherical velocity model and corrected for the Moho depth variation and relief. All the sources were relocated. As the result of simultaneous iterative inversion we get 3-D P and S velocity anomalies in the crust and uppermost mantle, Moho depth and corrected source parameters. The resulting images fit well with the existing tectonic elements in the study area. In the crust, a narrow P and S low-velocity anomaly marks the position of the Dead Sea Transform (DST) that is interpreted as sediments in the shallower layer and a zone of fractured and deformed rocks in the middle and lower crust. There is a narrow (50 km wide) band of thickening of the crust along the DST in the Arava valley between the Red Sea and Dead Sea and some 100 km north of the Dead Sea. This zone may be associated with the minimum of the lithospheric strength and, therefore, explain the location of the DST in the Arava Valley. The velocity anomalies under the crust and the map of the Moho depth clearly distinguish the oceanic (Levant basin) and continental types of crust (Asia Minor, Zagros, Cyprus and Eratosthenes Mount). Verification of the results takes an important part in this study. Inversions of different starting models and independent processing of data subsets show high robustness of the results. Synthetic tests clearly show the limits of the resolving power of the inversion with the existing data set.

[1] In this study we present the new tomographic code ANITA which provides 3-D anisotropic P and isotropic S velocity distribution based on P and S traveltimes from local seismicity. For the P anisotropic model, we determine four parameters for each parameterization cell. This represents an orthorhombic anisotropy with one predefined direction oriented vertically. 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 nonlinear iterative inversion procedure is similar to that used in the LOTOS code. We have implemented this algorithm for the updated data set of central Java, part of which was previously used for the isotropic inversion. It was obtained that the crustal and uppermost mantle velocity structure beneath central Java is strongly anisotropic with 7-10% of maximal difference between slow and fast velocity in different directions. In the forearc (area between southern coast and volcanoes), the structure of both isotropic and anisotropic structure is strongly heterogeneous. Variety of anisotropy orientations and highly contrasted velocity patterns can be explained by a complex block structure of the crust. Beneath volcanoes we observe faster velocities in vertical direction, which is probably an indicator for vertically oriented structures (channels, dykes). In the crust beneath the middle part of central Java, north to Merapi and Lawu volcanoes, we observe a large and very intense anomaly with a velocity decrease of up to 30% and 35% for P and S models, respectively. Inside this anomaly E-W orientation of fast velocity takes place, probably caused by regional extension stress regime. In a vertical section we observe faster horizontal velocities inside this anomaly that might be explained by layering of sediments and/or penetration of quasi-horizontal lenses with molten magma. In the mantle, trench parallel anisotropy is observed throughout the study area. Such anisotropy in the slab entrained corner flow may be due to presence of B-type olivine having predominant axis parallel to the shear direction, which appears in conditions of high water or/and melting content.

SUMMARY Seismic and volcanic activities in Central Java, Indonesia, the area of interest of this study, are directly or indirectly related to the subduction of the Indo‐Australian plate. In the framework of the MERapi AMphibious EXperiments (MERAMEX), a network consisting of about 130 seismographic stations was installed onshore and offshore in Central Java and operated for more than 150 days. In addition, 3‐D active seismic experiments were carried out offshore. In this paper, we present the results of processing combined active and passive seismic data, which contain traveltimes from 292 local earthquakes and additional airgun shots along three offshore profiles. The inversion was performed using the updated LOTOS‐06 code that allows processing for active and passive source data. The joint inversion of the active and passive data set considerably improves the resolution of the upper crust, especially in the offshore area in comparison to only passive data. The inversion results are verified using a series of synthetic tests. The resulting images show an exceptionally strong low‐velocity anomaly (−30 per cent) in the backarc crust northward of the active volcanoes. In the upper mantle beneath the volcanoes, we observe a low‐velocity anomaly inclined towards the slab, which probably reflects the paths of fluids and partially melted materials in the mantle wedge. The crust in the forearc appears to be strongly heterogeneous. The onshore part consists of two high‐velocity blocks separated by a narrow low‐velocity anomaly, which can be interpreted as a weakened contact zone between two rigid crustal bodies. The recent Java Mw= 6.3 earthquake (2006/05/26‐UTC) occurred at the lower edge of this zone. Its focal strike slip mechanism is consistent with the orientation of this contact.

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.

No abstract