TS7.8

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. 4 years of teleseismic earthquake data were processed automatically (without human intervention), from 724 temporary and permanent broadband stations of the AlpArray deployment including ocean-bottom seismometers. We have obtained an objective image of anisotropic structure in and around the Alpine region, at a spatial resolution that is unprecedented. 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. 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 western and eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. We discuss earlier concepts of Alpine geodynamics in the light of these new observational constraints. 

How to cite: Bokelmann, G., Hein, G., Kolinsky, P., Bianchi, I., and Working Group, A.: Shear-Wave Splitting in the Alpine Region, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5377, https://doi.org/10.5194/egusphere-egu21-5377, 2021.

In the frame of the AlpArray project we analyse teleseismic data from permanent and temporary stations of the greater Alpine region to study seismic discontinuities in the entire lithosphere. We use broadband S-to-P converted signals from below the seismic stations. In order to avoid sidelobes, no deconvolution or filtering is applied and S arrival times are used as reference. We show a number of north-south and east-west profiles through the greater Alpine area. The Moho signals are always seen very clearly, and also negative velocity gradients below the Moho are visible in a number of profiles. The subducting European Moho is visible in the Eastern Alps west of 13.5°E (the eastern edge of the Tauern Window) and reaches there about 60km depth at 47°N. East of about 13.5°E, the image of the Moho changes completely. No south dipping European Moho is found anymore, instead the Moho is shallowing towards the Pannonian Basin. This suggests severe post-nappe emplacement modifications east of about 13.5°E, most probably associated with delamination of the mantle lithosphere within the formerly subducting European slab, i.e. mantle that separated from the crustal parts of the Alpine-West Carpathian orogen during the last ca. 20 Ma when the Pannonian basin formed and the ALCAPA block underwent its E-directed lateral extrusion.

Ratschbacher, L., Frisch, W., Linzer, H.-G. and Merle, O. (1991) Lateral extrusion in the Eastern Alps, Part 2: Structural analysis. Tectonics, vol.10, No.2, 257-271.

How to cite: Kind, R., Schmid, S., Yuan, X., and Heit, B.: The seismic structure of the lithosphere in the greater Alpine area from S-to-P converted waves, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-442, https://doi.org/10.5194/egusphere-egu21-442, 2021.

We present potential scenarios of the European and Adriatic plates’ collision that formed the Alps and the neighbouring mountain belts. Our results are based on teleseismic body-wave data from the AlpArray-EASI complementary experiment (2014-2015, Hetényi et al., Tectonophysics 2018) and the AlpArray Seismic Network (Hetényi et al., Surv. Geophys. 2018).  Tomography of seismic velocities in the upper mantle along a ca. 200 km broad and 540 km long north-south transect images steady southward thickening of the lithosphere beneath the Bohemian Massif  and northward dipping East-Alpine lithospheric keel. Thanks to the dense spacing of the AlpArray Seismic Network stations and high-quality data, the high-resolution tomography resolves for the first time two sub-parallel down-going high-velocity heterogeneities beneath the Eastern Alps, instead of a single, thick anomaly. The southern heterogeneity, which we relate to the subducted Adriatic plate, is more distinct than the northern one, which loses its connection with the shallow parts. Moreover, amplitudes and size of this heterogeneity decrease in cross-sections perpendicular to the strike of the Alps when shifting towards the Central Alps. The presented collision scenarios consider the smaller northern heterogeneity as (1) a remnant of a delaminated early phase subduction of the European plate with the reversed polarity relative to that in the Western Alps, (2) a piece of continental and oceanic lithosphere together, or, (3) a fragment of a quite extended lithosphere margin foundering in a preceding phase of the Adriatic subduction.

How to cite: Plomerova, J., Zlebcikova, H., Hetenyi, G., Vecsey, L., Babuska, V., and Working Group, A.: European-Adriatic plate collision in teleseismic tomography of the Eastern Alps, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10724, https://doi.org/10.5194/egusphere-egu21-10724, 2021.

We perform a teleseismic P-wave travel time tomography to examine geometry and slab structure of the upper mantle beneath the Alpine orogen. Vertical component data of the extraordinary dense seismic network AlpArray are used which were recorded at over 600 temporary and permanent broadband stations deployed by 24 different European institutions in the greater Alpine region, reaching from the Massif Central to the Pannonian Basin and from the Po plain to the river Main. Mantle phases of 370 teleseismic events between 2015 and 2019 of magnitude 5.5 and higher are evaluated automatically for direct and core diffracted P arrivals using a combination of higher-order statistics picking algorithms and signal cross correlation. The resulting database contains over 170.000 highly accurate absolute P picks that were manually revised for each event. The travel time residuals exhibit very consistent and reproducible spatial patterns, already pointing at high velocity slabs in the mantle.

For predicting P-wave travel times we consider a large computational box encompassing the Alpine region up to a depth of 600 km within which we allow 3D-variations of P-wave velocity. To account for influences of the strongly heterogeneous crust that cannot be resolved with teleseismic data, we integrate a complex three-dimensional crustal model directly into our model. Outside the box we assume a spherically symmetric earth and apply the Tau-P method to calculate travel times and ray paths. These are injected at the boundaries of the regional box and continued using the fast marching method (Rawlinson et al. 2005). We invert differences between observed and predicted traveltimes for P-wave velocities inside the box. Velocity is discretized on a regular grid with a spacing of about 25x25x15 km. The misfit reduction reaches values of over 80% depending on damping and smoothing parameters.

The resulting model shows several steeply dipping high velocity anomalies following the Alpine arc. The most prominent structure stretches from the western Alps into the Apennines mountain range reaching depths of over 500 km. Two further anomalies of high complexity extending down to a depth of 300 km are located below the central and eastern Alps, both being detached from the lithosphere and separated by a clear gap below the western part of the Tauern window. The central anomaly shows mainly southwards dipping, whereas the eastern anomaly is mainly dipping to the northeast. We compare our results to former studies, confirming lateral positions of the anomalies. However, the new results can benefit from the superior resolution capabilities of the dense AlpArray seismic network, providing more accurate insights into depth extent, dip angle and directions. We perform various general, as well as purpose-built resolution tests, to verify the capabilities of our setup to resolve slab gaps as well as different possible slab dipping directions.

How to cite: Paffrath, M. and Friederich, W. and the the AlpArray Working Group: Teleseismic P-wave travel time tomography of the Alpine upper mantle using AlpArray seismic network data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8429, https://doi.org/10.5194/egusphere-egu21-8429, 2021.

