Arrival angles of teleseismic fundamental mode Rayleigh waves across the AlpArray

The dense AlpArray network allows studying seismic wave propagation with high spatial resolution. Here we introduce an array approach to measure arrival angles of teleseismic Rayleigh waves. The approach combines the advantages of phase correlation as in the two-station method with array beamforming to obtain the phase-velocity vector. 20 earthquakes from the first two years of the AlpArray project are selected, and spatial patterns of arrival-angle deviations across the AlpArray are shown in maps, depending on period and earthquake location. The cause of these intriguing spatial patterns is discussed. A simple wave-propagation modelling example using an isolated anomaly and a Gaussian beam solution suggests that much of the complexity can be explained as a result of wave interference after passing a structural anomaly along the wave paths. This indicates that arrival-angle information constitutes useful additional information on the Earth structure, beyond what is currently used in inversions

Ambient-noise tomography of the wider Vienna Basin region

We present a new 3-D shear-velocity model for the top 30 km of the crust in the wider Vienna Basin region based on surface waves extracted from ambient-noise cross-correlations. We use continuous seismic records of 63 broad-band stations of the AlpArray project to retrieve interstation Green’s functions from ambient-noise cross-correlations in the period range from 5 to 25 s. From these Green’s functions, we measure Rayleigh group traveltimes, utilizing all four components of the cross-correlation tensor, which are associated with Rayleigh waves (ZZ, RR, RZ and ZR), to exploit multiple measurements per station pair. A set of selection criteria is applied to ensure that we use high-quality recordings of fundamental Rayleigh modes. We regionalize the interstation group velocities in a 5 km × 5 km grid with an average path density of ∼20 paths per cell. From the resulting group-velocity maps, we extract local 1-D dispersion curves for each cell and invert all cells independently to retrieve the crustal shear-velocity structure of the study area. The resulting model provides a previously unachieved lateral resolution of seismic velocities in the region of ∼15 km. As major features, we image the Vienna Basin and Little Hungarian Plain as low-velocity anomalies, and the Bohemian Massif with high velocities. The edges of these features are marked with prominent velocity contrasts correlated with faults, such as the Alpine Front and Vienna Basin transfer fault system. The observed structures correlate well with surface geology, gravitational anomalies and the few known crystalline basement depths from boreholes. For depths larger than those reached by boreholes, the new model allows new insight into the complex structure of the Vienna Basin and surrounding areas, including deep low-velocity zones, which we image with previously unachieved detail. This model may be used in the future to interpret the deeper structures and tectonic evolution of the wider Vienna Basin region, evaluate natural resources, model wave propagation and improve earthquake locations, among others

Shear-wave velocity structure beneath the Dinarides from the inversion of Rayleigh-wave dispersion

Highlights • Rayleigh-wave phase velocity in the wider Dinarides region using the two-station method. • Uppermost mantle shear-wave velocity model of the Dinarides-Adriatic Sea region. • Velocity model reveals a robust high-velocity anomaly present under the whole Dinarides. • High-velocity anomaly reaches depth of 160 km in the northern Dinarides to more than 200 km under southern Dinarides. • New structural model incorporating delamination as one of the processes controlling the continental collision in the Dinarides. The interaction between the Adriatic microplate (Adria) and Eurasia is the main driving factor in the central Mediterranean tectonics. Their interplay has shaped the geodynamics of the whole region and formed several mountain belts including Alps, Dinarides and Apennines. Among these, Dinarides are the least investigated and little is known about the underlying geodynamic processes. There are numerous open questions about the current state of interaction between Adria and Eurasia under the Dinaric domain. One of the most interesting is the nature of lithospheric underthrusting of Adriatic plate, e.g. length of the slab or varying slab disposition along the orogen. Previous investigations have found a low-velocity zone in the uppermost mantle under the northern-central Dinarides which was interpreted as a slab gap. Conversely, several newer studies have indicated the presence of the continuous slab under the Dinarides with no trace of the low velocity zone. Thus, to investigate the Dinaric mantle structure further, we use regional-to-teleseismic surface-wave records from 98 seismic stations in the wider Dinarides region to create a 3D shear-wave velocity model. More precisely, a two-station method is used to extract Rayleigh-wave phase velocity while tomography and 1D inversion of the phase velocity are employed to map the depth dependent shear-wave velocity. Resulting velocity model reveals a robust high-velocity anomaly present under the whole Dinarides, reaching the depths of 160 km in the north to more than 200 km under southern Dinarides. These results do not agree with most of the previous investigations and show continuous underthrusting of the Adriatic lithosphere under Europe along the whole Dinaric region. The geometry of the down-going slab varies from the deeper slab in the north and south to the shallower underthrusting in the center. On-top of both north and south slabs there is a low-velocity wedge indicating lithospheric delamination which could explain the 200 km deep high-velocity body existing under the southern Dinarides

