EDP Sciences logo
Submit your paper
Search
Menu
Home All issues Volume 195 (2024) BSGF - Earth Sci. Bull., 195 (2024) 7 Full HTML
Issue
BSGF - Earth Sci. Bull.
Volume 195, 2024
Special Issue Some applications of LA-ICP-MS U-Pb geochronology: A tribute to Jean Louis Paquette
Article Number 7
Number of page(s) 18
DOI https://doi.org/10.1051/bsgf/2024003
Published online 23 May 2024

© J. Berger et al., Published by EDP Sciences 2024

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

The West European Variscan belt (Fig. 1) is the result of the convergence between two supercontinents, Laurussia and Gondwana, ultimately leading to their collision starting in the Early Carboniferous. This convergence is chiefly accommodated by the subduction of the Rheic Ocean, which separated the Avalonian terrains from the Gondwanan margin. However, the number, palaeogeography, lifetime and even the existence of other oceanic domains are poorly constrained and these uncertainties lead to various geodynamic models of pre-Carboniferous convergence and subductions (compare Kroner and Romer, 2013 and Franke et al., 2017; or for the western branch, Lardeaux et al., 2014 and Ballèvre et al., 2014). These models (see Vanderhaeghe et al., 2020 for a recent review) can be subdivided into those involving (i) one single ocean (generally the Rheic, see von Raumer and Stampfli, 2008; Kroner and Romer, 2013), (ii) several (up to three) oceanic basins (Matte, 2001; Franke et al., 2017), (iii) monocyclic continuous subduction zones (see Matte 1986, 2001, 2007; Ballèvre et al., 2014) or (iv) polycyclic subduction-collision models involving northward subduction of an oceanic lithosphere until the Silurian-Lower Devonian (generating the so-called eo-variscan event/suture) and southward subduction of another oceanic domain during Middle Devonian to Lowermost Carboniferous times (Faure et al., 2005, 2008, 2009; Lardeaux et al., 2014). More specifically, the identification of the subducted ocean(s) located along the southern Variscan domain is debated. It is referred to as the Galicia, (South) Brittany and Massif central ocean (or a combination of these appellations; Matte, 1986; Matte et al., 2001, 2007), the Rheic ocean (von Raumer and Stampfli, 2008; Nance et al., 2010; Kroner and Romer, 2013), or a branch of the Rheic ocean (Paris and Robardet, 1998) or the Saxo-thuringian ocean (Franke et al., 2017).

Considering that paleo-geographic and paleo-kinematic data are less robust for Paleozoic orogens compared to Cenozoic ones (see for example Jolivet et al., 2021), the evolution of pre-orogenic subductions is strongly based on the sedimentological, petrological and tectonic records acquired by either high pressure metamorphic units or supra-subduction complexes. One of the key arguments to favor/support a particular model is the age of the high pressure event(s). Recent dating performed on high pressure rocks (see Paquette et al., 2017; Lotout et al., 2018, 2020; de Hoÿm de Marien et al., 2023) suggest that eo-variscan 430-390 Ma ages obtained in the past for the HP metamorphism in the French Massif Central (see Faure et al., 2009) are artifacts and that most if not all the eclogites and related HP units could be related to a single Upper Devonian-Lowermost Carboniferous high pressure event as it is generally recognized in the Armorican massif (Bosse et al., 2000, 2005; Ballèvre et al., 2014). There is less ambiguity on the age of supra-subduction magmatism spanning from 379 to 355 Ma for calc-alkaline plutons (Bernard-Griffith et al., 1985; Cuney et al., 1993; Bertrand et al., 2001; Pin and Paquette, 2002; Fig. 1) and Middle Devonian to Lowermost Carboniferous for basaltic and felsic lavas with arc and back-arc affinities (Sieder and Ohnenstetter, 1986; Pin and Paquette, 1997; Pin and Paquette, 2002; Fig. 1). Published zircon U-Pb ages indicates that calc-alkaline plutons in the French Massif Central emplaced within about 30 Myr near the Devonian-Carboniferous transition (Bernard-Griffiths et al., 1985; Pin and Paquette, 2002). However, this time span is sometimes interpreted as covering two successive geological events considering that the plutonic rocks are locally sheared and recrystallized (Faure et al., 1997; Bellot and Roig, 2007). These authors suggest that igneous crystallization of diorites and tonalites occurred between 380−360 Ma while the 360−355 Ma ages represent deformation and recrystallization during the so-called D2 tectono-metamorphic event (Faure et al., 2009), despite the fact that original geochronological studies interpreted all the ages as dating igneous emplacement.

This contribution focuses on a suite of calc-alkaline plutons called the Tonalitic Line of Limousin and grossly aligned on a NW-SE strike from Vendée (SE Amorican massif) to SW French Massif central. We will discuss the nature of Late Devonian magmatism in Western French Massif central (Fig. 1) in the light of new ages and petrological data combined to published results. We will also discuss the implications for the pre-collisional evolution of the West European Variscan belt considering that the calc alkaline plutons could represent a single short-lived magmatic pulse preceding the stacking of the main litho-tectonic units.

thumbnail Fig. 1

Map of the Variscan belt in France showing the main Devonian-Early Carboniferous magmatic markers for subduction (Variscan zones modified after Ballèvre et al., 2014).

2 Geological setting

The French Massif central belongs to the internal units of the Moldanubian zone of the Western Variscan belt (Faure et al., 2009; Lardeaux et al., 2014; Martinez-Catalan et al., 2021; Fig. 1). This belt formed after Cambro-Ordovician rifting of the Northern Gondwanan margin and opening of (an) oceanic domain(s), an event is coeval with the formation of large volume of predominantly felsic magmatic rocks now forming Furongian-Ordovician orthogneiss in metamorphic units. The Devonian convergence, accommodated by subduction between Laurussia and Gondwana, is evidenced by numerous eclogites whose age have been re-evaluated recently at Middle Devonian to Lowermost Carboniferous times (Lotout et al., 2018, 2020; Benmammar, 2021; de Hoÿm de Marien et al., 2023). Final collision of continental blocks occurred during the Lower Carboniferous and triggered a massive phase of syn- to late-orogenic, dominantly silicic magmatism derived from crustal and mantle sources and lasting from about 360 Ma to 300 Ma (Vanderhaeghe et al., 2020).

Middle to Upper Devonian igneous rocks supposedly emplaced in a supra-subduction setting (Fig. 1) form isolated massifs in the internal part of the Variscan belt and especially in the French Massif Central, the neighboring South Armorican zone, the Vosges (Skrzypek et al., 2012) and Paleozoic crystalline massifs of the Alps (Bertrand et al., 2000; Guillot and Ménot, 2009; Guillot et al., 2012; Bergomi et al., 2017). In the Armorican massif, Devonian subalkaline basaltic to rhyolitic rocks occur in low grade Upper allochthonous units in Vendée (basalt and andesite of La Meilleraie interpreted either as arc or back arc magmatism: Thiéblemont and Cabanis, 1986; Ducassou et al., 2011) and at St Georges sur Loire (undated but assumed to be Silurian or Devonian; Cartier and Faure, 2004; Ducassou et al., 2011). In the French Massif central Fammenian to Tournaisian supra-subduction zone magmatism is documented in Morvan (basalt and andesites of the Somme Unit; Pin and Paquette, 2002) and in the so-called Brévenne ophiolite (bimodal volcanic and plutonic basaltic-rhyolitic associations of the Brévenne Unit; Sieder and Ohnenstetter, 1986; Pin and Paquette, 1997; Leloix et al., 1999). However, the most striking expression of this Middle to Late Devonian magmatic episode is the Tonalite Belt of Limousin (Didier and Lameyre, 1971) underlined by diorite and tonalitic plutons extending from Vendée to southwestern French Massif Central (Bernard Griffith et al., 1985; Pfeiffer, 1986; Shaw et al., 1993; Cuney et al., 2001a, 2001b; Pin and Paquette, 2002; Figs. 1 and 2a ). The present study focuses on Western Massif Central plutons (Fig. 2a) that intruded either into the Upper Gneiss Unit (UGU) made of migmatitic paragneiss with numerous lenses of retrogressed eclogites (Pfeiffer 1986) or into the uppermost amphibolite facies Leyme unit for the Figeac and St Céré massifs representing the southernmost occurrences of these plutonic associations (Pfeiffer, 1985, 1986; Duguet et al., 2007). Pin and Paquette (2002) also described two similar igneous complexes in the Northern French Massif Central (Aydat and La Marche plutons). It is worth to note that diorite-tonalite plutons are never found in the greenschist facies metapelites of the Parautochtonous nor in the amphibolite facies orthogneiss and metapelites (with very rare eclogitic lenses) of the Lower Gneiss Unit (LGU).

A representative example of the Tonalitic belt in Limousin is the St Féréole-Tulle massif (Fig. 2b Peiffer, 1986). It has a curved morphology parallel to the tectonic contact between the LGU and UGU (Floc’h, 1983; Bellot and Roig, 2007) and intrudes the migmatitic paragneiss and retrogressed eclogites belonging to the UGU (Fig. 2b). Pfeiffer (1986) describe different lithologies from hectometer to kilometer-size hornblendite and hornblende gabbro clusters to diorite and tonalite forming the main mass of plutons. The textural fabrics vary from igneous with euhedral plagioclase, subhedral amphibole and interstitial quartz and biotite to a deformed porphyroclastic gneiss recrystallized under solid state. Metamorphic foliations of gneissic diorites are subparallel to those measured in host gneisses (Pfeiffer, 1986), i.e. fabrics that are related to the D2 deformation event (Bellot and Roig, 2007; Faure et al., 2009). The bulk structure of the magmatic body resembles an unrooted laccolith transposed into the foliation of host paragneisses (Pfeiffer, 1986).