Noise cross-correlations provide a good azimuthal coverage, limited only by the distribution of noise sources and the layout of the stations used. It is therefore a promising method to constrain azimuthal anisotropy. As noise cross-correlations consist mainly of surface waves, they are especially sensitive to the crust and provide good depth constraints, as opposed to SKS-splitting data that are more sensitive to the upper mantle. We use the AlpArray network as well as stations from permanent networks all across Europe to perform time-domain beamforming on noise cross-correlations. The extent and density of the AlpArray network allows us to obtain reliable measurements all across the Alps. We divide the area in smaller zones using all stations outside the zone as sources and all stations inside as a sub-array for beamforming. This allows us to estimate the quality of our measurements in a region where strong lateral heterogeneities make measurements challenging, by estimating the magnitude of bias due to heterogeneities using the cos(theta) amplitude and evaluating uncertainties with bootstrap. This way, we measure Rayleigh wave azimuthal anisotropy in several period bands between 15 s and 60 s period. Inversion of dispersion curves in specific areas allows us to constrain the depth of the observed anisotropy. The results are broadly similar to results from SKS-splitting as they are generally parallel to the mountain belt. However, we observe lower anisotropy at short periods (40 seconds and less) in the Alps themselves than in surrounding regions. We also observe several structures in the crust that are not observed with SKS-splitting data. The most striking is a strong and spatially coherent NE-oriented anisotropy to the NW of the Alps that is possibly related to Variscan inheritance (at 40 seconds and less, in the upper and lower crust).  In the Northern Apennines, we observe anisotropy perpendicular to the belt at 30 s period (middle crust) that correlates well with an area of strong radial anisotropy recently observed by Alder et al (in review) at 30 km depth. 

How to cite: Soergel, D., Pedersen, H., Bodin, T., Paul, A., and Stehly, L.: Imaging azimuthal anisotropy in the alpine crust from ambient noise beamforming., EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13514, https://doi.org/10.5194/egusphere-egu21-13514, 2021.

We have successfully derived a new 3-D high resolution shear wave velocity model of the crust and uppermost mantle of a large part of W-Europe from transdimensional ambient-noise tomography. This model is intended to contribute to the development of the first 3-D crustal-scale integrated geophysical-geological model of the W-Alps to deepen understanding of orogenesis and its relationship to mantle dynamics.

We used an exceptional dataset of 4 years of vertical-component, daily seismic noise records (2015 - 2019) of more than 950 permanent broadband seismic stations located in and around the Greater Alpine region, complemented by 490 temporary stations from the AlpArray sea-land seismic network and 110 stations from Cifalps dense deployments.

We firstly performed a 2-D data-driven transdimensional travel time inversion for group velocity maps from 4 to 150 s (Bodin & Sambridge, 2009). The data noise level was treated as a parameter of the inversion problem, and determined within a Hierarchical Bayes method. We used Fast Marching Eikonal solver (Rawlinson & Sambridge, 2005) jointly with the reversible jump algorithm to update raypath geometry during inversion. In the inversion of group velocity maps for shear-wave velocity, we set up a new formulation of the approach proposed by Lu et al (2018) by including group velocity uncertainties. Posterior probability distributions on Vs and interfaces were estimated by exploring a set of 130 millions synthetic 4-layer 1-D Vs models that allow for low-velocity zones. The obtained probabilistic model was refined using a linearized inversion. For the ocean-bottom seismometers of the Ligurian-Provencal basin, we applied a specific processing to clean daily noise signals from instrumental and oceanic noises (Crawford & Webb, 2000) and adapted the inversion for Vs to include the water column.

Our Vs model evidences strong variations of the crustal structure along strike, particulary in the subduction complex. The European crust includes lower crustal low-velocity zones and a Moho jump of ~8-12 km beneath the W-boundary of the external crystalline massifs. We observe a deep LVZ structure (50 - 80 km) in the prolongation of the European continental subduction beneath the Ivrea body. The striking fit between the receiver functions ccp migrated section across the Cifalps profile and this new Vs model validate its reliability.

How to cite: Nouibat, A., Stehly, L., Paul, A., Brossier, R., Bodin, T., Schwartz, S., and Working Group, A.: First step towards an integrated geophysical-geological model of the W-Alps: A new Vs model from transdimensional ambient-noise tomography, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3197, https://doi.org/10.5194/egusphere-egu21-3197, 2021.

The CIFALPS receiver-function (RF) profile in the southwestern Alps provided the first seismological evidence of continental subduction in the Alps, with the detection of waves converted on the European Moho at 75-80 km depth beneath the western edge of the Po basin (Zhao et al., 2015). To complement the CIFALPS profile and enhance our knowledge of the lithospheric structure of the Western Alps, we installed CIFALPS2, a temporary network of 55 broadband seismic stations that operated for ~14 months (2018-2019) across the North-Western Alps (Zhao et al., 2018). The CIFALPS2 line runs from the Eastern Massif Central to the Ligurian coast, across the Mont-Blanc and Gran Paradiso massifs and the Ligurian Alps. Seismic stations were installed along a quasi-linear profile with a spacing of 7-10 km.

We will show 2 receiver-function CCP (common-conversion point) depth-migrated sections along the CIFALPS2 profile, the first one across the Alps, and the second one across the Ligurian Alps and the Po basin. The time-to-depth migration of RF data is based on the new 3-D Vs model of the Greater Alpine region derived by Nouibat et al. (2021) using transdimensional ambient noise tomography on a large dataset including the AlpArray seismic network. Depth sections across the Vs model are also useful for interpreting the RF CCP sections as they have striking similarities.

The images of the lithospheric structure of the NW Alps along CIFALPS2 are surprisingly different from those of the SW Alps along CIFALPS. The deepest P-to-S converted phases on the European Moho are detected at 60-65 km depth beneath the Ivrea-Verbano zone, that is 15 km less than on CIFALPS. The negative polarity converted phase interpreted as the base of the Ivrea body mantle flake on the CIFALPS section is still visible on CIFALPS2, but with a lower amplitude. The RF section confirms the existence of a jump of the European Moho of ~10 km amplitude in less than 10 km distance, which is located within a few km from the western boundary of the Mont Blanc external crystalline massif. All these observations are confirmed by the Vs model that also displays a less deep continental subduction than on CIFALPS, weaker S-wave velocities in the Ivrea body wedge, and the jump of the European Moho.