Crustal Thinning From Orogen to Back-Arc Basin: The Structure of the Pannonian Basin Region Revealed by P-to-S Converted Seismic Waves

We present the results of P-to-S receiver function analysis to improve the 3D image of the sedimentary layer, the upper crust, and lower crust in the Pannonian Basin area. The Pannonian Basin hosts deep sedimentary depocentres superimposed on a complex basement structure and it is surrounded by mountain belts. We processed waveforms from 221 three-component broadband seismological stations. As a result of the dense station coverage, we were able to achieve so far unprecedented spatial resolution in determining the velocity structure of the crust. We applied a three-fold quality control process; the first two being applied to the observed waveforms and the third to the calculated radial receiver functions. This work is the first comprehensive receiver function study of the entire region. To prepare the inversions, we performed station-wise H-Vp/Vs grid search, as well as Common Conversion Point migration. Our main focus was then the S-wave velocity structure of the area, which we determined by the Neighborhood Algorithm inversion method at each station, where data were sub-divided into back-azimuthal bundles based on similar Ps delay times. The 1D, nonlinear inversions provided the depth of the discontinuities, shear-wave velocities and Vp/Vs ratios of each layer per bundle, and we calculated uncertainty values for each of these parameters. We then developed a 3D interpolation method based on natural neighbor interpolation to obtain the 3D crustal structure from the local inversion results. We present the sedimentary thickness map, the first Conrad depth map and an improved, detailed Moho map, as well as the first upper and lower crustal thickness maps obtained from receiver function analysis. The velocity jump across the Conrad discontinuity is estimated at less than 0.2 km/s over most of the investigated area. We also compare the new Moho map from our approach to simple grid search results and prior knowledge from other techniques. Our Moho depth map presents local variations in the investigated area: the crust-mantle boundary is at 20–26 km beneath the sedimentary basins, while it is situated deeper below the Apuseni Mountains, Transdanubian and North Hungarian Ranges (28–33 km), and it is the deepest beneath the Eastern Alps and the Southern Carpathians (40–45 km). These values reflect well the Neogene evolution of the region, such as crustal thinning of the Pannonian Basin and orogenic thickening in the neighboring mountain belts

Evaluation of a Semi-Airborne Electromagnetic Survey Based on a Multicopter Aircraft System

The semi-airborne electromagnetic (EM) method has the potential to reach deeper exploration depths than purely airborne EM approaches. The concept of the method is to deploy high-power transmitters on the ground, which excite subsurface currents and induce strong magnetic fields, and to measure the corresponding EM fields with a passive airborne receiver instrument. Following recent conceptual developments of the semi-airborne EM technique deployed on helicopters, we performed a 10 km2 semi-airborne EM survey near Münster (Germany) based on a multicopter aircraft system. For this purpose, horizontal electric dipole (HED) transmitters were installed in the survey area and were surveyed individually. Magnetic transfer functions were determined and a model of the conductivity of the study area was derived. Despite restrictions such as low payload capacity and multicopter-related EM noise, we were able to estimate spatially and spectrally consistent transfer functions of high quality up to a distance of 2 km from the respective transmitter. Our results could be validated with independent results from a magnetotelluric and a direct current sounding. The study demonstrates that an unmanned aircraft system (UAS) is suitable for semi-airborne EM application and that such a system can be beneficial where ground-based methods and manned techniques become impractical

Evaluation of a Semi-Airborne Electromagnetic Survey Based on a Multicopter Aircraft System

