Anizotropní tomografie svrchního pláště pod Evropou

Název: Anizotropní tomografie svrchního pláště pod Evropou Autor: Helena Žlebčíková Katedra: Katedra geofyziky, Matematicko-fyzikální fakulta, Univerzita Karlova Školící pracoviště: Geofyzikální ústav Akademie věd České republiky (GFÚ AV ČR) Vedoucí: RNDr. Jaroslava Plomerová, DrSc., GFÚ AV ČR Konzultanti: RNDr. Vladislav Babuška, DrSc., GFÚ AV ČR RNDr. Luděk Vecsey, Ph.D., GFÚ AV ČR Abstrakt: Výzkum seismické anizotropie kontinentální plášťové litosféry odvozené ze společné inverze/interpretace směrových variací odchylek v časech šíření teleseismických vln P a parametrů štěpení vln SKS naznačuje, že orientaci os symetrie anizotropie je potřeba uvažovat obecně ve 3D. Mnohé tomografické studie nicméně anizotropii objemových vln zanedbávají zcela nebo se omezují pouze na azimutální nebo radiální anizotropii. Proto jsme vyvinuli kód AniTomo pro sdruženou anizotropní-izotropní tomografii svrchního pláště. Kód AniTomo modeluje 3D rozložení anizotropie a perturbací izotropních rychlostí vln P ve svrchním plášti inverzí relativních odchylek v časech šíření teleseismických vln P. Předpokladem je slabá anizotropie s hexagonální symetrií. Kód připouští oba typy hexagonální symetrie, tj. s "rychlou" osou a a "pomalou" rovinou (a,c) nebo s "pomalou" osou b a "rychlou" rovinou (a,c). Navíc osa symetrie může být...Title: Anisotropic tomography of the European upper mantle Author: Helena Žlebčíková Department: Department of Geophysics, Faculty of Mathematics and Physics, Charles University Training institution: Institute of Geophysics of the Czech Academy of Sciences (IG CAS) Supervisor: RNDr. Jaroslava Plomerová, DrSc., IG CAS Consultants: RNDr. Vladislav Babuška, DrSc., IG CAS RNDr. Luděk Vecsey, Ph.D., IG CAS Abstract: Large-scale seismic anisotropy of the continental mantle lithosphere derived from joint inversion/interpretation of directional variations of P-wave travel-time residuals and SKS-wave splitting calls for orientation of the symmetry axes to be treated generally in 3D. Nevertheless, most of the tomography studies neglect the anisotropy of the body waves completely or they are limited to either azimuthal or radial anisotropy. Therefore, we have developed a code called AniTomo for coupled anisotropic-isotropic travel-time tomography of the upper mantle. The novel code allows inversion of relative travel-time residuals of teleseismic P waves simultaneously for 3D distribution of P-wave isotropic- velocity perturbations and anisotropy of the upper mantle. We assume weak anisotropy of hexagonal symmetry with either the 'high-velocity' a axis or the 'low-velocity' b axis. The symmetry axis is allowed to be...Katedra geofyzikyDepartment of GeophysicsFaculty of Mathematics and PhysicsMatematicko-fyzikální fakult

Two subduction-related heterogeneities beneath the Eastern Alps and the Bohemian Massif imaged by high-resolution P-wave tomography

