European Crust Research Articles

Rayleigh wave attenuation and phase velocity maps of the greater Alpine region from ambient noise

We use seismic ambient noise data from 724 publicly available broadband seismic stations across central Europe to create detailed phase velocity and attenuation maps of Rayleigh waves, focusing on short periods down to 3 s. We interpret these maps in terms of the underlying physical processes relevant to the nature of continental crust. Through a regionalized interpretation based on tectonic settings, we highlight the significant role of fluid-filled fractures in the attenuation of surface waves. Our findings indicate a close connection between the time elapsed since the last tectonic activity in the European crust and the attenuation coefficient values. Additionally, we observe a pronounced decrease in attenuation coefficient values at periods below 6 s. The anti-correlation between attenuation coefficient and phase velocity in recently active tectonic regions suggests that fluid-filled fractures are likely the dominant factor governing seismic attenuation in the European crust.

Structural Evolution, Exhumation Rates, and Rheology of the European Crust During Alpine Collision: Constraints From the Rotondo Granite—Gotthard Nappe

AbstractThe rheology of crystalline units controls the large‐scale deformation geometry and dynamics of collisional orogens. Defining a time‐constrained rheological evolution of such units may help unravel the details of collisional dynamics. Here, we integrate field analysis, pseudosection calculations and in situ garnet U–Pb and mica Rb–Sr geochronology to define the structural and rheological evolution of the Rotondo granite (Gotthard nappe, Central Alps). We identify a sequence of four (D1–D4) deformation stages. Pre‐collisional D1 brittle faults developed before Alpine peak metamorphism, which occurred at 34–20 Ma (U–Pb garnet ages) at 590 ± 25°C and 0.9 ± 0.1 GPa. The reactivation of D1 structures controlled the rheological evolution, from D2 reverse mylonitic shearing at amphibolite facies (520 ± 40°C and 0.8 ± 0.1 GPa) at 18–20 Ma (white mica Rb–Sr ages), to strike‐slip, brittle‐ductile shearing at greenschist‐facies D3 (395 ± 25°C and 0.4 ± 0.1 GPa) at 14–15 Ma (white mica and biotite Rb–Sr ages), and then to D4 strike‐slip faulting at shallow conditions. Although highly misoriented for the Alpine collisional stress orientation, D1 brittle structures controlled the localization of D2 ductile mylonites accommodating fast (∼3 mm/yr) exhumation rates due to their weak shear strength (<10 MPa). This structural and rheological evolution is common across External Crystalline Massifs (e.g., Aar, Mont Blanc), suggesting that the European upper crust was extremely weak during Alpine collision, its strength controlled by weak ductile shear zones localized on pre‐collisional deformation structures, that in turn controlled localized exhumation at the scale of the orogen.

Ambient‐Noise Wave‐Equation Tomography of the Alps and Ligurian‐Provence Basin

AbstractTaking benefit of the AlpArray temporary network and permanent networks in W‐Europe, we construct a 3‐D onshore‐offshore velocity model of the crust and upper mantle using ambient‐noise wave‐equation tomography. We use a frequency‐dependent phase traveltime misfit function in an iterative procedure to refine a recent 3‐D Vs model computed from a Bayesian two‐step ambient noise tomography. Observed waveforms consist of vertical‐component noise correlations from 600 broadband stations in the Alps and surroundings, including ocean‐bottom seismometers in the Ligurian sea. We perform 3‐D inversion in the 5–85 s period range. In the long‐period band (20–85 s), an elastic approximation is considered, while in the 5–20 s band, we account for the effect of water layer in the Ligurian sea by applying a fluid‐solid coupling for acoustic‐elastic waveform simulations. The resulting Vs model enhances the shape and contrast of velocity structures, accounting for 3‐D and finite‐frequency effects. It emphasizes the deep sediments of the Ligurian‐Provence basin and focuses on the low‐velocity anomalies of the crust in the W‐Alps. We obtain a high‐resolution Moho depth map covering the Alps and Ligurian sea. In the W‐Alps, this map confirms the deepening of the European crust following the subduction beneath Adria and the existence of major structures such as the Moho jump beneath the external crystalline massifs and shallow depths associated with the Ivrea Body. It provides further constraints on the deep structure beneath the Ligurian‐Provence basin, regarding the lateral and along‐strike crustal‐thickness variations from the oceanic domain to the conjugate margins.