Available geochronological data (U-Pb zircon by ID-TIMS or SHRIMP) for the emplacement of diorites and tonalites range from 379 to 355 Ma (Bernard Griffiths et al., 1985; Cuney et al., 1993; Bertrand et al., 2001; Pin and Paquette, 2002). All authors converge to interpret these dates as representing igneous crystallization ages.

thumbnail Fig. 2

Geological maps showing the location of calc-alkaline plutons and their relation with tectono-metamorphic units. (a) Simplified geological map of the Limousin area adapted from published BRGM maps, the position of studied plutons is indicated together with samples labels. (b) Zoom on the St-Féréole-Tulle massif located within UGU but near the LGU/UGU contact. Foliation trajectories from Bellot and Roig (2007).

3 Sampling and analytical methods

Seven massifs were sampled from the central Limousin to the Lot Valley in southwestern Massif Central (Tab. 1, Fig. 2a), including, from NW to SE (Pfeiffer, 1987 and Fig. 2): Burgnac (sample labels: BIA, FRO), Saint Jean de Ligoure (SJL, SHB, SPL), Saint Julien le Vendômois (SAJ), Tulle-Saint Féréole (JAM, TUL, MAR), Beaulieu-sur-Dordogne (DBD, BEA), Saint-Céré (STC) and Capdenac (FIG, DCAP). They belong to the Upper Gneiss Unit, except for St-Céré and Capdenac which were ascribed to the amphibolite-facies Leyme unit structurally above the UGU.

Mineral chemical data (Tab. S1) were acquired with a SXFive Cameca electron microprobe hosted by the “Centre de microcaractérisation Raimond Castaing” in Toulouse, France. The beam voltage and current were set at 15 kV and 10 nA, respectively and calibration was done with a set of natural silicate and oxide standards. Major and trace-elements on pulverized bulk rocks were analyzed at the Central Analytical Facilities (CAF) of the Stellenbosch University, South Africa. Rock powders were dried, mixed with lithium metaborate and fused to produce glass beads. Major elements were measured by XRF while trace-elements were analyzed by LA-ICP-MS on the same glass beads. Relative standard deviations between measured and reference values for international standards (BE-N, BHVO, JG-1) analyzed by XRF is around 20% for oxides <0.1 wt%, below 5% for oxides <5 wt% and below 1% for oxides >10 wt%. Relative accuracy on trace-elements (monitored with secondary standard BHVO) is typically around 0.2−10%, the highest values representing either some specific elements (Nb, Ta, Pb) or low concentrations (<1 μg/g).

Nd-Sr isotopic data was obtained at GET-OMP according to the following procedure: 100 mg of whole rock powder was weighed in a teflon beaker and dissolved in a mixture HF/HNO3/HCl 2:1:1. Nd and Sr were extracted from the matrix using a combination of Sr-Spec, Thru-spec and Ln-Spec resins. An equivalent of 500 ng Sr and 150 ng Nd were run on a Thermo Scientific Triton+ mass spectrometer. NBS987 and La Jolla isotopic standards were regularly run during the measurements. Standard reproducibility for 143Nd/144Nd ratio is 0.510858 ± 18 (n = 50) for La Jolla and 0.710248 ± 35 (n = 70) for 87Sr/86Sr ratio of NBS987. However, part of the samples (TUL1, DBD, SAJ12A, SJL, AUB1) were run following Li et al. (2010) procedure, i.e. with Sm isobaric corrections. For those, Sm-corrected 145Nd/144Nd ratios are indicated and Sm-doped JNdi standards were run, they give a 143Nd/144Nd ratio of 0.512108 ± 8 and 145Nd/144Nd ratio of 0.348406 ± 5. Typical procedural blanks are 50 pg for Nd and 150 pg for Sr.

Zircons were separated with routine magnetic and heavy-liquid methods and mounted in epoxy. They were imaged with cold cathode optical cathodoluminescence at the University of Mons using a Cambridge Image Technology model 8200 Mk5 system. U-Pb analyses were done at Géosciences Montpellier, France (platform AETE-ISO of the OSU OREME) with a Teledyne G2 excimer laser probe coupled to a ThermoFinnigan Element XR high-resolution ICP-MS. The pre-ablation blank signal was acquired during 15s followed by about 30 s of ablation (spot size of 25 μm, repetition rate of 4 Hz, fluence of 6 J/cm2) and signal acquisition (see Bruguier et al., 2017 for further details). Pb/U and Pb/Pb isotopic ratios were calibrated against primary reference zircon G91500 (Wiedenbeck et al., 1995). Secondary standards GJ-1 (slightly discordant with an upper intercept TIMS age at 608.5 ± 1.5 Ma; Jackson et al., 2004) and Plesovice (concordant TIMS age at 337.1 ± 0.4 Ma; Sláma et al., 2008) were analyzed every five unknowns and yielded a concordia age of 604 ± 3 Ma (n = 15) and 339 ± 2 Ma (n = 13), respectively. Every zircon sample mount was analyzed with 25 laser spots. Only analyses with regular ablation spectra (i.e. similar in trend to standards, with gently decreasing counts per second and gently increasing Pb/U isotopic ratios during ablation) were selected for further processing. Irregular spectra characterized by abrupt drop(s) or increase(s) in signal (probably due to ablation on fractures or mineral inclusions) were systematically discarded. Raw data were reduced with the software Glitter. For spots yielding regular ablation spectra, the signal was integrated from about 5 s after laser ignition to about 3 s before laser shut off. Reduced and calibrated results were plotted and calculated with IsoplotR (Vermeesh, 2018). Two uncertainties are quoted with Concordia ages, following the recommendation of Horstwood et al. (2016). The first one is based on weighted mean statistic and is given by IsoplotR along with MSWD for concordance and equivalence. The second one (written into brackets) integrates the excess variance (1.6%) calculated on the basis of long term analysis of a secondary reference material (here GJ1) in the laboratory, including the uncertainties on primary reference material (here 91500) and also the errors on decay constants. Note that the obtained dates will be discussed with the second uncertainty, this allow better comparison of results obtained in different laboratories.

Table 1

Location, name, structure and modal proportion of samples used in this study.

4 Petrography, mineral chemistry and thermobarometry

Hornblendites have been sampled in the Tulle massif (Tab. 1, Fig. 2, samples TUL3 and JAM10). They consist of subhedral coarse (up to 1 cm wide) pargasite (98−99 vol.%) showing minute brown lamellar Schiller exsolutions of Fe-Ti oxides with interstitial anhedral altered plagioclase (∼1 vol.%) and ilmenite often included in amphibole (Fig. 3a, Tab. 1). The latter has high Mg# (60−63), low Si (6.05−6.30 apfu) and variable Ti (0.11−0.22 apfu) (Fig. 4). All plagioclase grains are altered into epidote-albite associations and their composition could not be determined.

Hornblende-gabbro samples in the Tulle massif (TUL1 and JAM11) have 31−41 vol.% plagioclase, 54−60 vol.% amphibole and 5−9 vol.% Fe-Ti oxides (Tab. 1). Hornblende form coarse (2−3 mm wide) subhedral grains with Fe-Ti oxide inclusions and plagioclase is present as a mosaic of submillimetric euhedral altered crystals (Fig. 3b). Pargasite has lower Mg# (50−55), Si (5.90−6.25 apfu) and Ti (0.08−0.15 apfu) compared to hornblendites and plagioclase is calcic (An85-92).