We present high-resolution tomographic images of the upper mantle beneath the Eastern Alps and the adjacent Bohemian Massif (BM) in the north based on recordings from the AlpArray-EASI and AlpArray seismic networks. The tomography locates the Alpine high-velocity perturbations between the Periadriatic Lineament and the Northern Alpine Front. The northward-dipping lithosphere is imaged down to ∼ 200–250 km of depth, without signs of delamination. The small amount of crustal shortening compared to that in the Western Alps and the bimodal character of the positive perturbations with a separation beneath the Tauern Window indicate a dual source of the velocity heterogeneity, most probably formed by a mixture of a fragment of detached European plate and the Adriatic plate subductions. A detached high-velocity heterogeneity, sub-parallel to and distinct from the Eastern Alps heterogeneity, is imaged at ∼ 100–200 km beneath the southern part of the BM. We associate this anomaly with the western end of a SW–NE-striking heterogeneity beneath the south-eastern part of the BM, imaged in models of larger extent. The strike, parallel with the Moldanubian–Brunovistulian mantle–lithosphere boundary in the BM and with the westernmost part of the Carpathian front, leads us to consider potential scenarios relating the heterogeneity to (1) a remnant of the delaminated European plate, (2) a piece of continental-and-oceanic lithosphere mixture related to the building of the BM, particularly to the closure of the old Rheic ocean during the MD–BV collision, or (3) a lithospheric fragment going through to the NW between the Eastern Alps and Western Carpathians fronts in a preceding subduction phase. The study is dedicated to our outstanding and respected colleague Vladislav Babuška, who coined innovative views on the European lithosphere and died on 30 March 2021

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

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

Anisotropic tomography of the European upper mantle

Title: Anisotropic tomography of the European upper mantle Author: Helena Žlebčíková Department: Department of Geophysics, Faculty of Mathematics and Physics, Charles University Training institution: Institute of Geophysics of the Czech Academy of Sciences (IG CAS) Supervisor: RNDr. Jaroslava Plomerová, DrSc., IG CAS Consultants: RNDr. Vladislav Babuška, DrSc., IG CAS RNDr. Luděk Vecsey, Ph.D., IG CAS Abstract: Large-scale seismic anisotropy of the continental mantle lithosphere derived from joint inversion/interpretation of directional variations of P-wave travel-time residuals and SKS-wave splitting calls for orientation of the symmetry axes to be treated generally in 3D. Nevertheless, most of the tomography studies neglect the anisotropy of the body waves completely or they are limited to either azimuthal or radial anisotropy. Therefore, we have developed a code called AniTomo for coupled anisotropic-isotropic travel-time tomography of the upper mantle. The novel code allows inversion of relative travel-time residuals of teleseismic P waves simultaneously for 3D distribution of P-wave isotropic- velocity perturbations and anisotropy of the upper mantle. We assume weak anisotropy of hexagonal symmetry with either the 'high-velocity' a axis or the 'low-velocity' b axis. The symmetry axis is allowed to be..

Anizotropní tomografie svrchního pláště pod Evropou

Název: Anizotropní tomografie svrchního pláště pod Evropou Autor: Helena Žlebčíková Katedra: Katedra geofyziky, Matematicko-fyzikální fakulta, Univerzita Karlova Školící pracoviště: Geofyzikální ústav Akademie věd České republiky (GFÚ AV ČR) Vedoucí: RNDr. Jaroslava Plomerová, DrSc., GFÚ AV ČR Konzultanti: RNDr. Vladislav Babuška, DrSc., GFÚ AV ČR RNDr. Luděk Vecsey, Ph.D., GFÚ AV ČR Abstrakt: Výzkum seismické anizotropie kontinentální plášťové litosféry odvozené ze společné inverze/interpretace směrových variací odchylek v časech šíření teleseismických vln P a parametrů štěpení vln SKS naznačuje, že orientaci os symetrie anizotropie je potřeba uvažovat obecně ve 3D. Mnohé tomografické studie nicméně anizotropii objemových vln zanedbávají zcela nebo se omezují pouze na azimutální nebo radiální anizotropii. Proto jsme vyvinuli kód AniTomo pro sdruženou anizotropní-izotropní tomografii svrchního pláště. Kód AniTomo modeluje 3D rozložení anizotropie a perturbací izotropních rychlostí vln P ve svrchním plášti inverzí relativních odchylek v časech šíření teleseismických vln P. Předpokladem je slabá anizotropie s hexagonální symetrií. Kód připouští oba typy hexagonální symetrie, tj. s "rychlou" osou a a "pomalou" rovinou (a,c) nebo s "pomalou" osou b a "rychlou" rovinou (a,c). Navíc osa symetrie může být...Title: Anisotropic tomography of the European upper mantle Author: Helena Žlebčíková Department: Department of Geophysics, Faculty of Mathematics and Physics, Charles University Training institution: Institute of Geophysics of the Czech Academy of Sciences (IG CAS) Supervisor: RNDr. Jaroslava Plomerová, DrSc., IG CAS Consultants: RNDr. Vladislav Babuška, DrSc., IG CAS RNDr. Luděk Vecsey, Ph.D., IG CAS Abstract: Large-scale seismic anisotropy of the continental mantle lithosphere derived from joint inversion/interpretation of directional variations of P-wave travel-time residuals and SKS-wave splitting calls for orientation of the symmetry axes to be treated generally in 3D. Nevertheless, most of the tomography studies neglect the anisotropy of the body waves completely or they are limited to either azimuthal or radial anisotropy. Therefore, we have developed a code called AniTomo for coupled anisotropic-isotropic travel-time tomography of the upper mantle. The novel code allows inversion of relative travel-time residuals of teleseismic P waves simultaneously for 3D distribution of P-wave isotropic- velocity perturbations and anisotropy of the upper mantle. We assume weak anisotropy of hexagonal symmetry with either the 'high-velocity' a axis or the 'low-velocity' b axis. The symmetry axis is allowed to be...Katedra geofyzikyDepartment of GeophysicsMatematicko-fyzikální fakultaFaculty of Mathematics and Physic