The Moho beneath the Ligurian Alps is detected at 25-30 km depth both on the RF and on the Vs depth sections. Moving northwards, this Ligurian Moho is separated from the Adriatic Moho by a puzzling S-dipping set of P-to-S converted waves with negative polarity. The crust of the Ligurian Alps is characterized by a set of north-dipping negative-polarity converted waves at 10 to 20 km depth beneath the Valosio massif, which is a small internal crystalline massif of (U)HP metamorphic rocks located north of Voltri. The similarity of this set of negative-polarity conversions to the one observed beneath the Dora Maira massif on the CIFALPS profile suggests that it may be a relic of the Alpine structure overprinted by the opening of the Ligurian sea.

How to cite: Paul, A., Nouibat, A., Zhao, L., Solarino, S., Schwartz, S., Malusà, M., Stehly, L., Aubert, C., Dumont, T., Eva, E., Guillot, S., Pondrelli, S., Salimbeni, S., and Working Group, A.: Striking differences in lithospheric structure between the north- and south-western Alps: insights from receiver functions along the Cifalps profiles and a new Vs model, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9391, https://doi.org/10.5194/egusphere-egu21-9391, 2021.

The Alpine orogeny, which was formed by subduction of the European plate beneath the Adria plate, is considered as one of the world’s foremost natural laboratories for the study of orogenic processes. In contrast to other mountain belts, the Western Alpine belt is curved and affected by three-dimensional effects. Due to differences in stress distribution and rheological properties of crustal rocks, the Moho geometry and crustal structure along different sections differ, in particular in the vicinity of the continental subduction complex.

To better understand the configuration of continental subduction along a profile that crosscuts the North-Western Alps, we combine receiver function analysis with computation of synthetic receiver functions and gravity anomaly modeling to precise the subduction structures and estimate a crustal 2D shear wave velocity and density model. Seismic data come from the CIFALPS2 (China-Italy-France Alps seismic survey) temporary experiment, which operated from 2018 to 2020. We use a 2D hybrid waveform simulation method (Zhao et al., 2008) that is reliable and efficient and has a better response to 2D structures compared to conventional 1D waveform inversion methods, in particular for the dipping Moho interface of the subduction complex. We compute synthetic receiver functions for a large set of models compatible with surface geology data, which are then processed to obtain synthetic CCP depth-migrated stacks. Furthermore, we model the Bouguer gravity data along the same profile to obtain preferred density distribution. The nature of rocks in the subduction complex can be inferred from our synthetical models.

Compared to the results of the CIFALPS profile in the Central Western Alps (Zhao et al., 2015), the subduction along the CIFALPS2 profile has a shallower dip angle, which is a significant difference between the two sections. As for velocity and density models, the two sections have a high velocity and high-density wedge in the subduction complex. We argue that the reason for the difference in crustal structures between the two sections may be related to the difference in stress distribution.

Zhao, L., et al. (2008). "A two-dimensional hybrid method for modeling seismic wave propagation in anisotropic media." Journal of Geophysical Research 113(B12).

Zhao, L., et al. (2015). "First seismic evidence for continental subduction beneath the Western Alps." Geology 43(9): 815-818.

How to cite: Mao, Y., Zhao, L., Paul, A., Solarino, S., Aubert, C., Dumont, T., Eva, E., Guillot, S., Malusà, M. G., Pondrelli, S., Salimbeni, S., and Schwartz, S.: Configuration of continental subduction beneath the Western Alps: results using forward modeling of receiver functions and gravity data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8117, https://doi.org/10.5194/egusphere-egu21-8117, 2021.

In the framework of the SeismoTeCH project, which aims at advancing our understanding of seismotectonic processes in Switzerland, we present the first 3-D attenuation model of the upper crust for the Central Alps and their northern foreland. The 3-D inversions derive the quality factor Q (1/attenuation) using path attenuation t^∗ observations for 4,192 distributed earthquakes recorded on permanent and temporary stations, including both velocity and acceleration records for the period 2002-2019. We followed a procedure of gradational inversions, in which a series of inversions with increasingly grid complexity are performed, with the goal of obtaining a useful Q model everywhere despite the varied data distribution. The Qs and Qp results show large-scale features in the upper crust, which are consistent with a recently improved high-resolution velocity models of the same region and serve to refine the interpretations of crustal structures from Vp and Vp/Vs. For example, the foreland region of southern Germany and northern Switzerland show a low Q crustal block bounded by high Q regions in the uppermost layer between -2.5 and 2.0 km depth. This markedly correlates with the overlying surface geology, where low Q areas coincide with the Molasse Basin, and the transition between low and high Q regions outline the geological boundary between the Molasse and the Mesozoic sediments towards north and the Alpine front to the south. At depths ranging between 2.0 - 6.5 km, low Q is imaged along the Rhone valley in the Valais in southwest Switzerland. This region presents the transition between the Centrals and Western Alps and hosts the presently seismically most active fault zones. As the attenuation of fractured areas is enhanced by fluids, low Q values may relate here to distributed microfractures that produce greater fracture connectivity and permeability in a relatively higher strain-rate zone. These geophysical constraints seem to support crustal scale fluid flow along fracture networks as manifest by the prominent occurrence of hot springs in this area. On the other hand, the moderate-to-high Qs and Qp (400-800) along with low Vp/Vs ratio and high Vs observed in the external Aar Massif could be indicative of metamorphic processes leading to different Vp/Vs ratios compared to the basement in the northern foreland (Black Forest Massif), and possibly image the continuation of the massif 20-30 km further to the northeast. In combination with recently developed Vp and Vs velocity models, the developed 3-D attenuation models provide additional constraints in terms of composition and physical properties of the uppermost crust of the central Alps as well as crucial input for next generation seismic hazard models of Switzerland, allowing for a more realistic prediction of earthquake related ground motions.

How to cite: Lanza, F., Diehl, T., Eberhart-Phillips, D., Herwegh, M., Fäh, D., and Wiemer, S.: 3-D Qp and Qs Seismic Attenuation for the Central Alps and their Foreland, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8028, https://doi.org/10.5194/egusphere-egu21-8028, 2021.

The Alpine orogen and the Apennine system are part of the complex tectonic setting in the Mediterranean Sea caused by the convergence between Africa and Eurasia. Between 30 Ma and 15 Ma the Calabrian subduction retreated in a southeast direction pulling Corsica and Sardinia away from the Eurasian landmass. In this extensional setting, the Ligurian Sea was formed as a back-arc basin. The rifting jumped 15 Ma ago in the Tyrrhenian Sea leaving Corsica and Sardinia in a stable position relative to Eurasia.