Other samples were classified as diorite, quartz diorite and tonalite according to their normalized volumetric amount of quartz in the quartz-plagioclase subsystem (<5%, between 5 and 20% and >20%, respectively) but they generally share similar textural features (Figs. 3c and 3d). They show variable amounts of amphibole (from 2 to 20 vol%), plagioclase (48−60 vol.%), biotite (10−25 vol.%) quartz (2−21 vol.%) with minor epidote (0−4 vol.%), apatite (<1.5 vol.%) and Fe-Ti oxides (<4 vol.%) (Tab. 1). In samples with perfectly preserved igneous textures (BEA, FIG, DBD, SAJ12A, SAJ13, STC; Fig. 3c) plagioclase is millimetric, tabular euhedral and zoned. It also occurs as small euhedral inclusions into amphibole and biotite. Amphibole is coarse (1−3 mm) but systematically interstitial into the network of plagioclase tabs and contains numerous inclusions of Ti-magnetite. Biotite forms large anhedral flakes with amphibole and rare epidote inclusions. The latter is commonly associated with biotite, apatite, quartz and zircon in interstitial clusters (particularly in tonalites). Localized deformation and recrystallization is observed in some samples (BIA, FRO, MAR1, MAR2, SHB, SJL; SPL Fig. 3d). The texture locally becomes porphyroclastic with elongated grains of amphibole and biotite surrounding partly recrystallized porphyroclasts of plagioclase. Plagioclase recrystallized along its rims, locally exhibiting a mortar texture and some large grains show undulatory extinction with development of subgrains, but igneous concentric zoning is still preserved. Undeformed zones are preserved with similar textures to those described previously. DCAP sampled in the active quarry of the Capdenac massif is a porphyroclastic gneiss with almost no relic of the original igneous texture except for large zoned plagioclase grains. Amphibole in quartz diorite is on average slightly more magnesian (Mg#: 58−47) and Si rich (6.35−6.70 apfu) compared to amphibole in tonalites (Mg#: 41−53; 6.20−6.75 Si pfu). Similar variations are observed for biotite (Mg#: 54−63 in quartz diorite and 49−54 in tonalites). There is no systematic differences in composition between amphibole from deformed versus undeformed samples. However, a small fringe of actinote can develop around hornblende in some deformed samples. Plagioclase is on average more calcic in quartz diorite (An39-48) than in tonalites (An32-45) (Fig. 4).

The amphibole-plagioclase barometer of Molina et al. (2015) and the amphibole-only thermometer of Liao et al. (2021) were used to estimate pressure and temperature of crystallization. They can be applied on a wide range of bulk and mineral compositions (although not for hornblendites in which plagioclase is totally altered) and have the following calibration uncertainties: ±35 °C for Liao et al. (2021) and ±1.5 to ±2.3 kbar for Molina et al. (2015). About 100 amphibole-plagioclase pairs were used for thermobarometry and results range from 4.5 to 8.5 kbar and 640−840 °C. When looking at averages per sample (Fig. 4), most calculated pressures fall within 5−8 kbar corresponding to middle or lower crustal depths (19−30 km) and 700−800 °C, i.e. temperature above the water-saturated tonalitic solidus (Schmidt and Thompson, 1996). Considering the large uncertainties on pressure estimations all the calculated pressures for diorites and tonalites are equivalent within errors. There are therefore no systematic differences between samples affected by solid-state deformation/recrystallization (MAR1, SJL) and those with igneous textures (FIG, SAJ12A, DBD, STC) or when comparing those found in the so called eo-variscan Leyme units (FIG, STC) to plutons from the Upper Gneiss Unit (MAR1, SJL, SAJ12A, DBD).

thumbnail Fig.3

Representative microphotographs of hornblendite TUL3 (a), hornblende gabbro TUL1 (b), hornblende-biotite quartz diorite SAJ12A (c) and hornblende-biotite tonalite MAR1 (d). Despite post-magmatic alteration, igneous textures are generally preserved (a, b, c) but a weak foliation can develop locally (d). Note the interstitial epidote-apatite-zircon-biotite cluster in sample MAR1 (d).
thumbnail Fig.4

Mineral composition and thermobarometry. a) Composition and classification of amphiboles (a); composition of plagioclases (b) and average P-T results per sample calculated with amphibole only and amphibole-plagioclase thermobarometers (see text) (c). The error bars represent the 1σ standard deviation of individual pressures and temperatures calculated for each sample.

5 Geochemistry