The semi-airborne electromagnetic (EM) method has the potential to reach deeper exploration depths than purely airborne EM approaches. The concept of the method is to deploy high-power transmitters on the ground, which excite subsurface currents and induce strong magnetic fields, and to measure the corresponding EM fields with a passive airborne receiver instrument. Following recent conceptual developments of the semi-airborne EM technique deployed on helicopters, we performed a 10 km2 semi-airborne EM survey near Münster (Germany) based on a multicopter aircraft system. For this purpose, horizontal electric dipole (HED) transmitters were installed in the survey area and were surveyed individually. Magnetic transfer functions were determined and a model of the conductivity of the study area was derived. Despite restrictions such as low payload capacity and multicopter-related EM noise, we were able to estimate spatially and spectrally consistent transfer functions of high quality up to a distance of 2 km from the respective transmitter. Our results could be validated with independent results from a magnetotelluric and a direct current sounding. The study demonstrates that an unmanned aircraft system (UAS) is suitable for semi-airborne EM application and that such a system can be beneficial where ground-based methods and manned techniques become impractical

Transversely isotropic lower crust of Variscan central Europe imaged by ambient noise tomography of the Bohemian Massif

The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons

Shear wave splitting in the Alpine region

To constrain seismic anisotropy under and around the Alps in Europe, we study SKS shear wave splitting from the region densely covered by the AlpArray seismic network. We apply a technique based on measuring the splitting intensity, constraining well both the fast orientation and the splitting delay. Four years of teleseismic earthquake data were processed, from 723 temporary and permanent broad-band stations of the AlpArray deployment including ocean-bottom seismometers, providing a spatial coverage that is unprecedented. The technique is applied automatically (without human intervention), and it thus provides a reproducible image of anisotropic structure in and around the Alpine region. As in earlier studies, we observe a coherent rotation of fast axes in the western part of the Alpine chain, and a region of homogeneous fast orientation in the Central Alps. The spatial variation of splitting delay times is particularly interesting though. On one hand, there is a clear positive correlation with Alpine topography, suggesting that part of the seismic anisotropy (deformation) is caused by the Alpine orogeny. On the other hand, anisotropic strength around the mountain chain shows a distinct contrast between the Western and Eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. The new observational constraints, especially the splitting delay, provide new information on Alpine geodynamics

High-Resolution Crustal S-wave Velocity Model and Moho Geometry Beneath the Southeastern Alps: New Insights From the SWATH-D Experiment

We compiled a dataset of continuous recordings from the temporary and permanent seismic networks to compute the high-resolution 3D S-wave velocity model of the Southeastern Alps, the western part of the external Dinarides, and the Friuli and Venetian plains through ambient noise tomography. Part of the dataset is recorded by the SWATH-D temporary network and permanent networks in Italy, Austria, Slovenia and Croatia between October 2017 and July 2018. We computed 4050 vertical component cross-correlations to obtain the empirical Rayleigh wave Green’s functions. The dataset is complemented by adopting 1804 high-quality correlograms from other studies. The fast-marching method for 2D surface wave tomography is applied to the phase velocity dispersion curves in the 2–30 s period band. The resulting local dispersion curves are inverted for 1D S-wave velocity profiles using the non-perturbational and perturbational inversion methods. We assembled the 1D S-wave velocity profiles into a pseudo-3D S-wave velocity model from the surface down to 60 km depth. A range of iso-velocities, representing the crystalline basement depth and the crustal thickness, are determined. We found the average depth over the 2.8–3.0 and 4.1–4.3 km/s iso-velocity ranges to be reasonable representations of the crystalline basement and Moho depths, respectively. The basement depth map shows that the shallower crystalline basement beneath the Schio-Vicenza fault highlights the boundary between the deeper Venetian and Friuli plains to the east and the Po-plain to the west. The estimated Moho depth map displays a thickened crust along the boundary between the Friuli plain and the external Dinarides. It also reveals a N-S narrow corridor of crustal thinning to the east of the junction of Giudicarie and Periadriatic lines, which was not reported by other seismic imaging studies. This corridor of shallower Moho is located beneath the surface outcrop of the Permian magmatic rocks and seems to be connected to the continuation of the Permian magmatism to the deep-seated crust. We compared the shallow crustal velocities and the hypocentral location of the earthquakes in the Southern foothills of the Alps. It revealed that the seismicity mainly occurs in the S-wave velocity range between ∼3.1 and ∼3.6 km/s