Within the framework of the AlpArray research initiative a long-term seismological experiment was conducted in the Ligurian Sea to investigate the lithospheric structure and the seismicity in the Ligurian basin. The passive seismic network consisted of 29 broad-band ocean bottom stations from Germany and France. It was in operation between June 2017 and February 2018.

Two seismicity clusters occurred in the centre of the Ligurian Basin. The 18 earthquakes are located in the lower crust and in the upper-most mantle at depths between 10 km and 16 km. Re-location was performed only using picks from the OBS in the centre of the Ligurian Sea to avoid artifacts from the complex 3D velocity structure of the basin. Mantle refractions Pn and Sn have apparent velocities of 8.2 km/s and 4.7 km/s. The low Vp-Vs-ratio of 1.72 indicates a more brittle behaviour of the mantle material.

Fault plane solutions were determined for four events using also the data of land stations in southern France, Corsica, Sardinia and northern Italy. The focal mechanisms are thrust faulting. Fault planes strike in a NE-SW direction, coinciding with the alignment of the events and the direction of the basin axis.

We interprete the two earthquake clusters related to the inversion of the Ligurian Basin where the basin’s centre is under compression and stresses are taken up by reactivated faults in the crust and uppermost mantle. The compressional forces could be caused by the convergence of Africa and Europe. In general, observations of earthquakes in continental mantle lithosphere are rare and they reveal on the one hand a strengthening of the crust and uppermost mantle during rifting and on the other hand they support the interpretation that rifting failed in the northern Ligurian Basin.

How to cite: Thorwart, M., Dannowski, A., Grevemeyer, I., Lange, D., Kopp, H., Petersen, F., Crawford, W., and Paul, A. and the AlpArray Working Group: Local Seismicity indicates basin inversion in the Ligurian Sea, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6411, https://doi.org/10.5194/egusphere-egu21-6411, 2021.

How to cite: Sonnet, M., Labrousse, L., Bascou, J., and Plunder, A.: Geophysical signature of the alpine slab: Field analogues and direct models, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-16332, https://doi.org/10.5194/egusphere-egu21-16332, 2021.

The Alps mountains and its forelands consist of a heterogeneous lithosphere, comprised of a multitude of tectonic blocks from different tectonic provinces with different thermo-physical properties. Patterns of seismicity distribution are also observed to vary significantly throughout the region. However, the relationship between seismicity and lithospheric heterogeneity has been often overlooked in previous studies. We present an overview of recent results that have attempted to address these questions through the use of integrated 3D modelling techniques, thereby including: (i) a gravity and seismic data constrained, 3D, density structural model of the lithosphere; (ii) a 3D thermal model constrained against available wellbore temperature data; and,  (iii) a 3D rheological model of the long-term lithospheric strength and effective viscosities. Our models support the existence of a first-order correlation between the distribution of seismicity (laterally and with depth) and the strength of the lithosphere, with the former being clustered mainly within weaker domains. Beneath the Alps, observed upper-crustal level (i.e., unimodal) seismicity correlates with a weaker lithosphere where plate strength is controlled by the thick crustal root. Whereas in the southern foreland, weaker zones are found preferentially around the stronger Adriatic indenter while in the northern foreland they are located in the crust beneath the the Upper Rhine Graben (URG). We found that this correlation is primarily controlled by resolved thermal gradients and is a function of the tectonic inheritance setting (e.g., UGR), crustal architecture (e.g., thickness of sediments, upper and lower crust) and LAB depth. Sediment thickness and topographic effects controls the shallow thermal filed (0 – 10 km) whereas the deeper thermal field is controlled by the thickness of felsic upper crust (higher radiogenic heat contribution), the mafic lower crust (less radiogenic heat contribution) and basal thermal boundary condition from LAB depth. Seismicity is bounded by specific isotherms, 450 ^oC in the crust and < 600 ^oC in the mantle, except in regions where slabs are imaged by seismic tomography models. This is in contrast to the recent proposition that convergence velocity is a first-order factor controlling seismicity in an orogen rather than its architecture. Fast convergence rates (e.g., Himalayas) have been related to the subduction of the cold crust to deeper crustal depths thereby leading to a deepening of the brittle  domain and to a bimodal (i.e., upper and lower crust) seismicity character. In contrast, slow convergence (e.g., Alps) is thought to lead to a hotter ductile lower crust thus limiting brittle deformation within the upper crust. We therefore end our contribution by opening a discussion on the relative role of convergence rates and lithospheric heterogeneities, inherited and/or developed during orogenesis, in controlling the seismicity. In doing so we carry out a comparison between observed seismicity and lithospheric architecture in the other mountain ranges of the western Alpine-Himalayan collision zone where  convergence velocities are of a similar order of magnitudes as Alps, i.e., the Betics, the Pyrenees and the Apennines but where seismicity is observed to occur both at upper and lower crustal levels.

How to cite: Kumar, A., Spooner, C., Scheck-Wenderoth, M., and Cacace, M.: How seismicity relates to lithospheric heterogeneity in the Alps, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11552, https://doi.org/10.5194/egusphere-egu21-11552, 2021.

Geodetic data play a crucial role in the detection of surface deformation related to active tectonic processes. The present study aimed to investigate the Northeastern Italian sector, characterized by a convergent regime due to the NNW-ward motion of the Adria microplate towards the Eurasian plate, at a rate of ~ 2mm/yr. N-S shortening is accommodated by fold and thrust systems in the Alpine chain and buried below the Friuli-Venetian plain sediments. We used InSAR and GNSS data respectively in 2015-2019 and 2000-2020 time interval to estimate the surface kinematics and deformation pattern of the area. We processed the SAR images acquired by the European satellites Sentinel 1A/B from ascending and descending tracks by using the Stanford Method for Persistent Scatterers (StaMPS). A post-processing of the resulting Line-Of-Sight (LOS) deformation time series was carried out by applying a spatial-temporal filter and calibrating using the velocities provided by GNSS stations. Finally, the post-processed ascending and descending LOS measurements were combined to solve for the vertical and horizontal (east-west) deformation components. We observed a positive vertical signal toward the Alps, in the northern region of Veneto and Friuli-Venezia Giulia. Moreover, we observed a significant negative vertical signal located in the plain and in the coastal zones due to the subsidence that strongly affects these areas. Horizontal velocities with rate of 1-2 mm/yr are observed close to main tectonic structures, especially in the eastern and the northwestern sector of the study area, where GNSS data reveal higher shortening rate.