Rocks belonging to the Tonalite belt of Limousin form a typical low-Ti, normal to high-K calc-alkaline suite (Tab. S2, Fig. 5). Hornblendites and hornblende gabbros are cumulates characterized by low SiO2 (40−43 wt%), high Al2O3 (16−23 wt%) and high MgO (8−12.5 wt%) (Fig. 6). They display the typical REE profile of amphibole segregates with a slight depletion in LREE compared to HREE ((La/Yb)N : 0.59−1.10) accompanied by a weak negative Eu anomaly (Eu/Eu*: 0.68−0.95, with Eu* = (Sm*Gd)0.5 (Fig. 7). They are depleted in very incompatible elements (especially in Th and U relative to Nb) compared to quartz diorites and tonalites, but share similar concentrations in Y and HREE (Fig. 7). Diorite (50 wt% SiO2; Mg#: 49), quartz diorites (52−54 wt% SiO2; Mg#: 49−54) and tonalites (54.5−60.7 wt% SiO2; Mg#: 44−53) form near linear trends in Harker diagrams (Fig. 6) marked by a decrease in TiO2, Al2O3, CaO and a slight increase in K2O with increasing differentiation. In most cases, when these sublinear trends are extrapolated toward low SiO2 values, they reach a composition of cumulate hornblende gabbros (Fig. 6). In REE and multi-element diagrams (Fig. 7), they are enriched in very incompatible elements (LREE, Rb, Ba, Th, U) with negative anomalies in Nb-Ta (Nb/Nb*: 0.12-0.67 with Nb*: (Th*La)0.5) and Ti, features recalling typical basic to intermediate arc lavas (Fig. 7). On average, the LREE enrichment and the magnitude of the Nb and Ti negative anomalies increases from diorite to tonalites (Fig. 7). The amplitude of the negative Eu anomaly also gently increases from diorite (Eu/Eu*: 0.94) to quartz-diorites (Eu/Eu*: 0.83−1.02) and tonalites (Eu/Eu*: 0.71−0.86).

The Sr-Nd isotopic composition of 10 representative samples (Fig. 8, Tab. S3) lies between 0.7050−0.7061 for initial 87Sr/86Sr and +0.2 to +2.0 for εNd. These values are similar to those reported by Shaw et al. (1993) for the Limousin Tonalite Belt but have slightly lower εNd compared to diorite and tonalites from the Northern and Eastern Massif Central (Pin and Paquette, 2002) (Fig. 8). Two UGU paragneisses have also been analyzed to test crustal assimilation, they plot between εNd −1 to −2 and 87Sr/86Sr 0.7095−0.7115 (Fig. 8).

thumbnail Fig.5

Classification of Late Devonian − Early Carboniferous calc-alkaline plutonic rocks from the French Massif Central and the Armorican Massif in the K2O-SiO2 diagram (a) and the AFM triangle plot (b). Literature data (white, black and grey symbols) are from Pfeiffer (1986); Shaw et al., (1993); Cuney et al., (2001) and Pin and Paquette (2002).
thumbnail Fig.6

Harker diagrams showing major element composition of calc-alkaline plutons in the French Massif Central and the Armorican massif. Literature data (white, black and grey symbols) are from 1986 (Pfeiffer); Shaw et al. (1993); Cuney et al., (2001) and Pin and Paquette (2002). Note the near linear trends drawn by diorites and tonalites extrapolating toward the composition of hornblende gabbros.
thumbnail Fig.7

Chondrite-normalized REE patterns and primitive mantle normalized multi-element diagrams for the calc-alkaline plutonic rocks investigated in this study. The compositions of volcanic rocks from the South Aegean arc are extracted from the Georoc database. Normalization values are from McDonough and Sun (1995).
thumbnail Fig.8

Sr-Nd isotopic diagram for calc-alkaline lavas/plutons from the French Massif Central (this study, Shaw et al., 1993; Cuney et al., 2001; Pin and Paquette, 2002) and for Brévenne back-arc magmatism (Pin and Paquette, 1997). Also shown, the composition of the Depleted Mantle (DM) at 360 Ma (calculated with the data of Salters and Stracke, 2004), the probable signature of the mantle wedge around 360 Ma as attested by calc-alkaline magmatism arc like magmatism. The green dotted mixing hyperbola represent mixing between Depleted Mantle and a UGU paragneiss to model mixing of the mantle source with sediments buried by subduction while the black mixing dotted curve simulates contamination of a depleted basalt by a UGU paragneiss to reproduce continental assimilation. Values in % indicate the relative mass of sedimentary component in the mix.

6 Zircon U-Pb geochronology

Three samples were selected for U-Pb dating on the basis of their geographical position and the amount of deformation/recrystallization (Fig. 9, Tab. S4). SJL is a tonalite from the Saint Jean de Ligoure massif located in northwest Limousin. This sample shows evidences for localized deformation and recrystallization but no macroscopic planar foliation. MAR1 is from the central Tulle-St Féréole massif. It is more affected by deformation than SJL as it shows a faint macroscopic gneissic fabric. FIG comes from the southernmost Capdenac massif and has perfectly preserved igneous textures at macroscopic and microscopic scales.

SJL yielded euhedral zircon crystals with elongation ratios around 1:3 and very regular oscillatory zoning (Fig. 9a). However, the grains were fractured with numerous minute inclusions of apatite (green on CL images), quartz (generally orange red) and epidote (black). Very few grains had a small (<10 μm wide) resorbed bright core that could not be analyzed. On the 25 spots made on the zircons, 12 were disturbed by the presence of inclusions and fractures generating very irregular ablation spectra (with brutal drop or jumps of the signal for Pb, Th and U isotopes). The remaining, undisturbed, 13 analyses give a concordant date of 359 ± 2 (± 5) Ma (MSWD: 0.97; Fig. 9b) and high Th and U contents (180−560 ppm Th and 390−1200 ppm U) translating into Th/U ratios between 0.43 and 0.62. Zircons from MAR1 are all euhedral, coarse (up to 400 μm long) and either elongated (1:5) or stocky (1:2). A lot of grains display a less luminescent core with sector (and sometimes faint oscillatory) zoning but the bright mantle zones have concentric oscillatory zoning (Fig. 9a). 20 spots yielded regular ablation spectra and very similar ellipses with a concordia date of 359 ± 3 (± 5) Ma (MSWD: 0.35) (Fig. 9c). Th/U ratios are high (0.43−0.78) but Th and U content are rather low (30−100 and 53−150 ppm respectively) compared to SJL zircons. Note that spots within darker cores yielded the same age that those made in mantle zones with oscillatory zoning (Fig. 9a). FIG had euhedral, often elongated (1:4 to 1:5) zircons but smaller than in the two other samples (<200 μm in length). On CL images (Fig. 9a), they generally display darker core with subtle oscillatory zoning and bright mantle/rims but no resorption textures between the zones. On the 25 spots made on the mount, 18 yielded regular ablation spectra with a concordia date of 360 ± 2 (± 5) Ma (MSWD: 1.1) (Fig. 9d). Th and U concentration were around 110−660 and 300−720 ppm respectively yielding Th/U ratios between 0.33 and 0.96. Similar ages were recovered on dark cores and on bright rims (Fig. 9a).

thumbnail Fig.9

Results of in situ LA-ICP-MS U-Pb zircon dating. (a) Cathodoluminescence images of representative zircon grains. Scale bars represent a length pf 20 μm and the quoted ages in Ma are calculated from 206Pb/238U ratios. (b, c, d) Concordia plots for tonalites SJL, MAR and FIG, respectively. MSWD and p-value are for concordance and equivalence. The two uncertainties are reported as 2σ.

7 Discussion

7.1 Differentiation of a mantle-derived hydrous calc-alkaline suite at mid to lower crustal levels

The diorites and tonalites from the Western French Massif Central form a medium- to high-K calc-alkaline suite (Fig. 5). Occasional granodiorites are also mentioned (Pfeiffer, 1986) but none were analyzed in this study. Diorites, quartz-diorites and tonalites form nearly linear trends in major element variation diagrams (Fig. 6). The composition of these plutonic rocks can be interpreted as melts and not crystal segregates, as confirmed by the absence of positive Eu and Sr anomalies despite high modal plagioclase content. The liquid line of descent formed by diorites and tonalites in major element variation diagrams is controlled by the fractionation of hornblendites and hornblende gabbros assemblages (Fig. 6). The most primitive diorites and quartz-diorites have 50−54 SiO2 wt%, Mg# around 49−54 and high alumina content, up to 20 wt% Al2O3 corresponding in composition to high alumina basalts and basaltic andesites found in volcanic arcs (Sisson and Grove, 1993). However, the quite low Mg#, compared to melt in equilibrium with mantle peridotite suggest that the most primitive samples were already slightly fractionated, as also confirmed by their quite low Ni and Cr contents (below 50 and 180 ppm respectively). The ubiquity of hornblende and titano-magnetite, even in the most primitive samples, with highly calcic plagioclase (up to 92 mol% An) in the hornblende gabbros is consistent with strongly hydrous oxidized parental magmas, typical of calc-alkaline suites found in continental magmatic arcs (Grove and Baker, 1984; Sisson and Grove, 1993). The trace element pattern of all samples except cumulates, also mimics volcanic rocks from calc-alkaline continental arcs (Fig. 7) with their typical Nb-Ta-Ti negative anomalies and their enrichment in fluid mobile elements (Rb, Ba, Sr, Pb). Moreover, near- but supra-solidus igneous epidote formed in interstitial clusters with biotite, zircon and apatite, also characterize hydrous, oxidized calc-alkaline magmas (Schmidt and Thompson, 1996; Schmidt and Poli, 2004). Considering that crystallization pressures determined with amphibole-plagioclase barometry are between 5 and 8 kbar for the French Massif Central diorites and tonalites (Fig. 4), supra-solidus late crystallization epidote probably indicate oxygen fugacities close to the NNO buffer (Schmidt and Thompson, 1996) during crystallization of tonalites.

Review of published (Shaw et al., 1993) and new isotopic data (Fig. 8) highlight the remarkable homogeneity of initial Sr and Nd isotopic compositions for the Limousin massifs (87Sr/86Sr: 0.7040−0.7061 and εNd: −1 to +2). There is no systematic variations of Sr-Nd isotopic composition with the lithological nature and SiO2 content (Fig. 8). Tonalites indeed have similar compositions to diorites and to primitive hornblende gabbro cumulates. Testing crustal assimilation of UGU paragneiss hosting the diorites-tonalites plutons by depleted Brévenne back arc basalt, representing a contemporaneous depleted mantle component, to explain the isotopic signatures of the calc-alkaline plutonic suite yield unrealistic levels of contamination (20−60%; Fig. 8) to reproduce the composition of Limousin diorites and tonalites. Instead, calculating contamination of a 360 Myr-old depleted MORB mantle (using values of Salters and Stracke, 2004 recalculated at 360 Ma) by sediments buried by subduction (UGU paragneiss were taken as the 360 Ma sedimentary components) yields results around 2% (Fig. 8). Such level of crustal input into mantle wedge is in agreement with estimates made for recent primitive arc magmas (Chauvel et al., 2009). Therefore crustal assimilation did not play a major role during differentiation of the Limousin diorite-tonalite suite. The isotopic compositions of the parental magmas are explained by melting of a mantle wedge contaminated by sediment-derived components buried by the slab with only a minor contribution of crustal assimilation as already suggested by various authors (Shaw et al., 1993; Cuney et al., 2001a; Pin and Paquette, 2002). Calc-alkaline lavas and plutons located in Northern and Eastern French Massif Central (FMC) have slightly different values (εNd: +1.5 to +5) probably representing a less contaminated mantle source or parental magma (Pin and Paquette, 2002). Globally, as proposed by Shaw et al. (1993), the focal zone represented by Late Devonian calc-alkaline plutons and lavas in the Sr-Nd diagram (Fig. 8) is the signature of the variously enriched Variscan mantle wedge under the French Massif Central around 360 Ma. Late Devonian magmatism with depleted isotopic signatures (εNdi up to +8 and 87Sr/86Sri down to 0.7040) characterizes some basic and felsic volcanic rocks of the Brévenne area, interpreted as a former continental back-arc basin (Pin and Paquette, 1997; Fig. 8). The dichotomy of enriched signatures for tonalite and diorites similar to arc magmas versus depleted signatures at the Brévenne back arc is consistent with geochemical studies on modern arc-back arc systems (Pearce et al., 2005; Pearce and Stern, 2006) describing a depleted mantle source at the back arc and an enriched mantle source below the arc.

7.2 Short-lived Late Devonian-Early Carboniferous subduction in the French Massif Central?

New U-Pb zircon dating of tonalites from Limousin yield dates ranging between 365 and 354 Ma (uncertainties included, Fig. 9). Considering that zircons are euhedral, show typical oscillatory growth zoning of igneous zircons (Corfu, 2003; Fig. 9a) and have high Th/U ratios, we interpret these dates as igneous crystallization ages. The results obtained on three different plutons are equivalent within uncertainties, although they are variously deformed and recrystallized. No correlation between geographical position and age of emplacement can be assumed based on these data.

Published zircon ages (ID-TIMS and SHRIMP; Bernard-Griffiths et al., 1985; Bertrand et al., 2001; Pin and Paquette, 2002) on calc-alkaline plutons across the French Massif Central all converge to a magmatic event ranging from Late Devonian to Lowermost Carboniferous (379−355 Ma). Recalculation of all these published ages using IsoplotR (Vermeesch, 2018) to allow better comparison yielded very similar ages compared to published ones (Fig. 10), except for the strongly discordant data of Bernard-Griffiths et al. (1985). The latter made an assumption on a recent Pb Loss and obtained two anchored discordia built only with two data points each that yielded upper intercepts of 379 ± 19 Ma and 355 ± 2 Ma, interpreted as igneous crystallization ages. Different upper intercept ages were obtained after recalculation with IsoplotR without assuming recent lead loss: 358 ± 28 Ma and 353 ± 27 Ma. Considering the very large uncertainties and the fact that each discordia is based on two analytical points only, the ages published by Bernard-Griffiths et al. (1985) have been ignored in the following discussion. Bertrand et al. (2001) dated a diorite from the northwestern FMC (Fig. 1) by ID-TIMS and SHRIMP. Recalculated TIMS data are discordant but 5 zircon fractions out of 6 are perfectly aligned and yield a precise and statistically robust age of 359 ± 3 Ma (MSWD: 1.1). SHRIMP data are more dispersed (MSWD: 6) but equivalent within uncertainties (356 ± 7 Ma); TIMS results should thus be preferred. Combining these reviewed data, with our new results (365−354 Ma, uncertainties included) and the precise data from Pin and Paquette (2002) for the North and Eastern FMC diorites (368−359 Ma, uncertainties included), it can be suggested that active margin magmatism in the French Massif central is rather short lived (368−354 Ma). Interpretation of the apparent spread of ages into two successive events that are igneous crystallization around 380 Ma and metamorphic recrystallization around 360 Ma as done by Faure et al. (1997) and Bellot and Roig (2007) is therefore irrelevant. The new and the published precise U-Pb zircon ages should be interpreted as representing one single igneous event.

Similar expressions of Devono-Carboniferous calc-alkaline magmatism are also preserved as Fammenian lavas in Morvan (372−359; Pin and Paquette, 2002), in the Armorican massif (373 +6/−11 Ma; Cuney et al., 1993) and the Variscan basement massifs of Western Alps (371−356 Ma; Bertrand et al., 2000; Guillot et al., 2012; Bergomi et al., 2017). The ID-TIMS data of Cuney et al. (1993) could not be recalculated and checked because measured isotopic ratios were not published in a table. The discordia is traced with 3 points not perfectly aligned but with one falling on the concordia chord. Extension, sedimentation and magmatism related to a back-arc or arc setting locally forming incomplete and dismembered “ophiolites” (Figs. 1 and 10) or basaltic lava flows interbedded with Devonian sediments are also preserved in the Armorican massif (but the igneous products are not dated; Thieblemont and Cabanis, 1994; Ducassou et al., 2011), the Eastern French Massif Central (371-357 Ma, uncertainties included; Sider et Ohnenstetter, 1986; Ohnenstetter et Sider, 1988 ; Feybesse et al., 1988 ; Pin et Paquette, 1997) and in the Vosges (372 ± 18 Ma; Skrzypek et al., 2012).

Subduction activity in the French Massif central is moreover supported by markers of HP metamorphism (Faure et al., 2009; Lardeaux et al., 2014). The polycyclic view of the Variscan belt in FMC (Faure et al., 1997) suggests that an early subduction event took place during a Late Silurian eo-variscan event. This is supported by dating of eclogites and related rocks yielding ages ranging from 430 to 390 Ma (Paquette et al., 1995; Gebauer et al., 1981; Ducrot et al., 1983; Berger et al., 2010). It is interesting to note that re-investigations of ID-TIMS zircon eo-variscan ages obtained by Paquette et al. (1995) on southern FMC eclogites with the in-situ LA-ICP-MS method (Paquette et al., 2017) did not reproduced the eo-variscan Early Devonian ages. Moreover, recent geochronological investigations on FMC eclogites (including new investigations on previously dated occurrences) involving in situ U-Pb dating, Lu-Hf and Sm-Nd garnet dating (Lotout et al., 2018, 2020; Benmammar, 2021; de Hoÿm de Marien et al., 2023) all yielded Middle Devonian to Lowermost Carboniferous ages (Fig. 10). Eclogites from the Upper Allochthonous units (Upper Gneiss Unit, Leptyno-amphibolite complex and Najac HP-LT unit) have Givetian to Frasnian U-Pb zircon ages (387 to ∼360 Ma; Lotout et al., 2018, 2022; Benmammar, 2021; De Hoÿm de Marien et al., 2023). Finally, it must be emphasized that the sedimentary record in the South Armorican zone (including the presence of concordant 400 Ma inherited magmatic zircons) suggests that relief variations controlled by an active subduction started as early as Lower Devonian (Ducassou et al., 2014).

Magmatic and metamorphic markers for subduction in the French Massif Central all range between 387 and 356 Ma (Fig. 10), suggesting at least ca 31 Myr of active subduction. This is short compared to orogenic systems characterized by long-lived subduction such as the Himalaya where at least 100 Myr of arc magmatism is recorded before collision (Burg and Bouilhol, 2019). This range is however comparable to the Western Alps, which is characterized by dispersed, short lived near-collision arc-like magmatism and 40 Myr long record of HP to UHP metamorphism (McCarthy et al., 2020). The calc-alkaline lavas and plutons in French Massif Central do not underline a transcrustal magmatic arc formed during a protracted subduction event; they are dispersed and isolated, suggesting that magmatic inputs were low and sparse. According to the metamorphic and igneous records, it can be suggested that subduction in the French Massif Central was short lived and thus probably consumed a small oceanic basin and/or hyper-extended continental margin comparable to the Alpine example. It is worth to note that the presence of post-nappes Visean calc-alkaline magmatism is sometimes taken as an evidence for renewed subduction after the main Tournaisian collision event (Pin and Paquette, 2002). The same authors also suggest that this magmatism can form during the post-collisional period without active subduction, by remelting of the mantle previously enriched by the Devonian subduction.

thumbnail Fig.10

(a) Compilation of geochronological data obtained on petrological markers of Middle to Late Devonian subduction in the French Massif Central. U-Pb zircon on calc-alkaline magmatism are from Bertrand et al. (2001), Pin and Paquette (2002) and this study; age of Brévenne magmatism is from Pin and Paquette (1997), recent ages U-Pb zircon and Sm-Nd/Lu-Hf garnet ages obtained on eclogites are from Lotout et al. (2018, 2020) and de Hoÿm de Marien et al. (2023), UGU anatexis as dated by U-Th-Pb monazite chemical age is from Melleton et al. (2009) and the curve indicating the relative igneous activity is from Chopin et al. (2023) traced with data compiled by Vanderhaeghe et al. (2020). (b) Tentative snapshot of the geodynamic situation around 360 Ma for the Western French Massif Central, with emphasis on diorites-tonalites, eclogites and migmatitic rocks belonging to the UGU (see text for explanations).

7.3 Some implications for the Variscan evolution of the Western French Massif Central

Three key implications can be proposed for the Variscan evolution of Western French Massif Central on the basis of the study of dioritic and tonalitic plutons.

(i) These igneous bodies, emplaced in supra-subduction position, intruded UGU migmatites and retrogressed UGU eclogite at 19−30 km depth between 365 and 354 Ma in Limousin (Fig. 2). Therefore, UGU Limousin eclogites were subducted to about 60 km depth (Bellot and Roig, 2007), exhumed to crustal depths and accreted to the upper plate of the subduction system before 365 Ma while subduction was still ongoing (Fig. 10). A process of relamination of subducted crust has recently been proposed to explain the syn-subduction transfer from the lower plate to the upper plate of UGU eclogites (De Hoÿm de Marien et al., 2023) but several other processes can be considered such as slab roll back (Faccenna and Brun, 2008) also supporting formation of the synchronous Brévenne back arc.

(ii) The migmatitic paragneisses hosting the Limousin eclogites were dated at 362 ± 4 Ma chemical U-Th-Pb monazite age; Melleton et al., 2009). Anatexis occurred around 6−10 kbar and 700−800 °C (Bellot and Roig, 2007), i.e. grossly at the same conditions than diorites and tonalites crystallization (Fig. 4). Partial melting in the UGU paragneisses is thus concomitant with calc-alkaline plutonism, meaning that the metasediments were partially molten during the intrusion of tonalites and diorites (Fig. 10).

(iii) The calc-alkaline plutons are found in the UGU with two occurrences in the amphibolite facies Leyme unit Duguet et al., 2007; see also Pfeiffer, 1986). The geological map of the St Féréole-Tulle massif (Fig. 2) and observations made by Pfeiffer (1986) clearly show that the shear zone separating LGU and UGU cuts across and deforms the igneous body while the latter is never found within LGU. These basic observations exclude the hypothesis of an early tectonic superposition of LGU and UGU during an Early Devonian eo-variscan event as proposed by Ledru et al. (1989), Faure et al. (2009) and Lardeaux (2014). Specifically, a Middle Devonian thrusting of eclogites-bearing UGU and LGU onto the low grade Parautochthonous unit in Limousin as proposed by Bellot and Roig (2007) is unlikely as Late Devonian calc-alkaline plutons are found in the upper allochthonous units only. Tectonic transport of UGU onto LGU and the Parautochtonous Unit consequently occurred after 365−354 Ma in the Western French Massif Central. According to the geochronological record of Variscan igneous activity in French Massif Central (Vanderhaeghe et al., 2020; Chopin et al., 2023; Fig. 10), the Late Devonian-Lowermost Carboniferous active margin magmatism was followed and partly coeval with crustal melting and emplacement of the single but huge granitic Guéret massif (359−351 Ma; Cartannaz et al., 2007; Marcoux et al., 2021) that is included within UGU (Gébelin et al., 2006). The main phase of syn- to late-Variscan magmatism ignited in the Visean with both calc-alkaline and anatectic granites cutting across the mylonitic contacts between tectono-metamorphic units (Vanderhaeghe et al., 2020). This igneous transition could probably mark the climax of Variscan collision, encompassing the stacking of units near the Tournaisian-Visean boundary.