Seismic structure of the Eastern European crust and upper mantle from probabilistic ambient noise tomography

The Eastern European lithosphere is a natural laboratory to study continental formation and evolution through time, comprising Archean continental remnants, Proterozoic rifts and belts, and younger accreted terranes. We investigate the seismic structure of the East European Craton (EEC) crust and uppermost mantle, and the transition from Precambrian to Phanerozoic Europe across the Trans European Suture Zone (TESZ) using probabilistic transdimensional ambient noise tomography. We cross-correlate noise recorded at broadband seismic stations from Eastern, Northern, and Central Europe, remove earthquake signals using continuous wavelet transform, and extract Rayleigh wave phase velocity dispersion curves. We invert these for the highest resolution shear wave velocity model of the Eastern European lithosphere to date, using Markov chain Monte Carlo Bayesian inversion. Our shear wave velocity model exhibits spatial correlation with major tectonic units and bears similarities with active seismic survey profiles in terms of seismic velocity patterns and main discontinuities. The crust thickens across the TESZ boundary and the mantle is seismically faster than beneath younger terranes, consistent with a less dense Precambrian lithosphere in the EEC. The crust and lithosphere beneath the Pannonian region is hyper-extended but the adjacent Transylvanian basin crust shows significant heterogeneity. The Precambrian building blocks of the EEC exhibit contrasting seismic fabrics. The Baltic orogens of Fennoscandia are underlain by uniform crust with a flat Moho, while Sarmatia shows alternating high and low velocity layers and a regional south-dipping crustal boundary from beneath the Ukrainian Shield towards the Crimean Peninsula. The last observation supports a geodynamic style driven by horizontal rather than vertical tectonics, with fundamental implications for the formation and evolution of early continents.

Investigating the Eastern Alpine–Dinaric transition with teleseismic receiver functions: Evidence for subducted European crust

The tectonic structure of the Eastern Alps is heavily debated with successive geophysical studies that are unable to resolve areas of ambiguity (e.g., the presence of a switch in subduction polarity and differing crustal models). In order to better understand this area, we produce a high resolution Moho map of the Eastern Alps based on a dense seismic broadband array deployment. Moho depths were derived from joint analysis of receiver function images of direct conversions and multiple reflections for both the SV (radial) and SH (transverse) components, which enables us to map overlapping and inclined discontinuities. We observe the European Moho to be underlying the Adriatic Moho from the west up to the eastern edge of the Tauern Window. East of the Tauern Window, a sharp transition from underthrusting European to a flat and thinned crust associated with Pannonian extension tectonics occurs, which is underthrust by both European crust in the north and by Adriatic crust in the south. The Adriatic lithosphere underthrusts northward below the Southern Alps and becomes steeper and deeper towards the Dinarides where it dips towards the north-east. Our results suggest that the steep high velocity region in the mantle below the Eastern Alps, observed in tomographic studies, is likely to be of European origin.

Genesis of the Eastern Adamello Plutons (Northern Italy): Inferences for the Alpine Geodynamics

The Corno Alto–Monte Ospedale magmatic complex crops out at the eastern border of the Adamello batholith, west of the South Giudicarie Fault (NE Italy). This complex includes tonalites, trondhjemites, granodiorites, granites and diorites exhibiting an unfoliated structure suggesting passive intrusion under extensional-to-transtensional conditions. Major, minor elements, REE and isotopic analyses and geochemical and thermodynamic modelling have been performed to reconstruct the genesis of this complex. Geochemical analyses unravel a marked heterogeneity with a lack of intermediate terms. Samples from different crust sections were considered as possible contaminants of a parental melt, with the European crust of the Serre basement delivering the best fit. The results of the thermodynamic modelling show that crustal melts were produced in the lower crust. Results of the geochemical modelling display how Corno Alto felsic rocks are not reproduced by fractional crystallization nor by partial melting alone: their compositions are intermediate between anatectic melts and melts produced by fractional crystallization. The tectonic scenario which favored the intrusion of this complex was characterized by extensional faults, active in the Southalpine domain during Eocene. This extensional scenario is related to the subduction of the Alpine Tethys in the Eastern Alps starting at Late Cretaceous time.