How to cite: Areggi, G., Merryman Boncori, J. P., Pezzo, G., Serpelloni, E., and Bonini, L.: Surface deformation analysis in Northeast Italy by using PS-InSAR and GNSS data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12724, https://doi.org/10.5194/egusphere-egu21-12724, 2021.

The Alps have an overall East-West orientation, which changes radically in their western termination, where they rotate southward into a N-S strike, and then eastward into an E-W strike, forming the arc of the Western Alps. This arc is commonly inferred to have formed during collision, due to indentation of the Adriatic plate into the European continental margin. Several models attempted to provide a kinematic explanation for the formation of this arched, lateral end of the Alps. Indeed, the radial nature of the transport directions observed along the arc of the Western Alps cannot be explained by a classic convergence model.
For more than 50 years the formation of this arc was been associated to westward-directed indentation of Adria, accommodated along East-West oriented strike-slip faults, a sinistral one in the South of the arc and a dextral one in the North. The dextral one correspond to the Insubric Fault. The sinistral strike-slip zone, inferred to be localized along the «Stura corridor» (Piedmont, Italy) would correspond to a displacement of 100 to 150 km according to palaeogeographical, and geometric analyses. However, field evidence is scarce and barely documented in the literature.
Vertical axis rotations of the Adriatic indenter also inferred to be syn-collisional could have influenced the acquisition of the morphology of the arc. Paleomagnetic analyses carried out in the Internal Zone and in the Po plain suggest a southward increading amount of counter-clockwise rotation of the Adriatic plate and the Internal Zone, varying from 20°-25° in the North to nearly 120° in the South.
Dextral shear zones possibly accommodating this rotation in some conceptual models is observed in several places below the Penninic Front and affect the Argentera massif to the south. However, the measured displacement quantities do not appear to be equivalent to those induced by such rotations.
The present study aims to constrain the kinematic evolution of the arc of the Western Alps through a multidisciplinary approach. The first aspect of this project is the structural analysis of the area (Stura corridor) inferred to accommodate large sinistral displacements allowing for the westward indentation of the Adriatic indenter. We discuss the general lack of field evidence supporting sinistral strike-slip movements, in contrast to large-scale compilation of structures suggesting the possible occurrence of such displacement. The second part consists of a palaeomagnetic study, in which new data are integred with a compilation of already existing data. This compilation shows that several parts of the arc in the External Zone did not suffer any Cenozoic rotations, hence suggesting that a proto-arc already axisted at the onset collision, as suggested by independent evidence of some paleogeographic reconstruction. Finally, 2D and 3D thermo-mechanical modeling in using the pTatin3D code is used to test which structural (geometrical), and rheological parameters affected the first-order morphology of the Western Alpin arc and its kinematics. The synthesis of these different approaches allows us to propose a new model explaining the kinematics and the mechanisms of formation of the Western Alps arc.

How to cite: Brunsmann, Q., Rosenberg, C., Bellahsen, N., and Le Pourhiet, L.: The arc of the Western Alps : understanding its kinematics and formation mechanisms based on new structural and paleomagnetic datas, and thermo-mechanical modelling., EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-16529, https://doi.org/10.5194/egusphere-egu21-16529, 2021.

The Austroalpine Ötztal Nappe shows pervasive Eoalpine and local Variscan high-pressure metamorphism and deformation in its southeastern end, which obscure prior structures. We used magmatic and detrital zircon U-Pb dating by laser ablation ICP-MS to identify the precursor units of the Ötztal Nappe and the relationships among them.

Magmatic protolith dating of granitoid othogneisses in the Ötztal basement yielded Ordovician ages (450 – 470 Ma). The zircons of the Ordovician magmatism are important markers in the detrital zircon record. The paragneisses of the Ötztal basement, in which the Ordovician granitoids intruded, show no Ordovician zircons. The partly calcareous metasediments of the Schneeberg Complex and the Laas Series record some Ordovician detrital zircons. While the Schneeberg Complex is in tectonic contact to the Ötztal Nappe (Klug & Froitzheim, subm.), the Laas Series is the post-Ordovician sedimentary cover of the Ötztal basement. A Permo-Triassic basal metasandstone of the Brenner Mesozoic shows next to a strong Ordovician zircon age population some Variscan and Permo-Triassic zircons.

Zircon dating allowed to identify pre-Ordovician basement with Ordovician intrusions covered by post-Ordovician-pre-Variscan and Permo-Mesozoic sediments. This supports the concept of a non-tectonic unity in the southeastern Ötztal Nappe outside of the Schneeberg Complex.

How to cite: Klug, L., Froitzheim, N., Tomaschek, F., and Lagos, M.: Ordovician zircons as detrital markers in the Ötztal Nappe (Austroalpine, Italy), EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2080, https://doi.org/10.5194/egusphere-egu21-2080, 2021.

Despite significant recent progress in the understanding and quantification of the parameters controlling deformation modes in carbonate multilayers within fold-and-thrust belts, the details of early deformation and faulting during the initial stages of large-scale thrusting remain poorly documented and understood. Aiming to narrow this knowledge gap, we have chosen to study the relatively low-strain carbonate multilayer footwall of the Belluno Thrust (BT), one of the most external and S-vergent thrusts of the eastern Southern Alps (Italy). The BT footwall is composed of a c. 600 m thick Meso-Cenozoic multilayer succession of shallow water carbonate and pelagic sedimentary units characterized by strong mineralogical heterogeneity, with calcite (32-98%), sheet silicates (1-27%), and quartz (1-37%) as principal components. Its structural framework reflects cumulative strain due to multiple deformation events and is defined by the superposition of different structures such as i) south-verging asymmetric folds, ii) faulted folds, cut by slip planes with centimetric to metric throw, iii) SC-C’ fabrics in the marly layers, and iv) cataclastic domains.  Structures recording the early shortening increments are generally well preserved mesoscopic upright folds. Asymmetric folds with gently N-dipping backlimbs and steeply S-dipping (or even overturned N-dipping) forelimbs, record further shortening of the early upright and symmetrical folds. Strain is strongly partitioned within the marly layers, with discrete faults commonly defined by multiple slip surfaces forming duplex geometries and SC-C’ fabrics and exploiting millimetric to centimetric marly beds as detachment layers. Thrusts and diffuse reverse faults not associated with any cataclasite localise along the backlimbs of the asymmetric folds, suggesting dominant layer-parallel shortening. Cataclasites develop instead along the thrust surfaces that cut across the steeply dipping (locally even overturned) forelimbs, where cataclastic flow becomes the dominant deformation mechanism. On the vertical forelimbs, cataclasis and strain localisation are commonly associated with veins, which contributed to harden the rock system.  