8 Conclusions

Diorites and tonalites in Western Massif Central formed after partial melting of a hydrous mantle wedge and fractional crystallization in the crust controlled by segregation of plagioclase and amphibole. New U-Pb zircon ages support emplacement in the continental upper plate of an active subduction context around 365−354 Ma (uncertainties included). The plutons crystallized at middle to lower crustal depths (19−30 km) and intruded migmatites and exhumed retrogressed eclogites of the Upper Gneiss Unit but also amphibolite facies paragneiss of the so-called Leyme unit. Field data, published observations and geological maps indicate that calc-alkaline plutons emplaced before the final tectonic piling of the UGU onto LGU but after the exhumation of UGU eclogites.

The magmatic pulse represented by calc-alkaline plutons across the different Variscan massifs of France (371−356 Ma) is partly synchronous with Late Devonian-Early Carboniferous U-Pb zircon ages recently obtained for HP metamorphism in the French Massif Central (387−357 Ma). The available magmatic and metamorphic record suggest that the subduction was short lived (∼30 Myr), suggesting either slow convergence rates and/or subduction of a small oceanic basin. Subduction and exhumation of UGU eclogites started before 365 Ma and active margin magmatism as represented by the diorites and tonalites argue that subduction was still ongoing between 365−354 Ma in the Western French Massif Central. All these results are against the existence of an eo-variscan event leading to an early superposition of LGU and UGU during Middle Devonian.

Further geochronological studies should focus on LGU eclogites that are not yet precisely dated and on the record of detrital zircons in clastic and volcaniclastic Devonian to Lowermost Carboniferous sediments. The latter may have preserved igneous and metamorphic Devonian detrital zircons witnessing former subduction-related units that are not currently cropping out.

Supplementary material

Table S1: Composition of amphibole and plagioclase plotted in Figure 4 and used for thermobarometry.

Table S2: Bulk rock major (wt%) and trace (?g/g) element contents of samples investigated in this study.

Table S3: Sr-Nd isotopic composition of representative samples, including host paragneisses.

Table S4: Results of LA-ICP-MS U-Pb dating.

Access here

Acknowledgments

We would like to thank Fabienne Barrère de Parseval for making the thin sections and Philippe de Parseval for his help with microprobe analysis. Guy and Maryse Chantepie are warmly acknowledged for showing us two gorgeous hornblendite outcrops. This work received financial support from the Université Paul Sabatier de Toulouse (IDEX nouveaux entrants) from a CNRS INSU TELLUS SYSTER grant. Nathan Cogné and one anonymous reviewer provided insightful comments to improve the manuscript. Valérie Bosse and Laurent Jolivet are thanked for their efficient editorial handling.

We would also like to acknowledge late Jean-Louis Paquette for his enormous and remarkable work, which has led to major advances in our understanding of Variscan evolution.