Lithospheric transdimensional ambient-noise tomography of W-Europe: implications for crustal-scale geometry of the W-Alps

SUMMARYA full understanding of the dynamics of mountain ranges such as the Alps requires the integration of available geological and geophysical knowledge into a lithospheric-scale 3-D geological model. As a first stage in the construction of this geo-model, we derive a new 3-D shear wave velocity model of the Alpine region, with a spatial resolution of a few tens of kilometres, making it possible to compare with geological maps. We use four years of continuous vertical-component seismic noise records to compute noise correlations between more than 950 permanent broad-band stations complemented by ∼600 temporary stations from the AlpArray sea-land seismic network and the Cifalps and EASI linear arrays. A specific pre-processing is applied to records of ocean–bottom seismometers in the Liguro-Provençal basin to clean them from instrumental and oceanic noises. We first perform a 2-D transdimensional inversion of the traveltimes of Rayleigh waves to compute group-velocity maps from 4 to $150\, \mathrm{ s}$. The data noise level treated as an unknown parameter is determined with a Hierarchical Bayes method. A Fast Marching Eikonal solver is used to update ray path geometries during the inversion. We use next the group-velocity maps and their uncertainties to derive a 3-D probabilistic Vs model. The probability distributions of Vs at depth and the probability of presence of an interface are estimated at each location by exploring a set of 130 million synthetic four-layer 1-D Vs models. The obtained probabilistic model is refined using a linearized inversion. Throughout the inversion for Vs, we include the water column where necessary. Our Vs model highlights strong along-strike changes of the lithospheric structure, particularly in the subduction complex between the European and Adriatic plates. In the South-Western Alps, our model confirms the existence of a low-velocity structure at $50-80\, \mathrm{ km}$ depth in the continuation of the European continental crust beneath the subduction wedge. This deep low-velocity anomaly progressively disappears towards the North-Western and Central Alps. The European crust includes lower crustal low-velocity zones and a Moho jump of $\sim \, 8-12$ km beneath the western boundary of the External Crystalline Massifs of the North-Western Alps. The striking fit between our Vs model and the receiver function migrated depth section along the Cifalps profile documents the reliability of the Vs model. In light of this reliability and with the aim to building a 3-D geological model, we re-examine the geological structures highlighted along the Cifalps profile.

Orogen‐Parallel Migration of Exhumation in the Eastern Aar Massif Revealed by Low‐T Thermochronometry

AbstractNew and published (U‐Th)/He data on zircon, apatite, and zircon fission track ages constrain the thermal overprint and cooling history of the eastern Aar Massif, Switzerland. The timing and pattern of cooling is in agreement with independent kinematic and age constraints from exposed shear zones. This suggests that the cooling ages mainly reflect exhumation and that long‐term exhumation‐dynamics were mainly controlled by crustal‐scale tectonic processes. Results of a statistical inverse model reveal significant diachrony in the timing of exhumation in the along‐strike direction. Maximum exhumation rates (1 mm/yr) were initially located in the central Aar Massif (from 22 to 10 Ma), then gradually migrated to the east between 10 Ma and present, while the central Aar Massif continued to exhume at slower rates (0.5 mm/yr). The diachrony in the timing of exhumation may be explained by lateral variations in the inherited thickness or the density of the accreted European crust. We attribute the increase in exhumation rates between 2 Ma and present to enhanced glacial erosion. Nevertheless, the post 2 Ma exhumation pattern reflects a continuation of noncylindrical massif “growth” in the eastward orogen‐parallel direction. This indicates that—although at slow rates—thick‐skinned and buoyancy‐driven compressional deformation, likely enhanced by the presence of easily erodible flysch units at the surface, might still be ongoing especially in the eastern Aar Massif. Noncylindrical massif‐growth is likely to also affect other External Crystalline Massifs or orogens, but may be overlooked because studies often focus on single orogen‐perpendicular transects.