Based on our systematic observations, we propose that deformation progressively evolved from folding and layer-parallel shortening (initial phases) to faulting and cataclasis (final phases) as a function of the dynamic interplay of the following factors: i) the geometrical relationships between fault orientation, fold attitude (forelimb and backlimb domains) and stress field, ii) the lithotype, which we conveniently account for by referring to the ratio between the cumulative thickness of the outcrop marly layers and the total measured stratigraphic thickness, iii) the involvement of fluids during deformation, iv) the mineral assemblage of the involved layers and v) the geometric framework of the domain localising strain with respect to the principal stress axes orientation. We conclude that these parameters play a major role in guiding strain localisation and partitioning during continuous shortening within fold-and-thrust belts. They also govern the transition from overall aseismic creep to coseismic rupturing at the scale of mesoscopic structures and, possibly, of the entire belt.

How to cite: Zuccari, C., Viola, G., Vignaroli, G., and Aldega, L.: What steers deformation within fold-and-thrust belts: insights from the carbonate multilayer footwall of the Belluno Thrust, Italian eastern Southern Alps, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1068, https://doi.org/10.5194/egusphere-egu21-1068, 2021.

The Dinarides fold and thrust belt on the Balkan Peninsula is the result of the long-lasting convergence between the Adriatic and Eurasian plates since the Mid-Jurassic. Late Jurassic obduction of ophiolites, Early Cretaceous composite nappe stacking, and subsequent continent-continent collision in the latest Cretaceous resulted in folding and thrusting that in the most external part of the Dinarides took place during the Middle Eocene – Oligocene. This extensive last phase of substantial crustal shortening and thickening was associated with flexural foreland basin deposition, resulting in Eo- to early Oligocene syntectonic units. These rocks and older Mesozoic carbonate platform units now form the mountain chain of the external Dinarides. So far, the driving mechanism behind the rock uplift was unknown and it was not clear when the present-day topography formed. Here we show that horizontal marine terraces preserved at elevations of up to 600 m in the external Dinarides are crucial to answer these questions.

We extracted horizontal surfaces, river incision profiles, and the Adriatic and Black Sea catchments from a digital elevation model (DEM). The extracted horizontal surfaces are interpreted as marine terraces because they are degradational, locally preserved in a staircase morphology, neither bedding- nor fault-related, and located close to the present-day Adriatic shoreline. The marine terraces stretch c. 600 km along-strike the entire Dinaric coastal region. Their spatial correlation agrees with the position of a reported positive P-wave tomography anomaly beneath the Dinarides. This up to 180 km deep anomaly correlates also with the thinnest part of the Adriatic lithosphere and the Adriatic-Black Sea drainage divide. The orogen-perpendicular river incisions profiles reveal a symmetric river incision pattern on both sides of the drainage divide. The mean amount of the river incision is equivalent to the mean elevation of the documented marine terraces. All results point to an orogen-wide surface uplift of the Dinarides.

Based on the geological record this post-collisional uplift event can be relatively dated to Oligocene-Miocene (28-17 Ma) and seems to be broadly contemporaneous with the emplacement of igneous rocks with mantle affinity (33-22 Ma) in the internal Dinarides. Previously published geophysical and petrological, as well as the new geomorphological data presented here suggest that the post-collisional reorganization of the Dinarides is attributed to an Oligocene-Miocene mantle delamination event, which results in uplift event affecting the entire Dinarides. We also show that no significant, orogen-scale deformation affected the uplifted Dinarides after the Early Miocene.

How to cite: Balling, P., Grützner, C., Tomljenović, B., Spakman, W., and Ustaszewski, K.: Post‑collisional mantle delamination in the Dinarides implied from staircases of Oligo‑Miocene uplifted marine terraces, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13148, https://doi.org/10.5194/egusphere-egu21-13148, 2021.

Mt. Ivanščica is one of inselbergs in the Internal Dinarides (NW Croatia) in the transitional area with Southern Alps. In this area, NNW-verging Dinaric structures are overprinted by S-verging Alpine structures. Mt. Ivanščica is composed of Mesozoic shallow to deep-marine sedimentary succession of the passive continental margin of Adriatic plate, which was facing the Neotethys ocean, overthrust by ophiolitic mélange. Here, we aim to present new preliminary structural data from pelagic sediments of Ivanščica Mt. in attempt to better understand tectonic history of this part of Internal Dinarides.

Mesozoic succession of Mt. Ivanščica is composed of Triassic clastic, volcanic and carbonate rocks overlain by Upper Triassic to Lower Jurassic shallow-marine carbonates. These are overlain by Jurassic pelagic carbonates and cherts followed by Tithonian−Valanginian pelagic “Aptychus Limestones”. The uppermost part of this succession is composed of Lower Cretaceous Oštrc fm., which conformably overlies the “Aptychus Limestone”. The Oštrc fm. is characterized by turbidites with ophiolitic detritus and represents syn-orogenic deposits presumed as formed in a front of advancing ophiolitic nappe(s).

The focus of our investigation is primarily on structural characteristics of the “Aptichus Limestones” and the Oštrc fm. The character of the contact between the “Aptychus Limestones” and underlying Upper Triassic to Lower Jurassic carbonates is still uncertain. According to Šimunić et al. (1982) “Aptychus limestones” unconformably overlays Triassic carbonates in periclinal geometry, while Babić (1974) suggests continuous condensed pelagic sedimentation throughout the Jurassic. In contrast with previous observations and interpretations, our observations suggest a tectonic contact, characterized by significantly different orientation of bedding and locally marked by fault gauge (clay) seams.

Structural analysis shows numerous gentle to open asymmetric folds in the “Aptychus Limestones” and closed chevron folds and isoclinal folds in overlaying Oštrc fm. Chevron folds and open to gentle asymmetric folds indicate NW vergence in present day orientation with fold axis parallel to the strike of the contact with underlying unit. Although different in shape and size, these folds are likely formed during the same tectonic event while their geometry is controlled by differences in rheological properties. Isoclinal folds occurring exclusively at the contact with ophiolitic mélange are characterized by E-W oriented fold axis and S dipping axial surfaces which is in a contrast with aforementioned folds. Thus, we assume that these folds originated from another, presumably older tectonic event. Bedding in Triassic dolomites uniformly dips towards the SE. Local occurrence of condensed pelagic limestones and radiolarian cherts is interpreted, as rheologically weak horizon ideal to form a décollement that, at least locally, could be interpreted to mark a thrust fault.