Anisotropic tomography of the European upper mantle

Title: Anisotropic tomography of the European upper mantle Author: Helena Žlebčíková Department: Department of Geophysics, Faculty of Mathematics and Physics, Charles University Training institution: Institute of Geophysics of the Czech Academy of Sciences (IG CAS) Supervisor: RNDr. Jaroslava Plomerová, DrSc., IG CAS Consultants: RNDr. Vladislav Babuška, DrSc., IG CAS RNDr. Luděk Vecsey, Ph.D., IG CAS Abstract: Large-scale seismic anisotropy of the continental mantle lithosphere derived from joint inversion/interpretation of directional variations of P-wave travel-time residuals and SKS-wave splitting calls for orientation of the symmetry axes to be treated generally in 3D. Nevertheless, most of the tomography studies neglect the anisotropy of the body waves completely or they are limited to either azimuthal or radial anisotropy. Therefore, we have developed a code called AniTomo for coupled anisotropic-isotropic travel-time tomography of the upper mantle. The novel code allows inversion of relative travel-time residuals of teleseismic P waves simultaneously for 3D distribution of P-wave isotropic- velocity perturbations and anisotropy of the upper mantle. We assume weak anisotropy of hexagonal symmetry with either the 'high-velocity' a axis or the 'low-velocity' b axis. The symmetry axis is allowed to be..

Anisotropic tomography of the European upper mantle

Title: Anisotropic tomography of the European upper mantle Author: Helena Žlebčíková Department: Department of Geophysics, Faculty of Mathematics and Physics, Charles University Training institution: Institute of Geophysics of the Czech Academy of Sciences (IG CAS) Supervisor: RNDr. Jaroslava Plomerová, DrSc., IG CAS Consultants: RNDr. Vladislav Babuška, DrSc., IG CAS RNDr. Luděk Vecsey, Ph.D., IG CAS Abstract: Large-scale seismic anisotropy of the continental mantle lithosphere derived from joint inversion/interpretation of directional variations of P-wave travel-time residuals and SKS-wave splitting calls for orientation of the symmetry axes to be treated generally in 3D. Nevertheless, most of the tomography studies neglect the anisotropy of the body waves completely or they are limited to either azimuthal or radial anisotropy. Therefore, we have developed a code called AniTomo for coupled anisotropic-isotropic travel-time tomography of the upper mantle. The novel code allows inversion of relative travel-time residuals of teleseismic P waves simultaneously for 3D distribution of P-wave isotropic- velocity perturbations and anisotropy of the upper mantle. We assume weak anisotropy of hexagonal symmetry with either the 'high-velocity' a axis or the 'low-velocity' b axis. The symmetry axis is allowed to be..

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