The SWATH-D Seismological Network in the Eastern Alps

Abstract The SWATH-D experiment involved the deployment of a dense temporary broadband seismic network in the Eastern Alps. Its primary purpose was enhanced seismic imaging of the crust and crust–mantle transition, as well as improved constraints on local event locations and focal mechanisms in a complex part of the Alpine orogen. The study region is a key area of the Alps, where European crust in the north is juxtaposed and partially interwoven with Adriatic crust in the south, and a significant jump in the Moho depth was observed by the 2002 TRANSALP north–south profile. Here, a flip in subduction polarity has been suggested to occur. This dense network encompasses 163 stations and complements the larger-scale sparser AlpArray seismic network. The nominal station spacing in SWATH-D is 15 km in a high alpine, yet densely populated and industrialized region. We present here the challenges resulting from operating a large broadband network under these conditions and summarize how we addressed them, including the way we planned, deployed, maintained, and operated the stations in the field. Finally, we present some recommendations based on our experiences.

Evidence for radial anisotropy in the lower crust of the Apennines from Bayesian ambient noise tomography in Europe

SUMMARYProbing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (>30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH > VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region.

The 3D thermal field across the Alpine orogen and its forelands and the relation to seismicity

Abstract Temperature exerts a first order control on rock strength, principally via thermally activated creep deformation and on the distribution at depth of the brittle-ductile transition zone. The latter can be regarded as the lower bound to the seismogenic zone, thereby controlling the spatial distribution of seismicity within a lithospheric plate. As such, models of the crustal thermal field are important to understand the localisation of seismicity. Here we relate results from 3D simulations of the steady state thermal field of the Alpine orogen and its forelands to the distribution of seismicity in this seismically active area of Central Europe. The model takes into account how the crustal heterogeneity of the region effects thermal properties and is validated with a dataset of wellbore temperatures. We find that the Adriatic crust appears more mafic, through its radiogenic heat values (1.30E-06 W/m3) and maximum temperature of seismicity (600 °C), than the European crust (1.3–2.6E-06 W/m3 and 450 °C). We also show that at depths of

The Maggia nappe: an extruding sheath fold basement nappe in the Lepontine gneiss dome of the Central Alps

The Lepontine gneiss dome represents a unique region of the Central Alps where Oligocene–Miocene amphibolite facies grade rocks and fold nappes of the deepest tectonic level of the Alpine orogenic belt are exposed in a tectonic window. The Cenozoic structures of the Maggia nappe reveals a giant tens of kilometre-scale tubular fold structure that cross-cut through the surrounding lower Penninic nappes from its root situated in the southern steep-belt of the Alps near Bellinzona. The Mesozoic sedimentary cover of the Maggia nappe is typical for the Helvetic stratigraphic domain. The age of formation of the lower Penninic fold-nappes by ductile detachment of the upper European crust during its underthrusting below the higher Penninic and Austroalpine orogenic lid and Adriatic indenter was estimated between 40 and 30 Ma. Maximal pressures of 8–9 kbars and temperatures of 600–700 °C were attended during and after the nappe emplacement some 30–22 Ma ago. The Maggia and surrounding nappes are crosscut by the isograds of the Barrovian regional metamorphism.

Refining the thermal structure of the European lithosphere by inversion of subsurface temperature data

Abstract We present an updated thermal model of the European lithosphere based on a new stochastic modeling work flow. We developed this work flow to estimate subsurface temperatures from site- to regional-scale, up to the depth of the lithosphere-asthenosphere boundary (LAB). Our model is composed of four layers, consisting of sediments, upper crust, lower crust and lithospheric mantle. We asigned thermal properties, including radiogenic heat production and temperature- and pressure-dependent bulk thermal conductivity on the base of broad-scale lithological variation within the European crust. We corrected thermal properties with a 1D steady-state temperature approximation, assuming only vertical heat flow. Using these corrected thermal properties, we calculated the 3D thermal field with a conjugate-gradient method, assuming fixed temperatures at the surface and at the base of the lithosphere. To obtain more robust results for our thermal model, we applied data assimilation, aiming at consistency between temperature and heat flow observations and tectonic model predictions. We used an Ensemble Smoother Multiple with Data Assimilation (ES-MDA) method to update prior estimates of thermal properties and the thermal field with temperature data. We calibrated our European thermal model with available regional thermal models due to the present lack of a unified dataset with public borehole temperature measurements. A large dichotomy is observed in the model along the Trans-European Suture Zone (TESZ). Northeast of the TESZ, geothermal gradients up to 10 km depth are mainly below 20 °C km−1. Southwest of the TESZ, gradients range from 20 °C km−1 near the Adriatic coast to more than 50 °C km−1 in volcanically and tectonically active regions. We show that the large scale thermal structure is locally and regionally perturbed by non-conductive heat transfer and affected by transient effects.

High-resolution surface wave tomography of the European crust and uppermost mantle from ambient seismic noise

Taking advantage of the large number of seismic stations installed in Europe, in particular in the greater Alpine region with the AlpArray experiment, we derive a new high-resolution 3-D shear wave velocity model of the European crust and uppermost mantle from ambient-noise tomography. The correlation of up to 4 yr of continuous vertical-component seismic recordings from 1293 broad-band stations (10 • W-35 • E, 30 • N-75 • N) provides Rayleigh wave group velocity dispersion data in the period band 5-150 s at more than 0.8 million virtual source-receiver pairs. 2-D Rayleigh wave group velocity maps are estimated using adaptive parametrization to accommodate the strong heterogeneity of path coverage. A probabilistic 3-D shear wave velocity model, including probability densities for the depth of layer boundaries and S-wave velocity values, is obtained by nonlinear Bayesian inversion. A weighted average of the probabilistic model is then used as starting model for the linear inversion step, providing the final V s model. The resulting S-wave velocity model and Moho depth are validated by comparison with previous geophysical studies. Although surface wave tomography is weakly sensitive to layer boundaries, vertical cross-sections through our V s model and the associated probability of the presence of interfaces display striking similarities with reference controlled-source seismology (CSS) and receiver function sections across the Alpine belt. Our model even provides new structural information such as an ∼8 km Moho jump along the CSS ECORS-CROP profile that was not imaged by the reflection data due to poor penetration across a heterogeneous upper crust. Our probabilistic and final shear wave velocity models have the potential to become new reference models of the European crust, both for crustal structure probing and geophysical studies including waveform modelling or full-waveform inversion.

A three-dimensional model of the Pyrenees and their foreland basins from geological and gravimetric data

We construct a three-dimensional geological model of the Pyrenees and their foreland basins with the Geomodeller. This model, which accounts for different sources of geological and geophysical informations, covers the whole Pyrenees, from the Atlantic Ocean to the Mediterranean Sea, and from the Iberian range to the Massif Central, down to 70 km depth. We model the geological structure with a stratigraphic column composed of a superposition of layers representing the mantle, lower, middle, and upper crusts. The sedimentary basins are described by two layers which allow us to make the distinction between Mesozoic and Cenozoic sediments, which are characterized by markedly different densities and seismic velocities. Since the Pyrenees result from the convergence between the Iberian and European plates, we ascribe to each plate its own stratigraphic column in order to be able to model the imbrication of Iberian and European crusts along this fossile plate boundary. We also introduce two additional units which describe the orogenic prism and the water column in the Bay of Biscay and in the Mediterranean Sea. The last ingredient is a unit that represents bodies of shallow exhumed and partly serpentinized lithospheric mantle, which are assumed to produce the positive Bouguer gravity anomalies in the North Pyrenean Zone. A first 3D model is built using only the geological information coming from geological maps, drill-holes, and seismic sections. We use the potential field method implemented in Geomodeller to interpolate these geological data. This model is then refined in order to better explain the observed Bouguer anomalies by adding new constraints on the main crustal interfaces. The final model explains the observed Bouguer anomalies with a standard deviation less than 3.4 mGal, and reveals anomalous deep structures beneath the eastern Pyrenees.

The Cenerian orogeny (early Paleozoic) from the perspective of the Alpine region

In the Alps, relicts of pre-Variscan basement are composed of metagreywackes and metapelites (partly migmatic) with intercalated amphibolites and sheets of Cambro–Ordovician peraluminous metagranitoids. Such gneiss terranes are the result of an orogenic type, which was globally widespread in early Paleozoic times. It caused the formation of several 100 km wide cratonized subduction–accretion complexes (SACs) hosting peraluminous arcs at the periphery of Gondwana. “Cenerian orogeny” is a newly suggested term for these early Paleozoic events, which culminate in the Ordovician. The justification for a separate name is given by three characteristics, which are significantly different compared to the Cadomian, Caledonian and Variscan orogenies: the age, the paleogeographic position and the tectonic setting. Other parts of the southern and central European crust might also have been generated by the cratonization of peri-Gondwanan SACs during the Cenerian orogeny.

Pre-Cadomian to late-Variscan odyssey of the eastern Massif Central, France: Formation of the West European crust in a nutshell

The East Massif Central (EMC), France, is part of the internal zone of the Variscan belt where late Carboniferous crustal melting and orogenic collapse have largely obliterated the pre- to early-Variscan geological record. Nevertheless, parts of this history can be reconstructed by using in-situ U-Th-Pb-Lu-Hf isotopic data of texturally well-defined zircon grains from different lithological units. All the main rock units commonly described in the EMC are present in the area of Tournon and include meta-sedimentary and meta-igneous rocks of the Upper Gneiss Unit (UGU) and of the Lower Gneiss Unit (LGU), as well as cross-cutting Variscan granitoid dikes and a heterogeneous granite coring the major Velay dome. Herein we demonstrate that the UGU and the LGU have markedly distinct zircon records. The results of this study are consistent with deposition of the protoliths of the paragneisses within a back-arc basin that was located adjacent to the Arabian-Nubian shield and/or the Saharan Metacraton during the late Ediacaran and collected detritus from the Gondwana continent. At ~545Ma some of these sedimentary rocks were affected by a first melting event that formed the protoliths of the LGU orthogneisses, those of which subsequently remelted at ca. 308Ma to form the Velay granite-migmatite dome. Protoliths of the UGU result mainly from a bimodal rift-related magmatism at ~480Ma, corresponding to melting of the Ediacaran sediments and depleted mantle. Zircon rims from the UGU additionally provide evidence for a metamorphic/migmatitic overprint during the Lower Carboniferous (~350–340Ma). Finally, several generations of granite dikes of which inherited zircons display characteristics of both the UGU and the LGU were protractedly emplaced from ~322Ma to ~308Ma, the youngest of which being coeval with the formation of the Velay dome. Our data further show that the vast majority of crustal material ultimately involved in the Variscan orogeny, which forms the present-day basement in the EMC, was derived from a sedimentary mixture of various components from the Gondwana continent deposited in Ediacaran times, with no evidence for the involvement of an older autochthonous crust.

The Iberian Plate: myth or reality?

The plate tectonics theory generally leads us to consider that Iberia was an independent plate separated from Europe by the North Pyrenean Fault (NPF). The NPF has been commonly interpreted as a transform fault associated with a huge counterclockwise transverse and rotational movement that allowed the opening of the Bay of Biscay and the relative eastward motion of Iberia during the Mesozoic. According to some interpretations, this movement may have generated an interplate gap several hundreds of km wide, which led to the creation of an oceanic crust during the Late Jurassic and Early Cretaceous. However, field studies recently carried out in the Pyrenees do not support these interpretations. The North Pyrenean Fault (NPF) of Tertiary age is observed in the central and eastern Pyrenees, where pioneering researchers defined it as separating the North Pyrenean Zone from the Axial Zone. However, this fault cannot be identified in the western part of the range to the west of the Ossau valley. Consequently, the geodynamic evolution of Iberia has always been dependent on Europe, especially during the failed oceanic rifting in the Mid-Cretaceous. Indeed, during this period, a central zone of crustal thinning occupied by turbiditic basins separated the European from the Iberian continental crust, with a very localized mantle exhumation found only in the Mauleon basin. Therefore, far from being an interplate range, the Pyrenees can neither be considered as an intraplate unit. We can define this orogenic belt as resulting from the “Tertiary tectonic inversion of a Mid-Cretaceous rift system”. According to this new interpretation, Iberia would not have been an isolated plate but represented an unstable, outlying part of Europe. Rather than displaying the features of a rigid lithospheric unit with well-defined boundaries, Iberia grouped together different crustal blocks undergoing specific movements at particular times. During the Mesozoic, normal, reverse or strike-slip displacements along their boundary faults generated several en échelon basins (Bilbao, Logroño, Soria and Maestrazgo in Spain; Parentis, Arzacq and the North Pyrenean Flysch Trough in France) whose diachronous development accounts very well for the opening of the Bay of Biscay and the relative W-E sinistral movement of the Iberian crust with respect to the European crust. In this way, we note a more marked paleogeographic shift in the western Basque region than in eastern Catalonia. Therefore, the Pyrenees were not generated by an “interplate collision”, but merely reflect the confrontation of two distended but continuous continental domains during the Tertiary, leading to the incipient underthrusting of the Iberian crust under the European crust. This zone of confrontation/convergence does not correspond in any way to the NPF, but rather comprises a complex imbricated structure revealed by a jump in the Moho, both in the Pyrenees and in the Cantabrian belt, which has recently become known as the “North Iberian Fault”.

The thrust zone of the Ligurian Penninic basal contact (Monte Fronté, Ligurian Alps, Italy)

ABSTRACTThis paper presents a geological map and interpretative cross-sections, which illustrate the structure of a major fault zone of the Alps, that is, the Ligurian segment of the Penninic Basal Contact (PBC), which led to the emplacement of the orogenic wedge onto the European crust. Chaotic deposits, whose origin is debated, characterize the footwall of this complex thrust zone. The lower part of the sequence containing chaotic deposits consists of low deformed paraconglomerates and exotic megablocks embedded with turbidites. Conversely, the highly deformed upper part of the sequence englobes fragments of substratum-derived succession and is bounded by thrust planes. The nature of these chaotic deposits suggests an origin by gravitational processes related to the unstable front of the advancing wedge associated with offscraping of tectonic slices during the thrusting of allochtonous nappes along the PBC thrust zone.

Role of crustal assimilation and basement compositions in the petrogenesis of differentiated intraplate volcanic rocks: a case study from the Siebengebirge Volcanic Field, Germany

The Siebengebirge Volcanic Field (SVF) in western Germany is part of the Cenozoic Central European Volcanic Province. Amongst these volcanic fields, the relatively small SVF comprises the entire range from silica-undersaturated mafic lavas to both silica-undersaturated and silica-saturated differentiated lavas. Owing to this circumstance, the SVF represents a valuable study area representative of intraplate volcanism in Europe. Compositions of the felsic lavas can shed some new light on differentiation of intraplate magmas and on the extent and composition of potential crustal assimilation processes. In this study, we provide detailed petrographic and geochemical data for various differentiated SVF lavas, including major and trace element concentrations as well as Sr–Nd–Hf–Pb isotope compositions. Samples include tephriphonolites, latites, and trachytes with SiO2 contents ranging between 53 and 66 wt%. If compared to previously published compositions of mafic SVF lavas, relatively unradiogenic 143Nd/144Nd and 176Hf/177Hf coupled with radiogenic 87Sr/86Sr and 207Pb/204Pb lead to the interpretation that the differentiated volcanic rocks have assimilated significant amounts of lower crustal mafic granulites like the ones found as xenoliths in the nearby Eifel volcanic field. These crustal contaminants should possess unradiogenic 143Nd/144Nd and 176Hf/177Hf, radiogenic 87Sr/86Sr, and highly radiogenic 207Pb/204Pb compositions requiring the presence of ancient components in the central European lower crust that are not sampled on the surface. Using energy-constrained assimilation–fractional crystallisation (EC-AFC) model calculations, differentiation of the SVF lithologies can be modelled by approximately 39–47 % fractional crystallisation and 6–15 % crustal assimilation. Notably, the transition from silica-undersaturated to silica-saturated compositions of many felsic lavas in the SVF that is difficult to account for in closed-system models is also well explained by such amounts of crustal assimilation.