Formation of isoclinal folds in the Oštrc fm. and the tectonic contact with ophiolitic mélange is preliminarily attributed to the Aptian-Albian nappe stacking known from the Internal Dinarides. In addition, we assume that the pelagic succession of the “Aptychus Limestones” together with the overlying Oštrc fm. and the ophiolitic mélange are thrusted over the Upper Triassic to Liassic carbonates sometime later, possibly during the final stage of Neotethys closure in the Internal Dinarides.

How to cite: Vukovski, M., Kukoč, D., and Tomljenović, B.: An overview on Structural position of Mesozoic succession of distal Adriatic continental margin on Ivanščica Mt. (NW Croatia), EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5818, https://doi.org/10.5194/egusphere-egu21-5818, 2021.

The Bükk and Uppony Hills (NE Hungary) are two adjacent structural units with correlations to the Northern Dinarides and Inner Western Carpathians (ALCAPA), respectively. These two units are separated by the Nekézseny Fault, which may therefore be considered as a presently displaced segment of the Dinaric-ALCAPA contact zone (Schmid et al. 2008). Along this contact zone, the Bükk-type Permo-Mesozoic formations are thrust over the Paleozoic and Senonian formations of the Uppony Unit (Schréter 1945). Despite of the Nekézseny Fault being a terrain boundary, its structural evolution has not been studied in details. Preliminary structural data suggested multiple faulting events between the latest Senonian and early Miocene (Fodor et al. 2005), however, the initial age of the contact zone has remained uncertain.

In this study a detailed structural analysis was carried out in order to understand the deformation geometry, kinematics and the timing of movements along the Nekézseny Fault. Our preliminary results show that the Nekézseny Fault developed in response to NW-SE shortening. Low-angle fractures within individual pebbles suggest an early (latest Cretaceous or early Paleogene) age for the NW-SE shortening, as pebble fracturing is limited to the early stage of diagenesis and requires soft or semi-consolidated fine-grained matrix.

The top-to-the-NW emplacement of the Bükk over the Uppony Unit was accompanied by the folding of the Senonian conglomerate in the footwall, where a large, almost isoclinal recumbent fold developed due to the estimated several km of displacement along the main contact zone. Despite of the similarity in the shortening directions, the top-to-the-NW shortening certainly post-dates the penetrative S-SE-vergent contractional structures present throughout the Bükk Hills, that are related to the latest Jurassic to Early Cretaceous nappe stacking and subsequent shortening (Csontos 1999). Microtectonic analysis of the Nekézseny Fault Zone proved that the main contact zone is a strongly distorted cataclastite zone, which suggests a late-stage low-temperature deformation. Similarly younger semi-ductile or low-temperature contractional structures (e.g. kink folds) were recognized in several parts of the Bükk Unit, all of which were dated tentatively to the late Cretaceous (Flórián-Szabó & Csontos 2002, Juhász 2020, McIntosh 2014, Koroknai et al. 2008, Scherman 2018). Our observations indicate that the top-to-the-NW displacement was much more extensive than previously thought and incorporated large part of the Bükk Unit. This shows that the top-to-the-NW displacement represents an important deformation phase, which should be integrated into the Mesozoic structural evolution of the Alpine-Dinaric area.

This study was supported by the research founds NKFIH OTKA 113013 and 134873, the ÚNKP-17-2 and ÚNKP-20-3 New National Excellence Program of the Ministry of Human Capacities.

References:
Csontos (1999): Bulletin of the Hungarian Geological Society 129, 4, 611-651.
Fodor et al. (2005): GeoLines 19: 141-161.
Juhász (2020): TDK thesis, ELTE, Budapest.
Koroknai et al. (2008): Journal of Structural Geology 30, 159-176.
McIntosh (2014): PhD thesis, University of Debrecen, Debrecen.
Scherman (2018): MSc thesis, ELTE, Budapest.
Schmid et al. (2008): Swiss Journal of Geosciences 101: 139-183.
Schréter (1945): Annual Report of the Geological Institute of Hungary, 1941-42: 197-237.

How to cite: Oravecz, É., Juhász, D., Götz, A., Kövér, S., Scherman, B., and Fodor, L.: Structural evolution of the Nekézseny Fault – a displaced segment of the Dinaric-ALCAPA contact zone in NE Hungary (Bükk and Uppony Hills), EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11499, https://doi.org/10.5194/egusphere-egu21-11499, 2021.

The Carpathians fold-and-thrust belt results from oblique collision of ALCAPA and Tisza-Dacia plates with the eastern European margin. It formed during the Oligocene and Miocene, propagating laterally from NW to SE as clearly demonstrated by balanced-cross sections (Nakapelyukh et al., 2017; Castellucio et al., 2016; Merten et al., 2010). The coeval development of the foreland basin (Roure et al., 1993) is revealed by an axial transport system that prograded from NW to SE, ultimately supplying sediments to the Black Sea (de Leeuw et al., 2020). However, lacking a regional synthesis and integration of thermochronology data, lateral propagation of exhumation in the orogen has not been demonstrated yet.

 We reconstruct the exhumation history of the entire Carpathians from the Oligocene onwards and link it with the development of the Carpathians foreland basin (CFB) using a source-to-sink approach. We compiled more than 500 apatite and zircon fission-track and (U-Th)/He ages from the literature. This comprehensive database was separated by region (Western, Eastern, and South-Eastern Carpathians) and by tectonic domain (as defined in Schmid et al., 2008). This partitioning allows for the inversion of large datasets, reflects the tectonic complexity of the belt, and avoids spurious spatial correlations (Schildgen et al., 2018). The thermochronology data was inverted using Pecube (Braun et al., 2012) to constrain exhumation rates in a Bayesian approach. We thus obtain estimates of exhumation rates through time along the belt (with their uncertainty) and convert these into bulk  sediment fluxes over time, permitting tracking of sediment routing from the eroding belt to the CFB. Ultimately, these data will be used to unravel deeper geodynamics, including the possible effects of slab detachment on the evolution of the belt and its foreland basin.

Key words: Low-temperature thermochronology, Carpathians, exhumation, source to sink, Pecube inversions.

How to cite: Roger, M., van der Beek, P., de Leeuw, A., and Husson, L.: Diachronous exhumation of the Carpathians from low-temperature thermochronology, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2648, https://doi.org/10.5194/egusphere-egu21-2648, 2021.