References

  • Ballèvre M, Martínez Catalán JR, López-Carmona A, Pitra P, Abati J, Fernández RD, et al. 2014. Correlation of the nappe stack in the Ibero-Armorican arc across the Bay of Biscay: a joint French-Spanish project. Geol Soc Lond Spec Publ 405: 77–113. [CrossRef] [Google Scholar]
  • Bellot J-P., Roig J-Y. 2007. Episodic exhumation of HP rocks inferred from structural data and PT paths from the southwestern Massif Central (Variscan belt, France). J Struct Geol 29: 1538–1557. [CrossRef] [Google Scholar]
  • Benmammar A, Berger J, Triantafyllou A, Duchene S, Bendaoud A, Baele J-M, et al. 2020. Pressure-temperature conditions and significance of Upper Devonian eclogite and amphibolite facies metamorphisms in southern French Massif central. BSGF-Earth Sci Bull 191. 1–18. [Google Scholar]
  • Benmammar A. 2021. Comment expliquer la dualité thermique des subductions pré-orogéniques ? : Exemple du métamorphisme de haute-pression dans le Massif central français. Unpublished PhD Thesis, Université Paul Sabatier de Toulouse. 212pp (available online). [Google Scholar]
  • Berger J, Féménias O, Ohnenstetter D, Bruguier O, Plissart G, Mercier J-CC , et al. 2010. New occurrence of UHP eclogites in Limousin (French Massif Central): age, tectonic setting and fluid-rock interactions. Lithos 118: 365–382. [CrossRef] [Google Scholar]
  • Bergomi MA, Dal Piaz GV, Malusà MG, Monopoli B, Tunesi A. 2017. The Grand St Bernard-Briançonnais Nappe System and the Paleozoic Inheritance of the Western Alps Unraveled by Zircon U-Pb Dating. Tectonics 36: 2950–2972. [CrossRef] [Google Scholar]
  • Bernard-Griffiths J, Gebauer D, Grunenfelder M, Piboule M. 1985. The tonalite belt of Limousin (French Central Massif); U-Pb zircon ages and geotectonic implications. Bull Soc Géol Fr 1: 523–529. [CrossRef] [Google Scholar]
  • Bertrand JM, Pidgeon R, Leterrier J, Guillot F, Gasquet D, Gattiglio M. 2000. SHRIMP and IDTIMS U-Pb zircon ages of the pre-Alpine basement in the Internal Western Alps (Savoy and Piemont). Schweiz Mineral Petrogr Mitt 80: 225–248. [Google Scholar]
  • Bertrand JM, Leterrier J, Cuney M, Brouand M, Stussi JM, Delaperrière E, et al. 2001. Géochronologie U-Pb sur zircons de granitoïdes du Confolentais, du massif de Charroux-Civray (seuil du Poitou) et de Vendée. Géol Fr 1: 167–189. [Google Scholar]
  • Bosse V, Feraud G, Ruffet G, Ballèvre M, Peucat J-J., De Jong K. 2000. Late Devonian subduction and early‐orogenic exhumation of eclogite‐facies rocks from the Champtoceaux Complex (Variscan belt, France). Geol J 35: 297–325. [CrossRef] [Google Scholar]
  • Bosse V, Féraud G, Ballevre M, Peucat J-J., Corsini M. 2005. Rb-Sr and 40Ar/39Ar ages in blueschists from the Ile de Groix (Armorican Massif, France): implications for closure mechanisms in isotopic systems. Chem Geol 220: 21–45. [CrossRef] [Google Scholar]
  • Bruguier O, Bosch D, Caby R, Vitale-Brovarone A, Fernandez L, Hammor D, et al. 2017. Age of UHP metamorphism in the Western Mediterranean: insight from rutile and minute zircon inclusions in a diamond-bearing garnet megacryst (Edough Massif, NE Algeria). Earth Planet Sci Let 474: 215–225. [CrossRef] [Google Scholar]
  • Burg J-P., Bouilhol P. 2018. Timeline of the South Tibet − Himalayan belt: the geochronological record of subduction, collision, and underthrusting from zircon and monazite U-Pb ages. Can J Earth Sci 56: 1318–1332. [Google Scholar]
  • Cartannaz C, Rolin P, Cocherie A, Marquer D, Legendre O, Fanning CM, et al. 2007. Characterization of wrench tectonics from dating of syn- to post-magmatism in the north-western French Massif Central. Int J Earth Sci 96: 271–287. [CrossRef] [Google Scholar]
  • Cartier C, Faure M. 2004. The Saint-Georges-sur-Loire olistostrome, a key zone to understand the Gondwana-Armorica boundary in the Variscan belt (Southern Brittany, France). Int J Earth Sci 93: 945–958. [CrossRef] [Google Scholar]
  • Chauvel C, Marini J-C., Plank T, Ludden JN. 2009. Hf-Nd input flux in the Izu-Mariana subduction zone and recycling of subducted material in the mantle. Geochem Geophys Geosys 10: Q01001. [CrossRef] [Google Scholar]
  • Chopin F, Leprêtre R, El Houicha M, Tabaud A-S., Schulmann K, Míková J, et al. 2023. U-Pb geochronology of Variscan granitoids from the Moroccan Meseta (Northwest Africa): Tectonic implications. Gondwana Res 117: 274–294. [CrossRef] [Google Scholar]
  • Corfu F, Hanchar JM, Hoskin PWO, Kinny P. 2003. Atlas of Zircon Textures. Rev Mineral Geochem 53: 469–500. [CrossRef] [Google Scholar]
  • Cuney M, Stussi J-M., Brouand M, Dautel D, Michard A, Gros Y, et al. 1993. Géochimie et géochronologie U/Pb des diorites quartziques du Tallud et de Moncoutant: nouveaux arguments pour une extension de la Ligne Tonalitique Limousine en Vendée. C R Acad Sci Série 2316: 1383–1390. [Google Scholar]
  • Cuney M, Brouand M, Stussi JM, Virlogeux D. 2001 a Le complexe plutonique de Charroux-Civray (Vienne): témoin du magmatisme infra-carbonifère dans le segment occidental de la chaîne varisque européenne. Géol Fr 1-2: 143–166. [Google Scholar]
  • Cuney M, Brouand M, Stussi J-M. 2001b. Le magmatisme hercynien en Vendée. Corrélations avec le socle du Poitou et l’ouest du Massif central français. Géol Fr 1-2: 117–142. [Google Scholar]
  • de Hoÿm de Marien L, Pitra P, Poujol M, Cogné N, Cagnard F, Le Bayon B. 2023. Complex geochronological record of an emblematic Variscan eclogite (Haut-Allier, French Massif Central). J Metam Geol xx 1–29. [Google Scholar]
  • Didier J. Lameyre J. 1971. Les roches granitiques du Massif Central. In: Symposium J Jung Plein Air Service, Clermont-Ferrand, 133–155. [Google Scholar]
  • Ducassou C, Ballèvre M, Lardeux H, Robin C. 2011. Evidence for pre-orogenic, Early Devonian rifting in the Variscan belt: stratigraphy and structure of the Palaeozoic cover of the Mauges Unit (Upper Allochthon, Armorican massif, France). Int J Earth Sci 100: 1451–1475. [CrossRef] [Google Scholar]
  • Ducrot J, Lancelot JR, Marchand J. 1983. Datation U-Pb sur zircons de l’éclogite de la Borie (Haut-Allier, France) et conséquences sur l’évolution anté-hercynienne de l’Europe occidentale. Earth Planet Sci Lett 62: 385–394. [CrossRef] [Google Scholar]
  • Duguet M, Le Breton N, Faure M. 2007. P-T paths reconstruction of a collisional event: The example of the Thiviers-Payzac Unit in the Variscan French Massif Central. Lithos 98: 210–232. [CrossRef] [Google Scholar]
  • Faccenna C, Brun J-P. 2008. Exhumation of high-pressure rocks driven by slab rollback. Earth Planet Sci Lett 272: 1–7. [CrossRef] [Google Scholar]
  • Faure M, Leloix C, Roig J-Y. 1997. L’évolution polycyclique de la chaîne hercynienne. Bull Soc Géol Fr 168: 695–705. [Google Scholar]
  • Faure M, Mézème EB, Duguet M, Cartier C, Talbot J-Y. 2005. Paleozoic tectonic evolution of medio-Europa from the example of the French Massif Central and Massif Armoricain. J Virtual Explor 19: 1–25. [CrossRef] [Google Scholar]
  • Faure M, Bé Mézème E, Cocherie A, Rossi P, Chemenda A, Boutelier D. 2008. Devonian geodynamic evolution of the Variscan Belt, insights from the French Massif Central and Massif Armoricain. Tectonics 27: TC2005. [CrossRef] [Google Scholar]
  • Faure M, Lardeaux J-M., Ledru P. 2009. A review of the pre-Permian geology of the Variscan French Massif Central. C R Geosci 341: 202–213. [CrossRef] [Google Scholar]
  • Faure M, Mezeme EB, Cocherie A, Melleton J, Rossi P. 2009b. The South Millevaches Middle Carboniferous crustal melting and its place in the French Variscan belt. Bull Soc Géol Fr 180: 473–481. [CrossRef] [Google Scholar]
  • Faure M, Li X-H., Lin W. 2017. The northwest-directed “Bretonian phase” in the French Variscan Belt (Massif Central and Massif Armoricain): a consequence of the Early Carboniferous Gondwana-Laurussia collision. C R Géosci 349: 126–136. [CrossRef] [Google Scholar]
  • Feybesse JL, Lardeaux JM, Johan V, Tegyey M, Dufour E, Lemiere B. 1988. La série de la Brévenne (Massif central français): une unité dévonienne charriée sur le complexe métamorphique des Monts du Lyonnais à la fin de la collision varisque. C R Acad Sci Série 2307: 991–996. [Google Scholar]
  • Floc’h, JP, 1983. La série métamorphique du Limousin central: une traverse de la branche ligérienne de l’orogène varisque, de l’Aquitaine à la zone broyée d’Argentat (Massif central français). Unpublished thesis, University of Limoges, France.xx [Google Scholar]
  • Franke W, Cocks LRM, Torsvik TH. 2017. The Palaeozoic Variscan oceans revisited. Gondwana Res 48: 257–284. [CrossRef] [Google Scholar]
  • Gebauer D, Bernard-Griffiths J, Gruenenfelder M. 1981. U-Pb zircon and monazite dating of a mafic-ultramafic complex and its country rocks; example, Sauviat-sur-Vige, French Central Massif. Contrib Mineral Petrol 76: 292–300. [CrossRef] [Google Scholar]
  • Gébelin A, Martelet G, Chen Y, Brunel M, Faure M. 2006. Structure of late Variscan Millevaches leucogranite massif in the French Massif Central: AMS and gravity modelling results. J Struct Geol 28: 148–169. [CrossRef] [Google Scholar]
  • Grove TL, Baker MB. 1984. Phase equilibrium controls on the tholeiitic versus calc-alkaline differentiation trends. J Geophys Res Solid Earth 89: 3253–3274. [CrossRef] [Google Scholar]
  • Guillot F, Bertrand J-M., Bussy F, Lanari P, Cosma L, Pin C. 2012. Early Variscan I-type pluton in the pre-Alpine basement of the Western Alps: The ca. 360Ma Cogne diorite (NW-Italy). Lithos 153: 94–107. [CrossRef] [Google Scholar]
  • Guillot S, Menot RP. 2009. Paleozoic evolution of the External Crystalline Massifs of the Western Alps. C R Geosci 341: 253–265. [CrossRef] [Google Scholar]
  • Horstwood MSA, Košler J, Gehrels G, Jackson SE, McLean NM, Paton C, et al. 2016. Community-Derived Standards for LA-ICP-MS U-(Th-)Pb Geochronology − Uncertainty Propagation, Age Interpretation and Data Reporting. Geostand Geoanal Res 40: 311–332. [CrossRef] [Google Scholar]
  • Jolivet L, Baudin T, Calassou S, Chevrot S, Ford M, Issautier B, et al. 2021. Geodynamic evolution of a wide plate boundary in the Western Mediterranean, near-field versus far-field interactions. BSGF − Earth Sci Bull 192: 48. [CrossRef] [EDP Sciences] [Google Scholar]
  • Kroner U, Romer RL. 2013. Two plates − Many subduction zones: the Variscan orogeny reconsidered. Gondwana Res 24: 298–329. [CrossRef] [Google Scholar]
  • Lardeaux J-M. 2014. Deciphering orogeny: a metamorphic perspective Examples from European Alpine and Variscan belts: Part II: Variscan metamorphism in the French Massif Central − a review. Bull Soc Geol Fr 185: 281–310. [CrossRef] [Google Scholar]
  • Lardeaux JM, Schulmann K, Faure M, Janoušek V, Lexa O, Skrzypek E, et al. 2014. The moldanubian zone in the French Massif Central, Vosges/Schwarzwald and Bohemian Massif revisited: differences and similarities. Geol Soc Lond Spec Publ 405: 7–44. [CrossRef] [Google Scholar]
  • Ledru P, Lardeaux J-M., Santallier D, Autran A, Quenardel JM, Floc’h JP, et al. 1989. Où sont les nappes dans le Massif central français ? Bull Soc Géol Fr 605–618. [CrossRef] [Google Scholar]
  • Leloix C, Faure M, Feybesse J-L. 1999. Hercynian polyphase tectonics in the northeast French Massif Central: the closure of the Brévenne Devonian-Dinantian rift. Int J Earth Sci 88: 409–421. [CrossRef] [Google Scholar]
  • Li C-F., Li X-H., Li Q-L., Guo J-H., Li X-H., Yang Y-H. 2012. Rapid and precise determination of Sr and Nd isotopic ratios in geological samples from the same filament loading by thermal ionization mass spectrometry employing a single-step separation scheme. Anal Chim Acta 727: 54–60. [CrossRef] [Google Scholar]
  • Liao Y, Wei C, Rehman HU. 2021. Titanium in calcium amphibole: Behavior and thermometry. Am Mineral 106: 180–191. [CrossRef] [Google Scholar]
  • Lotout C, Pitra P, Poujol M, Anczkiewicz R, Van Den Driessche J. 2018. Timing and duration of Variscan high-pressure metamorphism in the French Massif Central: a multimethod geochronological study from the Najac Massif. Lithos 308: 381–394. [CrossRef] [Google Scholar]
  • Lotout C, Poujol M, Pitra P, Anczkiewicz R, Van Den Driessche J. 2020. From burial to exhumation: emplacement and metamorphism of mafic eclogitic terranes constrained through multimethod petrochronology, case study from the Lévézou Massif (French Massif Central, Variscan Belt). J Petrol 61: egaa 046. [Google Scholar]
  • Marcoux É, Barré B, Pichavant M, Poujol M. 2021. Âge et genèse de la coupole granitique à métaux rares (Sn, Li, Nb-Ta, W) de Montebras (Creuse, Massif central français). BSGF − Earth Sci Bull 192: 16. [CrossRef] [EDP Sciences] [Google Scholar]
  • Martínez Catalán JR, Schulmann K, Ghienne J-F. 2021. The Mid-Variscan Allochthon: Keys from correlation, partial retrodeformation and plate-tectonic reconstruction to unlock the geometry of a non-cylindrical belt. Earth-Sci Rev 220: 103700 [CrossRef] [Google Scholar]
  • Matte P. 1986. Tectonics and plate tectonics model for the Variscan belt of Europe. Tectonophysics 126: 329–374. [CrossRef] [Google Scholar]
  • Matte P. 2001. The Variscan collage and orogeny ( 480–290 Ma) and the tectonic definition of the Armorica microplate: a review. Terra Nova 13: 122–128. [CrossRef] [Google Scholar]
  • Matte P. 2007. Variscan thrust nappes, detachments, and strike-slip faults in the French Massif Central: Interpretation of the lineations. Geol Soc Am Mem 200: 391–402. [Google Scholar]
  • McCarthy A, Tugend J, Mohn G, Candioti L, Chelle-Michou C, Arculus R, et al. 2020. A case of Ampferer-type subduction and consequences for the Alps and the Pyrenees. Am J Sci 320: 313–372. [CrossRef] [Google Scholar]
  • McDonough WF, Sun SS. 1995. The composition of the Earth. Chem Geol 120: 223–253. [CrossRef] [Google Scholar]
  • Melleton J, Faure M, Cocherie A. 2009. Monazite U-Th/Pb chemical dating of the Early Carboniferous syn-kinematic MP/MT metamorphism in the Variscan French Massif Central. Bull Soc Géol Fr 180: 283–292. [CrossRef] [Google Scholar]
  • Molina JF, Moreno JA, Castro A, Rodríguez C, Fershtater GB. 2015. Calcic amphibole thermobarometry in metamorphic and igneous rocks: new calibrations based on plagioclase/amphibole Al-Si partitioning and amphibole/liquid Mg partitioning. Lithos 232: 286–305. [CrossRef] [Google Scholar]
  • Nance RD, Gutiérrez-Alonso G, Keppie JD, Linnemann U, Murphy JB, Quesada C, Strachan RA, Woodcock, NH. 2010. Evolution of the Rheic ocean. Gondwana Res 17: 194–222. [CrossRef] [Google Scholar]
  • Ohnenstetter M, Sider H. 1988. Contraintes géochimiques apportées par le magmatisme sur le développement du bassin marginal ensialique du Beaujolais au Dévonien. Bull Soc Géol Fr 4: 499–510. [CrossRef] [Google Scholar]
  • Paquette J-L., Monchoux P, Couturier M. 1995. Geochemical and isotopic study of a norite-eclogite transition in the European Variscan belt: implications for U/Pb zircon systematics in metabasic rocks. Geochim Cosmochim Acta 59: 1611–1622. [CrossRef] [Google Scholar]
  • Paquette J-L., Ballèvre M, Peucat J-J., Cornen G. 2017. From opening to subduction of an oceanic domain constrained by LA-ICP-MS U-Pb zircon dating (Variscan belt, Southern Armorican Massif, France). Lithos 294: 418–437. [CrossRef] [Google Scholar]
  • Paris F, Robardet M. 1990. Early Palaeozoic palaeobiogeography of the Variscan regions. Tectonophysics 177: 193–213. [CrossRef] [Google Scholar]
  • Pearce JA, Stern RJ, Bloomer SH, Fryer P. 2005. Geochemical mapping of the Mariana arc-basin system: Implications for the nature and distribution of subduction components. Geochem Geophys Geosys 6: Q07006. [CrossRef] [Google Scholar]
  • Pearce J, Stern B. 2006. Origin of Back-Arc Basin Magmas: Trace Element and Isotope Perspectives. Geophys Monog Ser 166: 63–86. [Google Scholar]
  • Pfeiffer M-T. 1985. La ligne tonalitique du Limousin : sa composition, sa signification, son implication dans la structuration varisque. Unpublished PhD Thesis, Université de Limoges, France. [Google Scholar]
  • Pfeiffer M-T. 1986. La signification de la ligne tonalitique du Limousin. Son implication dans la structuration varisque du Massif Central français. C R Acad Sci Série 2303: 305–310. [Google Scholar]
  • Pin C, Paquette J-L. 1997. A mantle-derived bimodal suite in the Hercynian Belt: Nd isotope and trace element evidence for a subduction-related rift origin of the Late Devonian Brévenne metavolcanics, Massif Central (France). Contrib Mineral Petrol 129: 222–238. [CrossRef] [Google Scholar]
  • Pin C, Paquette JL. 2002. Sr-Nd isotope and trace element evidence for a Late Devonian active margin in northern Massif-Central (France). Geodin Acta 15: 63–77. [Google Scholar]
  • Pochon A, Poujol M, Gloaguen E, Branquet Y, Cagnard F, Gumiaux C, et al. 2016. U-Pb LA-ICP-MS dating of apatite in mafic rocks: Evidence for a major magmatic event at the Devonian-Carboniferous boundary in the Armorican Massif (France). Am Mineral 101: 2430–2442. [CrossRef] [Google Scholar]
  • Salters VJM, Stracke A. 2004. Composition of the depleted mantle. Geochem Geophys Geosys 5: Q05B07. [Google Scholar]
  • Schmidt MW, Thompson AB. 1996. Epidote in calcalkaline magmas; an experimental study of stability, phase relationships, and the role of epidote in magmatic evolution. Am Mineral 81: 462–474. [CrossRef] [Google Scholar]
  • Schmidt MW, Poli S. 2004. Magmatic epidote. Rev Mineral Geochem 56: 399–430. [CrossRef] [Google Scholar]
  • Shaw A, Downes H, Thirlwall MF. 1993. The quartz-diorites of Limousin: elemental and isotopic evidence for Devono-Carboniferous subduction in the Hercynian belt of the French Massif Central. Chem Geol 107: 1–18. [CrossRef] [Google Scholar]
  • Sider H, Ohnenstetter M. 1986. Field and petrological evidence for the development of an ensialic marginal basin related to the Hercynian orogeny in the Massif Central, France. Geol Rundsch 75: 421–443. [CrossRef] [Google Scholar]
  • Sisson TW, Grove TL. 1993. Experimental investigations of the role of H2O in calc-alkaline differentiation and subduction zone magmatism. Contrib Mineral Petrol 113: 143–166. [CrossRef] [Google Scholar]
  • Skrzypek E, Tabaud A-S., Edel J-B., Schulmann K, Cocherie A, Guerrot C, et al. 2012. The significance of Late Devonian ophiolites in the Variscan orogen: a record from the Vosges Klippen Belt. Int J Earth Sci 101: 951–972. [CrossRef] [Google Scholar]
  • Sláma J, Košler J, Condon DJ, Crowley JL, Gerdes A, Hanchar JM, et al. 2008. Plešovice zircon — A new natural reference material for U–Pb and Hf isotopic microanalysis. Chem. Geol. 249: 1-35. [CrossRef] [Google Scholar]
  • Thiéblemont D, Cabanis B. 1986. Découverte d’une association de volcanites d’arcs et de basaltes de type “MORB” dans la formation paléo-volcanique silurienne de la Meilleraie, Vendée France. C R Acad Sci Paris Série 2232: 641–646 [Google Scholar]
  • Thiéblemont D, Cabanis B. 1994. Magmatic Record of the Geodynamic Evolution of Brittany and Vendée During Paleozoic Times: Trace Element Constraints, in: Chantraine J, Rolet J, Santallier DS, Piqué A, Keppie JD (Eds.), Pre-Mesozoic Geology in France and Related Areas Springer Berlin Heidelberg, Berlin, Heidelberg xx 220–230. [CrossRef] [Google Scholar]
  • Vanderhaeghe O, Laurent O, Gardien V, Moyen J-F., Gébelin A, Chelle-Michou C, et al. 2020. Flow of partially molten crust controlling construction, growth and collapse of the Variscan orogenic belt: the geologic record of the French Massif Central. BSGF − Earth Sci Bull 191: 25. [CrossRef] [EDP Sciences] [Google Scholar]
  • Vermeesch P. 2018. IsoplotR: A free and open tool box for geochronology. Geosci Front 9: 1479–1493. [CrossRef] [Google Scholar]
  • Wiedenbeck M, Alle P, Corfu F, Griffin WL, Meier M, Oberli F et al., 1995. 3 natural zircon standards for U-Th-Pb, Lu-Hf, trace-element and REE analyses. Geostand Newsl 19: 1–23. [CrossRef] [Google Scholar]