Within-plate migration of alkaline basaltic centers is generally related to translation of the lithosphere with a hot spot above a largely stationary mantle plume. Here, we report, for the first time, a hitherto unrecognized migration of small-sized Late Miocene to Early Quaternary alkali-basaltic volcanic centers at the transition from Eastern Alps to the Pannonian Basin. The volcanic centers migrated along the South Burgenland High in a regular sequence over a 95 km distance from NNE to SSW between 11 Ma and 1.7 Ma. The basaltic magmatism was also associated with thinning of the crust and lithosphere and an increase of the thermal gradient as previous studies testify. In detail, three stages of migration are recognized as follows: (i) Stage 1 with a 55 km SSW-toward shift of volcanism between 11 and 5 Ma; (ii) Stage 2 with a 35 km-WSW shift and enlargement of the distribution of volcanic centers between ca. 5 and 3.5 Ma; and (iii) Stage 3 with a S-directed shift of ca. 25 km between 3.5 and 1.7 Ma. We propose that this pattern and mechanism of migration of the volcanic centers along the South Burgenland High resulted from thermally induced progressive thinning of the lithosphere over a mantle plume underneath the ALCAPA block, which was moving from SSW to the NNE between 11 and 1.7 Ma, interrupted by a marked eastward shift between 5 and 3.5 Ma (Stage 2). The Stage 1 shift of alkali-basaltic volcanism can be also observed in the Little Hungarian Plain and South Slovakia volcanic fields, here rotated in a sinistral wrench corridor, and the Stage 2 in the Balaton-Bakony of the Pannonian Basin (north of Lake Balaton). The migration of volcanic centers correlates with orientation and timing of regional shortening phases and inversion of the Alpine-Carpathian-Pannonian system although the overall amount of shortening remains uncertain. The convergence of the Alpine-Carpathian-Pannonian-Dinaric system is driven by the northward motion of the Adriatic microplate. Previous balancing of the shortening of Eastern Alps resulted in an approximately constant convergence rate of ca. 10 mm/yr since 20 Ma ago although the system changed to overall compression since ca. 6 Ma and lower shortening rates can be expected. Interestingly, the reasonable estimated migration rates of volcanic centers are all in a similar order of magnitude between 6.5 and 13.8 mm/yr as the mentioned 10 mm/yr. This is nearly an order of magnitude larger than present-day rates of 1 – 2 mm/yr with E to NE-directed motion, but similar to the global hotspot-fixed reference frame, which implies a ca. northeast-directed motion of >10 mm/yr. Although this could represent an explanation of observed migration of volcanic centers, we exclude it because of the short tracks (maximum 95 km) and small volumes of volcanic products implying rather local plumes in the upper mantle.

Acknowledgements: This work has been supported by the Austrian Science Fund (FWF), grant no. M-1343.

How to cite: Neubauer, F. and Cao, S.: Migration of Late Miocene to Quaternary alkaline magmatism at the Alpine-Pannonian transition area: Significance for deformation rates and decoupling of lithospheric levels, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2015, https://doi.org/10.5194/egusphere-egu21-2015, 2021.

The Pannonian Basin is a continental extensional basin system with various depocentres within the Alpine–Carpathian–Dinaridic orogenic belt. Along the western basin margin, exhumation along the Rechnitz, Pohorje, Kozjak, and Baján detachments resulted in cooling of diverse crustal segments of the Alpine nappe stack (Koralpe-Wölz and Penninic nappes); the process is constrained by variable thermochronological data between ~25–23 to ~15 Ma. Rapid subsidence in supradetachment sub-basins indicates the onset of sedimentation in the late Early Miocene (Ottnangian? or Karpatian, from ~19 or 17.2 Ma). In addition to extensional structures, strike-slip faults mostly accommodated differential extension between domains marked by large low-angle normal faults. Branches of the Mid-Hungarian Shear Zone (MHZ) also played the role of transfer faults, although shear-zones perpendicular to extension also occurred locally.

During this period, the distal margin of the large tilted block in the hanging wall of the detachment system, the pre-Miocene rocks of the Transdanubian Range (TR) experienced surface exposure, karstification, and terrestrial sedimentation. The situation changed after ~15–14.5 Ma when faulting, subsidence, and basin formation shifted north-eastward. Migration of normal faulting resulted in fault-controlled basin subsidence within the TR which lasted until ~8 Ma.

3D thermo-mechanical lithospheric and basin-scale numerical models predict similar spatial migration of the depocenters from the orogenic margin towards the basin center. The reason for this migration is found in the interaction of deep Earth and surface processes. A lithospheric and smaller crustal-scale weak zones inherited from a preceding orogenic structure localize initial deformation, while their redistribution controls asymmetric extension accompanied by the upraising of the asthenopshere and flexure of the lithosphere. Models suggest ~4–5 Myr delay of the onset of sedimentation after the onset of crustal extension and ~150–200 km of shift in depocenters during ~12 Myr. These modeling results agree well with our robust structural and chronological data on basin migration.

Simultaneously with or shortly after depocenter migration, the southern part of the former rift system, mostly near the MHZ, underwent ~N–S shortening; the basin fill was folded and the boundary normal faults were inverted. The style of deformation changed from pure contraction to transpression. The Baján detachment could be slightly folded, although its synformal shape could also be considered a detachment corrugation. Deformation was dated to ~15–14 Ma (middle Badenian) in certain sub-basins while in other sub-basins deformation seems to be continuous throughout the late Middle Miocene from ~15 Ma to ~11.6 Ma.

Another contractional pulse occurred in the earliest Late Miocene, between ~11.6 and ~9.7 Ma while the western part of the TR was still affected by extensional faulting and subsidence. All these contractional deformations can be linked to the much larger fold-and-thrust belt that extends from the Southern and Julian Alps through the Sava folds region in Slovenia. Contraction is still active, as indicated by recent earthquakes in Croatia.

Mol Ltd. largely supported the research. The research is supported by the scientific grant NKFI OTKA 134873 and the Slovenian Research Agency (research core funding No. P1-0195).

How to cite: Fodor, L., Balázs, A., Csillag, G., Dunkl, I., Héja, G., Kelemen, P., Kövér, S., Németh, A., Nyíri, D., Oravecz, É., Selmeczi, I., Trajanova, M., Vrabec, M., and Vrabec, M.: Migration of deformation, subsidence, and basin formation in the SW Pannonian Basin (central Europe) and the change to contractional deformation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13300, https://doi.org/10.5194/egusphere-egu21-13300, 2021.