Cite this article as: Berger J, Beau-Hurdebourcq L, Serrano J, Benoit M, Grégoire M, Benmammar A, Duchene S, Bruguier O, Baele J-M. 2024. Short-lived active margin magmatism preceding Variscan collision in the Western French Massif CentralCourte pulsation du magmatisme de marge active avant la collision continentale varisque dans le Massif central français occidental, BSGF - Earth Sciences Bulletin 195: 7.

All Tables

Table 1

Location, name, structure and modal proportion of samples used in this study.

In the text

All Figures

thumbnail Fig. 1

Map of the Variscan belt in France showing the main Devonian-Early Carboniferous magmatic markers for subduction (Variscan zones modified after Ballèvre et al., 2014).
In the text
thumbnail Fig. 2

Geological maps showing the location of calc-alkaline plutons and their relation with tectono-metamorphic units. (a) Simplified geological map of the Limousin area adapted from published BRGM maps, the position of studied plutons is indicated together with samples labels. (b) Zoom on the St-Féréole-Tulle massif located within UGU but near the LGU/UGU contact. Foliation trajectories from Bellot and Roig (2007).
In the text
thumbnail Fig.3

Representative microphotographs of hornblendite TUL3 (a), hornblende gabbro TUL1 (b), hornblende-biotite quartz diorite SAJ12A (c) and hornblende-biotite tonalite MAR1 (d). Despite post-magmatic alteration, igneous textures are generally preserved (a, b, c) but a weak foliation can develop locally (d). Note the interstitial epidote-apatite-zircon-biotite cluster in sample MAR1 (d).
In the text
thumbnail Fig.4

Mineral composition and thermobarometry. a) Composition and classification of amphiboles (a); composition of plagioclases (b) and average P-T results per sample calculated with amphibole only and amphibole-plagioclase thermobarometers (see text) (c). The error bars represent the 1σ standard deviation of individual pressures and temperatures calculated for each sample.
In the text
thumbnail Fig.5

Classification of Late Devonian − Early Carboniferous calc-alkaline plutonic rocks from the French Massif Central and the Armorican Massif in the K2O-SiO2 diagram (a) and the AFM triangle plot (b). Literature data (white, black and grey symbols) are from Pfeiffer (1986); Shaw et al., (1993); Cuney et al., (2001) and Pin and Paquette (2002).
In the text
thumbnail Fig.6

Harker diagrams showing major element composition of calc-alkaline plutons in the French Massif Central and the Armorican massif. Literature data (white, black and grey symbols) are from 1986 (Pfeiffer); Shaw et al. (1993); Cuney et al., (2001) and Pin and Paquette (2002). Note the near linear trends drawn by diorites and tonalites extrapolating toward the composition of hornblende gabbros.
In the text
thumbnail Fig.7

Chondrite-normalized REE patterns and primitive mantle normalized multi-element diagrams for the calc-alkaline plutonic rocks investigated in this study. The compositions of volcanic rocks from the South Aegean arc are extracted from the Georoc database. Normalization values are from McDonough and Sun (1995).
In the text
thumbnail Fig.8

Sr-Nd isotopic diagram for calc-alkaline lavas/plutons from the French Massif Central (this study, Shaw et al., 1993; Cuney et al., 2001; Pin and Paquette, 2002) and for Brévenne back-arc magmatism (Pin and Paquette, 1997). Also shown, the composition of the Depleted Mantle (DM) at 360 Ma (calculated with the data of Salters and Stracke, 2004), the probable signature of the mantle wedge around 360 Ma as attested by calc-alkaline magmatism arc like magmatism. The green dotted mixing hyperbola represent mixing between Depleted Mantle and a UGU paragneiss to model mixing of the mantle source with sediments buried by subduction while the black mixing dotted curve simulates contamination of a depleted basalt by a UGU paragneiss to reproduce continental assimilation. Values in % indicate the relative mass of sedimentary component in the mix.
In the text
thumbnail Fig.9

Results of in situ LA-ICP-MS U-Pb zircon dating. (a) Cathodoluminescence images of representative zircon grains. Scale bars represent a length pf 20 μm and the quoted ages in Ma are calculated from 206Pb/238U ratios. (b, c, d) Concordia plots for tonalites SJL, MAR and FIG, respectively. MSWD and p-value are for concordance and equivalence. The two uncertainties are reported as 2σ.
In the text
thumbnail Fig.10

(a) Compilation of geochronological data obtained on petrological markers of Middle to Late Devonian subduction in the French Massif Central. U-Pb zircon on calc-alkaline magmatism are from Bertrand et al. (2001), Pin and Paquette (2002) and this study; age of Brévenne magmatism is from Pin and Paquette (1997), recent ages U-Pb zircon and Sm-Nd/Lu-Hf garnet ages obtained on eclogites are from Lotout et al. (2018, 2020) and de Hoÿm de Marien et al. (2023), UGU anatexis as dated by U-Th-Pb monazite chemical age is from Melleton et al. (2009) and the curve indicating the relative igneous activity is from Chopin et al. (2023) traced with data compiled by Vanderhaeghe et al. (2020). (b) Tentative snapshot of the geodynamic situation around 360 Ma for the Western French Massif Central, with emphasis on diorites-tonalites, eclogites and migmatitic rocks belonging to the UGU (see text for explanations).
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.