Chaîne varisque
Open Access
BSGF - Earth Sci. Bull.
Volume 191, 2020
Chaîne varisque
Numéro d'article 25
Nombre de pages 56
Publié en ligne 23 septembre 2020

© O. Vanderhaeghe et al., Published by EDP Sciences 2020

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Migmatites and granites are the main constituents of the continental crust and their petrogenesis and emplacement are intimately linked to orogenic evolution (Thompson and Connolly, 1995; Sawyer, 1998; Brown, 2001; Foster et al., 2001; Vanderhaeghe, 2009; Sawyer et al., 2011; Weinberg, 2016; Závada et al., 2018). Various heat sources have been proposed to cause high-temperature metamorphism and partial melting of orogenic roots comprising an increase in radioactive heat production of the thickened crust, an increase in the basal heat flux associated with delamination of the lithospheric mantle and heat advection through the emplacement of mantle-derived magmas (Houseman et al., 1981; England and Thompson, 1984; Henk et al., 2000; Annen and Sparks, 2002; Vanderhaeghe et al., 2003; Vanderhaeghe and Duchêne, 2010; Ueda et al., 2012). In turn, partial melting has a profound impact on the rheology at the scale of the rock and at the one of the entire crust (Vigneresse et al., 1996; Brown and Solar, 1998; Solar et al., 1998; Rosenberg, 2001; Vanderhaeghe and Teyssier, 2001a, 2001b; Gébelin et al., 2006; Schulmann et al., 2008; Vanderhaeghe, 2009), which is expressed by an intimate link between deformation and melt/solid segregation (Sawyer, 1994; Brown and Rushmer, 1997; Brown and Solar, 1998; Weinberg and Searle, 1998; Vanderhaeghe, 1999; Hasalová et al., 2008; Hasalova et al., 2011; Weinberg et al., 2013). Partial melting potentially triggers mechanical decoupling of the subducted continental crust from the downgoing slab as it has been proposed for example in the Norwegian Caledonides or the Variscan Bohemian Massif (Labrousse et al., 2011; Gordon et al., 2016; Závada et al., 2018). Partial melting has also been identified as the key parameter controlling lateral flow of the deep root of orogenic belts through the activation of vertical shear zones (Solar et al., 1998; Weinberg and Mark, 2008) leading to the formation of orogenic plateaux (Vanderhaeghe and Teyssier, 2001a, 2001b; Vanderhaeghe et al., 2003; Gerbault et al., 2005; Cagnard et al., 2006; Chardon et al., 2009). Finally, the ubiquitous presence of large domes cored by migmatites in the exhumed roots of orogenic belts (Whitney et al., 2004; Vanderhaeghe, 2009) as well as the spatial-temporal correlations between the emplacement of granitic laccoliths and the activation of low-angle detachments (Lister and Baldwin, 1993; Vanderhaeghe, 1999; Searle et al., 2009; Whitney et al., 2013) suggest that the presence of a partially molten crust and the migration of granitic melts controls the behaviour of the orogenic crust during orogenic gravitational collapse.

The Variscan belt of Western Europe (Fig. 1) has long been recognized as particularly rich in migmatites and granitoids (e.g. Zwart, 1967) and is thus the perfect target to investigate the impact of partial melting and magmatism on orogenic evolution. The large spread of radiochronological ages obtained on migmatites and magmatic rocks (from c. 390 Ma to c. 290 Ma) and the wide variety of petrological and geochemical characteristics of the magmatic rocks indicate that the 100 Ma long tectonic history of the Variscan belt has been punctuated by the emplacement of magmas implying the contribution of both the crust and the mantle (Letterrier, 1978; Cuney et al., 1990; Solgadi et al., 2007; Bussien et al., 2008; Couzinié et al., 2014; von Raumer et al., 2014; Laurent et al., 2017). In addition to the identification of the sources and geodynamic context of these magmas, arises the question of the impact of low-viscosity and low-density silicate melts on the dynamic evolution of the Variscan orogen.

Despite the significant exposure of migmatites and granitoids, their implication on the tectonic evolution of the Variscan belt has not been fully explored, to the noticeable exceptions of some papers on the French Massif Central (Malavieille et al., 1990; Burg and Vanderhaeghe, 1993; Costa and Rey, 1995; Vanderhaeghe et al., 1999) and Central Iberia, the Vosges, and Bohemia (Henk, 2000; Henk et al., 2000; Schulmann et al., 2008, 2014; Lardeaux et al., 2014; Rubio Pascual et al., 2016). The goal of this paper is to discuss the impact of partial melting and magmatism on the tectonic evolution of the Variscan orogenic belt of Western Europe based on a synthesis of structural, petrological, geochemical, geochronological and sedimentological data available for the Variscan basement of the French Massif Central. This region offers a unique section through the Variscan crust that recorded, from Silurian to Permian, a prolonged history of burial and exhumation associated with the construction and destruction of the belt, respectively. The tectonic evolution is particularly marked by the generation of migmatites under HP, MP and LP metamorphic conditions and by varied plutonic rocks emplaced from the middle Devonian to the early Permian. This history is complemented by the P-T-t record of lower crustal xenoliths brought back to the surface by Cenozoic volcanoes and by unmetamorphosed to low-grade volcano-sedimentary series deposited from the middle Devonian (Givetian) to the Permian that constrain the topographic evolution of the belt.

In this paper, we propose a new geodynamic model for the generation and flow of partially molten rocks and magmas during orogenic evolution from construction by tectonic accretion of subducted continental units followed by lateral growth of the orogenic belt associated with construction of an orogenic plateau and eventually to gravitational collapse. This geodynamic model is also nourished by new data recently published in companion papers comprising (i) Lu-Hf tracing of igneous and detrital zircon in gneisses and plutonic rocks (Chelle-Michou et al., 2017; Couzinié et al., 2017, 2019), (ii) U-Pb dating of zircon and monazite by LA-ICP-MS on the Carboniferous plutonic rocks (Chelle-Michou et al., 2017; Laurent et al., 2017) and (iii) a detailed petrogenetic model for these granitoids (Moyen et al., 2017).

thumbnail Fig. 1

Tectonic map of the Variscan belt in Western Europe. Continental terranes, Avalonia, Saxo-Thuringia, Armorica-Barandia, Brunia, are separated by ophiolitic sutures, namely the Rheic suture and the Medio-European suture, to the north and south of Armorica respectively. The internal zone of the Variscan belt comprises high-grade nappes overlying a parautochtonous unit belonging to the northern margin of Gondwana.

2 Geology of the French Massif Central: a window through the Variscan belt

2.1 The Variscan belt: continental blocks, oceanic sutures, allochthonous terranes and paleogeographic reconstructions

The Variscan belt has been first defined as a post-Cambrian and pre-Permian mobile belt based on analysis of stratigraphic unconformities, structures and nappes (Suess, 1883; Bertrand, 1887). It is extending from East Asia to the tip of South America running through Central Europe and along the edges of North America (Matte, 2001). Paleomagnetic data indicate that the Variscan belt formed as a consequence of the convergence between Laurussia (Laurentia + Baltica) and Gondwana resulting in the Pangea supercontinent (Scotese and McKerrow, 1990; Unrug, 1997; Tait et al., 2000). However, the number of oceanic sutures and the former sizes of the oceanic basins are discussed as developed below.

Pioneer work correlating data from the eastern (Bohemia) and western (Iberia) terminations of the belt led to (i) the definition of the main geological-tectonic zones, (ii) the identification of the principal Paleozoic tectonic events based on the relationship between structures and Ediacaran to Carboniferous sedimentary deposits in the external zone and along the foreland, and (iii) the recognition of high-grade nappes in the internal zone (Stille, 1924; Kossmat, 1927; Gaertner, 1937; Demay, 1948; Gèze, 1949). The relationship between sedimentation and deformation, best exposed along its external domains, allowed the identification of continental blocks such as Avalonia, Armorica, Saxo-Thuringia, Barrandia and Brunia, all of which preserve Cambrian unconformities more or less affected by Variscan deformation and metamorphism (Matte, 1986, 1991; Franke and Engel, 1986; Franke, 1989, 2000; Kroner and Romer, 2013). These continental blocks are separated by oceanic sutures marked by ophiolitic assemblages or mélanges. All authors agree on the Rheic suture, also called the Lizard-Rheno-Hercynian suture, that corresponds to the former Rheic Ocean, to the south of Avalonia and to the north of Saxo-Thuringia and Armorica (Fig. 1) (Matte, 1991; Franke, 2000; Ballèvre et al., 2014; Franke et al., 2017). The nature of the terranes and the presence of sutures south of Armorica is however debated. Despite uncertainties, we favor the existence of multiple sutures (designated as “secondary sutures” on Fig. 1) based on (i) the presence of ultramafic and mafic high-pressure rocks of different ages at different structural levels (Girardeau et al., 1986; Dubuisson et al., 1989; Faure et al., 1997; Bosse et al., 2000; Ballèvre et al., 2009; Berger et al., 2006, 2010a; Lardeaux, 2014; Lardeaux et al., 2014; Lotout et al., 2018), and (ii) the occurrence of remnants of subordinate Devonian rift and/or oceanic basins (Sider and Ohnenstetter, 1986; Skrzypek et al., 2012).

In the French Massif Central, until the middle of the 20th century, the prevailing model attributed the high-grade granitic-gneisses to a crystalline basement and the paragneisses and schists to a sedimentary cover deposited from the Neoproterozoic throughout the Paleozoic (Jung, 1953; Roques, 1971). In the absence of geochronological data, metamorphism was considered as polycyclic, the basement being metamorphosed during the Neoproterozoic, and then, together with its sedimentary cover, during the Caledonian and the Variscan orogenies (Forestier, 1961; Chenevoy and Ravier, 1971; Roques, 1971). This view has been profoundly modified first by the results of absolute dating demonstrating ubiquitous Variscan reworking (Gebauer et al., 1981; Pin and Lancelot, 1982; Rolin et al., 1982; Duthou et al., 1994). Another turning point was to interpret the felsic/mafic gneisses with a tholeiitic signature, defined as the Leptynite-Amphibolite Complex (LAC), as remnants of ophiolites (Forestier, 1961; Briand and Piboule, 1979; Cabanis et al., 1983; Maillet et al., 1984; Mercier et al., 1985; Piboule and Briand, 1985; Briand et al., 1988; Dubuisson et al., 1989; Pin, 1990). In the Armorican Massif, ophiolites define the Medio-European suture also referred to as the Galicia-South Brittany suture, or the Eo-Variscan suture (Bernard-Griffiths and Cornichet, 1985; Hanmer, 1977; Faure et al., 1997; Ballèvre et al., 2009). This suture corresponds to the former Medio-European Ocean (also designated as the Galicia-Massif Central Ocean or the Paleotethys Ocean) that is located to the south of the Armorica-Barrandia continental block (Matte, 1986, 2001; Pin, 1990; Stampfli et al., 2013) (Fig. 1).

Based on this analysis, most authors agree that the Variscan belt of Western Europe results from tectonic accretion of the Avalonia and Armorica ribbon-shaped continental terranes and the closure of intervening oceanic basins, marked by (i) the Iapetus-Tornquist suture north of Avalonia, (ii) the Rheic/Rheno-Hercynian suture between Avalonia and Armorica (Matte, 2001; von Raumer et al., 2003) and (iii) the Medio-European suture between Armorica and Gondwana (Matte, 1986, 2001; Pin, 1990; Stampfli et al., 2013) (Figs. 2A and 2B). All models invoke closure of the Iapetus Ocean during the Ordovician concomitantly to the rapid opening of the Rheic Ocean (e.g. Hamilton and Murphy, 2004) by reactivation of a Neoproterozoic suture along the northern Gondwana margin (Linnemann et al., 2007). The driving force for this oceanisation has been attributed either to slab pull from the northward subducting Iapetus (Murphy et al., 2006) or to roll-back of a southward subducting slab beneath the Gondwana margin (Martínez Catalán et al., 2009). Paleogeographic reconstruction differ for the early Ordovician period with regard to the position of Armorica that has implications on the existence and size of the Medio-European Ocean (Tait et al., 2000; Matte, 2001; Martínez Catalán et al., 2007, 2009; Faure et al., 2008; Ballèvre et al., 2009; Shail and Leveridge, 2009; Nance et al., 2010; Torsvik et al., 2012; Kroner and Romer, 2013; Stampfli et al., 2013; Lardeaux, 2014; Lardeaux et al., 2014; Schulmann et al., 2014; Skrzypek et al., 2014; Edel et al., 2018). Indeed, sedimentologic and paleontological data lead to paleogeographic reconstructions indicating a position close to the South Pole for Gondwana together with Avalonia and Armorica, and between the tropics and the equator for Laurentia and Baltica (Paris and Robardet, 1990; Robardet et al., 1993; Fortey and Cocks, 2003; Robardet, 2003). This is consistent with paleomagnetic data indicating an Iapetus Ocean at least 3000 km wide (Tait et al., 2000; Cocks and Torsvik, 2002; Hamilton and Murphy, 2004). On the other hand, some paleomagnetic data indicate that Armorica remained attached to the Gondwana margin until the Devonian (Kössler et al., 1996), which is consistent with the continuity of benthic faunas from Armorica to the northern margin of Gondwana (Robardet, 2003) while other paleomagnetic data indicate that Armorica moved toward Laurussia during the Ordovician (Tait et al., 2000), implying closure of the Rheic Ocean and opening of a 2000–3000 km wide Medio-European Ocean (Tait et al., 1997; Cocks and Torsvik, 2006; Shaw and Johnston, 2016). As a result, two paleogeographic reconstructions clash, one favoring a large Medio-European Ocean between Gondwana and a continental ribbon including Armorica (Matte, 2001; von Raumer et al., 2003; Stampfli and Borel, 2004; Domeier and Torsvik, 2014; Domeier, 2016), while the other considers at most a small immature rift in this region and infer that the large oceanic realm is the Rheic Ocean (Martínez Catalán et al., 2007; Nance et al., 2010; Kroner and Romer, 2013) (Figs. 2C and 2D). On the other hand, tectonic reconstructions based on geological data implying multiple rifts and/or oceanic basins in between Armorica and Gondwana (Girardeau et al., 1986; Dubuisson et al., 1989; Faure et al., 1997; Lardeaux, 2014; Lardeaux et al., 2014), correspond to an intermediate proposition.

thumbnail Fig. 2

Geodynamic-tectonic models and paleogeographic reconstructions for evolution of the Variscan belt of Western Europe. A. The double subduction model (modified after Matte, 1986, 1991, 2001). B. The polycyclic model (modified after Faure et al., 1997, 2002, Faure et al., 2009a, 2009b and Lardeaux et al., 2014). C. Single (Rheic) Ocean model of Paris and Robardet (1990); Martínez Catalán et al. (2001); Nance et al. (2010). D. Mutliple Oceans (Rheic, Medio-European, …) model of Tait et al. (1997); Matte (2001); Stampfli and Borel (2004); Domeier and Torsvik (2014).

2.2 The main lithologic-tectonic units of the French Massif Central

The French Massif Central is one of the largest exposures of the Variscan belt of Western Europe. Its geology is synthesized with reference to 1:1 000 000 scale map (Chantraine et al., 1996, 2003) and to the lithologic-tectonic units defined on the basis of their distinct lithological, structural, metamorphic and geochronological record (Burg and Matte, 1978; Matte, 1986; Ledru et al., 1989; Quenardel et al., 1991; Faure et al., 2009a; Lardeaux, 2014). According to these authors, the French Massif Central is made of (i) low- to high-grade metamorphic nappes with a Devonian to Carboniferous tectonic record, (ii) late Devonian to late Carboniferous plutonic rocks, (iii) late Devonian to Carboniferous volcanic and carbonate to clastic sedimentary rocks affected by low-grade metamorphism, and (iv) unconformable late Carboniferous and Permian detrital sediments (Fig. 3). The French Massif Central is subdivided into a western part and an eastern part by the NNE-SSW trending Sillon Houiller Fault (Arthaud and Matte, 1975; Feybesse, 1981) (Fig. 3). The latter crosscuts a set of NW-SE to NNW-SSE trending dextral shear zones that is particularly well-developed in the western part of the French Massif Central and connects to the South Armorican shear zone (Lerouge and Quenardel, 1988; Gébelin et al., 2007). The nappe pile is also cross-cut by low-angle detachements (Malavieille et al., 1990; Burg et al., 1993; Faure, 1995; Gardien et al., 1997; Bellot, 2007).

thumbnail Fig. 3

Geological map of the French Massif Central (modified after Chantraine et al., 2003). Metamorphic rocks comprise (i) a low-grade metasedimentary sequence attributed to the Lower Paleozoic (Cambrian to Lower Carboniferous), (ii) micaschists and paragneisses of uncertain age but considered as Neoproterozoic to Lower Paleozoic deposits, (iii) orthogneisses with a Cambrian or Ordovician age. The Middle Allochthonous Unit, the Parautochtonous Unit and the Lower Gneiss Unit are made of these metasedimentary rocks and orthogneisses. The Upper Gneiss Unit is characterized by relics of granulite facies metamorphism and by enclaves-boudins of mafic and ultramafic rocks affected by eclogite facies metamorphism, designated as the Leptynite-Amphibolite Complex (LAC). The color scheme of magmatic rocks is indicative of their age. Granitoids are distinguished according to their petrologic types in Figure 6.

2.2.1 Metamorphic nappes

The interpretation of the LAC as a suture marking the boundary between an Upper Gneiss Unit (UGU) and a Lower Gneiss Unit (LGU) led to the reinterpretation of the structure of the FMC in terms of nappe stacking (Burg and Matte, 1978; Ledru et al., 1989) but also to comparisons between the Variscan belt of Western Europe and the Himalaya-Tibet orogen (Dewey and Burke, 1973; Mattauer and Etchecopar, 1976; Burg and Matte, 1978; Autran and Cogné, 1980; Bard et al., 1980; Matte, 1986, 1991). Since its first description, the definition of the LAC has fluctuated and appears to cover a variety of rocks, with emphasis either on the bimodal magmatic association or on the high-pressure metamorphism affecting the mafic and ultramafic rocks (see discussion in Santallier et al., 1988). In addition, the LAC has been recognized at different structural positions and characterized by various metamorphic conditions and ages that led to the addition of a Middle Allochthonous Unit (MAU) in an intermediate position between the UGU and LGU (Girardeau et al., 1986; Dubuisson et al., 1989; Berger et al., 2010a, 2010b; Lotout et al., 2018). The proposed superposition of the UGU (and of the MAU) over the LGU relies on contrasted metamorphic record delineating an inverted metamorphic gradient and locally on the identification of tectonic contacts (Faure et al., 1979; Burg et al., 1984). Based on these characteristics, the main units of the metamorphic nappe stack are, from top to bottom (Tab. 1):

  • The Upper Gneiss Unit (UGU): It is made of garnet/cordierite-bearing diatexites associated with metatexites derived from orthogneisses and paragneisses with relics of granulite facies mineral parageneses. It has recorded a typical maximum pressure of 10 kbar for a temperature up to 900 °C retrogressed into amphibolite facies and greenschist facies pointing to decompression and cooling (Audren et al., 1987; Lardeaux et al., 2001; Schulz et al., 2001; Bellot and Roig, 2007; Schulz, 2009) (Fig. 4). Geochronological data are consistent with a Devonian to early Carboniferous age for these HP granulite facies migmatites (Duthou et al., 1981, 1994; Lafon, 1986; Schulz, 2014; Do Couto et al., 2016). The Upper Gneiss Unit contains boudins and enclaves of the Leptynite Amphibolite Complex (LAC) that have preserved HT eclogitic relics first discovered in the Haut-Allier and in the Rouergue (Forestier, 1961; Lasnier, 1968; Nicollet, 1977) but then found in most if not all of the LAC enclosed in the UGU (Burg and Matte, 1978; Gardien and Lardeaux, 1991; Mercier et al., 1991a, 1991b; Gardien, 1993). These eclogites have typically recorded PT conditions of c. 15 kbar for c. 750 °C (Santallier, 1981; Le Breton et al., 1986; Godard, 1990; Bellot and Roig, 2007) but more extreme conditions above 28 kbar are reported in the Lyonnais (Lardeaux et al., 2001) (Fig. 4). The few available geochronological data for the French Massif Central have served to propose a late Silurian to Devonian age for this HT eclogite facies metamorphism (Pin and Lancelot, 1982; Ducrot et al., 1983; Paquette et al., 1995; Do Couto et al., 2016).

    thumbnail Fig. 4

    Synthetic P-T-t paths of the Upper Gneiss Unit, Middle Allochthonous Unit, and of the Lower Gneiss Unit. The Upper Gneiss Unit is made of granulitic migmatites enclosing mafic and ultramafic enclaves of the Leptynite-Amphibolite Complex, whereas the Middle Allochthonous is not migmatitic. A. The Leptynite-Amphibolite Complex contains UHP and eclogitic relics dated at the transition from late Silurian to late Devonian as indicated by the numbers in the circles (c. 432 to c. 377 Ma). B. These rocks are retrogressed into the granulite facies dated from the middle to the late Devonian (c. 384 to c. 360 Ma) in the migmatitic gneisses of the Upper Gneiss Unit. C. The Lower Gneiss Unit comprises rare HP/LT rocks in the Limousin dated at c. 376 Ma. Rocks of the Lower Gneiss Unit have also recorded early Carboniferous (348–320 Ma) and late Carboniferous (304–302 Ma) LP amphibolite facies metamorphism.

  • The Middle Allochthonous Unit (MAU): Identified in Limousin, it consists mainly of micaschists associated with an ophiolitic assemblage that is similar to the Leptynite Amphibolite Complex. In contrast to the Upper Gneiss Unit, the LAC of the MAU is marked by relics of LT eclogite facies metamorphism (Girardeau et al., 1986; Dubuisson et al., 1989; Berger et al., 2010a, 2010b). The maximum pressure recorded by a garnet-kyanite assemblage is 29 kbar but the temperature did not exceed 660 °C consistent with subduction to 100 km depth (Berger et al., 2010a) (Fig. 4). The pressure peak in Najac is 18 kbar and the temperature only reached 600 °C (Lotout et al., 2018). These rocks have then been retrogressed into amphibolite facies. Published geochronological data for this LT eclogite facies metamorphism point to a late Silurian age in Central Limousin (Berger et al., 2010a) and middle Devonian in Najac (Lotout et al., 2018).

  • The Lower Gneiss Unit (LGU): It is made of micaschists, paragneisses and felsic orthogneisses with minor carbonates and mafic rocks. The LGU is affected by greenschist to amphibolite facies metamorphism depicting a Carboniferous Barrovian MP/MT gradient reaching partial melting locally masked by retrogression into greenschist facies (Nicollet, 1978; Burg et al., 1984, Burg et al., 1989a, 1989b; Schulz et al., 2001; Gébelin et al., 2009; Gardien et al., 2011) (Fig. 4).

  • The Para-Autochtonous Unit (PAU): It is made of micaschists, paragneisses and orthogneisses affected by pervasive Carboniferous greenschist facies metamorphism locally preserving relics of amphibolite facies (Bellot and Roig, 2007) (Fig. 4).

These high-grade nappes, identified throughout the French Massif Central, are structurally overlain by Metamorphic Upper Units (MUU) with limited lateral extent:

2.2.2 Protoliths age for the high-grade rocks, basement of the Devonian to Carboniferous volcanic and sedimentary units

The first proposed ages for the different lithological units exposed in the French Massif Central relied on relative chronology based on stratigraphy, biostratigraphy and cross-cutting relationships (Jung, 1953; Roques, 1971). The intensity of metamorphism and deformation affecting rocks forming the UGU, MAU, and LGU precludes the identification of fossils in high-grade metasedimentary rocks and the recognition of initial magmatic-sedimentary contacts. Thus, protoliths age determination relying mostly on analogies with stratigraphic ages determined in other regions are subject to caution. Recent age determinations based on radiometric dating of inherited zircon grains (distinguished from Variscan metamorphic zircon) provide (i) a maximum age (the youngest age obtained on inherited grains) for the deposition of the sediments protoliths of paragneisses, and (ii) a crystallization age for magmatic rocks, protoliths of orthogneisses, and/or amphibolites and other mafic metaplutonics. Ages spreading from the late Neoproterozoic (Ediacaran) to the early Paleozoic (Cambrian to Ordovician) have been attributed to most of the metasedimentary sequences (Franke, 2000; Chantraine et al., 2003; Linnemann et al., 2007). Protoliths ages of metamorphic rocks (Pin and Duthou, 1990; R’Kha Chaham et al., 1990; Duthou et al., 1994; Melleton et al., 2010; Chelle-Michou et al., 2017; Couzinié et al., 2017, 2019; Lotout et al., 2018) (Fig. 5) indicate (i) an Ediacaran to Ordovician age for metasediments deposited along the former margin of the Gondwana continent, (ii) Ediacaran, Cambrian and Ordovician ages for orthogneisses, and (iii) a Cambrian to Ordovician age for the mafic and ultramafic rocks (Tab. 2). According to these data, the only identified Proterozoic basement in the French Massif Central corresponds to the paragneisses and orthogneisses with Ediacaran protoliths.

thumbnail Fig. 5

Geochronological constraints on the emplacement and deposition ages of pre-Variscan rocks. For orthogneisses, the emplacement age corresponds to the Rb-Sr whole rock isochron age (circles) or the the youngest U-Pb, U-Th-Pb age obtained on magmatic zircon or monazite. For paragneisses, a maximum deposition age is provided by the youngest inherited age. Cambrian and Ordovician orthogneisses do not appear associated with a given nappe but are found throughout the French Massif Central.

2.2.3 Migmatites and granulites

In addition to the HP granulite facies migmatites that are part of the UGU, MP and LP migmatites, developed to the expense of late Neoproterozoic to early Paleozoic paragneisses and orthogneisses, are exposed at the lowest structural level of the nappe pile in the Limousin, Velay and Montagne Noire (Williamson et al., 1996; Downes et al., 1997; Ledru et al., 2001; Gébelin et al., 2009; Barbey et al., 2015; Villaros et al., 2018). In the Limousin, PT conditions of 850–750 °C and 5–6 kbar have been estimated for the cordierite-bearing migmatitic gneisses and associated lenses of leucosomes of the Millevaches massif in the Limousin dated at c. 315 Ma (Gébelin et al., 2009). In the Montagne Noire, migmatites coring the Caroux dome have recorded P-T conditions of 6 kbar and 720 °C (Rabin et al., 2015). They contain magmatic zircon with U-Pb ages spreading from c. 330 to c. 300 Ma (Franke et al., 2011; Faure et al., 2014; Roger et al., 2015). In the Velay dome, migmatites display contrasting P-T conditions and ages as a function of their structural position. Along the margin of the dome, migmatites have recorded partial melting with biotite remaining stable at a pressure over 5 kbar for a temperature over 750 °C (Montel et al., 1992). These migmatites contain monazite that yield microprobe U-Pb ages ranging from c. 329 to 314 Ma (Mougeot et al., 1997; Cocherie et al., 2005; Bé Mézémé et al., 2006). In contrast, migmatites coring the Velay dome have recorded fluid-absent biotite melting in the cordierite stability field with P-T estimates of 2–5 kbar and 760–850 °C (Montel et al., 1992; Barbey et al., 1999, 2015; Villaros et al., 2018). These migmatites display a Rb-Sr whole rock age of 298 ± 8 Ma (Rb–Sr whole-rock, Caen-Vachette et al., 1981) and contain monazite dated at 301 ± 5 Ma by U-Pb (Mougeot et al., 1997).

The nature of the lower crust underlying these migmatites is documented by xenoliths included in late Variscan plutonic rocks and in Cenozoic basaltic lavas. Biotite-sillimanite enclaves incorporated in diorite intrusive in the southern margin of the Velay dome have recorded P-T conditions of 8–10 kbar for 700–800 °C (Montel, 1985). The main mineral paragenesis of felsic and mafic granulites xenoliths in Cenozoic lavas indicate a pressure of 8–10 kbar for a temperature of 700–800 °C followed by isothermal decompression at 5–6 kbar (Leyreloup, 1974; Downes and Leyreloup, 1986). A few metaigneous xenoliths have preserved relics pointing to a pressure as high as 14 kbar for a temperature of c. 900 °C. These granulites contain zircon grains that yield U-Pb ages ranging from 320 to 280 Ma (Downes et al., 1991). The geochemical signatures of the granulites with a metasedimentary protolith indicate a refractory or residual character while the mafic and felsic metaigneous granulites display characteristics of a calc-alkaline liquid and of a cumulate, respectively (Dupuy et al., 1979; Downes et al., 1990; Vielzeuf and Vidal, 2012). The representation of the deep part of the cross sections (Fig. 8) is in part based on these characteristics.

2.2.4 Late Devonian to Carboniferous plutonic rocks

The c. 70 to 80 Ma-long magmatic activity recorded by plutonic and minor volcanic rocks of the French Massif Central is characterized by an extreme variety of petrological types and geochemical signatures. In this section, we describe these petrologic and geochemical characteristics and discuss their significance integrating their spatial distribution and their ages of emplacement (Tab. 3). The petrology and geochemistry of granitoids, as described in the 1/1 000 000 geological map of France (Chantraine et al., 2003) follows the classification by Barbarin (1999). Four types are distinguished (Fig. 6):

  • Calc-alkaline granitoids equivalent to the amphibole-rich calc-alkaline granitoids (ACG) of Barbarin (1999) include tonalites, granodiorites and granites. In addition to amphibole, they may contain pyroxene, and frequently include enclaves and small volumes of gabbro to diorites (i.e. more mafic terms of the same series). They are essentially similar to arc granitoids typically attributed to fractional crystallization of an enriched mafic magma.

  • Peraluminous granites and leucogranites correspond to muscovite or muscovite + biotite bearing granites sensu-stricto equivalent to the muscovite-bearing peraluminous granitoids (MPG) of Barbarin (1999). They are attributed to relatively cold (< 850 °C, water present- or muscovite breakdown) melting of metasediments (Gardien et al., 1995; Villaros et al., 2018) and might represent the first melt extracted from the partially molten source at the onset of partial melting.

  • Peraluminous granites to granodiorites, corresponding to the cordierite-bearing peraluminous granitoids (CPG) of Barbarin (1999). They probably relate to relatively hot (≥ 850 °C, biotite breakdown) melting of a continental felsic source (ortho or paragneisses) (Gardien et al., 1995; Barbey et al., 1999; Villaros et al., 2018). They might represent a partially molten source with a high-melt fraction (i.e. diatexites) implying inefficient melt and/or magma extraction but possibly also some solid settling.

  • High-K sub-alkaline granitoids, i.e. K-feldspar porphyritic calc-alkaline granitoids (KCG) of Barbarin (1999). They are porphyritic granites to granodiorites, commonly amphibole-bearing, and they typically contain abundant accessory minerals such as apatite and titanite. They contain micro-granular mafic enclaves, and are associated with intermediate plutonic rocks (diorites, to tonalites, to monzodiorites) of similar, “vaugneritic” (see below) composition to which they are probably petrogenetically related (Moyen et al., 2017).

In addition, granites and migmatites are associated with mafic but potassic plutonic rocks locally known as “vaugnerites” (Michon, 1987; Sabatier, 1991). Vaugnerites range from diorites to syenites and consist of amphibole, biotite, clinopyroxene, plagioclase, rare orthopyroxene, interstitial K-feldspar and quartz (Sabatier, 1991). Vaugnerites are K-, LILE- and LREE-rich mafic to intermediate rocks, pointing to an origin by partial melting of a mantle source enriched by the addition of crustal components, probably during earlier subduction (Rapp et al., 2010; Couzinié et al., 2014, 2016). They form the most mafic components of the KCG suites. Interestingly, they are undistinguishable from CPG and MPG from an isotopic point of view (Sr, Nd or Hf) (Williamson et al., 1992; Couzinié et al., 2016; Moyen et al., 2017).

The study of vaugnerites and their counterparts in other orogenic settings worldwide indicate that their characteristics unlikely result from the crustal contamination of basaltic magma on their way to the surface and are rather primarily inherited from a mantle source enriched by crustal components (Williams, 2004; Laurent et al., 2011, 2014; Prelević et al., 2012; Campbell et al., 2014; Couzinié et al., 2016). It should be mentioned that vaugnerites are similar to the “durbachites” described in other parts of the Variscan belt (Sabatier, 1991; von Raumer et al., 2014). Such rock types are present in most if not all orogenic settings elsewhere in the world where they are called “appinites”, “redwitzite”, “high Sr-Ba granitoids” (Fowler et al., 2001) or “sanukitoids” in an Archean context (Martin et al., 2005; Heilimo et al., 2010). In the FMC, the isotopic similarity between vaugnerites and crust-derived granites shows that the crustal component was derived from the local crust, probably during (continental) subduction shortly prior to melting (Couzinié et al., 2016; Moyen et al., 2017). A similar model is proposed for the c. 345 Ma old KCG granitoids (the Blatna suite) of the Bohemian Massif (Janousek and Holub, 1997; Janoušek et al., 2004), suggesting that KCG may derive from similar petrogenetic processes as vaugnerites and typically represent their differentiated products. However, this is not so clear for the FMC where KCG may also derive from interactions between vaugnerites and melts derived from the local crust, such as CPG and MPG (Solgadi et al., 2007; Laurent et al., 2017; Moyen et al., 2017). In either case, KCG are genetically linked to the vaugnerites so they are classified together with them in the following as “mantle-derived” granitoids.

thumbnail Fig. 6

Petrology of granitoids of the French Massif Central. Plutons are distinguished on the basis of the petrology of their dominant facies according to the classification of Barbarin (1999). ACG-type (calc-alkaline granites) are interpreted as arc magmas originated in an Andean-type continental margin by partial melting of an enriched mantle contaminated by the crust of the upper plate and/or mixed with crutal magmas. MPG-type (or muscovite-bearing peraluminous granites) are attributed to muscovite dehydration- or water present melting of a dominantly metasedimentary source. CPG-type (or cordierite-bearing peraluminous granites) are attributed to biotite dehydration melting of orthogneisses. KCG-type (or K-rich calc-alkaline granites) typically contain abundant micromafic enclaves and are attributed to mixing between magmas generated by partial melting of the crust and a magma generated by partial melting of an enriched lithospheric mantle represented by Mg-K diorites (the so-called vaugnerites).

2.2.5 Carboniferous volcanic and sedimentary deposits

In the northern part of the French Massif Central, cystalline rocks are locally capped by middle Carboniferous (Visean: c. 347–325 Ma) volcanic and marine deposits, as exemplified in the Sioule, Morvan and Brévenne regions, implying that they were exhumed and below sea-level at this time (Bertaux et al., 1993; Franke, 2014). In the southern part of the French Massif Central, Devonian to mid-Carboniferous (Visean) carbonates and turbidites are unconformably deposited on plutonic rocks and Ediacaran to Ordovician orthogneisses and paragneisses, in an underfilled foreland basin at the front of a propagating, low-grade fold-and-thrust system (Franke and Engel, 1986; Souquet et al., 2003). Namurian (c. 330–325 Ma) olistoliths and probably Westphalian (c. 325–304 Ma) turbidites and coarse conglomerate deposits attest for the erosion of a growing mountain belt to the north (Engel et al., 1978, 1981; Engel, 1984).

Clastic sediments associated with minor volcanics of late Carboniferous (Stephanian) and Permian ages are unconformably deposited on top of the crystalline rocks in extensional basins associated with strike-slip shear zones marking the waning stages of the Variscan orogeny (Ménard and Molnar, 1988; Van Den Driessche and Brun, 1992; Becq-Giraudon, 1993; Becq-Giraudon et al., 1996). Presence of coal in the Stephanian basins attests for a high-geothermal gradient (Copard et al., 2000), which could be attributed to the juxtaposition of the sediments to high-grade rocks freshly exhumed along low-angle detachments. Permian basins are typically wider than Stephanian ones suggesting progressive peneplanation of the Variscan topography at the end of the Carboniferous.

In summary, the superposition of the UGU over the MAU, LGU and PAU define an inverted metamorphic gradient with HP granulite facies migmatites overlying amphibolite facies paragneisses and orthogneisses with locally preserved LT eclogites. The structural, petrological and geochronological record of these nappes document a diachronous history of burial, exhumation and emplacement spreading from the late Silurian to the Devonian. The top of the nappe pile is locally marked by Metamorphic Upper Units (MUU). The lowest structural level is composed of MP to LP migmatites exposed in the Limousin, and in the Velay and Montagne Noire domes. The nappe pile is dissected by strike-slip shear zones and low-angle detachments. It is intruded by plutonic rocks with ages ranging from late Devonian to late Carboniferous-Permian. These crystalline rocks are capped by volcanic and sedimentary rocks deposited in intramontane and foreland basins with ages also ranging from late Devonian to Permian. The detailed structural relationships between these different lithological units and the timing of geological events is further presented in the next section.

2.3 Architecture and P-T-t record of nappes of the western and eastern French Massif Central

As stated above, the identification of the LAC as a suture led to the revival, after the pioneer proposition of Demay (1948), of the nappe concept in the French Massif Central and several models have been proposed ever since. The early tectonic model (Burg and Matte, 1978; Matte, 1986) highlights three distinct nappes cored by the migmatites of the UGU; namely from north to south, the Sioule, the Haut Allier – Marvejols, and the Rouergue nappes (Fig. 7). In this model each outcropping zone of the UGU and associated LAC corresponds to a locally rooted nappe that, in turn delineates a suture and thus a former oceanic basin. In contrast, subsequent models (Faure et al., 2009a; Matte, 1991, 2001; Lardeaux, 2014) propose that these three nappes form only one with the LAC representing a single ocean rooted beneath the Paris Basin. Consequently, the discontinuous outcrops of the UGU with enclaves of LAC are interpreted as klippes.

The high-grade nappes are characterized by a penetrative composite foliation resulting from superimposed structures and metamorphic parageneses that is typically parallel to the tectonic contacts. This foliation is dominantly shallow dipping (Burg and Matte, 1978; Matte, 1986; Faure et al., 2009a) but is locally steeply dipping such as in the Livradois or in the Monts du Lyonnais (Feybesse et al., 1988; Lardeaux and Dufour, 1987; Gardien et al., 1990, 2011). Away from strike slip shear zones, the lineation associated with this composite foliation is dominantly E-W to WNW-ESE trending in the western part of the French Massif Central and in the Sioule region but is dominantly N-S to NNE-SSW trending in the eastern part (Faure et al., 2009a). Moreover, this foliation is in place affected by regional upright folding such as in the Limousin (Burg and Matte, 1978; Girardeau et al., 1986; Matte, 1986). The structure of the nappe pile is blurred by numerous granitic plutons and original contacts are reworked by thrusts, strike-slip shear zones, low-angle detachments and high-angle normal faults. This complex structural record and the scarcity of outcrops impedes the identification of the original tectonic contacts between the nappes in most places. Nevertheless, in order to discuss the tectonic evolution of potentially distinct nappes, in the following sections, we review available data regarding the structural position and P-T-t record of the lithological-tectonic units presented above distinguishing the Western and Eastern parts of the French Massif Central separated by the Sillon Houiller Fault (Figs. 810).

thumbnail Fig. 7

Cross sections depicting the previously proposed nappe structure. A. Multiple sutures-nappe model (modified after Burg and Matte, 1978). B. Single suture model associated with a basement duplex structure (Matte, 1991). C. Single suture model associated with a stack of three nappes and an “unknown Proterozoic basement” (Faure et al., 2009a, 2009b).

thumbnail Fig. 8

Cross sections of the French Massif Central. The location of the cross sections (A-A’, B-B’, C-C’, and D-D’) is indicated on the geological map Figure 2. Same legend as Figure 3 with the addition of a granulitic lower crust intruded by mantle-derived mafic magmas. The upper part of the sections is constrained by field observations. The shaded lower part of the sections is less constrained and is based on scarce geophysical data that allow the prolongation of some structures at depth and on exposed sections of the Variscan lower crust in the Southern Alps (Ivrea Zone) and in Calabria.

thumbnail Fig. 9

P-T-t constraints on the metamorphic history of the Upper Gneiss Unit and of the Middle Allochthonous Unit. The Leptynite-Amphibolite Complex is characterized by HT eclogite facies metamorphism retrogressed into granulite facies by isothermal decompression and then into amphibolite facies by a decrease in temperature. Granulite facies to amphibolite facies metamorphism is also recorded in the migmatitic gneisses hosting the LAC. The PT path is depicted by the white arrows and numbers correspond to radiometric ages.

thumbnail Fig. 10

P-T-t constraints on the metamorphic history of the Lower Gneiss Unit. The P-T-t paths of orthogneisses and paragneisses of the LGU indicate first an increase in temperature followed by isobaric cooling. This is particularly well illustrated by the P-T-t path of the south Velay which is characterized by a HT/LP gradient. The PT path is depicted by the white arrows and numbers correspond to radiometric ages.

2.3.1 Western French Massif Central (W-FMC)

In the northern part of the W-FMC, the Aigurande region is exposing metamorphic rocks unconformably overlain by Mesozoic sediments of the Paris Basin to the north and is delimited by the La Marche shear zone to the south (Fig. 8 cross section AA’). Metamorphic rocks present a polyphased structural and metamorphic history associated with an inverted metamorphic gradient characterized by the superposition, from top to bottom, of the UGU, the LGU and PAU (Quenardel and Rolin, 1984; Faure et al., 1990). At the top of the nappe pile, the UGU is dominated by diatexites and metatexitic orthogneisses and paragneisses with rare quartzites. These rocks display a dominant garnet-sillimanite-cordierite mineral paragenesis with relics of kyanite that attest for HP partial melting followed by retrogression during decompression. The amphibolite facies foliation of these migmatites is associated with top-to-the SE sense of shear criteria considered to record nappe emplacement (Faure et al., 1990). The boundary between the UGU and the LGU is marked by boudins of eclogitic amphibolites and of ultramafics attributed to the LAC. Amphibolites that are part of the LAC, yield 40Ar/39Ar dates on amphibole at 389 ± 8 Ma interpreted as the age of amphibolite facies metamorphism (Boutin and Montigny, 1993). A mylonitic shear zone underlines the contact between the migmatitic units of the LGU and micaschists attributed to the PAU (Faure et al., 1990). This contact was first interpreted as a thrust responsible for burial of the PAU beneath the LGU (Quenardel and Rolin, 1984). In contrast, retrogression of the garnet-biotite dominant foliation of the micaschists associated with top-to-the NE kinematic criteria has been reinterpreted as reflecting exhumation of the PAU during a period of regional extension estimated at 325–312 Ma based on syntectonic leucogranites emplacement (Faure et al., 1990). It is noteworthy that these micaschists of the PAU contain relics of kyanite and staurolite attesting for a so far undated MP/MT amphibolite facies metamorphism before retrogression into greenschist facies.

To the north of the FMC, the crystalline basement beneath the Paris Basin has been sampled in the Couy deep borehole down to 3500 m (Fig. 1). Granulite facies migmatites of the UGU are associated with metabasites attributed to the LAC with a Cambrian-Ordovician protolith as constrained by a Sm-Nd isochron of 494 ± 17 Ma and a U-Pb zircon date of 497 ± 13 Ma (Pagel et al., 1992). These rocks display a NE-SW trending foliation steeply dipping to the SE. The granulite facies mineral paragenesis yield a pressure ranging from 15 to 9 kbar and a temperature from 900 to 650 °C. Retrogression into the amphibolite facies is recorded at c. 6 kbar for c. 600 °C (Burg et al., 1989a, 1989b; Ballèvre and Balé, 1992). Amphibolites yield an Rb-Sr isochron of 387 ± 2 Ma interpreted as dating high-grade metamorphism. 40Ar/39Ar ages on amphibole and biotite from granulite facies amphibolites and orthogneisses, respectively range from c. 385 Ma to c. 379 Ma and point to rapid cooling below 300 °C before the end of the Devonian (Costa and Maluski, 1988). Accordingly, the P-T-t record of the UGU sampled in the Couy deep borehole is similar to the one of the UGU exposed in the Aigurande region and these rocks are attributed to the same nappe that will be designated as the northern nappe in the following.

South of the Aigurande region, the northern Limousin region exposes a slightly different nappe package. The UGU is mainly made of migmatitic paragneisses with ubiquitous relics of HP metamorphism expressed as eclogites and numerous garnet amphibolite boudins (Le Breton et al., 1986). Below, the LGU exposed in the core of the Thaurion and Meuzac antiforms (Fig. 8, section AA’) is dominated by amphibolite facies orthogneisses with some paragneisses. In addition to the high-grade units recognized in the Aigurande region, several authors have identified a Middle Allochthonous Unit (MAU) stacked in between the UGU and LGU (Girardeau et al., 1986; Dubuisson et al., 1989; Berger et al., 2010a, 2010b). The UGU contains eclogite facies metabasites attributed to the LAC yielding U-Pb zircon ages pointing to a crystallization of their protolith between c. 489 and c. 475 Ma. The eclogitic UHP event is dated at 412 ± 10 Ma followed by a resetting potentially linked to partial melting and retrogression into the granulite facies at 382 ± 7 Ma (Berger et al., 2010a), which is also consistent with previous whole rock Rb-Sr isochrons from c. 385 to c. 375 Ma (Duthou, 1978; Duthou et al., 1994), U-Pb on zircon at 383 ± 5 Ma (Lafon, 1986) and U-Th-Pb on monazite between c. 378 and c. 374 Ma (Faure et al., 2008). The relationship between this nappe pile and the one exposed in the Aigurande region is not clearly identified but the differences displayed by their cooling histories point to distinct exhumation histories. In the following, the nappe exposed in the northern Limousin will be designated as the Central nappe.

In the southern Limousin the nappe pile is deformed in a serie of upright folds designated as the Uzerche synform and the Tulle antiform (Ledru et al., 1989) (Fig. 8, section AA’). The UGU is dominated by granulite facies migmatitic paragneisses with eclogitic mafic boudins yielding peak metamorphic conditions at c. 15 kbar and c. 750 °C (Santallier, 1981; Bellot and Roig, 2007). In a lower structural position, migmatitic paragneiss and orthogneiss that contain boudins/enclaves of eclogites and garnet-spinel peridotites with peak P-T conditions at c. 15 kbar and c. 700 °C similar to the ones of the UGU (Santallier, 1981; Ledru et al., 1989; Bellot and Roig, 2007). These eclogitic boudins have been first attributed to the LGU but could be part of the MAU according to their structural position and P-T record. Structurally below these rocks, migmatitic paragneisses display a syn-migmatitic foliation underlined by cordierite-sillimanite-bearing leucosomes that document a retrograde P-T path from c. 9.5 kbar at c. 850 °C to c. 6 kbar at c. 600 °C and locally to c. 3.5 kbar at c. 550 °C (Bellot and Roig, 2007). These migmatites yield a U-Pb age on zircon at 382 ± 5 Ma interpreted as the age of partial melting (Lafon, 1986) and a variety of 40Ar/39Ar ages on micas as well as U-Th-Pb ages on monazite spreading from c. 350 Ma to c. 315 Ma that might reflect progressive exhumation and cooling, or partial reseting owing to recrystallization, during the Carboniferous (Costa, 1992; Gebelin et al., 2004; Melleton et al., 2009). The high-grade nappes are overlain by the Thiviers-Payzac Unit and by the Génis Unit (Guillot and Doubinger, 1971; Guillot and Lefevre, 1975; Santallier, 1981; Bellot and Roig, 2007) that are affected by a prograde Barrovian metamorphic gradient with peak conditions at c. 9 kbar for c. 750 °C at the lowest structural level. The Thiviers-Payzac Unit is intruded by pre-tectonic calc-alkaline dolerite dykes dated at 363 ± 10 Ma by K-Ar on whole rock (Bellot and Roig, 2007). The presence of high-grade rocks of the UGU sandwiched in between lower-grade rocks is interpreted to correspond to vertical extrusion of the high-pressure UGU into an orogenic wedge affected by Barrovian metamorphism (Bellot and Roig, 2007), as envisioned in models of Chemenda et al. (1996) or Escher and Beaumont (1997). According to these data, the South Limousin UGU is characterized by a younger cooling history than the North Limousin one, which is used to define what will be referred as the southern nappe in the following.

The oldest Variscan plutonic rocks identified in the Western part of the French Massif Central are late Devonian and display a variety of petrological and geochemical signatures. Small plutons emplaced into the UGU (ACG type of Barbarin 1999), with a composition ranging from gabbro to granodiorite and signatures from calc-alkaline to tholeiitic, define a broadly linear trend referred to as the “Limousin tonalite line” (Peiffer, 1986; Cuney et al., 1990, 1993). These rocks first yielded TIMS U-Pb zircon ages of 379 ± 19 Ma and 355 ± 2 Ma (Bernard-Griffiths et al., 1985) and have since provided more precise ages of 365 ± 3 Ma and 360 ± 1 Ma by the same method (Pin and Paquette, 2002). In contrast, the contemporaneous Guéret pluton (Turpin et al., 1990) is a cordierite-bearing peraluminous granite (CPG in the nomenclature of Barbarin, 1999) dated in the late Devonian (Berthier et al., 1979). It forms a c. 1 km thick laccolith overlying the cordierite-bearing migmatitic gneisses of the UGU (Lameyre et al., 1988; Dupis et al., 1990; Gébelin et al., 2006). These plutonic rocks are cross-cut by a system of E-W to NW-SE trending dextral shear zones including the La Marche, Courtine and Pradines shear zones that started their activity at about 350 Ma as attested by 40Ar/39Ar on syntectonic biotite in the Aigurande plateau (Gébelin et al., 2007). These shear zones control the syntectonic emplacement of CPG and MPG plutons from 345 to 310 Ma (Lafon and Respaut, 1988; Lerouge and Quenardel, 1988; Alexandrov et al., 2001; Gébelin et al., 2007, 2009; Rolin et al., 2009, 2014), bounded at their roof by detachment zones (Faure and Pons, 1991; Gébelin et al., 2007).

Leucogranites are interpreted to be generated by partial melting of a metasedimentary middle crust (Cuney et al., 1990; Williamson et al., 1996; Moyen et al., 2017; Villaros et al., 2018). This is consistent with their structural position relative to their host rocks, as exemplified by the Millevaches laccolith rooted into cordierite-bearing migmatitic paragneisses that have been affected by partial melting under MP/MT amphibolite facies conditions (6 kbar for 850 °C) as attested by leucosomes localized in strike-slip shear bands (Gébelin, 2004; Gébelin et al., 2006, 2009). This high-grade metamorphism and partial melting has been dated at 315 ± 4 Ma and 316 ± 2 Ma, coeval with the syntectonic emplacement of leucogranites in the Pradines dextral strike-slip shear zone at 313 ± 4 Ma (Gébelin et al., 2009). The dextral strike-slip shear zones are cross-cut by high-angle normal faults and low-angle detachments such as the NE-SW trending Nantiat and Bussière shear zones and the N-S trending Argentat shear zone (Fig. 3) (Gébelin et al., 2007, 2009). The footwall of these detachment zones represented sites of strong (meteoric) fluid-rock-deformation interactions during the late Carboniferous (Dusséaux, 2019; Dusséaux et al., 2019). Exhumation along these low-angle mylonite zones is constrainted to be middle (Visean) to late Carboniferous (Stephanian) in age, by argon thermochronology (Alexandrov et al., 2000; Roig et al., 2002; Gébelin, 2004; Gébelin et al., 2007; Rolin et al., 2014). The activity of these shallow dipping detachments and strike-slip shear zones, which are part of the Sillon Houiller Fault, controlled the deposition of late Carboniferous coal-bearing sediments in small extensional and pull-apart basins (Feybesse, 1981; Thiéry et al., 2009). The onset of deposition in these basins during the Visean is confirmed by a fireclay dated at 332 ± 4 Ma by U-Pb TIMS on zircon in northern Limousin (Bruguier et al., 1998). The Permian clastic sedimentary deposits of the Brive basin unconformably overlie these late Carboniferous deposits and mark the end of the Variscan tectonic activity in this region.

Although scarce, available geophysical data for the W-FMC provide some constraints on the deep crustal structure. Gravity data indicate that most plutons are laccoliths with an average thickness of about 1 km but that can locally reach up to 3 km (Gébelin et al., 2006; Joly et al., 2008, 2009). As a complement, seismic data allow to identifiy the prolongation of surface structures at depth (Bitri et al., 1999). In the upper crust, the main feature is that most reflectors appear to match the projection of the dominantly shallow-dipping high-grade fabric parallel to lithological-tectonic contacts identified at the surface. These reflectors are only crosscut and offset by high-angle faults and low-angle detachments. For example, high-angle faults separating the North and South Limousin and affecting high-grade rocks coring the Meuzac antiform, are offsetting reflectors marking the contacts between the UGU, MAU, LGU and PAU and are rooted into a high reflectivity zone at 10 km depth. Similarly, the Argentat shear zone corresponds to a several km thick zone of reflectors shallowly dipping to the West rooting into a reflective zone at about 10 km depth. Accordingly, the Argentat shear zone might be interpreted as the breakaway zone of a low-angle detachment rooted in the brittle-ductile transition. These high reflectivity zones are overlying a c. 10 km thick seismically transparent middle crustal zone with a relative low density that has been interpreted to be composed of granitic material (Bitri et al., 1999). Alternatively, this middle crust, coring low amplitude (few km) and long wavelength (tens of km) dome-shaped structures, could be composed by migmatites, as it has been proposed for similar structures detected beneath low-angle detachments in the South Armorican Massif (Bitri et al., 2010). Beneath this transparent middle crustal zone, the lower crust is typically marked by its high reflectivity from c. 20 km down to the Moho at c. 30 km depth (Bitri et al., 1999). These characteristics are used to constrain the deep part of the cross sections (Fig. 8) beyond information provided by surface outcrops.

2.3.2 Eastern French Massif Central (E-FMC)

The continuity of the nappes from the W-FMC to the E-FMC across the Sillon Houiller Fault is not easily established. Nevertheless, as for the W-FMC, the P-T-t record of high-grade rocks of the E-FMC points to the diachronous emplacement of several distinct nappes.

In the northern part of the E-FMC, the position of the late Devonian Brévenne back-arc, to the south of the Morvan arc represented by the late Devonian Somme Unit, has been used to infer a southward subduction of the Rheic Ocean along a continental active margin (Faure et al., 1997, 2008). Sparse outcrops of retrogressed eclogites, serpentinised peridotites and amphibolites affected by HP metamorphism are characteristic of the UGU (Gardien et al., 1988; Godard, 1990). The timing of exhumation of these high-grade rocks and their structural relationships with the Devonian volcanic-sedimentary sequences are ill-defined. However, it has been proposed that they were already exhumed before the middle Devonian at the time of formation of the Morvan arc and Brévenne back-arc (Faure et al., 1997, 2008). In that case, the UGU exposed in the Morvan region might correspond to the northern nappe described in the W-FMC. Because of these uncertainties, this part of the French Massif Central is not represented on cross sections of Figure 8.

In the Sioule area (Fig. 8, section BB’), the Variscan basement shows a dominant shallow-dipping composite foliation deformed in a broad antiform cored by granitic plutons and displaying an inverted metamorphic gradient (Faure et al., 1993, 2002). The top of the nappe pile is made of cordierite-bearing diatexites and migmatitic orthogneisses and paragneisses that have recorded isothermal decompression from 12–13 kbar to 2–3 kbar at 650–700 °C. Both lithologies display a composite foliation bearing a NE-SW trending lineation (Audren et al., 1987; Schulz et al., 2001; Schulz, 2009). Serpentinite boudins and granulitic relics allow to assign these rocks to the UGU (Ravier and Chenevoy, 1979). Metamorphic monazite with a U-Th-Pb EPMA age at 416 ± 15 Ma is attributed to HP metamorphism (Do Couto et al., 2016) but retrogression under amphibolite facies has not been dated. The contact with the lower-grade underlying micaschists attributed to the PAU has first been interpreted as a thrust (Burg and Matte, 1978; Ledru et al., 1989) but has been then attributed to an extensional detachment reflecting exhumation of the PAU during Visean regional extension dated at 337–336 Ma by Ar thermochronology on mica and amphibole of syntectonic granites (Faure et al., 1993, 2002). Despite uncertainties on the timing of exhumation of the UGU in the Morvan and Sioule regions, we propose to consider that they are part of the same northern nappe exhumed, at least partly, during the late Devonian as previously proposed by Faure et al. (Faure et al., 1997, 2008). The southern boundary of the Sioule-Morvan high-grade nappe is marked by a dextral strike-slip corridor (Hermitage shear zone) crossing the Forez and Brévenne regions and localizing the emplacement of syn-tectonic MPG plutons dated from 332 ± 2 to 326 ± 2 Ma (U-Pb on zircon, Laurent et al., 2017). These high-grade metamorphic and plutonic rocks are capped by Visean undeformed volcanics, volcaniclastics and granophyres represented by the “tufs anthracifères” series dated at 336 ± 5 Ma by a Rb-Sr isochron (Carrat and Zimmermann, 1984; Sider and Ohnenstetter, 1986; Delfour et al., 1989; Leloix et al., 1999; Faure et al., 2002). A similar age of 332 ± 2 Ma has been obtained by U-Pb TIMS on zircon from a rhyolite sampled in the Decazeville basin (Bruguier et al., 1998). These Visean volcanic rocks are locally associated with marine deposits, which has been used to propose that the high-grade rocks of the nappe pile were exhumed but remained below sea-level at this time (Franke, 2014).

South of the Hermitage shear zone, several fragments of the UGU, composed of migmatitic paragneisses and orthogneisses grading from metatexites to diatexites, are exposed in the Combrailles, Cézallier, Artense, Livradois, Truyère, Haut-Allier, Lyonnais and Vivarais regions (Feybesse et al., 1988; Lardeaux and Dufour, 1987; Gardien, 1993; Gardien and Lardeaux, 1991; Mercier et al., 1991a, 1991b; Gardien et al., 2011) (Fig. 8, sections BB’, CC’, DD’). We propose here that these litho-tectonic units are part of the central nappe. The migmatitic gneisses of the UGU yield whole rock Rb-Sr isochrons ranging from 385 to 375 Ma in the Lyonnais (Duthou et al., 1981, 1994) and a monazite age of 360 ± 4 Ma by U-Th-Pb EMPA in the Livradois (Gardien et al., 2011; Vanderhaeghe et al., 2013). These ages are tentatively interpreted to record the transition from an early stage of HP granulite facies (at least 10 kbar for c. 800 °C) followed by decompression to 6 kbar. In the northern root zone of this nappe, exposed in the Lyonnais and Livradois, the NE-SW trending regional foliation of the UGU migmatites is steeply dipping to the north (Lardeaux and Dufour, 1987; Feybesse et al., 1988; Gardien et al., 1990, 2011) (Fig. 8, sections CC’, DD’). North of the Livradois, these migmatites display a penetrative C/S fabric consistent with a top to the south sense of shear (Koné, 1985; Gardien et al., 2011; Vanderhaeghe et al., 2013). To the north of the Lyonnais, the late Devonian volcanic-sedimentary series of the Brévenne unit are characterized by upright folds with an NNE-SSW trending axial planar schistosity under greenschist facies metamorphism (Fig. 8 section DD’). The contact between the Brévenne Unit and the UGU of the Lyonnais is marked by transposition of previous structures into NE-SW trending shallow-dipping shear zone (Feybesse et al., 1988), also delineated by a syntectonic granite displaying a C/S fabric consistent with a top to the NW sense of shear (Feybesse et al., 1988; Leloix et al., 1999). Deformation of this granite is dated at 339–337 Ma by 40Ar/39Ar on recrystallized muscovite (Faure et al., 2002), which is consistent with both Rb-Sr whole rock isochrons obtained by Gay et al. (1981) on the syntectonic granite and by Vialette (1973) on genetically linked hypovolcanics (Leloix et al., 1999; Faure et al., 2002), giving ages at 339 ± 8 and 336 ± 5 Ma respectively. These data suggest that at least part of the exhumation of the UGU in this region occurred during the early Carboniferous and was associated with top to the NW shearing. In the Livradois (Fig. 8, section CC’), the southern part of the migmatitic UGU nappe is cut by dextral-reverse transcurrent shear zones and associated syntectonic peraluminous granodiorite and leucogranite plutons (CPG and MPG) dated by TIMS U-Pb on zircon at 315 ± 4 and 311 ± 18 Ma, respectively (Solgadi et al., 2007; Gardien et al., 2011; Vanderhaeghe et al., 2013), which provides a maximum age for the exhumation of these rocks. The geochronological data obtained on the HP migmatites of the UGU in the Lyonnais and Livradois are significantly younger than the ones obtained in the Sioule area and we propose to attribute these rocks to the central nappe.

In the Haut-Allier area (Fig. 8, section CC’), the foliation of the UGU flattens to the south and delineates a dome cored by migmatites attributed to the LGU (Burg and Matte, 1978; Gardien et al., 2011). In the Artense, Truyere, and Marvejols regions, amphibolite facies paragneisses attributed to the LGU and greenschist facies micaschists attributed to the PAU define an inverted metamorphic gradient with respect to the overlying UGU (Ledru et al., 1989; Mercier et al., 1991a). Amphibolites of the UGU contain thin tonalitic to trondhjemitic layers interpreted to reflect high-pressure partial melting (Nicollet and Leyreloup, 1978; Pin and Lancelot, 1982). This high-pressure metamorphism has been dated by zircon TIMS U-Pb analyses at 432 + 20/− 10 Ma in the Haut Allier (Ducrot et al., 1983), 415 ± 6 Ma in Marvejols (Pin and Lancelot, 1982) and at 413 ± 23 and 408 ± 7 by Pb-Pb on zircon and a Sm-Nd whole rock and garnet isochron respectively in the Rouergue (Paquette et al., 1995). Retrograde amphibolite facies metamorphism is recorded by lower intercepts defined by discordant U-Pb zircon data at 345–340 Ma obtained on paragneisses and amphibolites (Pin and Lancelot, 1982). The UGU-LGU contact is intruded by the Margeride laccolithic composite pluton (Couturié et al., 1979; Couturié and Caen-Vachette, 1980; Talbot et al., 2005). The main porphyritic monzogranite yields U-Pb zircon dates of 334 ± 9 Ma (Respaut, 1984) and of 313 ± 2 Ma (Laurent et al., 2017) and a U-Pb on monazite of 314 ± 3 Ma (Pin, 1979). The leucogranitic facies cross cutting the monzogranite is dated by various methods from 311 to 298 (Couturié and Caen-Vachette, 1980; Lafon and Respaut, 1988; Monié et al., 2000; Laurent et al., 2017). In the Vivarais (Fig. 8, section DD’), several klippes of UGU made of migmatites including amphibolite enclaves overly paragneisses and orthogneisses of the LGU, outlining an inverted metamorphic gradient (Gardien, 1993). Rocks of the UGU have recorded HT eclogitic to HP granulitic metamorphic conditions up to 15 kbar at c. 800 °C, followed by retrogression at 5 kbar at c. 550 °C (Gardien, 1993). The HP migmatitic gneisses and amphibolites of the UGU contain zircon with metamorphic rims dated by U-Pb at 351.5 ± 3.0 Ma and 343.5 ± 2.6 Ma, respectively and interpreted to represent retrogression of the eclogite into amphibolite facies after tectonic accretion during continental collision (Chelle-Michou et al., 2017). The underlying paragneisses display low-pressure granulite facies metamorphism associated with widespread partial melting dated at ca. 308 Ma (Chelle-Michou et al., 2017). These zones in the Vivarais, Marvejols and Rouergue, with an inverted metamorphic gradient might correspond to the southern tip of the UGU central nappe.

To the south of the Rouergue, the Najac eclogites have recorded a pressure of 15–20 kbar and a temperature of 560–630 °C (Lotout et al., 2018). Eclogite-facies metamorphism is dated by U-Pb on zircon at 385.5 ± 2.3 Ma, by Lu-Hf on garnet cores at 382.8 ± 1.0 Ma, and by a Sm-Nd amphibole-garnet-whole rock isochron of 376.7 ± 3.3 Ma. Subsequent exhumation and cooling below c. 500 °C are constrained by an apatite U-Pb age at 369 ± 13 Ma. The maximum temperature reached by these eclogites is significantly lower than the ones recorded by eclogites enclosed in the granulite facies migmatites of the UGU. Moreoever, geochronological data available for the Najac clogites point to burial coeval with the UGU central nappe exposed in the Haut-Allier, Marvejols and Rouergue, but also for a much older exhumation. These data suggest that the Najac eclogites are part of a nappe equivalent to the MAU, which has a different P-T-t history than the UGU.

The Velay and Montagne Noire migmatite domes correspond to the lowermost exposed structural level of the E-FMC (Figs. 4 and 8 section CC’, DD’). These migmatites originated by partial melting of paragneisses and orthogneisses under MP to LP conditions at the bottom of the nappe pile from c. 320 to c. 300 Ma (Montel et al., 1992; Williamson et al., 1996; Downes et al., 1997; Barbey et al., 1999, 2015; Ledru et al., 2001; Villaros et al., 2018) and contain dismembered enclaves of vaugnerites (Michon, 1987; Sabatier, 1991; Ledru et al., 2001; Couzinié et al., 2014). The 100 km-wide, first-order Velay dome is delineated by the foliation of the surrounding paragneisses and orthogneisses (Lagarde et al., 1994). The synmigmatitic foliation defines subdomes with an average diameter of 10–20 km and characterized by a radial distribution of the HT mineral lineation (Ledru et al., 2001). Accordingly, the Velay dome is a crustal-scale structure interpreted to represent the exhumed middle part of the orogenic crust. The western boundary of the dome is subvertical and its southern part is overturned toward the south (Burg and Vanderhaeghe, 1993; Vanderhaeghe et al., 1999; Ledru et al., 2001). Along the eastern side of the dome, in the Vivarais region, the symigmatitic foliation is shallow-dipping and small klippen of the HP migmatites belonging to the UGU central nappe, as described above, are exposed above LGU grading into MP to LP migmatites (Gay et al., 1982; Gardien, 1993; Gardien and Lardeaux, 1991; Chelle-Michou et al., 2017). The superposition of the HP migmatites attributed to the UGU on top of the LGU is inferred to represent an inverted metamorphic gradient preserved from the time of the nappe emplacement. In turn, the downsection transition from mid-amphibolite facies metamorphism documented in the LGU just beneath the UGU to MP and LP migmatites of the Velay dome corresponds to a normal metamorphic gradient marked by an important increase in temperature with depth.

The northern part of the Velay dome is delimited by the Pilat mylonitic low-angle detachment intersected by high angle cataclastic normal faults (Malavieille et al., 1990; Gardien et al., 1997). Heterogeneous granite (i.e. diatexites) in the core of the dome shows a relatively tight cluster of U-Pb ages on zircon and monazite, in the range of 307 ± 2 to 301 ± 5 Ma (Mougeot et al., 1997; Couzinié et al., 2014; Chelle-Michou et al., 2017; Laurent et al., 2017). Magmatic zircon grains and rims with ages ranging from 340 to 300 Ma (Chelle-Michou et al., 2017; Laurent et al., 2017) point to a protracted history of zircon growth and/or dissolution-reprecipitation in these diatexites. To the north of the Velay the Gouffre d’Enfer syntectonic granite, emplaced within the Pilat detachment, yielded an unprecise whole rock Rb-Sr age of 322 ± 26 Ma (Caen-Vachette et al., 1984; Vitel, 1988). To the south of the Velay dome, a Rb-Sr isochron of 302 ± 4 Ma has been obtained for the Rocles syntectonic granite (Caen-Vachette et al., 1981) emplaced in a detachment that has been overturned later on (Burg and Vanderhaeghe, 1993; Vanderhaeghe et al., 1999; Bouilhol et al., 2006). U-Th-Pb ages by EMPA on monazite of 324 ± 4 Ma and 325 ± 5 Ma have also been obtained for the the Rocles granite (Bé Mézémé, 2005; Bé Mézémé et al., 2006), which overlap with a more reliable U-Pb zircon LA-ICPMS age of 320.3 ± 3.8 Ma (Couzinié, 2017). This suggests that while the Rocles granite was emplaced at c. 320 Ma, the Rb-Sr system was probably reset at 302 Ma. Migmatites and heterogeneous granite of the core of the dome are intruded by small plutons and dykes of CPG/MPG with sharp contacts designated as late-migmatitic (Montel and Abdelghaffar, 1993) dated by U-Th-Pb on monazite (LA-ICPMS) from 307 ± 2 to 297 ± 4 Ma (Didier et al., 2013). They are characterized by an aluminous and potassic signature and are rich in rounded enclaves of metapelites and of microdiorite suggesting a metapelitic source together with a contribution from a magma generated by partial melting of an enriched mantle (Villaros et al., 2018). Microgranites with similar chemical signatures have yielded identical U-Th-Pb EMPA ages on monazite at about 300 Ma (306 ± 12, 291 ± 9 Ma) but also at about 255 Ma (257 ± 8 and 252 ± 11 Ma), which points to either a Permian crystallisation or to hydrothermal perturbation at this time (Montel et al., 2002). High-angle normal faults rooting in the low-angle detachments, delimit the late Carboniferous St. Etienne, Jaujac and Alès extensional basins filled by coarse detrital sediments including clasts from the migmatites and granites interbedded with volcanic deposits (rhyolitic ash fall tuffs and accretionary lapilli) and coal-bearing sediments (Becq-Giraudon et al., 1996).

The Montagne Noire is characterized by a dome structure with a core of migmatitic gneisses and a mantle of low-grade metasedimentary rocks (Gèze, 1949; Arthaud et al., 1966; Ellenberger, 1967; Engel et al., 1978, 1981; Demange, 1980; Demange and Jamet, 1985). Migmatites and orthogneisses in the core of the dome display a prolate finite strain ellipsoid indicative of constriction with a subhorizontal long axis parallel to the axis of the elliptical shape dome (Echtler and Malavieille, 1990; Matte et al., 1998; Rabin et al., 2015). Migmatites coring the Caroux dome in the Montagne Noire have recorded P-T conditions of 6 kbar and 720 °C (Rabin et al., 2015) and they are juxtaposed to low-grade LP/HT micaschists (Thompson and Bard, 1982; Demange and Jamet, 1985). Magmatic zircon from the migmatites yield U-Pb ages spreading from c. 330 to c. 300 Ma (Franke et al., 2011; Faure et al., 2014; Roger et al., 2015). Migmatites contain mafic enclaves that have preserved relicts of eclogite facies metamorphism (Demange, 1980; Faure et al., 2014). Zircon grains from the eclogite facies rocks yield two age peaks, one at c. 360 Ma and the other at c. 315 Ma, with a few discordant ages pointing at an early Paleozoic heritage (Faure et al., 2014; Whitney et al., 2015). Faure et al. (2014) attribute the 360 Ma age to eclogite facies metamorphism whereas Whitney et al. (2015), based on REE signatures of the different zircon zones, argue that eclogite facies is recorded by the 315 Ma age. The Montalet syntectonic granite emplaced along the northern side of the Montagne Noire yield a U-Pb age on zircon of 294 Ma (Poilvet et al., 2011). This syntectonic granite is juxtaposed to the late Carboniferous Graissessac basin along a mylonitic to cataclastic detachment (Van Den Driessche and Brun, 1992). Many different models have been proposed for the Montagne Noire, the dome structure having been interpreted as (i) an anticline developed by horizontal shortening (Matte et al., 1998) (Fig. 7), (ii) a ductile layer exhumed in a pull-apart basin (Nicolas et al., 1977; Rey et al., 2011; Whitney et al., 2015), (iii) a diapir (Charles et al., 2009), (iv) an eroded antiformal stack (Malavieille, 2010) or (v) a Metamorphic Core Complex formed as a consequence of gravitational collapse of the Variscan belt during convergence (Echtler, 1990; Echtler and Malavieille, 1990; Aerden, 1998; Aerden and Malavieille, 1999) or during regional extension (Van Den Driessche and Brun, 1992).

3 Previous tectonic-geodynamic reconstructions and debated issues

In this section, we present the various tectonic-geodynamic reconstructions that were proposed at the scale of the Variscan belt with a focus on the French Massif Central, with the aim of identifying robust features but also to point at the discrepancies and shortcomings.

3.1 Monocyclic doubly-vergent orogen model

Early geodynamic reconstructions at the scale of Western Europe invoked a doubly-plunging subduction system based on the presence of the two major suture zones described above, namely the Rheic suture and the Medio-European suture, and opposite vergence of structures on the northern and southern sides of the orogenic belt (Matte, 1986, 1991, 2001) (Fig. 2A). According to the latest version of this model, Silurian subduction of Ordovician oceanic basins along the southern branch of the Variscan Belt was followed by Devonian continental subduction and Carboniferous collision between Laurussia and Gondwana (Matte, 2001). This model, mostly relying on a synthesis of structural, metamorphic and sedimentary data, provides a first order tectonic-geodynamic framework for the Variscan belt. In contrast, other geodynamic models assuming that Armorica remained attached to the north Gondwana margin do not include the closure of a Medio-European Ocean (e.g. Nance et al., 2010) and thus fail to account for the presence of ophiolitic assemblages to the south of Armorica. Still, the reconstructions proposed by Matte (1986, 1991, 2001) elude some key features including (i) the presence of a late Devonian arc and back-arc association in the Vosges (Skrzypek et al., 2012) and the Eastern Massif Central (Somme and Brévenne Units, see above); (ii) the mechanisms of exhumation of eclogites and migmatites, and, most of all, (iii) syn-orogenic partial melting and magmatism that are notably absent from the cross sections despite the vast exposed surface represented by migmatites and granites. Note also that in the scenario proposed by Matte (1986, 1991, 2001), the HP rocks of the UGU are issued from the southern margin of Armorica, i.e. the upper plate relative to Silurian-Devonian subduction. Accordingly, from their position, it is unclear how these rocks could have been buried to granulite facies and then exhumed back to the surface.

3.2 Polycyclic orogenic model

Faure et al. (1997) refined the geodynamic model elaborated by Matte (1986, 1991, 2001) in order to account for the late-Devonian magmatic arc and back-arc inferring a southward subduction of the Rheic Ocean based on the relative position of the Somme Unit to the north of the Brévenne Unit (Fig. 3). These authors propose a polcyclic model for the tectonic evolution of the Variscan belt in Western Europe with two distinct orogenic phases (Fig. 2B). The first orogenic phase is associated with northward oceanic subduction during the Silurian and Devonian leading to the closure of the Medio-European Ocean followed by continental subduction during the early Devonian and exhumation of the UGU before the late Devonian. This model provides an explanation for high-pressure metamorphism, which is attributed to burial of rocks from the Gondwana margin that first have been dragged with the downgoing plate and then decoupled from it allowing their exhumation. It should be noted that mafic and ultramafic rocks of the LAC are described by Faure et al. (2008) as part of the UGU but are not considered to represent an ophiolite, which is at odds with the interpretations of most authors. For example, to explain the intimate association of mafic and felsic HP rocks, Lardeaux (2014) proposes that the UGU and the LAC represent the formerly thinned continental margin of Gondwana (Fig. 2B). Irrespective of these differences, all these authors invoke a second orogenic phase associated with southward subduction of the Rheic Ocean during the late Devonian beneath an upper plate, made of the exhumed, HP UGU nappe. The latter undergoes extension as attested by the late Devonian Brévenne back-arc rift basin (Fig. 2B). In this model, partial melting of the UGU is attributed to retrogression owing to decompression during syn-orogenic exhumation whereas partial melting of the underthrusted LGU is interpreted to be caused by burial beneath the UGU (Fig. 2B). At the lithospheric scale, a high-velocity anomaly beneath the Paris Basin detected in tomographic models is proposed to represent a remnant of this southward subducting slab (Averbuch and Piromallo, 2012). As an alternative, Lardeaux et al. (2001) advocated that the late Devonian Brévenne back-arc basin was opened above a northward subducting slab (corresponding to the former Medio European Ocean) at the time of exhumation of the UGU and was closed as a consequence of collision between Armorica and Gondwana. Despite differences in the polarity of subduction, all these authors attribute Carboniferous deformation, metamorphism and magmatism in the French Massif Central to crustal thickening owing to post Devonian continental collision followed by delamination of the lithospheric root beneath the orogenic belt (Lardeaux et al., 2001; Faure et al., 2002; Lardeaux, 2014; Lardeaux et al., 2014). This model considers implicitely that the UGU represents a single nappe exhumed and emplaced above the LGU before the opening of the late Devonian Brévenne back-arc basin.

3.3 Collision versus syn-orogenic extension during the Carboniferous

The tectonic-geodynamic setting leading to the construction-evolution of the Variscan belt has also been actively debated for the Carboniferous period. The mid- to late-Carboniferous is marked by (i) regional-scale strike-slip shear zones and a NW-SE trending stretching lineation parallel to the belt throughout the French Massif Central (Arthaud and Matte, 1975; Mollier and Bouchez, 1982; Lardeaux and Dufour, 1987; Gébelin et al., 2007, 2009); (ii) the exhumation of metamorphic rocks and the emplacement of syntectonic plutons beneath low-angle detachements (Faure and Pons, 1991; Roig and Faure, 2000); and (iii) fold and thrust belts in the foreland of the Variscan belt as exemplified in the Ardennes (Sintubin et al., 2009) and in the Montagne Noire (Echtler, 1990). These features have been interpreted as reflecting either lateral escape of crustal blocks during progressive construction of the orogenic belt following continental collision (Arthaud and Matte, 1975; Matte, 1986; Burg et al., 1987; Lardeaux and Dufour, 1987; Gébelin et al., 2007), syn-orogenic extension of the thickened orogenic crust (Burg et al., 1993; Faure, 1995), or syn-convergent extension (Gébelin et al., 2009).

For the late Carboniferous-Permian period, gravitational collapse of the Variscan belt has been recognized through the development of a rift system superimposed on the thickened crust (Ménard and Molnar, 1988) and by the identification of low-angle detachments in the French Massif Central (Malavieille et al., 1990; Burg et al., 1993; Faure, 1995; Gardien et al., 1997; Vanderhaeghe et al., 1999; Bouilhol et al., 2006). For the same late Carboniferous period, the prevailing tectonic-geodynamic model for the nearby Pyrénées, south of the FMC, favours a context of transpression (Gleizes et al., 1997; Laumonier et al., 2010; Denèle et al., 2014; Cochelin et al., 2017) although models invoking regional extension have also been proposed (Wickham et al., 1987; Gibson, 1991). On the other hand, a synthesis at the scale of western Europe reveals that the late Devonian to early-mid Carboniferous (Visean) sedimentary record, preserved locally in the internal zone and in the foreland of the Variscan belt, is dominated by platform carbonates and pelagic deposits, which is not in favour of the presence of a Tibetan-type elevated orogenic plateau for this time period and questions the pertinence of the concept of orogenic collapse applied to the Variscan belt (Franke, 2014). Nevertheless, these Visean carbonate deposits have not been identified in the intermediate part of the French Massif Central, in between the Brévenne region and the Montagne Noire. The near absence of late Devonian to middle Carboniferous deposits in the W-FMC, contrasts with their abundance in the E-FMC and suggests either that it was at higher altitude at this time or that these deposits have been eroded. At last, as pointed by Franke et al. (2011), it should be noted that foreland deposits are underlain in the Montagne Noire and in the Pyrénées by rocks affected by HT/LP metamorphism and deformation during the Variscan orogeny, which are not typical characteristics of external zones of orogenic belts.

3.4 Impact of partial melting and magmatism on Variscan tectonics?

Despite the predominance of plutonic and migmatitic rocks in the French Massif Central, and more generally in the Variscan belt (Zwart, 1967), early tectonic models partly eluded the relationships between HT metamorphism, magmatism and orogenic evolution (Matte, 1986, 1991). During the 1980s, progress on experimental petrology and thermodynamic modelling (Wyllie, 1977; Vielzeuf and Holloway, 1988; Spear and Cheney, 1989; Gardien et al., 1995) has allowed to identify the P-T conditions of partial melting during orogenic evolution and discuss their tectonic-geodynamic significance (England and Thompson, 1984; Thompson and Connolly, 1995). In the French Massif Central, partial melting has been identified throughout the high-grade nappe pile and is associated with HP, MP and LP metamorphism. High-pressure granulite to amphibolite facies migmatites dated from the early to the middle Devonian are interpreted to record either partial melting during exhumation of the UGU (Faure et al., 1997, 2008) or to magmatic arc accretion (Lardeaux, 2014). Intermediate pressure migmatites typically associated with Barrovian metamorphism with ages spreading from the middle Devonian to the middle Carboniferous (from c. 370 to c. 310 Ma) are interpreted as reflecting continental collision marked by nappe emplacement (Montel et al., 1992; Ledru et al., 1994; Faure et al., 2008; Barbey et al., 2015) and/or a mantle heat supply owing to asthenospheric upwelling as a consequence of lithospheric delamination (Faure et al., 2002, 2009a). The last partial melting event is characterized by low-pressure migmatites and granites exposed in large domes is attributed to a temperature increase caused by the emplacement of mantle-derived magmas (vaugnerites) during late-orogenic collapse (Montel et al., 1992; Ledru et al., 1994; Barbey et al., 2015). Each of these migmatite types are therefore dominantly regarded as recording discrete, short-lived events of partial melting associated with prograde metamorphism. This is in conflict with the fact that the solubility of accessory minerals in silicate liquids during melting and crystallization is positively correlated with temperature (Harrison and Watson, 1983; Watson and Harrison, 1983; Boehnke et al., 2013). Consistently, thermodynamic modelling shows that zircon and monazite, the most common geochronometers in migmatites, preferentially crystallize during the retrograde (cooling) path and only seldom preserve a record of the prograde path (Kelsey et al., 2008; Yakymchuk and Brown, 2014).

The discrepancies beween tectonic and geodynamic models proposed for different parts of the Variscan belt of Western Europe illustrate the difficulty in tracking the continuity of terranes and sutures. They potentially point to a non-cylindrical structure marked by discontinuous ribbon-shaped continental terranes separated by immature oceanic basins. Furthermore, these discrepancies reflect differences in the interpretation of the geological record in terms of metamorphism, magmatism and deformation. Until now, the existing models have failed to integrate the impact of partial melting and of the generation of magmatic rocks on orogenic evolution. We hereby provide a novel synthetic model for the geodynamic-tectonic evolution of the Variscan belt integrating most of the existing data summarized in the previous paragraphs. In contrast with previous models, we emphasize the pivotal role of protracted partial melting and magmatism on the rheology of the orogenic crust and thus on the tectonic-geodynamic evolution of the Variscan belt.

4 New model for the geodynamic-tectonic evolution of the Variscan belt of Western Europe

4.1 Pre-Variscan configuration: the North Gondwana hyper-extended margin

As presented in Section 2.1, the relative position of the cratons, the numbers and sizes of continental terranes and oceanic basins, and the polarity of subductions are all actively debated owing to discrepancies between paleomagnetic and paleobiostratigraphic data, and to uncertainties in tectonic reconstructions. The ophiolitic assemblages with different protoliths and peak ages as well as contrasting peak temperatures identified in the French Massif Central, in the Armorican Massif, and in the Vosges (Hanmer, 1977; Bernard-Griffiths and Cornichet, 1985; Girardeau et al., 1986; Santallier et al., 1988; Burg et al., 1989a, 1989b; Dubuisson et al., 1989; Mercier et al., 1991a; Bosse et al., 2000; Lardeaux et al., 2001; Ballèvre et al., 2009; Berger et al., 2010a; Skrzypek et al., 2012; Lotout et al., 2018), might represent different sutures that might have corresponded to former rift and/or oceanic basins separating former crustal blocks. This scheme is consistent with the pre-Variscan record of the micaschists, paragneisses and orthogneisses forming the UGU, MAU, LGU and their relative para-autochthon exposed in the French Massif Central and the Pyrenées (Fig. 5, Tabs. 13). Namely, Ediacaran sediments intruded by Cambrian granitic plutons represent a Cadomian basement for the thick detrital sedimentary sequences of lower Paleozoic age deposited during the post-Pan-African dislocation of the Gondwana continent (Duthou et al., 1981, 1984; R’Kha Chaham et al., 1990; Alexandrov et al., 2001; Chelle-Michou et al., 2017; Couzinié et al., 2017, 2019). Other orthogneisses with intrusive contacts (dyke networks, contact metamorphism) of Ordovician age, represent former laccoliths emplaced during extension along the northern margin of the Gondwana continent (Lasnier, 1968; Bernard-Griffiths, 1975; Bernard-Griffiths et al., 1977; Duthou et al., 1981, 1984; R’Kha Chaham et al., 1990; Pin and Marini, 1993; Barbey et al., 2001; Deloule et al., 2002; Roger et al., 2004; Cocherie et al., 2005; Castiñeiras et al., 2008; Melleton et al., 2010; Lotout et al., 2017). The calc-alkaline to tholeiitic chemical signature of the LAC also suggests an emplacement of this bimodal magmatic suite in a continental to oceanic rift environment during the Cambrian-Ordovician (Pin and Lancelot, 1982; Bodinier et al., 1986; Briand et al., 1995; Chelle-Michou et al., 2017). Pre-Neoproterozoic rocks, such as the Archean to Paleoproterozoic Icartian gneisses exposed in the northern Armorican Massif (D’Lemos et al., 1990; Le Corre et al., 1991), have not been identified in the French Massif Central. As already mentioned by other authors (Bouchardon et al., 1989; Pin, 1990; Faure et al., 1997; Lardeaux, 2014; Franke et al., 2017), these data are consistent with the development of a thinned continental margin with one or several immature rifts and/or oceanic basins along the hyperextended northern continental margin of the Gondwana craton.

Accordingly, a two-stage geodynamic model is proposed for the pre-Variscan configuration of the French Massif Central. During the Cambrian (Figs. 16A and 17A), subduction of the Iapetus Ocean along the Gondwana margin was marked by the emplacement of calc-alkaline plutonic and volcanosedimentary rocks (D’Lemos et al., 1990; Le Corre et al., 1991). Crustal extension and opening of rifts at the rear of this active margin during the Ordovician might have been favoured by retreat of the Iapetus slab leading to the opening of the Medio-European Ocean (Figs. 16B and 17B). Such a hyper-extended margin might have spread over several hundred kilometres from the Ordovician to early Silurian and might conciliate the apparent contradiction between paleomagnetic reconstructions implying a distance of more than 2000 km between Gondwana and Laurussia during the Ordovician and paleontologic data pointing to small communicating basins. Accordingly, the Variscan orogenic belt is essentially made of a reworked continental crust made of a Cadomian basement affected by post Pan-African Ordovician hyper-extension, intruded by granitoids and covered by thick volcanic-sedimentary series. Reworking of this hyper-extended margin occurred as a consequence of convergence between Laurussia and Gondwana from the Silurian to the Carboniferous (Figs. 17CE).

Table 1

Geochronological constraints on emplacement or deposition of pre-Variscan rocks in the French Massif Central. 1 =  Melleton et al., 2010; 2 =  Alexandrov et al., 2001; 3 =  Alexandre, 2007; 4 =  Lafon, 1986; 5 =  Bernard-Griffiths, 1975; 6 =  Bernard-Griffiths et al., 1977; 7 =  Gebauer et al., 1981; 8 =  Berger et al., 2010b; 9 =  Lasserre et al., 1980; 10 =  Paquette et al., 1995; 11 =  Pin and Lancelot, 1978; 12 =  Pin and Lancelot, 1982; 13 =  Maurel et al., 2003; 14 =  Faure et al., 2017; 15 =  Lotout et al., 2017; 16 =  Ducrot et al., 1979; 17 =  Cocherie et al., 2005; 18 =  Pitra et al., 2012; 19 =  Roger et al., 2004; 20 =  Lescuyer and Cocherie, 1992; 21 =  Trap et al., 2017; 22 =  Couzinié et al., 2017; 23 =  Caen-Vachette, 1979; 24 =  Bé Mézémé et al., 2006; 25 =  R’Kha Chaham et al., 1990; 26 =  Couzinié et al., 2017; 27 =  Chelle-Michou et al., 2017; 28 =  Duthou et al., 1981.

Table 2

A. Geochronological data on Variscan magmatic rocks in the Western part of French Massif Central. 1 =  Boutin and Montigny, 1993; 2 =  Petitpierre and Duthou, 1980; 3 =  Rolin et al., 1982; 4 =  Gébelin et al., 2007; 5 =  Roig et al., 1996; 6 =  Choukroune et al., 1983; 7 =  Berthier et al., 1979; 8 =  Duthou, 1978; 9 =  Cartannaz et al., 2007a, 2007b; 10 =  Cartannaz, 2006; 11 =  Bé Mézémé, 2005; 12 =  Ducrot et al., 1983; 13 =  Berger et al., 2010a; 14 =  Duthou, 1978; 15 =  Rolin et al., 2009; 16 =  Bernard-Griffiths et al., 1977; 17 =  Lafon, 1986; 18 =  Faure et al., 2008; 19 =  Pin and Paquette, 2002; 20 =  Bernard-Griffiths et al., 1985; 21 =  Bertrand et al., 2001; 22 =  Thiéry, 2010; 23 =  Holliger et al., 1986; 24 =  Joly, 2007; 25 =  Alexandrov et al., 2000; 26 =  Lafon and Respaut, 1988; 27 =  Cuney et al., 2002; 28 =  Gébelin, 2004; 29 =  Roig et al., 2002; 30 =  Monié et al., 2000; 31 =  Faure et al., 2009b; 32 =  Thiéry et al., 2009; 33 =  Gébelin et al., 2009. B. Geochronological data on Variscan magmatic rocks in the Eastern part of French Massif Central: 34 =  Costa and Maluski, 1988; 35 =  Costa, 1990; 36 =  Hottin and Calvez, 1988; 37 =  Do Couto et al., 2016; 38 =  Faure et al., 2002; 39 =  Schulz, 2009; 40 =  Pin (unpublished) cited in Duthou et al., 1984; 41 =  Pin and Barbarin (unpublished) cited in Duthou et al., 1984; 42 =  Saint-Joanis, 1975; 43 =  Kosztolanyi, 1971; 44 =  Vialette (unpublished) cited in Duthou et al., 1984; 45 =  Laurent et al., 2017; 46 =  Cocherie, 2007; 47 =  Gardien et al., 2011; 48 =  Schulz, 2014; 49 =  Couturié et al., 1979; 50 =  Respaut, 1984; 51 =  Pin, 1979; 52 =  Isnard, 1996; 53 =  Lafon and Respaut, 1988; 54 =  Pin, 1981; 55 =  Pin and Lancelot, 1982; 56 =  Legendre et al., 2009; 57 =  Costa, 1989; 58 =  Paquette et al., 1995; 59 =  Pin, 1981; 60 =  Maluski and Monié, 1988; 61 =  Duguet, 2003; 62 =  Thiéry, 2010; 63 =  Delfour and Guerrot, 1997; 64 =  Choulet et al., 2012; 65 =  Pin and Paquette, 1997; 66 =  Faure et al., 2002; 67 =  Duthou et al., 1994; 68 =  Costa et al., 1993; 69 =  Gay et al., 1981; 70 =  Feybesse et al., 1995; 71 =  Duthou et al., 1998; 72 =  Caen-Vachette et al., 1984; 73 =  Gourgaud, 1973; 74 =  Cocherie, 2007; 75 =  Bé Mézémé et al., 2006; 76 =  Mougeot et al., 1997; 77 =  Bouilhol et al., 2006; 78 =  Bé Mézémé, 2005; 79 =  Couzinié et al., 2014; 80 = Costa unpublished cited in Malavieille et al., 1990; 81 =  Didier et al., 2013; 82 =  Batias and Duthou, 1979; 83 =  Briand et al., 2002; 84 =  Caron, 1994; 85 =  Doublier et al., 2006; 86 =  Monié et al., 2000; 87 =  Vialette et al., 1979; 88 =  Brichau et al., 2008; 89 =  François, 2009; 90 =  Vialette and Sabourdy, 1977; 91 =  Hamet and Mattauer, 1977; 92 =  Mialhe, 1980; 93 =  Chauvet et al., 2012; 94 =  Maluski et al., 1991; 95 =  Franke et al., 2011; 96 =  Doublier et al., 2015; 97 =  Whitney et al., 2015; 98 =  Faure et al., 2014; 99 =  Roger et al., 2015; 100 =  Faure et al., 2010; 101 =  Pitra et al., 2012; 102 =  Matte et al., 1998; 103 =  Franke et al., 2011; 104 =  Poilvet et al., 2011.

Table 3

Pressure-Temperature-time data constraining the evolution of metamorphic rocks of the French Massif Central. UGU = Uppper Gneiss Unit, LGU = Lower Gneiss Unit, PAU = Para-autochtonous Unit, EU = Upper Unit, GU/TPU/St SU = Thyviers Payzac- Genis Unit-St Savadour, PFTB = Paleozoic Fold Thrust Belt (Mt Noire). The pressure is expressed in kbar, the temperature in degrees Celius and the ages are given in Ma. In bold are U-Pb on zircon, in italic are Ar-Ar ages on micas, underlined ages are whole rock Rb-Sr ages and the doubly underlined are U-Th- Pb ages on monazite. 1 =  Costa and Maluski, 1988; 2 =  Burg et al., 1989a, 1989b; 3 =  Boutin and Montigny, 1993; 4 =  Berger et al., 2010a; 5 =  Berger et al., 2010b; 6 =  Santallier, 1981; 7 =  Ducrot et al., 1983; 8 =  Bellot and Roig, 2007; 9 =  Gébelin, 2004; 10 =  Costa, 1992; 11 =  Melleton et al., 2009; 12 =  Lafon, 1986; 13 =  Godard, 1990; 14 =  Audren et al., 1987; 15 =  Schulz et al., 2001; 16 =  Schulz, 2009; 17 =  Do Couto et al., 2016; 18 =  Delor et al., 1986; 19 =  Lotout et al., 2018; 20 =  Faure et al., 2008; 21 =  Delor et al., 1986; 22 =  Joanny et al., 1989; 23 =  Bodinier and Burg, 1981; 24 =  Burg et al., 1986; 25 =  Delor et al., 1987; 26 =  Burg and Leyreloup, 1989; 27 =  Costa, 1990; 28 =  Mercier et al., 1991a, 1991b; 29 =  Briand et al., 1988; 30 =  Lardeaux et al., 2001; 31 =  Dufour, 1985; 32 =  Pin and Lancelot, 1982; 33 =  Costa et al., 1993; 34 =  Gardien, 1990; 35 =  Gardien and Lardeaux, 1991; 36 =  Gardien, 1993; 37 =  Gardien et al., 2011; 38 =  Schulz, 2014; 39 =  Schulz et al., 1996; 40 =  Santallier 1981; 41 =  Bellot, 2001; 42 =  Gébelin et al., 2007; 43 =  Caen-Vachette et al., 1984; 44 =  Vitel, 1988; 45 =  Briand and Gay, 1978; 46 =  Briand, 1978; 47 =  Burg et al., 1984; 48 =  Costa, 1989; 49 =  Montel et al., 1992; 50 =  Barbey et al., 2015; 51 =  Monier, 1980; 52 =  Autran and Guillot, 1978; 53 =  Caron, 1994; 54 =  Arnaud et al., 2004; 55 =  Faure et al., 2010; 56 =  Faure et al., 2014; 57 =  Roger et al., 2015; 58 =  Maluski et al., 1991; 59 =  Matte et al., 1998; 60 =  Rabin et al., 2015.

4.2 Late Silurian to Devonian subduction and HP partial melting of terranes issued from the North Gondwana margin

As described in the previous sections, the oldest evidence for the onset of the Variscan orogenic cycle corresponds to the HP metamorphism identified in mafic and ultramafic rocks of the LAC (Fig. 4) pointing to subduction of these units (Lasnier, 1968, 1971; Pin and Vielzeuf, 1983; Bouchardon et al., 1989; Gardien et al., 1990). The presence of UHP mineral assemblages in the LAC of the Lyonnais (Lardeaux et al., 2001) and Limousin (Berger et al., 2010a, 2010b), argues for burial of these units to a depth of more than 100 km. Nevertheless, two types of eclogites are distinguished, namely (i) HT eclogites that are forming enclaves-boudins into migmatites of the UGU as exemplified in the Aigurande, north Limousin, Sioule, Livradois, Lyonnais, Rouergue regions (Pin and Vielzeuf, 1983; Dufour, 1985; Burg et al., 1989a, 1989b; Gardien, 1990; Gardien et al., 1990; Lardeaux et al., 2001; Faure et al., 2008); and (ii) LT eclogites that are hosted by micaschists not affected by partial melting as illustrated by the MAU exposed in south Limousin or in Najac (Berger et al., 2010a; Lotout et al., 2018). According to the few ­available geochronological data, subduction of these units is diachronous and occurred from the late Silurian to the Devonian (Fig. 18AC) (Pin and Lancelot, 1982; Ducrot et al., 1983; Paquette et al., 1995; Berger et al., 2010a; Do Couto et al., 2016; Lotout et al., 2018).

A similar nappe pile has been described in the South Armorican Massif, with the distinction of HT and LT eclogites (Ballèvre et al., 2009). The core of the Champtoceaux nappe is made of granulite facies migmatitic gneisses, with enclaves of ultramafic rocks with HT eclogite facies relics (Marchand, 1981; Ballèvre et al., 1987, 2002; Godard, 2001). These rocks overlay lower grade orthogneisses and micaschists that contain mafic boudins with LT eclogite facies relics. LT eclogites and blueschist facies rocks of the Ile de Groix are in a similar structural position as the MAU in the FMC and have been interpreted to represent the suture of the Galicia-South Brittany Lower Paleozoic Ocean. 40Ar/39Ar and Rb-Sr dating of blueschist facies rocks, yield c. 360 Ma for the metamorphic peak attributed to subduction of this ocean at the Devonian-Carboniferous transition. Greenschist facies retrogression of these rocks reflects cooling and exhumation dated in the early Carboniferous as constrainted by ages ranging from 355 to 340 Ma (Bosse et al., 2000, 2005; Paquette et al., 2017).

Mafic-ultramafic enclaves with relictual HT eclogitic facies metamorphism enclosed in granulite facies migmatites was also identified in other Variscan massifs. In the Vosges, garnet-lherzolite belonging to a subcontinental lithospheric mantle exhumed from more than 150 km are present in the varied gneiss unit dominated by granulite facies migmatitic gneiss (Altherr and Kalt, 1996; Brueckner and Medaris, 2000; O’Brien and Rötzler, 2003). The Bohemian massif also comprises large mafic and ultramafic bodies enclosed in high-grade gneisses (Kusbach et al., 2012). In the French Massif Central, ages interpreted to record eclogite facies peak metamorphism are typically 20 to 30 Ma older to the ones attributed to granulite facies metamorphism and partial melting of their host migmatitic paragneisses and orthogneisses (Duthou et al., 1981, 1994; Pin and Lancelot, 1982; Chelle-Michou et al., 2017).

The significance of this association and of this age gap is debated. The question is whether it represents (i) a pre-metamorphic association reflecting intrusion of mafic magmas into the crust or tectonic accretion of oceanic and continental terranes, (ii) a syn-metamorphic extrusion of mantle into the lower crust, (iii) or a subduction of the continental crust and mixing of mantle and crustal units. Tectonic-geodynamic models attributing the succession of eclogite and granulite facies metamorphism to oceanic subduction followed by continental collision (Matte, 1986; Girardeau et al., 1986; Dubuisson et al., 1989; Ledru et al., 1989; Lardeaux et al., 2001) provide an explanation for the HP metamorphic condition but do not account for the systematic incorporation of the LAC into the migmatitic UGU. Models invoking a vertical extrusion of the mantle into the partially molten granulitic lower crust account for the presence of mafic and ultramafic rocks into the migmatitic continental rocks but do not propose a driving force for such a process (Kusbach et al., 2012). A more appealing proposition is that pieces of suprasubduction mantle might be incorporated into the orogenic belt during relamination of orogenic crust by flow of partially molten crustal units decoupled from the downgoing slab (Lexa et al., 2011; Kusbach et al., 2015). Another option is that the association of remnants of an oceanic suture together with lithospheric mantle into granulite facies migmatitic gneisses, was achieved by mixing of subducted units into the mantle (O’Brien and Rötzler, 2003; Faure et al., 2008). In this latter scenario, granulite-facies partial melting of the UGU would be driven either by decompression of the UGU during syn-orogenic exhumation (Faure et al., 2008) or owing to thermal relaxation about 30 Ma after subduction (O’Brien and Rötzler, 2003). We concur to the latter proposition as it will be developed in the next section.

4.3 Middle Devonian to early Carboniferous syn-orogenic exhumation of the partially molten subducted crustal units with mantle enclaves

The geological record of the middle Devonian to early Carboniferous period is marked by the association of HT eclogites and HP granulite-facies migmatites of the UGU. These are only present south of the Nort-sur-Erdre Fault in the Armorican Massif and of the Lalaye-Lubine Fault in the Vosges Massif (Fig. 1), which delineate the main suture south of the Armorica-Barrandia continental block (Faure et al., 1997). North of this suture, the inprint of Variscan metamorphism is absent or limited. The typical retrogression and transposition of the granulite facies mineral assemblage and foliation of the UGU into an amphibolite facies foliation records isothermal decompression caused by rapid exhumation (Pin and Vielzeuf, 1983; Dufour, 1985; Burg et al., 1989a, 1989b; Gardien et al., 1990; Mercier et al., 1991a, 1991b; Lardeaux et al., 2001; Faure et al., 2008) (Fig. 4). The position of the exhumed high-grade nappes relative to the Nort-sur-Erde Fault and Lalaye-Lubine Fault implies that this suture was reworked as a steep syn-orogenic detachement allowing the exhumation of the continental units previously entrained in subduction and then decoupled from the slab (Figs. 18B and 18C). The inverted metamorphic gradient at the contact between the UGU and the LGU (Nicollet, 1978; Burg et al., 1984; Burg et al., 1989a, 1989b; Schulz et al., 2001) suggests that the exhumation of the UGU is associated with burial of the LGU. This is consistent with transposition of the granulite facies foliation of the UGU into an amphibolite facies synmigmatitic foliation (Forestier, 1961; Burg and Matte, 1978) and with the absence of HP/HT eclogitic relicts in the LGU. Thrusting of the UGU over the LGU is locally corroborated by kinematic data indicating a top-to-the SE sense of shear (Faure et al., 1979, 2009a; Burg et al., 1984). However, this contact is in many places reworked and retrogressed into a greenschist facies fabric associated with a top-to-the NW sense of shear attributed to regional extension leading to exhumation of the LGU and of the PAU (Faure et al., 1979, 2008, 2009a).

Available geochronological and stratigraphic data indicate that exhumation of the UGU started before the late Devonian with the emplacement of the northern nappe exposed in the Aigurande plateau, the Morvan area, and sampled in the Couy borehole (Costa and Maluski, 1988; Godard, 1990; Boutin and Montigny, 1993). The exhumation of the UGU forming the central and southern nappe (Limousin, Sioule-Combrailles, Livradois – Haut-Allier, Lyonnais) is bracketed between the late Devonian and the Visean and thus postdates the exhumation of the northern nappe (Pin and Lancelot, 1982; Melleton et al., 2009; Gardien et al., 2011; Do Couto et al., 2016; Chelle-Michou et al., 2017). Migmatitic paragneisses from the Limousin and the Sioule display a subhorizontal foliation bearing a NW-SE trending lineation mostly associated with top-to-the-NW kinematic criteria interpreted to accommodate syn-orogenic exhumation of the UGU (Bellot and Roig, 2007; Do Couto et al., 2016). In the Livradois, the foliation of the migmatites attributed to the UGU is steeply-dipping and is associated with dextral top-to-the SE kinematic criteria interpreted to record extrusion of the UGU (Gardien et al., 2011; Vanderhaeghe et al., 2013). Such a mechanism of vertical extrusion for the formation of high-grade nappes is consistent with the interpretation of the Champtoceaux nappe, in the Armorican Massif, as a fold nappe (Ballèvre et al., 2009) and has also been proposed to account for the structure of the varied gneiss unit in the Vosges (Skrzypek et al., 2014) and the Bohemian Massif (Schulmann et al., 2014).

Several lines of reasoning suggest that the multiple nappes model is more coherent with the hyperextended margin model made of thinned continental blocks separated by small immature oceanic basins than with a pre-Variscan configuration characterized by a single, large oceanic domain between Gondwana and Armorica. First, the inverted metamorphic isograds, requiring limited thermal relaxation after underthrusting, fits better with a model considering a system of several small nappes with at most a few tens of kilometres of lateral expansion, rather than the case of a single, several hundred kilometres long nappe. Indeed, thermal-mechanical modelling of orogenic evolution indicates that inverted metamorphic/thermal gradients are not sustainable for more than about 100 km equivalent to about 10 Ma with a convergence rate of 1 cm/yr (Huerta et al., 1996; Henry et al., 1997; Vanderhaeghe et al., 2003). Moreover, the multiple nappes model predicts diachronous HP metamorphism, which seems to be confirmed by the currently available geochronological data. Finally, the absence of large volumes of calc-alkaline magmatism associated with subduction initiation and steady-state subduction further suggests that forced subduction of several hyper-extended basins is a more likely scenario than the spontaneous subduction intiation of a mature oceanic crust (see McCarthy et al., 2018).

Following this rationale, following other authors, we propose that the construction of the Variscan orogenic belt was achieved by tectonic accretion of continental units that were previously subducted and then decoupled from the downgoing slab (Faure et al., 1997, 2008; Lardeaux et al., 2001) rather than by continental collision and indentation and thickening of the overriding plate as proposed in early tectonic reconstructions (Matte, 1986; Franke, 1989; Ledru et al., 1989). Again, available geochronological data for the UGU and LAC are consistent with partial melting at high-pressure (10–20 kbars) of paragneisses, orthogneisses and amphibolites c. 20 to 30 Ma after subduction (Nicollet and Leyreloup, 1978; Duthou et al., 1981, 1994; Pin and Lancelot, 1982; Chelle-Michou et al., 2017). Although such data are scarce and of unequal robustness from a nappe to another, this time span is broadly consistent with that required for thermal relaxation after burial (England and Thompson, 1984; Vanderhaeghe et al., 2003). Following models advocated for the tectonic accretion of HP mafic-ultramafic rocks enclosed in migmatitic units (O’Brien and Rötzler, 2003; Labrousse et al., 2011; Gordon et al., 2016; Závada et al., 2018), we propose that partial melting of subducted units forming the protolith of the UGU triggered their mechanical decoupling from the slab. In this scenario, mafic and ultramafic enclaves would represent fragments of previously subducted oceanic crust and lithospheric mantle entrained by the buoyant and low-viscosity, partially molten rocks, on their way back towards the surface. This process corresponds to syn-orogenic exhumation by vertical extrusion of the partially molten nappes. According to this proposition, the UGU represents the part of the thinned continental margin that has been subducted and then decoupled from the downgoing slab, while the LGU represents the part of the former continental margin that was underthrusted beneath the UGU during its syn-orogenic exhumation. The age gradient for retrogression and cooling of the UGU, ranging from c. 385 Ma in the north to c. 340 Ma in the south, is consistent with progressive exhumation of nappes successively decoupled from a north-vergent and southward-retreating subduction slab (Fig. 18). We defined here a northern, central and southern nappe but this proposition and the number of nappes and of repeated UGU/LGU alternations remains to be clarified. Another consequence of partial melting of subducted fertile continental units might be the percolation of felsic melts into the suprasubduction mantle wedge, contributing to its enrichment in incompatible elements. Such a process might explain the signature of mafic magmas emplaced during the Carboniferous, as described in the next section.

4.4 Late Devonian exhumation of nappes in between retreating slabs

The geological record of the late Devonian period (Figs. 11, 18B and 18C) is marked by a great diversity of information that is difficult to reconcile in a single geodynamic context. In the northern part of the French Massif Central, plutonic and volcanic rocks (ACG) exposed respectively in the Somme and Limousin regions have been interpreted as a continental magmatic arc (Pin and Paquette, 1997; Faure et al., 2008). The tholeiitic to calc-alkaline volcanic rocks of the Brévenne Unit have been attributed to a back-arc rift basin (Bébien, 1971; Pin and Lancelot, 1982; Sider and Ohnenstetter, 1986). Rocks of the same age and geochemical signatures exposed in the Central Bohemian Plutonic Complex have been interpreted as arc magmas originated in an Andean-type continental margin (Janousek and Holub, 1997; Janoušek et al., 2004; Schulmann et al., 2009). The positive εNd of these magmatic rocks points to a dominant juvenile contribution (Fig. 13A). However, structural, metamorphic and geochronologic record demonstrate that the late Devonian is also marked by (1) the exhumation of the central nappe in the Limousin and Livradois (Bellot and Roig, 2007; Melleton et al., 2009; Gardien et al., 2011; Do Couto et al., 2016); and (2) the emplacement of cordierite-bearing peraluminous granite laccoliths (CPG) (Berthier et al., 1979; Pin and Paquette, 1997; Bertrand et al., 2001; Cartannaz, 2006; Gébelin et al., 2009). The crust-derived Guéret-type plutons formed through melting of dominantly metasedimentary sources, possibly with the contribution of a mafic igneous lower crust (Downes and Duthou, 1988; Downes et al., 1997). These data imply that the end of opening of the Brévenne rift and the construction of the Morvan continental arc are coeval with partial melting of subducted continental units. Furthermore, this activity is contemporaneous to deposition of detrital sediments and carbonates in the underfilled southern foreland basin at the front of a propagating thrust system, currently exposed in the Pyrénées (Franke and Engel, 1986; Souquet et al., 2003).

To the south of Armorica, geochronological data obtained on blueschists exposed in the Ile de Groix and on low-temperature eclogites at the base of the Champtoceaux nappe indicate that high-pressure metamorphism attributed to subduction occurred at about 370–360 Ma (Bosse et al., 2000, 2005; Paquette et al., 2017). North of the Nort sur Erdre fault, the Saint Georges sur Loire unit, a late Silurian to early Devonian back-arc basin, affected by blueschist facies metamorphism, has also potentially been subducted at this time (Ledru et al., 1986; Cartier et al., 2001). Mid to late Devonian subduction, but with a southward dip, is likewise recorded north of the Medio-European suture, in the Saxo-Thuringian terrane to the north of the Vosges and of the Bohemian Massif (Schulmann et al., 2014; Skrzypek et al., 2014). In the Vosges Massif the “ligne des Klippes” represents an immature oceanic basin opened during the late Devonian and inverted in the Lower Carboniferous (Skrzypek et al., 2012). North of Armorica, in South West England, the Lizard suture represents an early Devonian immature oceanic domain obducted during the early Carboniferous (Floyd and Leveridge, 1987; Clark et al., 1998; Shail and Leveridge, 2009).

Accordingly, late Devonian extension in the internal zone of the Variscan belt occurred in an overall context of plate convergence as implied by subduction/exhumation of nappes, foreland propagation of deformation and deposition of clastic sediments. Therefore, two geodynamic scenari are evoked to account for this extension, namely (i) opening of a back-arc associated with northward subduction of the Medio European Ocean being contemporary with the southward subduction of the Rheic Ocean followed by continental collision (monocyclic model, Matte, 1986; Ledru et al., 1989; Lardeaux et al., 2001), or (ii) renewed extension of the Eo-Variscan orogenic belt associated with southward subduction of the Rheic Ocean after closure of the Medio European Ocean by northward subduction (polycyclic model, (Pin, 1990; Leloix et al., 1999; Faure et al., 1997, 2009a; Lardeaux, 2014; Lardeaux et al., 2014). The cooling ages obtained on the UGU of the Aigurande plateau indicates that at least some of the high-grade nappes were exhumed to the surface before late Devonian rifting and emplacement of the Brévenne volcanics, which is more consistent with the polycyclic model.

Alternatively, to resolve the apparent discrepancy between an overall context of plate convergence and local extension, we propose that late Devonian geodynamics was characterized by convergence accommodated by subduction but also by retreat of the northern (Rheic) and southern (Medio-European) slabs. This explains simultaneous closure of immature oceanic basins entrained in subduction, tectonic accretion of ribbon-shaped terranes and extension of the upper plate. Hence, the exhumation of previously subducted continental units might have occurred without substantial erosion owing to the opening of space toward the surface in between the retreating slabs as envisioned in the conceptual models of Vanderhaeghe and Duchêne (2010). In this scenario, the top to the NW kinematic criteria marking the contact between the Brévenne Unit and the UGU of the Monts du Lyonnais might have accommodated the last stage of syn-orogenic exhumation of the high-grade nappe as proposed for the exhumation of the UGU in the Limousin at the same period (Bellot and Roig, 2007).

However, this model does not fully account for lateral variations of the tectonic context, evidenced by exhumation of nappes and emplacement of peraluminous granites in the western French Massif Central coeval with opening of a continental to oceanic rift in its eastern part. Two directions are proposed for future reflexion. The first one would be to consider a context of oblique plate convergence, which might allow for opening of pull-apart basins along releasing bends and concomitant burial and exhumation of tectonic units along restraining bends. The second one would be to infer differential retreat between the western and eastern side of the French Massif Central that might have been accommodated by a precursor of the Sillon Houiller strike-slip fault (Fig. 3). Slab retreat potentially controls the tectonic-magmatic geological history of the French Massif Central throughout the Carboniferous as described in the following sections.

thumbnail Fig. 11

Late Devonian to Early Carboniferous magmatism. The western part of the French Massif Central is dominated by plutonic rocks of the ACG-type and CPG-type while the eastern part of the French Massif Central also comprises volcanics with a tholeiitic to calc-alkaline signature.

4.5 Carboniferous (c. 345–310 Ma) building of an orogenic plateau by lateral flow of the partially molten orogenic root

The Carboniferous (345–310 Ma) geological record of the FMC is marked in the northern part by E-W to NW-SE trending regional-scale dextral strike-slip shear zones and orogen-parallel stretching localizing the emplacement of mantle and crust-derived magmas and in the southern part by thrusts and folds in the Cévennes, Albigeois and Montagne Noire, coeval with deposition of clastic sediments in the foreland (Figs. 1, 3, 12, 18D) (Arthaud and Matte, 1975; Engel, 1984; Feist and Galtier, 1985; Arnaud and Burg, 1993; Faure et al., 1999). The strike-slip shear zones, exemplified in the Limousin (Faure and Pons, 1991; Mollier and Bouchez, 1982; Roig and Faure, 2000; Gébelin et al., 2007, 2009), in the Forez (Barbarin and Belin, 1982), in the Livradois (Gardien et al., 2011; Vanderhaeghe et al., 2013), and in the Lyonnais (Lardeaux and Dufour, 1987), merge to the northwest with the South Armorican shear zone (Lerouge and Quenardel, 1988; Carlier De Veslud et al., 2004; Gébelin et al., 2007; Rolin et al., 2009, 2014). As presented above, these strike-slip shear zones have been interpreted either to record a period of transpression (Lardeaux and Dufour, 1987; Gébelin et al., 2007) or to reflect the transition from collision to orogen-parallel extension (Ledru and Autran, 1987; Faure et al., 1993, 2009b; Roig and Faure, 2000).

After a 10 Ma gap from c. 355 to c. 345 Ma, magmatism proceeded with the emplacement of plutonic and volcanic rocks throughout the middle Carboniferous (Fig. 12). In Section 3.4, following previous authors (Williamson et al., 1996, 1997; Barbey et al., 2015; Laurent et al., 2017; Moyen et al., 2017), MPG and CPG (and their volcanic counterparts) are proposed to derive from a mixed crustal source comprising ortho- and paragneisses from the LGU. Partial melting of the orogenic root at this time is documented by MP migmatites dated at c. 315 Ma underlying the Millevaches granitic laccolith (Gébelin et al., 2009) and by MP migmatites mantling the Velay dome with ages ranging from c. 330 to c. 315 Ma (Mougeot et al., 1997; Cocherie et al., 2005; Bé Mézémé et al., 2006). Although volumetrically less abundant, mantle-derived magmatism is expressed almost exclusively in the East Massif Central (only a few occurences are described in Limousin) by the presence of widespread vaugnerites, forming decimeter- to kilometer-sized bodies intimately associated with the granites and migmatites. The dual geochemical character of vaugnerites (low SiO2, high Mg# together with elevated K, LILE and LREE contents), as well as their Sr-Nd-Pb-Hf isotope signatures (Turpin et al., 1988; Couzinié et al., 2016), is consistent with partial melting of a lithospheric mantle that was previously enriched in incompatible elements (Sabatier, 1991; Solgadi et al., 2007; von Raumer et al., 2014 Couzinié et al., 2016). The variety of the geochemical signatures of some granitoids and the fact that they contain micromafic enclaves indicate a mixture of crustal-derived and mantle-derived magmas, as illustrated in the Livradois (Solgadi et al., 2007). This model would be particularly relevant to explain the origin of the compositionally intermediate KCG and, to a lesser extent CPG. The Margeride granite, for example, shows an intermediate signature, interpreted as reflecting interaction with mafic melts (Williamson et al., 1992; Laurent et al., 2017). Collectively, granites exposed in the eastern French Massif Central show the implication of different sources that are progressively affected by partial melting as the temperature of the orogenic root increases (Fig. 14). Geochronological data indicate that (i) granites and vaugnerites emplaced together from c. 340 to c. 300 Ma and (ii) there is a progressive younging of U-Pb emplacement ages of both granites and vaugnerites from the North (Forez, Lyonnais) to the South (Cévennes) in the considered time period, pointing to the southward migration of a lithospheric-scale thermal anomaly resulting in both crust- and mantle melting (Laurent et al., 2017).

In the North-Eastern Massif Central, the extrusive equivalents of the granitoids record a drastic shift in composition during the middle Carboniferous. They change from mafic, low to medium-K calc-alkaline to felsic and high K (potassic rhyolites, and associated microgranites), broadly similar to MPG in composition. These rhyolites are middle to upper Visean in age as constrained by a U-Pb zircon age of 336.9 ± 3.2 Ma (Cartannaz et al., 2007a), identical within uncertainty to the age of the underlying granitoids (Laurent et al., 2017). Moreover, both the early, and the middle Carboniferous magmatic rocks do include some mantle-related components, with nevertheless clear differences in the nature of the mantle component between these two periods (Fig. 14A). Up to the early Carboniferous (before c. 345–340 Ma), the isotopic composition of the mantle-derived magmatic rocks (late Devonian Brévenne and lower Carboniferous series) are on the mantle array, consistent with an asthenospheric mantle source. On the other hand, the middle (and late) Carboniferous mantle-derived magmatic rocks (vaugnerites and lamprophyres) are shifted towards “crustal” compositions, reflecting the growing influence of a recycled crustal component in the mantle below the French Massif Central (Fig. 13). A similar and coeval evolution is described in the Bohemian Massif (Janousek and Holub, 1997; Janoušek et al., 2004). In the French Massif Central, there is not only a temporal evolution, but also a spatial distribution of mantle sources. While the enriched mantle is centered on the Velay complex, the non-enriched mantle domain is located in the Northern part of the region (Beaujolais and Morvan). Similarly, different mantle compositions are still featured today in mantle xenolith from Cenozoic volcanoes (Lenoir et al., 2000), displaying a comparable spatial distribution. Furthermore, the volumetrically more abundant Cenozoic magmatism occurs in the domain where the mantle is more enriched (and probably more fertile).

The preservation of platform carbonates of Visean age in the northern part of the French Massif Central and around the Montagne Noire, indicates that these regions were below sea level at this time. Other sedimentary sequences of this period corresponds to the erosion products of subcontemporaneous volcanics and plutonics traping organic matter in continental basins delimited by strike-slip shear zones or normal faults (Bertaux et al., 1993; Bruguier et al., 1998; Thiéry et al., 2009).

In order to reconcile all observations, we propose that exhumation of migmatites along strike-slip shear zones coeval with foreland propagation of the deformation front correspond to growth of an orogenic plateau by lateral flow of the partially molten orogenic root in a context of plate convergence associated with southward slab retreat (Figs. 17 and 18D). This plateau spreads from south of the Brévenne region to north of the Montagne Noire as constrained by the deposition of plateform carbonates. In this model, the partially molten orogenic root also flows toward the north beneath the former late Devonian Brévenne rift basin. Such horizontal flow, already proposed for the Vosges (Skrzypek et al., 2014) and the Bohemian Massifs (Lexa et al., 2011; Schulmann et al., 2014; Kusbach et al., 2015), explains the voluminous granitic magmatism observed at that time. The systematic occurrence of vaugnerites along with granitic magmas point to mantle melting that might be caused by decompression as a consequence of slab retreat.

thumbnail Fig. 12

Middle Carboniferous magmatism. Magmatic rocks with ages ranging from c. 345 to c. 310 Ma are widespread throughout the northern part of the French Massif Central indicating the presence of a partially molten source at depth during this period. Plutonic rocks display a variety of geochemical signatures encompassing MPG-type, CPG-type, and KCG-type pointing to the contribution of crustal and mantle sources.

thumbnail Fig. 13

Geochemistry of magmatic rocks of the French Massif Central. a) Two mantle sources in the Massif Central (calculated at 315 Ma, an average between the c. 335 Ma lavas and the ca 305 Ma lamprophyres and vaugnerites): note the clear difference between the pre-335 mafic magmas (lavas from the Brévenne Unit (Pin and Paquette, 1997) and diverse lavas from NE Massif Central (Pin and Paquette, 2002)), and the post-335 Ma lavas (lamprophyres: (Agranier 2001); enclaves in granites (Pin, 1990) and vaugnerites (Williamson et al., 1992); underplated mafic magmas (or cumulates) found as enclave in the Cenozoic Bournac volcano (Downes et al., 1990)). b) Change in the nature of the lavas, in Shand (1943) A/CNK vs. A/NK and Peccerillo and Taylor (1976) SiO2–K2O diagrams. Pre-335 Ma lavas are mafic and metaluminous, whereas post-335 Ma lavas are felsic, high-K and peraluminous, essentially similar to MPG granites forming at the same age.

thumbnail Fig. 14

Geochemical characteristics of magmatic rocks exposed in the East French Massif Central. Summary of the geochemical characteristics of E-FMC granitoids, in a A/CNK (molar Al2O3/CaO + Na2O + K2O diagram, Shand 1943) vs. FSMB ((FeOt + MgO) × (Sr(wt.%) + Ba(wt.%)) diagram (Laurent et al., 2014). This diagram separates granitoids related to different sources (Moyen et al., 2017), and shows that KCG are primarily related to the differentiation of vaugnerites (with minor crustal components occasionaly); MPG are related to a metasedimentary source; CPG are generated from a source dominated by orthogneisses, but with more common involvement of either a metasedimentary component (particularly pronounced in the Velay complex) or a mafic component (e.g. the Margeride granite).

4.6 Late Carboniferous to Permian (c. 305–295 Ma) gravitational collapse and exhumation of the partially molten root of the Variscan orogenic belt

The late Carboniferous (Figs. 15 and 18E) is marked by regional scale extension of the Variscan belt of Western Europe (Ménard and Molnar, 1988), which is particularly illustrated in the French Massif Central by the activation of low-angle detachments (Malavieille et al., 1990; Van Den Driessche and Brun, 1992; Burg et al., 1993; Faure, 1995; Gardien et al., 1997; Vanderhaeghe et al., 1999; Bouilhol et al., 2006) controlling the exhumation of migmatites in the core of domes such as the Velay in the central-east French Massif Central (Dupraz and Didier, 1988; Burg and Vanderhaeghe, 1993; Ledru et al., 2001) and the Montagne Noire in the southern French Massif Central (Gèze, 1949; Nicolas et al., 1977; Echtler and Malavieille, 1990; Van Den Driessche and Brun, 1992; Rabin et al., 2015; Trap et al., 2017) (Fig. 3).

The HT/LP metamorphic conditions recorded by the migmatites coring the Velay and Montagne Noire domes and by the lower crustal granulites have been classically interpreted in terms of a sudden increase in temperature at the end of the Carboniferous (Marignac et al., 1980; Dupraz and Didier, 1988; Lardeaux et al., 2001; Ledru et al., 2001; Couzinié et al., 2014; Barbey et al., 2015). However, such a catastrophic event is not required by the data. Indeed, as discussed in the previous section, (i) the oldest ages obtained on migmatites in the Limousin and around the Velay and Montagne Noire domes, and (ii) the coeval emplacement of granitic plutons and vaugnerites indicates that the crust and the mantle were already partially molten as early as 340 Ma ago and remained so throughout the Carboniferous (Laurent et al., 2017; Vanderhaeghe et al., 1999). A similar time range has been proposed for the duration of the partial melting event in the Variscan crustal segment exposed in the Ivrea-Verbano Zone (Guergouz et al., 2018) but also for the root of the Grenville orogenic belt (Turlin et al., 2018). The continuous magma extraction from a partially molten source at depth results in the emplacement of granitic plutons (CPG, MPG) at higher structural level in combination with the entrainment of source material, which fully explains their chemical characteristics (Villaros et al., 2018). Final crystallization of migmatites is constrained by the ages of late-migmatitic dikes at c. 297 Ma in the Velay (Montel et al., 2002; Didier et al., 2013) and by the emplacement of a syntectonic granite at c. 294 Ma in the Montagne Noire (Poilvet et al., 2011).

The size of the Velay and Montagne Noire domes, i.e. several tens of km in diameter, indicates that they are crustal-scale features (Van Den Driessche and Brun, 1992; Burg and Vanderhaeghe, 1993; Ledru et al., 2001). This view of a pervasively partially molten orogenic root is corroborated by the diversity of the protoliths (orthogneisses, paragneisses, amphibolites…) that have molten to form the migmatites (Downes and Duthou, 1988; Williamson et al., 1992; Downes et al., 1997; Ledru et al., 2001; Barbey et al., 2015; Rabin et al., 2015). Moreover, the fact that the domes are circumscribed by the foliation of the host paragneisses and orthogneisses (Lagarde et al., 1994; Echtler and Malavieille, 1990; Mattauer et al., 1996) and by the syn-melting foliation of the migmatites (Burg and Vanderhaeghe, 1993; Ledru et al., 2001; Rabin et al., 2015) indicates that the migmatites correspond to a mechanically coherent partially molten body. In the Velay, the synmigmatitic foliation delineates subdomes of about 10–20 km in diameter and bears a radially distributed HT mineral lineation pointing to the role of gravitational instabilities (Ledru et al., 2001; Vanderhaeghe, 2009). Such subdomes have been described in Naxos Aegean domain and interpreted as reflecting crustal-scale convection (Kruckenberg et al., 2011; Vanderhaeghe et al., 2018). In the Montagne Noire, migmatites in the core of the dome display a prolate finite strain ellipsoid indicative of constriction with a subhorizontal long axis parallel to the axis of the elliptical shape dome, which is coeval with exhumation of the core of the dome (Echtler and Malavieille, 1990; Van Den Driessche and Brun, 1992; Mattauer et al., 1996; Aerden and Malavieille, 1999; Charles et al., 2009; Rabin et al., 2015). Following this rationale, the development of crustal-scale migmatite domes at 305–295 Ma marks the time at which the long-lived, partially molten lower-middle crust cooled rapidly as a consequence of rapid exhumation owing to gravitational collapse. If crustal thinning was faster than thermal relaxation, isothermal decompression of the migmatites would result in an increase of the geothermal gradient as observed in the E-FMC (Montel et al., 1992; Barbey et al., 1999, 2015; Ledru et al., 2001).

To explain the protracted melting of the root of the Variscan orogenic crust, one may invoke the effects of an increase in radioactive heat production of the thickened crust, combined with the increase in mantle heat flux owing to the southwards removal of the lithospheric mantle slab underneath the FMC. Coeval crustal thickening and lithospheric mantle thinning is the best case scenario to produce a high geothermal gradient in the continental crust (Vanderhaeghe and Duchêne, 2010) that might last for several tens of Myrs (Ueda et al., 2012). This model is also supported by the nature and spatial/temporal evolution of granitoid and vaugnerite magmatism (Moyen et al., 2017; Laurent et al., 2017). In addition to the thermal input at the base of the orogenic crust, partial melting might also be enhanced by exhumation, as dehydration melting reactions are crossed during decompression (Thompson and Connolly, 1995).

In the external zone of the Variscan belt, specifically in the Montagne Noire area, migmatites are present beneath and deform a metasedimentary sequence affected by recumbent folds and low-grade metamorphism. This is paradoxal as this scheme does not match the classical model of thermal relaxation after nappe stacking (Franke et al., 2011). Moreover, the deposition age of the protolith of part of the metasedimentary sequence (Visean to Namurian) overlaps with the geochronological record of HT/LP metamorphism and granitic intrusion, spreading from c. 330 to c. 300 Ma. In order to resolve this paradox, we propose that sediment deposition, HT/LP metamorphism and granite emplacement were indeed coeval but occurred in laterally remote units that were subsequently juxtaposed owing to lateral horizontal flow of the partially molten orogenic root beneath the sedimentary rocks and their upper crustal basement (Figs. 18E and 19). In this scenario, the migmatites coring the Montagne Noire dome represent partially molten rocks that were located beneath the orogenic plateau since the early Carboniferous and have flown laterally from the internal to the external zone in the late Carboniferous. The migmatites and associated granites were then exhumed, cooled and juxtaposed to the metasedimentary nappes along a detachment. The Variscan basement exposed in the Axial Zone of the Pyrenees and in the North Pyrenean Massifs displays similar geological features (de Saint Blanquat, 1993; Gleizes et al., 1997; Cochelin et al., 2017) and the partially molten orogenic root might have flown toward the foreland beneath the current day Pyrenees providing an explanation for their peculiar Carboniferous-Permian structural and metamorphic record.

Comparison between the eastern and the western part of the FMC shows that large, late Carboniferous migmatite domes are only exposed in the eastern part of the FMC, which points to the Sillon Houiller as a major tectonic divide between these two parts of the FMC, at least during the late orogenic tectonic evolution. The Sillon Houiller corresponds to a subvertical sinistral strike-slip shear zone (Grolier and Letourneur, 1968; Arthaud and Matte, 1975) cross-cuting the Lower to Middle Carboniferous transcurrent strike-slip shear zones. It is in places sealed by Visean volcanics, and localizes the deposition of late Carboniferous to Permian clastic sediments accompanied with the emplacement of rhyolite at the onset of extension (Bonijoly and Castaing, 1984; Joly et al., 2008, 2009; Lapierre et al., 2008; Thiéry et al., 2009). Scarce U-Pb data indicate an onset of deposition in some basins as early as 330 Ma (Bosmoreau and Decazeville Basins), but most ash beds yield ages of 300–295 Ma that constrain a predominance of syntectonic intramontane basins (Jaujac, Bosmoreau, Alès, Bertholène, Graissessac and Roujan-Neffies basins) during the late Carboniferous (Bruguier et al., 1998, 2003). In the Livradois, cooling of the migmatites and granites is recorded by argon thermochronology on micas and K-feldspar between 307 and 300 Ma. It is associated with exhumation along a top-to the west low-angle detachment that controls the deposition of coal-bearing sediments in the Brassac basin (Gardien et al., 2011; Vanderhaeghe et al., 2013). The discordant contact of these deposits with the LGU (locally migmatitic) and UGU attests that these metamorphic units reached surface exposure at the end of Carboniferous and that exhumation and crystallization of the partially molten orogenic crust and of crustal melts was essentially terminated by 295 Ma.

In the western Massif Central, this period is marked by ductile-brittle detachments such as the Argentat fault that also accommodate a component of strike-slip displacement (Bellot and Roig, 2007). Accordingly, the activity of the Sillon Houiller straddles the transition from the orogenic plateau development to gravitational collapse of the Variscan belt. The limited lateral offset of the Sillon Houiller compared to its the length is consistent with its role as a transfer fault accommodating a larger amount of N-S extension in the eastern part of the French Massif Central relative to the western part (Burg et al., 1990). The presence of migmatite domes solely in the eastern part of the French Massif Central is consistent with this model. Currently, except for a slight variation beneath the Cenozoic rift, the continental crust displays a constant thickness of c. 30 km on both sides of the Sillon Houiller (Ziegler and Dèzes, 2006). This suggests that, at the time of orogenic gravitational collapse, the larger amount of surface extension in the eastern French Massif Central was compensated in the orogenic root by flow of the partially molten rocks from the northern and western French Massif Central toward the Velay and Montagne Noire dome in the southeastern part of the French Massif Central. Moreover, the presence of migmatites in the Montagne Noire and in the Pyrenees, beneath the foreland sedimentary sequence affected by recumbent folds and low-grade to HT/LP metamorphism, also suggests that the partially molten orogenic root has flown toward the foreland. This model provides a potential explanation for the enigmatic high geothermal gradient identified in the external zone and for propagation of crustal thickening in the foreland beneath supracrustal rocks. Our model is also compatible with the shallow northward plunge (0–30°) of the mineral and stretching lineation displayed along the Sillon Houiller (Grolier and Letourneur, 1968; Bonijoly and Castaing, 1984).

The P-T-t record of the lower crust documented by felsic and mafic granulitic xenoliths with a maximum pressure of 14 kbar and a temperature of 900 °C followed by near isothermal decompression (Leyreloup, 1974; Montel, 1985; Downes and Leyreloup, 1986) is consistent with a hot and thick orogenic crust that has then been affected by thinning. U-Pb ages on zircon from the granulites from 320 to 280 Ma cover the transition from crustal thickening to gravitational collapse (Downes and Leyreloup, 1986; Costa and Rey, 1995). The geochemical characteristics of the granulites suggest that they represent a mixture of (i) resisters (rocks that were not prone to melt), (ii) solid residues after melt extraction, (iii) cumulates (iv) and residual melts (Dupuy et al., 1979; Pin and Vielzeuf, 1983; Vielzeuf et al., 1990; Downes et al., 1997), which is complementary to the evidences for crustal-derived magmas emplaced at higher structural levels. Such features are also identified in exposed sections of the Variscan lower crust in the Ivrea Zone (Percival, 1992; Barboza et al., 1999; Bea and Montero, 1999; Schaltegger and Gebauer, 1999; Guergouz et al., 2018) or in Calabria (Schenk, 1980, 1981, 1989).

thumbnail Fig. 15

Late Carboniferous-Permian magmatism. Magmatic rocks with ages ranging from c. 305 to c. 295 Ma are localized along the Sillon Houiller and are present as plutons and in the core of large migmatite domes (Velay, Montagne Noire) in the eastern part of the French Massif Central. The presence of plutonic rocks to the south of the French Massif Central suggests that the partially molten source has migrated since the Lower Carboniferous from the internal to external zone of the Variscan belt. The combination of MPG-type, CPG-type and KCG-type granites emplaced during this period is consistent with the contribution of crustal and mantle sources and is interpreted as reflecting the impact of southward slab retreat that is more pronounced beneath the eastern part of the French Massif Central.

4.7 Pertinence of the proposed geodynamic-tectonic model compared to physical modelling of the dynamics of orogenic belts

In this section, we assess the significance of the geological record of rocks forming the Variscan belt exposed in the French Massif Central in terms of the thermal-mechanical evolution of orogenic belts as investigated by physical modelling.

Eclogites, which are the oldest metamorphic rocks identified in the French Massif Central, and are present at the highest structural level of the nappe pile, have recorded HP/LT conditions. Such conditions imply burial at more than 30 km depth and more rapidely than the effect of thermal diffusion, typically a subduction rate of more than 1 cm/yr (Huerta et al., 1996; Henry et al., 1997).

modelling of the thermal evolution of orogenic belts has shown that a time delay of 20 to 30 Myrs after crustal thickening is required for thermal diffusion and heat production through natural decay of radioactive isotopes to lead to significant partial melting of the orogenic root (England and Thompson, 1984; Thompson and Connolly, 1995). This 20 to 30 Myrs gap is consistent with the geochronological record of HT eclogites and granulites preserved in the UGU. In turn, the thermal impact of removal of the lithospheric mantle root is more rapid and dramatic in terms of increase in temperature in the crust (Houseman et al., 1981; Arnold et al., 2001). Thinning of the lithospheric mantle is so efficient that it has been proposed as a mechanism to account for HT/LP metamorphism in back-arc basins marked by a relatively thin crust (Collins, 2002). In the case of a convergent plate boundary marked by slab retreat and tectonic accretion, both the radioactive heat production and the mantle heat flux are increased concomitantly, which is the most favorable scenario for a hot orogenic root (Sandiford and Powell, 1990; Vanderhaeghe and Duchêne, 2010; Arnold et al., 2001). In such a case, thermal relaxation after removal of the lithospheric mantle occurs over about 100 Myrs (Ueda et al., 2012), which is roughly the duration of HT metamorphism and magmatism invoked in the geodynamic-tectonic model presented in this paper.

Exhumation of rocks entrained in subduction entails mechanical decoupling of these rocks from the downgoing slab, which in turns indicates that their buoyancy reached their mechanical strength (Chemenda et al., 1996; Escher and Beaumont, 1997; Warren et al., 2008). In turn, partial melting might be particularly efficient in decreasing the strength of buried rocks and thus favouring their decoupling from the downgoing slab and their exhumation.

The presence of a low-viscosity orogenic root is causing horizontal flow of the thickened crust (Artyushkov, 1973; Molnar and Lyon-Caen, 1988; Bird, 1991; Royden, 1996). This horizontal flow, driven by the gravity force associated with lateral variations of the weight of the crustal column, occurs preferentially toward a free boundary or a mechanically weak zone. In the case of an advancing plate boundary, i.e. indentation, horizontal flow occurs preferentially laterally and is associated with the activation of strike-slip shear zones (Royden, 1997; Cagnard et al., 2006). The low-viscosity layer might also flow in the direction of convergence, toward the foreland (Henk, 2000; Vanderhaeghe et al., 2003). The presence of a low-viscosity layer also impedes the maintenance of an irregular topographic surface (Artyushkov, 1973), which leads to the formation of an orogenic plateau (Molnar et al., 1993; Vanderhaeghe et al., 2003). If the zone of weak crust is maintained along its boundaries by stronger crustal sections, the orogenic plateau might be maintained (Vanderhaeghe et al., 2003; Cook and Royden, 2008). In the contrary, horizontal flow of the low-viscosity orogenic root will lead to redistribution of the orogenic crust until total decay of the gravity force (Rey et al., 2001). The style of extension is controlled by the rheology of the crust and thus its temperature (Buck, 1991; Rey et al., 2009).

Given its low density and low viscosity, a partially molten orogenic root is susceptible to develop diapiric Rayleigh-Taylor instabilities (Ramberg, 1968; Talbot, 1979; Perchuk et al., 1992; Weinberg and Schmeling, 1992; Cruden et al., 1995) and even convective instabilities (Weinberg, 1997; Vanderhaeghe et al., 2018).

5 Conclusion

The structural, petrological, geochemical and geochronological record of the French Massif Central provides an archive of the thermal-mechanical evolution of the Variscan belt in Western Europe and documents the impact of partial melting and magmatism during orogenic evolution, as summarized in the following and in Figures 18 and 19.

The pre-Variscan paleogeography (Figs. 16 and 17) is marked by hyper-extension of the northern margin of the Gondwana supercontinent, evidenced by the coeval emplacement of alkaline granitoids and of tholeiitic to calc-alkaline bimodal mafic-felsic magmas of the Leptynite Amphibolite Complex during Ordovician times (485–460 Ma). This setting is particularly favorable for trapping of voluminous detrital sediments and the emplacement of alkaline magmas that might represent the protoliths of the main portion of the metagreywackes, metapelites and orthogneisses forming the UGU, MAU, LGU and PAU.

The presence of HT eclogite facies mafic and ultramafic enclaves of the Leptynite Amphibolite Complex into the migmatitic HP granulites of the UGU indicates subduction of the immature oceanic crust (Fig. 18A) together with continental ribbons (Fig. 18B). The 20–30 Ma difference in age between the eclogite-facies metamorphism recorded by the LAC corresponds to the time required for thermal relaxation and partial melting. Partial melting potentially triggered decoupling of the UGU from the downgoing slab allowing for the syn-orogenic exhumation of the partially molten UGU entraining pieces of mafic and ultramafic rocks forming the LAC and representing previously subducted oceanic lithosphere and/or part of the suprasubduction mantle. The percolation into the suprasubduction mantle of felsic melt segregated from the partially molten UGU would have possibly contributed to its enrichment in HFSE, REE and LILE required to form the source of the later vaugnerites.

During the Devonian-Carboniferous transition, opening of a rift in the internal zone of the Variscan belt is accompanied by the emplacement of the c. 360 Ma ACG and low-K calc-alkaline to tholeiitic magmas of the Brévenne Unit and Limousin tonalitic line. This is coeval with thrust propagation in the external zone and deposition of detrital sediments in the foreland that we tentatively attribute to the outward propagation of the orogenic belt in a context of slab retreat (Fig. 18C). Such a context might have facilitated syn-orogenic exhumation of the buoyant and low-viscosity partially molten felsic units previously dragged with the subducting slab.

Regional scale transcurrent shear zones and foreland propagation of thrusts accommodate the development of an orogenic plateau by gravity-driven lateral flow of the partially molten orogenic root, owing to a major Carboniferous thermal anomaly of lithospheric extent. This is indeed associated with high-temperature/medium-pressure metamorphism and the emplacement of syntectonic plutons from 345 to 310 Ma (Fig. 18D). The diversity of the geochemical signatures of the magmatic rocks encompassing MPG, CPG and their volcanic equivalents, KCG and Mg-K vaugnerites indicates a contribution of both the crust and the mantle. The latter along with the southward younger emplacement ages of these magmatic rocks is consistent with progressive retreat toward the south of the northward plunging slab during the Carboniferous.

The late Carboniferous gravitational collapse of the Variscan orogenic crust (Fig. 18E) is accommodated by extension of the upper crust and by lateral flow and exhumation of the partially molten root. This is marked by the formation of crustal-scale domes cored by LP migmatites and heterogeneous granites coeval with the emplacement of syntectonic laccoliths in the footwall of low-angle detachments. These are complementary to a refractory granulitic lower crust formed by protracted high-temperature metamorphism, partial melting and melt/solid segregation.

According to this new model consistent with physical modelling, continuous and protracted presence of melt in the root of the orogenic crust plays a crucial role in the tectonic evolution of the Variscan belt by (i) triggering syn-orogenic exhumation of subducted continental units decoupled from the downgoing slab; (ii) controlling the formation and lateral development of an orogenic plateau; and finally, (iii) guiding the formation of metamorphic core complexes during orogenic gravitational collapse. Crustal melting starts with segregation of melts from the subducted oceanic and continental units in the Devonian. Subsequently during the Carboniferous, the emplacement of plutons and volcanics during the building of the orogenic plateau has a contribution from the partially molten crust and from the mantle. The maintenance of a partially molten crust for several tens of Myrs is probably favored by the combined effects of radioactive heat production and increasing mantle heat flux owing to removal of the lithospheric mantle slab. It ended with the extraction of differentiated magmas and crystallization of the collapsed partially molten orogenic root. The contrasting Carboniferous geological record between the Western and Eastern French Massif Central separated by the Sillon Houiller is consistent with a more pronounced slab retreat in the East toward the Paleotethys free boundary. The Eastern part of the French Massif Central is indeed characterized by (i) the abundance of Mg-K diorites (vaugnerites) and KCG-type granites indicating the contribution of mantle-derived magmas, and (ii) a widespread extension associated with the development of the Velay and Montagne Noire migmatite domes (Fig. 19).

thumbnail Fig. 16

Pre-Variscan geodynamic configuration. A. During the Cambrian, the Gondwana margin is marked by the emplacement of calc-alkaline magmas attributed to an enriched mantle source above a subducting slab. The size of the Iapetus is about 2000–3000 km wide. B. During the Ordovician, the Avalonia and Armorica continental ribbons are separated from the Gondwana margin. The Rheic Ocean corresponds to the future Rheic suture exposed in southern Great Britain. The Medio-European Ocean corresponds to the Leptynite-Amphibolite Complex forming boudins and enclaves in high-grade nappes of the Moldanubian allochthonous terrane (see Fig. 3). The tholeiitic to calc-alkaline signature of the LAC is interpreted as reflecting an emplacement of the magmatic protoliths in a back-arc or immature oceanic setting. Alkaline magmas intrusive in Ediacaran sedimentary sequences correspond to orthogneisses preserved in the LGU and PAU and are attributed to opening of a series of rifts leading to hyperextension of the Gondwana margin. These features are consistent with retreat of the southward plunging Iapetus slab that will eventually lead to tectonic accretion of Avalonia to the Laurussia craton.

thumbnail Fig. 17

Plate scale geodynamic reconstruction for the Paleozoic (modified after Domeier, 2016; Domeier and Torsvik, 2014; Matte, 2001). A. Cambrian 500 Ma (541–485 Ma). The Gondwana margin is at the South Pole, Laurentia is at the Equator and the Iapetus Ocean is at least 3000 km wide. The Gondwana margin is an active plate boundary marked by subduction of the Iapetus. B. Ordovician 470 Ma (485–444 Ma). The margin of Gondwana is marked by hyperextension resulting in the separation of Avalonia and Armorica and the opening of the Rheic and Medio-European Oceans. C. Silurian 430 Ma (444–419 Ma). The Iapetus Ocean has closed and Avalonia has been tectonically accreted to Laurussia. The Rheic Ocean started to subduct beneath the margin of Laurussia and Medio-European Ocean, between Armorica and Gondwana is at its maximum width. D. Devonian 380 Ma (419–359 Ma). Armorica is bounded by subduction zones with opposite vergence resulting in the formation of the Variscan belt. E. Carboniferous 330 Ma (359–299 Ma). The Variscan orogenic front progresses from the hinterland to the foreland in association with slab retreat.

thumbnail Fig. 18

Geodynamic-tectonic model of the Variscan bel exposed in the French Massif Central. A. Silurian (c. 420 Ma) subduction of the Medio-European Ocean. B. Middle Devonian (c. 385 Ma) subduction and partial melting of the hyper-extended northern continental margin of the Gondwana craton. C. Late Devonian (c. 365 Ma) slab retreat, opening of the Brévenne rift and syn-orogenic exhumation of high-grade partially molten nappes. D. Lower Carboniferous (c. 330 Ma) development of an orogenic plateau by lateral flow of partially molten orogenic root in a context of plate convergence and southward slab retreat. E. Late Carboniferous (c. 300 Ma) gravitational collapse of the Variscan belt in a context of southward slab delamination accommodated by lateral and upward flow of the partially molten orogenic root concomitant with brittle extension of the upper crust.

thumbnail Fig. 19

3D model for the late Carboniferous crustal and lithospheric scale structure of the Variscan belt beneath the French Massif Central. Gravitational collapse is accommodated by (i) lateral flow of the partially molten orogenic root from the western to the eastern side of the French Massif Central and from the internal to the external zone toward the south, and (ii) upward flow to form migmatite domes within Metamorphic Core Complexes. The Sillon Houiller accommodates differential slab retreat between the eastern and western parts of the French Massif Central.


This paper results from a long gestation process. First of all, Olivier Vanderhaeghe is very grateful to Patrick Ledru and Jean-Michel Caron for their initiation to the geology of the Massif Central Français and also to Jean-Pierre Burg for his encouragements in developing the structural analysis of migmatites. The authors address a special thank to Philippe Rossi, director of the French geological mapping program, and the BRGM for financial support at different stages. Arnaud Villaros acknowledges financial support from the Labex VOLTAIRE (ANR-10-LABX-100-01) and Jean-François Moyen, Simon Couzinié and Arnaud Villaros were supported by credits from the INSU 2016 Syster program.

A first version of this manuscript, submitted to an IJES special volume following the Variscan 2015 conference, has benefited from the stimulating viewpoint of Roberto Weinberg and from the critical opinion of an anonymous reviewer. As a result, the present version was considerably augmented and was itself significantly improved owing to the constructive review of Jacques Malavieille and from the challenging comments of an anonymous reviewer not sharing the same vision of this segment of the Variscan belt. The authors thank Olivier Averbuch and Laurent Jolivet for their careful and balanced editorial handing.


  • Aerden DGAM. 1998. Tectonic evolution of the Montagne Noire and a possible orogenic model for syncollisional exhumation of deep rocks, Variscan belt, France. Tectonics 17: 62–79. [CrossRef] [Google Scholar]
  • Aerden DG, Malavieille J. 1999. Origin of a large-scale fold nappe in the Montagne Noire Variscan belt, France. J. Struct. Geol. 21: 1321–1333. [CrossRef] [Google Scholar]
  • Alexandre P. 2007. U-Pb zircon SIMS ages from the French Massif Central and implication for the pre-Variscan tectonic evolution in Western Europe. Comptes Rendus Geosci. 339: 613–621. [Google Scholar]
  • Alexandrov P, Cheilletz A, Deloule É, Cuney M. 2000. 319 ± 7 Ma crystallization age for the Blond granite (northwest Limousin, French Massif Central) obtained by U/Pb ion-probe dating of zircons. Comptes Rendus Académie Sci. – Ser. IIA-Earth Planet. Sci. 330: 617–622. [Google Scholar]
  • Alexandrov P, Floc’h J-P, Cuney M, Cheilletz A., 2001. Datation U-Pb à la microsonde ionique des zircons de l’unité supérieure de Gneiss dans le Sud Limousin, Massif Central. Comptes Rendus Académie Sci. – Ser. IIA-Earth Planet. Sci. 332: 625–632. [Google Scholar]
  • Altherr R, Kalt A. 1996. Metamorphic evolution of ultrahigh-pressure garnet peridotites from the Variscan Vosges Mts. (France). Chem. Geol. 134: 27–47. [Google Scholar]
  • Annen C, Sparks RSJ. 2002. Effects of repetitive emplacement of basaltic intrusions on thermal evolution and melt generation in the crust. Earth Planet. Sci. Lett. 203: 937–955. [CrossRef] [Google Scholar]
  • Arnaud F, Burg J-P. 1993. Microstructures des mylonites schisteuses : cartographie des chevauchements varisques dans les Cévennes et détermination de leur cinématique. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 317: 1441–1447. [Google Scholar]
  • Arnaud F, Boullier AM, Burg J-P. 2004. Shear structures and microstructures in micaschists: the Variscan Cévennes duplex (French Massif Central). J. Struct. Geol. 26: 855–868. [CrossRef] [Google Scholar]
  • Arnold J, Jacoby WR, Schmeling H, Schott B. 2001. Continental collision and the dynamic and thermal evolution of the Variscan orogenic crustal root – Numerical models. J. Geodyn. 31: 273–291. [CrossRef] [Google Scholar]
  • Arthaud F, Matte P. 1975. Les décrochements tardi-hercyniens du sud-ouest de l’Europe. Géométrie et essai de reconstitution des conditions de la déformation. Tectonophysics 25: 139–171. [CrossRef] [Google Scholar]
  • Arthaud F, Mattauer M, Proust F. 1966. La structure et la microtectonique des nappes hercyniennes de la Montagne Noire. In: Colloque « Étages Tectoniques », À La Baconnière, NeuchâTel, 247 p. [Google Scholar]
  • Artyushkov EV. 1973. Stresses in the lithosphere caused by crustal thickness inhomogeneities. J. Geophys. Res. 78: 7675–7708. [CrossRef] [Google Scholar]
  • Audren C, Feybesse JL, Tegyey M, Triboulet C. 1987. Relations entre déformations et cristallisations et chemins « PTtd » des micaschistes polyphasés d’Echassières. Modèle d’évolution géodynamique. Géol Fr. 2: 43–45. [Google Scholar]
  • Autran A, Cogné J. 1980. La zone interne de l’orogène varisque dans l’Ouest de la France et sa place dans le développement de la chaîne hercynienne. Mém. BRGM 191–202. [Google Scholar]
  • Autran A, Guillot PL. 1978. L’évolution métamorphique du Limousin (Massif Central français) au Paléozoïque. Relation entre les cycles Calédonien et Varisque. In: La Chaîne Varisque d’Europe Moyenne et Occidentale, pp. 211–226. [Google Scholar]
  • Averbuch O, Piromallo C. 2012. Is there a remnant Variscan subducted slab in the mantle beneath the Paris basin? Implications for the late Variscan lithospheric delamination process and the Paris basin formation. Tectonophysics 558: 70–83. [CrossRef] [Google Scholar]
  • Ballèvre M, Balé P. 1992. Forage scientifique de Sancerre-Couy: tectonique et métamorphisme. Géologie Fr. 135–138. [Google Scholar]
  • Ballèvre M, Kienast JR, Paquette J-L. 1987. Le métamorphisme éclogitique dans la nappe hercynienne de Champtoceaux (Massif Armoricain). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 305: 127–131. [Google Scholar]
  • Ballèvre M, Capdevila R, Guerrot C, Peucat J-J. 2002. Discovery of an alkaline orthogneiss in the eclogite-bearing Cellier Unit (Champtoceaux Complex, Armorican Massif): a new witness of the Ordovician rifting. Comptes Rendus Geosci. 334: 303–311. [CrossRef] [Google Scholar]
  • Ballèvre M, Bosse V, Ducassou C, Pitra P. 2009. Palaeozoic history of the Armorican Massif: models for the tectonic evolution of the suture zones. Comptes Rendus Geosci. 341: 174–201. [CrossRef] [Google Scholar]
  • 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]
  • Barbarin B. 1999. A review of the relationships between granitoid types, their origins and their geodynamic environments. Lithos 46: 605–626. [CrossRef] [Google Scholar]
  • Barbarin B, Belin JM. 1982. Mise en évidence du cisaillement ductile hercynien « St Gervais-L’Hermitage » (Massif central français). CR Acad Sci Paris II 294: 1377–1380. [Google Scholar]
  • Barbey P, Marignac C, Montel JM, Macaudiere J, Gasquet D, Jabbori J. 1999. Cordierite growth textures and the conditions of genesis and emplacement of crustal granitic magmas: the Velay granite complex (Massif Central, France). J. Petrol. 40: 1425–1441. [CrossRef] [Google Scholar]
  • Barbey P, Cheilletz A, Laumonier B. 2001. The Canigou orthogneisses (Eastern Pyrenees, France, Spain): an Early Ordovician rapakivi granite laccolith and its contact aureole. Comptes Rendus Académie Sci. – Ser. II A-Earth Planet. Sci. 332: 129–136. [Google Scholar]
  • Barbey P, Villaros A, Marignac C, Montel J-M. 2015. Multiphase melting, magma emplacement and PT-time path in late-collisional context: the Velay example (Massif Central, France). Bull. Soc. Geol. Fr. 186: 93–116. [CrossRef] [Google Scholar]
  • Barboza SA, Bergantz GW, Brown M. 1999. Regional granulite facies metamorphism in the Ivrea zone: Is the Mafic Complex the smoking gun or a red herring? Geology 27: 447–450. [CrossRef] [Google Scholar]
  • Bard JP, Burg JP, Matte P, Ribeiro A. 1980. La chaine hercynienne d’Europe occidentale en termes de tectonique des plaques. Géologie Eur. 108: 233–46. [Google Scholar]
  • Batias P, Duthou JL. 1979. Age Viséen supérieur du granite porphyroïde de Vienne-Tournon (Massif Central français). In: Proc. 7e La Réun. Annu. Sci. Terre Lyon. [Google Scholar]
  • Bé Mézémé E. 2005. Contribution de la géochronologie U-Th-Pb sur monazite à la compréhension de la fusion crustale dans la chaîne Varisque française et implication géodynamique. Orléans. [Google Scholar]
  • Bé Mézémé E, Cocherie A, Faure M, Legendre O, Rossi P. 2006. Electron microprobe monazite geochronology of magmatic events: examples from Variscan migmatites and granitoids, Massif Central, France. Lithos 87: 276–288. [CrossRef] [Google Scholar]
  • Bea F, Montero P. 1999. Behavior of accessory phases and redistribution of Zr, REE, Y, Th, and U during metamorphism and partial melting of metapelites in the lower crust: an example from the Kinzigite Formation of Ivrea-Verbano, NW Italy. Geochim. Cosmochim. Acta 63: 1133–1153. [CrossRef] [Google Scholar]
  • Bébien J. 1971. Eléments nouveaux sur le volcanisme dévono-dinantien de l’extrémité sud-ouest du faisceau synclinal du Morvan. Cr Seances Acad Sci Paris 273: 466–8. [Google Scholar]
  • Becq-Giraudon J-F. 1993. Problèmes de la biostratigraphie dans le Paléozoïque supérieur continental (Stéphanien-Autunien) du Massif Central. Geodin. Acta 6: 219–224. [CrossRef] [Google Scholar]
  • Becq-Giraudon J-F, Montenat C, Van Den Driessche J. 1996. Hercynian high-altitude phenomena in the French Massif Central: tectonic implications. Palaeogeogr. Palaeoclimatol. Palaeoecol. 122: 227–241. [CrossRef] [Google Scholar]
  • Bellot J-P. 2001. La structure de la croûte varisque du Sud-Limousin (Massif central français) et ses relations avec les minéralisations aurifères tardi-orogéniques: Apport des données géologiques, gitologiques, géophysiques et de la modélisation 3D. Université de Montpellier II. [Google Scholar]
  • Bellot J-P. 2007. Extensional deformation assisted by mineralised fluids within the brittle-ductile transition: Insights from the southwestern Massif Central, France. J. Struct. Geol. 29: 225–240. [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]
  • Berger J, Femenias O, Mercier JCC, Demaiffe D. 2006. A Variscan slow-spreading ridge (MOR-LHOT) in Limousin (French Massif Central): magmatic evolution and tectonic setting inferred from mineral chemistry. Mineral. Mag. 70, 175–185. [CrossRef] [Google Scholar]
  • Berger J, Féménias O, Ohnenstetter D, Bruguier O, Plissart G, Mercier J-CC, et al. 2010a. New occurrence of UHP eclogites in Limousin (French Massif Central): age, tectonic setting and fluid-rock interactions. Lithos 118: 365–382. [CrossRef] [Google Scholar]
  • Berger J, Femenias O, Ohnenstetter D, Plissart G, Mercier J-C. 2010b. Origin and tectonic significance of corundum-kyanite-sapphirine amphibolites from the Variscan French Massif Central. J. Metamorph. Geol. 28: 341–360. [CrossRef] [Google Scholar]
  • Bernard-Griffiths J. 1975. Essai sur la signification des âges au strontium dans une série métamorphique: le Bas Limousin, Massif-Central français. Université de Clermont, Unité d’enseignement et de recherche de sciences. [Google Scholar]
  • Bernard-Griffiths J, Cornichet J. 1985. Origin of eclogites from South Brittany, France: a Sm/Nd isotopic and REE study. Chem. Geol. Isot. Geosci. Sect. 52: 185–201. [CrossRef] [Google Scholar]
  • Bernard-Griffiths J, Cantagrel J-M, Duthou J-L. 1977. Radiometric evidence for an Acadian tectonometamorphic event in western Massif Central Francais. Contrib. Mineral. Petrol. 61: 199–212. [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. Société Géologique Fr. 1: 523–529. [CrossRef] [Google Scholar]
  • Bertaux J, Becq-Giraudon J-F, Jacquemin H. 1993. Les bassins anthracifères de la région de Roanne (Loire, Massif central) marqueurs d’une tectonique active durant le Viséen. Géologie Fr. 4: 3–10. [Google Scholar]
  • Berthier F, Duthou JL, Roques M. 1979. Datation géochronologique Rb/Sr sur roches totales du granite de Guéret (Massif Central). Age fini-Dévonien de mise en place de l’un de ses facies types. Bull BRGM 1: 59–72. [Google Scholar]
  • Bertrand M. 1887. La chaine des Alpes et la formation du continent européen. Bull. Soc. Geol. Fr. 3: 440–442. [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éologie Fr. 1: 167–189. [Google Scholar]
  • Bird P. 1991. Lateral extrusion of lower crust from under high topography in the isostatic limit. J. Geophys. Res. 96: 10275. [CrossRef] [Google Scholar]
  • Bitri A, Truffert C, Bellot J-P, Bouchot V, Ledru P, Milési J-P, et al. 1999. Imagery of crustal-scale As-Au-Sb hydrothermal palaeofields in the Variscan belt: vertical seismic reflection (GeoFrance 3D: French Massif Central). Comptes Rendus Acad. Sci. Ser. II – A Earth Planet. Sci. 329: 771–777. [Google Scholar]
  • Bitri A, Brun J-P, Gapais D, Cagnard F, Gumiaux C, Chantraine J, et al. 2010. Deep reflection seismic imaging of the internal zone of the South Armorican Hercynian belt (western France)(ARMOR 2/Géofrance 3D Program). Comptes Rendus Geosci. 342: 448–452. [CrossRef] [Google Scholar]
  • Bodinier J-L, Burg J-P. 1981. Evolution métamorphique et tectonique des séries cristallophylliennes du Rouergue occidental : mise en évidence d’un chevauchement dans la région de Najac (Aveyron). Bull BRGM I: 315–339. [Google Scholar]
  • Bodinier J-L, Giraud A, Dupuy C, Leyreloup A, Dostal J. 1986. Caracterisation geochimique des metabasites associees a la suture meridionale hercynienne ; Massif Central francais et Chamrousse (Alpes). Bull. Société Géologique Fr. 2: 115–123. [CrossRef] [Google Scholar]
  • Boehnke P, Watson EB, Trail D, Harrison TM, Schmitt AK. 2013. Zircon saturation re-revisited. Chem. Geol. 351: 324–334. [CrossRef] [Google Scholar]
  • Bonijoly D, Castaing CH. 1984. Fracturation et genèse des bassins stéphaniens du Massif central français en régime compressif. Ann Soc Geol Nord 103: 187–199. [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]
  • Bouchardon J-L, Santallier D, Briand B, Ménot R-P, Piboule M. 1989. Eclogites in the French Palaeozoic Orogen: geodynamic significance. Tectonophysics 169: 317–332. [CrossRef] [Google Scholar]
  • Bouilhol P, Leyreloup AF, Delor C, Vauchez A, Monié P. 2006. Relationships between lower and upper crust tectonic during doming: the mylonitic southern edge of the Velay metamorphic core complex (Cévennes-French Massif Central). Geodin. Acta 19: 137–153. [CrossRef] [Google Scholar]
  • Boutin R, Montigny R. 1993. Datation 40Ar-39Ar des amphibolites du complexes leptynoamphibolique du Plateau d’Aigurande: collision varisque à 390 Ma dans le Nord-Ouest du Massif Central français. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 316: 1391–1398. [Google Scholar]
  • Briand B. 1978. Métamorphisme inverse et chevauchement de type himalayen dans la série de la vallée du Lot. CR Acad Sci Paris 286: 729–731. [Google Scholar]
  • Briand B, Gay M. 1978. La série inverse de Saint-Geniez-d’Olt : évolution métamorphique et structurale. Bull. BRGM 3: 167–186. [Google Scholar]
  • Briand B, Piboule M. 1979. Les métabasites de la série de Marvejols (Massif central) : témoins d’un magmatisme tholéïtique d’arrière-arc cambro-ordovicien ? Bull BRGM 2(I): 131–171. [Google Scholar]
  • Briand B, Piboule M, Bouchardon JL. 1988. Diversite geochimique des metabasites des groupes leptyno-amphiboliques du Rouergue et de Marvejols (Massif Central) ; origine et implications. Bull. Société Géologique Fr. 4: 489–498. [CrossRef] [Google Scholar]
  • Briand B, Bouchardon J-L, Ouali H, Piboule M, Capiez P. 1995. Geochemistry of bimodal amphibolitic – Felsic gneiss complexes from eastern Massif Central, France. Geol. Mag. 132: 321–337. [CrossRef] [Google Scholar]
  • Briand B, Duthou J-L, Guerrot C, Chenevoy M. 2002. Les « granites à tablettes d’orthose » du Vivarais, témoins d’un magmatisme post-épaississement d’âge Dinantien inférieur ; identification d’une unité géologique Nord-Ouest-Vivarais. Comptes Rendus Geosci. 334: 741–747. [CrossRef] [Google Scholar]
  • Brichau S, Respaut J-P, Monié P. 2008. New age constraints on emplacement of the Cévenol granitoids, South French Massif Central. Int. J. Earth Sci. 97: 725–738. [CrossRef] [Google Scholar]
  • Brown M. 2001. Orogeny, migmatites and leucogranites: a review. J. Earth Syst. Sci. 110: 313–336. [CrossRef] [Google Scholar]
  • Brown M, Rushmer T. 1997. The role of deformation in the movement of granitic melt: views from the laboratory and the field. Deform. – Enhanc. Fluid Transp. Earth’s Crust Mantle 8: 111–144. [Google Scholar]
  • Brown M, Solar GS. 1998. Granite ascent and emplacement during contractional deformation in convergent orogens. J. Struct. Geol. 20: 1365–1393. [CrossRef] [Google Scholar]
  • Brueckne HK, Medaris LG. 2000. A general model for the intrusion and evolution of “mantle” garnet peridotites in high-pressure and ultra-high-pressure metamorphic terranes. J. Metamorph. Geol. 18: 123–133. [CrossRef] [Google Scholar]
  • Bruguier O, Becq-Giraudon JF, Bosch D, Lancelot JR. 1998. Late Visean hidden basins in the internal zones of the Variscan belt: U-Pb zircon evidence from the French Massif Central. Geology 26: 627–630. [CrossRef] [Google Scholar]
  • Bruguier O, Becq-Giraudon JF, Champenois M, Deloule E, Ludden J, Mangin D. 2003. Application of in situ zircon geochronology and accessory phase chemistry to constraining basin development during post-collisional extension: a case study from the French Massif Central. Chem. Geol. 201: 319–336. [CrossRef] [Google Scholar]
  • Buck WR. 1991. Modes of continental lithospheric extension. J. Geophys. Res. Solid Earth 96: 20161–20178. [CrossRef] [Google Scholar]
  • Burg J-P, Leyreloup A. 1989. Métamorphisme granulitique de roches granitiques en Rouergue (Massif Central). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 309: 719–725. [Google Scholar]
  • Burg JP, Matte PJ. 1978. A cross section through the French Massif Central and the scope of its Variscan geodynamic evolution. Z. Dtsch. Geol. Ges. 429–460. [Google Scholar]
  • Burg J-P, Vanderhaeghe O. 1993. Structures and way-up criteria in migmatites, with application to the Velay dome (French Massif Central). J. Struct. Geol. 15: 1293–1301. [CrossRef] [Google Scholar]
  • Burg JP, Matte P, Leyreloup A, Marchand J. 1984. Inverted metamorphic zonation and large-scale thrusting in the Variscan Belt: an example in the French Massif Central. Geol. Soc. Lond. Spec. Publ. 14: 47–61. [CrossRef] [Google Scholar]
  • Burg J-P, Delor C, Leyreloup A. 1986. Le massif du Lévézou et les séries adjacentes du Rouergue oriental. Nouvelles données pétrographiques et structurales. Géologie Fr. 3: 229–272. [Google Scholar]
  • Burg J-P, Bale P, Brun J-P, Girardeau J. 1987. Stretching lineation and transport direction in the Ibero-Armorican arc during the siluro-devonian collision. Geodin. Acta 1: 71–87. [CrossRef] [Google Scholar]
  • Burg J-P, Castaing C, Chantraine J, Hottin A-M, Kienast J-R, Mégnien C, et al. 1989a. Les formations métamorphiques traversées par le sondage de Sancere-Couy (programme GPF). Nouveau jalon de la chaîne varisqe. Comptes Rendus Acad. Sci. 2: 1819–1824. [Google Scholar]
  • Burg JP, Leyreloup AF, Romney F, Delor CP. 1989b. Inverted metamorphic zonation and Variscan thrust tectonics in the Rouergue area (Massif Central, France): PTt record from mineral to regional scale. Geol. Soc. Lond. Spec. Publ. 43: 423–439. [CrossRef] [Google Scholar]
  • Burg J-P, Brun J-P, Van Den Driessche J. 1990. Le sillon houiller du Massif Central français: faille de transfert pendant l’amincissement crustal de la chaîne. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 311: 147–152. [Google Scholar]
  • Burg J-P, Van den Driessche J, Brun J-P. 1993. Syn-to post-thickening extension in the Variscan Belt of Western Europe: modes and structural consequences. Géologie Fr. 3: 33–51. [Google Scholar]
  • Bussien D, Bussy F, Masson H, Magna T, Rodionov N. 2008. Variscan lamprophyres in the Lower Penninic domain (Central Alps): age and tectonic significance. Bull. Société Géologique Fr. 179: 369–381. [CrossRef] [Google Scholar]
  • Cabanis B, Guillot PL, Santallier D, Jaffrezic H, Meyer G, Treuil M. 1983. Apport des elements-traces a l’etude geochimique des metabasites du Bas Limousin. Bull. Société Géologique Fr. 7: 563–574. [CrossRef] [Google Scholar]
  • Caen-Vachette M. 1979. Age Cambrien des rhyolites transformées en leptynites dans la série métamorphique du Pilat (Massif Central français). C R Hebd Séanc Acad Sci Paris 289: 997–1000. [Google Scholar]
  • Caen-Vachette M, Couturié JP, Didier J. 1981. Age westphalien du granite de Rocles (Cévennes, massif Central Français). Comptes Rendus L’Académie Sci. 293: 957–960. [Google Scholar]
  • Caen-Vachette M, Gay M, Peterlongo J-M, Pitiot P, Vitel G. 1984. Age radiométrique du granite syntectonique du gouffre d’Enfer et du métamorphisme hercynien dans la série de basse pression du Pilat (Massif Central Français). Comptes-Rendus Séances Académie Sci. Sér. 2 Mécanique-Phys. Chim. Sci. Univers Sci. Terre 299: 1201–1204. [Google Scholar]
  • Cagnard F, Durrieu N, Gapais D, Brun J-P, Ehlers C. 2006. Crustal thickening and lateral flow during compression of hot lithospheres, with particular reference to Precambrian times. Terra Nova 18: 72–78. [CrossRef] [Google Scholar]
  • Campbell IH, Stepanov AS, Liang H-Y, Allen CM, Norman MD, Zhang Y-Q, et al. 2014. The origin of shoshonites: new insights from the Tertiary high-potassium intrusions of eastern Tibet. Contrib. Mineral. Petrol. 167. [CrossRef] [Google Scholar]
  • Carlier De Veslud CL, Alexandre P, Cuney M, Ruffet G, Cheilletz A, Virlogeux D. 2004. Thermochronologie 40Ar/39Ar et evolution thermique des granitoides meso-varisques du complexe plutonique de Charroux-Civray (Seuil de Poitou)(4 fig.). Bull. Soc. Geol. Fr. 175: 147–156. [CrossRef] [Google Scholar]
  • Caron C. 1994. Les minéralisations Ph-Zn associées au paléozoique inférieur d’Europe méridionale. Traçage isotopique Pb-Pb des gites de l’iglesiente(SW Sardaigne) et des Cévennes et évolution du socle encaissant par la géochronologie U-Pb, 40Ar-39Ar et K-Ar. [Google Scholar]
  • Carrat HG, Zimmermann J-L. 1984. Ages K-Ar des roches volcaniques du Morvan. Comptes-Rendus Séances Académie Sci. Sér. 2 Mécanique-Phys. Chim. Sci. Univers Sci. Terre 299: 801–803. [Google Scholar]
  • Cartannaz C. 2006. Magmatismes et déformations polyphasés : exemple des massifs de Guéret et de millevaches (massif central français) : origine des magmas et contexte de mise en place. Besançon. [Google Scholar]
  • Cartannaz C, Rolin P, Cocherie A, Henry P, Rossy M. 2007a. Feuille Aubusson, 667 p. [Google Scholar]
  • Cartannaz C, Rolin P, Cocherie A, Marquer D, Legendre O, Fanning CM, et al. 2007b. 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. [Google Scholar]
  • Cartier C, Faure M, Lardeux H. 2001. The Hercynian orogeny in the South Armorican Massif (Saint‐Georges‐sur‐Loire Unit, Ligerian Domain, France): rifting and welding of continental stripes. Terra Nova 13: 143–149. [CrossRef] [Google Scholar]
  • Castiñeiras P, Navidad M, Liesa M, Carreras J, Casas JM. 2008. U-Pb zircon ages (SHRIMP) for Cadomian and Early Ordovician magmatism in the Eastern Pyrenees: new insights into the pre-Variscan evolution of the northern Gondwana margin. Tectonophysics 461: 228–239. [CrossRef] [Google Scholar]
  • Chantraine J, Autran A, Cavelier C, et al. 1996. Carte géologique de la France. [Google Scholar]
  • Chantraine J, Autran A, Cavelier C, Clozier L. 2003. Carte géologique de la France à l’échelle du millionième. Bureau de recherches géologiques et minières. [Google Scholar]
  • Chardon D, Gapais D, Cagnard F. 2009. Flow of ultra-hot orogens: a view from the Precambrian, clues for the Phanerozoic. Tectonophysics 477: 105–118. [CrossRef] [Google Scholar]
  • Charles N, Faure M, Chen Y. 2009. The Montagne Noire migmatitic dome emplacement (French Massif Central): new insights from petrofabric and AMS studies. J. Struct. Geol. 31: 1423–1440. [CrossRef] [Google Scholar]
  • Chauvet A, Volland-Tuduri N, Lerouge C, Bouchot V, Monié P, Charonnat X, et al. 2012. Geochronological and geochemical characterization of magmatic-hydrothermal events within the Southern Variscan external domain (Cévennes area, France). Int. J. Earth Sci. 101: 69–86. [CrossRef] [Google Scholar]
  • Chelle-Michou C, Laurent O, Moyen J-F, Block S, Paquette J-L, Couzinié S, et al. 2017. Pre-Cadomian to late-Variscan odyssey of the eastern Massif Central, France: formation of the West European crust in a nutshell. Gondwana Res. 46: 170–190. [Google Scholar]
  • Chemenda AI, Mattauer M, Bokun AN. 1996. Continental subduction and a mechanism for exhumation of high-pressure metamorphic rocks: new modelling and field data from Oman. Earth Planet. Sci. Lett. 143: 173–182. [CrossRef] [Google Scholar]
  • Chenevoy M, Ravier J. 1971. Caractères généraux des métamorphismes du Massif Central. In: Symposium J. Jung, Clermont-Ferrand, pp. 109–132. [Google Scholar]
  • Choukroune P, Gapais D, Matte P. 1983. Tectonique hercynienne et deformation cisaillante: La faille ductile senestre de la Marche (Massif Central Français). CR Acad Sci 296: 859–862. [Google Scholar]
  • Choulet F, Faure M, Fabbri O, Monié P. 2012. Relationships between magmatism and extension along the Autun-La Serre fault system in the Variscan Belt of the eastern French Massif Central. Int. J. Earth Sci. 101: 393–413. [CrossRef] [Google Scholar]
  • Clark AH, Scott DJ, Sandeman HA, Bromley AV, Farrar E. 1998. Siegenian generation of the Lizard ophiolite: U-Pb zircon age data for plagiogranite, Porthkerris, Cornwall. J. Geol. Soc. 155: 595–598. [CrossRef] [Google Scholar]
  • Cochelin B, Chardon D, Denèle Y, Gumiaux C, Le Bayon B. 2017. Vertical strain partitioning in hot Variscan crust: Syn-convergence escape of the Pyrenees in the Iberian-Armorican syntax. Bull. Société Géologique Fr. 188: 39. [CrossRef] [Google Scholar]
  • Cocherie A. 2007. Datations U-Pb (laser-ICPMS-MC) sur zircons et U-Th-Pb sur monazites de granitoïdes du Massif central (carte de Firminy). Rapport BRGM. [Google Scholar]
  • Cocherie A, Baudin T, Autran A, Guerrot C, Fanning CM, Laumonier B. 2005. U-Pb zircon (ID-TIMS and SHRIMP) evidence for the early ordovician intrusion of metagranites in the late Proterozoic Canaveilles Group of the Pyrenees and the Montagne Noire (France). Bull. Société Géologique Fr. 176: 269–282. [CrossRef] [Google Scholar]
  • Cocks LRM, Torsvik TH. 2002. Earth geography from 500 to 400 million years ago: a faunal and palaeomagnetic review. J. Geol. Soc. 159: 631–644. [Google Scholar]
  • Cocks LRM, Torsvik TH. 2006. European geography in a global context from the Vendian to the end of the Palaeozoic. Geol. Soc. Lond. Mem. 32: 83–95. [CrossRef] [Google Scholar]
  • Collins WJ. 2002. Hot orogens, tectonic switching, and creation of continental crust. Geology 30: 535.<0535:HOTSAC>2.0.CO;2. [CrossRef] [Google Scholar]
  • Cook KL, Royden LH. 2008. The role of crustal strength variations in shaping orogenic plateaus, with application to Tibet: Crustal strength and plateau deformation. J. Geophys. Res. Solid Earth 113. [Google Scholar]
  • Copard Y, Disnar J-R, Becq-Giraudon J-F, Boussafir M. 2000. Evidence and effects of fluid circulation on organic matter in intramontane coalfields (Massif Central, France). Int. J. Coal Geol. 44: 49–68. [CrossRef] [Google Scholar]
  • Costa S. 1989. Age radiométrique 39Ar–40Ar du métamorphisme des series du Lot et du charriage du groupe leptyno-amphibolique de Marvejols (MCF). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 309: 561–567. [Google Scholar]
  • Costa S. 1990. De la collision continentale à l’extension tardi-orogénique, 100 millions d’années d’histoire varisque dans le Massif central français‎. Une étude chronologique par la méthode 40Ar–39Ar. Université Montpellier 2 Sciences et Techniques du Languedoc. [Google Scholar]
  • Costa S. 1992. East-west diachronism of the collisional stageInthe french Massif Central Implications for the European Variscan Orogen. Geodin. Acta 5: 51–68. [CrossRef] [Google Scholar]
  • Costa S, Maluski H. 1988. Datations par la méthode 39Ar-40Ar de matériel magmatique et métamorphique paléozoïque provenant du forage de Couy-Sancerre (Cher, France). Programme GPF. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 306: 351–356. [Google Scholar]
  • Costa S, Rey P, 1995. Lower crustal rejuvenation and growth during post-thickening collapse: Insights from a crustal cross section through a Variscan metamorphic core complex. Geology 23: 905–908. [CrossRef] [Google Scholar]
  • Costa S, Maluski H, Lardeaux J-M. 1993. 40Ar/39Ar chronology of Variscan tectono-metamorphic events in an exhumed crustal nappe: the Monts du Lyonnais complex (Massif Central, France). Chem. Geol. 105: 339–359. [CrossRef] [Google Scholar]
  • Couturié JP, Caen-Vachette M. 1980. Age Westphalien des leucogranites recoupant le granite de la Margeride (Massif Central français). CR Acad Sci Paris 291: 43–45. [Google Scholar]
  • Couturié JP, Vachette-Caen M, Vialette Y. 1979. Age Namurien d’un laccolite granitique différencié par gravité : le granite de la Margeride (Massif Central français). CR Acad Sci Paris 289: 449–452. [Google Scholar]
  • Couzinié S. 2017. Evolution of the continental crust and significance of the zircon record, a case study from the French Massif Central. Universiteit van Stellenbosch. [Google Scholar]
  • Couzinié S, Moyen J-F, Villaros A, Paquette J-L, Scarrow JH, Marignac C. 2014. Temporal relationships between Mg-K mafic magmatism and catastrophic melting of the Variscan crust in the southern part of Velay Complex (Massif Central, France). J. Geosci. 69–86. [Google Scholar]
  • Couzinié S, Laurent O, Moyen J-F, Zeh A, Bouilhol P, Villaros A. 2016. Post-collisional magmatism: Crustal growth not identified by zircon Hf-O isotopes. Earth Planet. Sci. Lett. 456: 182–195. [Google Scholar]
  • Couzinié S, Laurent O, Poujol M, Mintrone M, Chelle-Michou C, Moyen J-F, et al. 2017. Cadomian S-type granites as basement rocks of the Variscan belt (Massif Central, France): Implications for the crustal evolution of the north Gondwana margin. Lithos 286: 16–34. [CrossRef] [Google Scholar]
  • Couzinié S, Laurent O, Chelle-Michou C, Bouilhol P, Paquette J-L, Gannoun A-M, et al. 2019. Detrital zircon U-Pb-Hf systematics of Ediacaran metasediments from the French Massif Central: Consequences for the crustal evolution of the north Gondwana margin. Precambrian Res. 324: 269–284. [Google Scholar]
  • Cruden AR, Koyi H, Schmeling H. 1995. Diapiric basal entrainment of mafic into felsic magma. Earth Planet. Sci. Lett. 131: 321–340. [CrossRef] [Google Scholar]
  • Cuney M, Friedrich M, Blumenfeld P, Bourguignon A, Boiron MC, Vigneresse JL, et al. 1990. Metallogenesis in the French part of the Variscan orogen. Part I: U preconcentrations in pre-Variscan and Variscan formations – a comparison with Sn, W and Au. Tectonophysics 177: 39–57. [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. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 316: 1383–1390. [Google Scholar]
  • Cuney M, Alexandrov P, de Veslud CLC, Cheilletz A, Raimbault L, Ruffet G, et al. 2002. The timing of W-Sn-rare metals mineral deposit formation in the Western Variscan chain in their orogenic setting: the case of the Limousin area (Massif Central, France). Geol. Soc. Lond. Spec. Publ. 204: 213–228. [CrossRef] [Google Scholar]
  • de Saint Blanquat M. 1993. La faille normale ductile du massif du Saint Barthélémy. Evolution hercynienne des massifs nord-pyrénéens catazonaux considérée du point de vue de leur histoire thermique. Geodin. Acta 6: 59–77. [CrossRef] [Google Scholar]
  • Delfour J, Guerrot C. 1997. Âge Viséen inférieur du microgranite de Picampoix (Nièvre). Contribution à l’étude du magmatisme carbonifère du Morvan. Géologie Fr. 2: 3–12. [Google Scholar]
  • Delfour J, Beurrier M, Tegyey M, Lemiere B, Kerrien Y, Mouterde R, et al. 1989. Carte géologique de la France (1/50 000), feuille de Tarare (697). Orléans : BRGM. [Google Scholar]
  • Delor C, Leyreloup A, Bodinier J-L, Burg J-P. 1986. Découverte d’éclogites à glaucophane dans la klippe de Najac (Massif Central, France): nouveaux témoins océaniques d’un stade haute pression dans la chaîne de collision varisque. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 302: 739–744. [Google Scholar]
  • Delor C, Burg J-P, Guiraud M, Leyreloup A. 1987. Les métapélites à phengite-chloritoïde-grenat-staurotide-disthène de la klippe de Najac-Carmaux : nouveaux marqueurs d’un métamorphisme de haute pression varisque en Rouergue occidental. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 305: 589–595. [Google Scholar]
  • Deloule E, Alexandrov P, Cheilletz A, Laumonier B, Barbey P. 2002. In-situ U-Pb zircon ages for Early Ordovician magmatism in the eastern Pyrenees, France: the Canigou orthogneisses. Int. J. Earth Sci. 91: 398–405. [Google Scholar]
  • Demange M. 1980. Arguments for allochemistry of the cordierite migmatites in the agout massif (Montagne-Noire, France). Comptes Rendus Hebd. Séances Académie Sci. Sér. D 291: 367–370. [Google Scholar]
  • Demange M, Jamet P. 1985. Le stade majeur du métamorphisme est de type moyenne pression sur le flanc sud de la Montagne Noire dans la région de Labastide-Rouairoux (Tarn, France). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 301: 603–606. [Google Scholar]
  • Demay A. 1948. Tectonique antéstéphanienne du Massif central. Impr. nationale. [Google Scholar]
  • Denèle Y, Laumonier B, Paquette J-L, Olivier P, Gleizes G, Barbey P. 2014. Timing of granite emplacement, crustal flow and gneiss dome formation in the Variscan segment of the Pyrenees. Geol. Soc. Lond. Spec. Publ. 405: 265–287. [Google Scholar]
  • Dewey JF, Burke KC. 1973. Tibetan, Variscan, and Precambrian basement reactivation: products of continental collision. J. Geol. 81: 683–692. [CrossRef] [Google Scholar]
  • Didier A, Bosse V, Boulvais P, Bouloton J, Paquette J-L, Montel J-M, et al. 2013. Disturbance versus preservation of U–Th–Pb ages in monazite during fluid–rock interaction: textural, chemical and isotopic in situ study in microgranites (Velay Dome, France). Contrib. Mineral. Petrol. 165: 1051–1072. [CrossRef] [Google Scholar]
  • D’Lemos RS, Strachan RA, Topley CG. 1990. The Cadomian orogeny in the North Armorican Massif: a brief review. Geol. Soc. Lond. Spec. Publ. 51: 3–12. [CrossRef] [Google Scholar]
  • Do Couto D, Faure M, Augier R, Cocherie A, Rossi P, Li X-H, et al. 2016. Monazite U–Th–Pb EPMA and zircon U–Pb SIMS chronological constraints on the tectonic, metamorphic, and thermal events in the inner part of the Variscan orogen, example from the Sioule series, French Massif Central. Int. J. Earth Sci. 105: 557–579. [Google Scholar]
  • Domeier M. 2016. A plate tectonic scenario for the Iapetus and Rheic oceans. Gondwana Res. 36: 275–295. [CrossRef] [Google Scholar]
  • Domeier M, Torsvik TH. 2014. Plate tectonics in the late Paleozoic. Geosci. Front. 5: 303–350. [CrossRef] [Google Scholar]
  • Doublier MP, Potel S, Wemmer K. 2006. Age and grade of metamorphism in the eastern Monts de Lacaune-implications for the collisional accretion in Variscan externides (French Massif Central). Geodin. Acta 19: 391–407. [CrossRef] [Google Scholar]
  • Doublier MP, Potel S, Wemmer K. 2015. The tectono-metamorphic evolution of the very low-grade hangingwall constrains two-stage gneiss dome formation in the Montagne Noire (Southern France). J. Metamorph. Geol. 33: 71–89. [CrossRef] [Google Scholar]
  • Downes H, Duthou J-L. 1988. Isotopic and trace-element arguments for the lower-crustal origin of Hercynian granitoids and pre-Hercynian orthogneisses, Massif Central (France). Chem. Geol. 68: 291–308. [Google Scholar]
  • Downes H, Leyreloup A. 1986. Granulitic xenoliths from the French Massif Central – Petrology, Sr and Nd isotope systematics and model age estimates. Geol. Soc. Lond. Spec. Publ. 24: 319–330. [CrossRef] [Google Scholar]
  • Downes H, Dupuy C, Leyreloup AF. 1990. Crustal evolution of the Hercynian belt of Western Europe: Evidence from lower-crustal granulitic xenoliths (French Massif Central). Chem. Geol. 83: 209–231. [CrossRef] [Google Scholar]
  • Downes H, Kempton PD, Briot D, Harmon RS, Leyreloup AF. 1991. Pb and O isotope systematics in granulite facies xenoliths, French Massif Central: implications for crustal processes. Earth Planet. Sci. Lett. 102: 342–357. [CrossRef] [Google Scholar]
  • Downes H, Shaw A, Williamson BJ, Thirlwall MF. 1997. Sr, Nd and Pb isotopic evidence for the lower crustal origin of Hercynian granodiorites and monzogranites, Massif Central, France. Chem. Geol. 140: 289–289. [CrossRef] [Google Scholar]
  • Dubuisson G, Mercier J-C, Girardeau J, Frison J-Y. 1989. Evidence for a lost ocean in Variscan terranes of the western Massif Central, France. Nature 337: 729. [CrossRef] [Google Scholar]
  • Ducrot J, Lancelot JR, Reille JL. 1979. Datation en Montagne Noire d’un temoin d’une phase majeure d’amincissement crustal caracteristique de l’Europe prevarisque. Bull. Société Géologique Fr. 7: 501–505. [CrossRef] [Google Scholar]
  • Ducrot J, Lancelot JR, Marchand J. 1983. Datation U-Pb sur zircons de l’eclogite de la Borie (Haut-Allier, France) et consequences sur l’evolution ante-hercynienne de l’Europe occidentale. Earth Planet. Sci. Lett. 62: 385–394. [CrossRef] [Google Scholar]
  • Dufour E. 1985. Granulite facies metamorphism and retrogressive evolution of the Monts du Lyonnais metabasites (Massif Central, France). Lithos 18: 97–113. [CrossRef] [Google Scholar]
  • Duguet M. 2003. Evolution tectono-métamorphique des unités de type Thiviers-Payzac dans la chaîne hercynienne française (Massif Central et Vendée). Université d’Orléans. [Google Scholar]
  • Dupis A, Robin G, Durandeau A, Lameyre J, Vauchelle L. 1990. Etude geophysique de l’extremite occidentale du granite de Gueret. Bull. Société Géologique Fr. 6: 683–691. [CrossRef] [Google Scholar]
  • Dupraz J, Didier J. 1988. Le complexe anatectique du Velay (Massif central français): structure d’ensemble et évolution géologique. Géologie Fr. 4: 73–88. [Google Scholar]
  • Dupuy C, Leyreloup A, Vernieres J. 1979. The lower continental crust of the Massif Central (Bournac, France) – With special references to REE, U and Th composition, evolution, heat-flow production. Phys. Chem. Earth 11: 401–415. [CrossRef] [Google Scholar]
  • Dusséaux C. 2019. Topographic reconstructions of the Variscan belt of Western Europe through the study of fossil hydrothermal systems. University of Pymouth. [Google Scholar]
  • Dusséaux C, Gébelin A, Boulvais P, Gardien V, Grimes S, Mulch A. 2019. Meteoric fluid‐rock interaction in Variscan shear zones. Terra Nova. [Google Scholar]
  • Duthou JL. 1978. Les granitoides du Haut Limousin (Massif central francais). Chronologie Rb-Sr de leur mise en place. Le thermometamorphisme carbonifere. Bull. Société Géologique Fr. 7: 229–235. [Google Scholar]
  • Duthou JL, Piboule M, Gay M, Dufour E. 1981. Rb-Sr Dating of orthogranulites from the Monts du Lyonnais (Massif Central, France). Comptes Rendus Académie Sci. Sér. II 292: 749–752. [Google Scholar]
  • Duthou JL, Cantagrel J, Didier J, Vialette Y. 1984. Palaeozoic granitoids from the French Massif Central: age and origin studied by 87Rb/87Sr system. Phys. Earth Planet. Inter. 35: 131–144. [Google Scholar]
  • Duthou JL, Chenevoy M, Gay M. 1994. Rb-Sr middle Devonian age of cordierite bearing migmatites from Lyonnais area (French Massif Central). Comptes Rendus Acad. Sci. Ser. 2 Sci. Terre Planetes 319: 791–796. [Google Scholar]
  • Duthou J-L, Chenevoy M, Gay M. 1998. Présence d’un magmatisme d’âge Viséen moyen dans le versant sud du massif du Pilat, Massif central oriental ; conséquences. Comptes Rendus Académie Sci. – Ser. IIA-Earth Planet. Sci. 327: 749–754. [Google Scholar]
  • Echtler H. 1990. Geometry and kinematics of recumbent folding and low-angle detachment in the Pardailhan nappe (Montagne Noire, Southern French Massif Central). Tectonophysics 177: 109–123. [CrossRef] [Google Scholar]
  • Echtler H, Malavieille J. 1990. Extensional tectonics, basement uplift and Stephano-Permian collapse basin in a late Variscan metamorphic core complex (Montagne Noire, Southern Massif Central). Tectonophysics 177: 125–138. [CrossRef] [Google Scholar]
  • Edel JB, Schulmann K, Lexa O, Lardeaux JM. 2018. Late Palaeozoic palaeomagnetic and tectonic constraints for amalgamation of Pangea supercontinent in the European Variscan belt. Earth Sci. Rev. 177: 589–612. [CrossRef] [Google Scholar]
  • Ellenberger F. 1967. Replis de micaschistes et tectonique d’infrastructure au sein du massif gneissique du Caroux (zone axiale de la Montagne Noire). Impr. Louis-Jean. [Google Scholar]
  • Engel W. 1984. Migration of folding and flysch sedimentation on the southern flank of the variscan belt (Montagne Noire, Mouthoumet massif, Pyrenees). Z. Dtsch. Geol. Ges. 279–292. [Google Scholar]
  • Engel W, Feist R, Franke W. 1978. Synorogenic gravitational transport in the Carboniferous of the Montagne Noire (S-France). Z. Dtsch. Geol. Ges. 461–472. [Google Scholar]
  • Engel W, Feist R, Franke W. 1981. Le Carbonifère anté-stéphanien de la Montagne Noire : rapports entre mise en place des nappes et sédimentation. Bureau de recherches géologiques et minières. [Google Scholar]
  • England PC, Thompson AB. 1984. Pressure–temperature–time paths of regional metamorphism I. Heat transfer during the evolution of regions of thickened continental crust. J. Petrol. 25: 894–928. [CrossRef] [Google Scholar]
  • Escher A, Beaumont C. 1997. Formation, burial and exhumation of basement nappes at crustal scale: a geometric model based on the Western Swiss-Italian Alps. J. Struct. Geol. 19: 955–974. [CrossRef] [Google Scholar]
  • Faure M. 1995. Late orogenic carboniferous extensions in the Variscan French Massif Central. Tectonics 14: 132–153. [CrossRef] [Google Scholar]
  • Faure M, Pons J. 1991. Crustal thinning recorded by the shape of the Namurian-Westphalian leucogranite in the Variscan belt of the northwest Massif Central, France. Geology 19: 730.<0730:CTRBTS>2.3.CO;2. [CrossRef] [Google Scholar]
  • Faure M, Pin C, Mailhé D. 1979. Les roches mylonitiques associées au charriage du groupe leptyno-amphibolique sur les schistes du Lot dans la région de Marvejols (Lozère). CR Acad Sci Paris 288, 167–170. [Google Scholar]
  • Faure M, Prost AE, Lasne E. 1990. Déformation ductile extensive d’âge Namuro-Westphalien dans le plateau d’Aigurande, Massif central français. Bull. Société Géologique Fr. 6: 189–197. [CrossRef] [Google Scholar]
  • Faure M, Grolier J, Pons J. 1993. Extensional ductile tectonics of the Sioule metamorphic series (Variscan French Massif Central). Geol. Rundsch. 82: 461–474. [CrossRef] [Google Scholar]
  • Faure M, Leloix C, Roig J-Y. 1997. L’évolution polycyclique de la chaîne hercynienne. Bull. Société Géologique Fr. 168: 695–705. [Google Scholar]
  • Faure M, Charonnat X, Chauvet A. 1999. Structural map and tectonic evolution of the Cevennes para-autochthonous domain of the Hercynian belt (French Massif Central). Comptes Rendus Acad. Sci. Ser. II – A Earth Planet. Sci. 6: 401–407. [Google Scholar]
  • Faure M, Monié P, Pin C, Maluski H, Leloix C. 2002. Late Viséan thermal event in the northern part of the French Massif Central: new 40 Ar/39 Ar and Rb-Sr isotopic constraints on the Hercynian syn-orogenic extension. Int. J. Earth Sci. 91: 53–75. [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. [CrossRef] [Google Scholar]
  • Faure M, Lardeaux J-M, Ledru P. 2009a. A review of the pre-Permian geology of the Variscan French Massif Central. Comptes Rendus Geosci. 341: 202–213. [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. Société Géologique Fr. 180: 473–481. [Google Scholar]
  • Faure M, Cocherie A, Mézème EB, Charles N, Rossi P. 2010. Middle Carboniferous crustal melting in the Variscan Belt: New insights from U–Th–Pbtot. monazite and U–Pb zircon ages of the Montagne Noire Axial Zone (southern French Massif Central). Gondwana Res. 18: 653–673. [CrossRef] [Google Scholar]
  • Faure M, Cocherie A, Gaché J, Esnault C, Guerrot C, Rossi P, et al. 2014. Middle Carboniferous intracontinental subduction in the Outer zone of the Variscan belt (Montagne Noire axial zone, French Massif Central): Multimethod geochronological approach of polyphase metamorphism. Geol. Soc. Lond. Spec. Publ. 405: 289–311. [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. Comptes Rendus Géoscience 349: 126–136. [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. Comptes Rendus Géoscience 349: 126–136. [CrossRef] [Google Scholar]
  • Feist R, Galtier J. 1985. Découverte de flores d’âge namurien probable dans le flysch à olistolites de Cabrières (Hérault). Implication sur la durée de la sédimentation synorogénique dans la Montagne Noire (France méridionale). Comptes-Rendus Séances Académie Sci. Sér. 2 Mécanique-Phys. Chim. Sci. Univers Sci. Terre 300: 207–212. [Google Scholar]
  • Feybesse J-L. 1981. Tectonique de la région de Laroquebrou (Cantal, Massif Central français). Rôle de la déformation ductile et évolution du Sillon Houiller. Clermont-Ferrand. [Google Scholar]
  • Feybesse J-L, Lardeaux J-M, Johan V, Tegyey M, Dufour E, Lemiere B, et al. 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. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 307, 991–996. [Google Scholar]
  • Feybesse JL, Couturié JP, Ledru P, Johan, V. 1995. Les granites de la Margeride, de Chambon-le-Château et de Saint-Christophe-d’Allier (Massif Central) : des laccolites synchrones des derniers stades de l’épaississement varisque. Géologie Fr. 1: 27–45. [Google Scholar]
  • Floyd PA, Leveridge BE. 1987. Tectonic environment of the Devonian Gramscatho basin, south Cornwall: framework mode and geochemical evidence from turbiditic sandstones. J. Geol. Soc. 144: 531–542. [CrossRef] [Google Scholar]
  • Forestier F-H. 1961. Métamorphisme hercynien et antéhercynien dans le bassin du Haut-Allier (Massif central français). Faculté des Sciences de l’Université de Clermont-Ferrand. [Google Scholar]
  • Fortey RA, Cocks LRM. 2003. Palaeontological evidence bearing on global Ordovician-Silurian continental reconstructions. Earth-Sci. Rev. 61: 245–307. [Google Scholar]
  • Foster DA, Schafer C, Fanning CM, Hyndman DW. 2001. Relationships between crustal partial melting, plutonism, orogeny, and exhumation: Idaho-Bitterroot batholith. Tectonophysics 342: 313–350. [CrossRef] [Google Scholar]
  • Fowler MB, Henney PJ, Darbyshire DPF, Greenwood PB. 2001. Petrogenesis of high Ba-Sr granites: the Rogart pluton, Sutherland. J. Geol. Soc. 158: 521–534. [CrossRef] [Google Scholar]
  • François T. 2009. Contraintes géochimiques et géochronologiques sur l’origine et la mise en place des granites du Mont Lozère (Master 2). Université Montpellier 2 Sciences et Techniques du Languedoc. [Google Scholar]
  • Franke W. 1989. Variscan plate tectonics in Central Europe – Current ideas and open questions. Tectonophysics 169: 221–228. [CrossRef] [Google Scholar]
  • Franke W. 2000. The mid-European segment of the Variscides: tectonostratigraphic units, terrane boundaries and plate tectonic evolution. Geol. Soc. Lond. Spec. Publ. 179: 35–61. [Google Scholar]
  • Franke W. 2014. Topography of the Variscan orogen in Europe: failed-not collapsed. Int. J. Earth Sci. 103: 1471–1499. [CrossRef] [Google Scholar]
  • Franke W, Engel W. 1986. Synorogenic sedimentation in the Variscan Belt of Europe. Bull. Société Géologique Fr. 2: 25–33. [CrossRef] [Google Scholar]
  • Franke W, Doublier MP, Klama K, Potel S, Wemmer K. 2011. Hot metamorphic core complex in a cold foreland. Int. J. Earth Sci. 100: 753–785. [CrossRef] [Google Scholar]
  • Franke W, Cocks LRM, Torsvik TH. 2017. The Palaeozoic Variscan oceans revisited. Gondwana Res. 48: 257–284. [Google Scholar]
  • Gaertner (Von) HR. 1937. Der bau des Französischen Zentralplateaus. Geol Rundsch 48–68. [CrossRef] [Google Scholar]
  • Gardien V. 1990. Garnet and staurolite as relictual phases within the Low-Pressure facies series of the Pilat Unit (French Massif Central) – A record of polyphase tectonometamorphic reequilibration. Comptes Rendus Académie Sci. Sér. II 310: 233–240. [Google Scholar]
  • Gardien V. 1993. High to medium Pressure relics in the eastern Vivarais series (Eastern part of the French Massif-Central). Comptes Rendus Académie Sci. Sér. II 316: 1247–1254. [Google Scholar]
  • Gardien V, Lardeaux J-M. 1991. Découverte d’éclogites dans la synforme de Maclas: extension de l’unité supérieure des gneiss à l’Est du Massif Central. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 312: 61–68. [Google Scholar]
  • Gardien V, Lardeaux JM, Misseri M. 1988. The Monts du Lyonnais peridotites (MCF) – A record of Paleozoic upper mantle subduction. Comptes Rendus Acad. Sci. 307: 1967–1972. [Google Scholar]
  • Gardien V, Tegyey M, Lardeaux JM, Misseri M, Dufour E. 1990. Crust‐mantle relationships in the French Variscan chain: the example of the Southern Monts du Lyonnais unit (eastern French Massif Central). J. Metamorph. Geol. 8: 477–492. [CrossRef] [Google Scholar]
  • Gardien V, Thompson AB, Grujic D, Ulmer P. 1995. Experimental melting of biotite + plagioclase + quartz ± muscovite assemblages and implications for crustal melting. J. Geophys. Res. Solid Earth 100: 15581–15591. [Google Scholar]
  • Gardien V, Lardeaux J-M, Ledru P, Allemand P, Guillot S. 1997. Metamorphism during late orogenic extension; insights from the French Variscan belt. Bull. Société Géologique Fr. 168: 271–286. [Google Scholar]
  • Gardien V, Vanderhaeghe O, Arnaud N, Cocherie A, Grange M, Lécuyer C. 2011. Thermal maturation and exhumation of a middle orogenic crust in the Livradois area (French Massif Central). Bull. Société Géologique Fr. 182: 5–24. [Google Scholar]
  • Gay M, Peterlongo JM, Caen-Vachette M. 1981. Age radiométrique des granites en massifs allongés et en feuillets minces syn-tectoniques dans les Monts du Lyonnais (Massif central français). CR Acad Sci Paris 293: 993–996. [Google Scholar]
  • Gay M, Briand B, Chenevoy M, Piboule M. 1982. Évolution structurale de la série métamorphique du Vivarais oriental (Massif central). Bull Bur Rech Geol Min Fr 3: 219–232. [Google Scholar]
  • Gebauer D, Bernard-Griffiths J, Grünenfelder M. 1981. U-Pb zircon and monazite dating of a mafic-ultramafic complex and its country rocks. Contrib. Mineral. Petrol. 76: 292–300. [Google Scholar]
  • Gébelin A. 2004. Déformation et mise en place des granites (360-300Ma) dans un segment de la chaîne Varisque (plateau de Millevaches, Massif Central). Université Montpellier II-Sciences et Techniques du Languedoc. [Google Scholar]
  • Gebelin A, Martelet G, Brunel M, Faure M, Rossi P. 2004. Late Hercynian leucogranites modelling as deduced from new gravity data: the example of the Millevaches massif (Massif Central, France). Bull. Société Géologique Fr. 175: 239–248. [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]
  • Gébelin A, Brunel M, Monié P, Faure M, Arnaud N. 2007. Transpressional tectonics and Carboniferous magmatism in the Limousin, Massif Central, France: Structural and 40Ar/39Ar investigations. Tectonics 26. [Google Scholar]
  • Gébelin A, Roger F, Brunel M. 2009. Syntectonic crustal melting and high-grade metamorphism in a transpressional regime, Variscan Massif Central, France. Tectonophysics 477: 229–243. [CrossRef] [Google Scholar]
  • Gerbault M, Martinod J, Hérail G. 2005. Possible orogeny-parallel lower crustal flow and thickening in the Central Andes. Tectonophysics 399: 59–72. [CrossRef] [Google Scholar]
  • Gèze B. 1949. Etude géologique de la Montagne Noire et des Cévennes méridionales. Société géologique de France. [Google Scholar]
  • Gibson RL. 1991. Hercynian low-pressure-high-temperature regional metamorphism and subhorizontal foliation development in the Canigou massif, Pyrenees, France – Evidence for crustal extension. Geology 19: 380–383. [CrossRef] [Google Scholar]
  • Girardeau J, Dubuisson G, Mercier J-CC. 1986. Cinématique de mise en place des ophiolites et nappes crystallophiliennes du Limousin, Ouest du Massif Central francais. Bull. Société Géologique Fr. 2: 849–860. [CrossRef] [Google Scholar]
  • Gleizes G, Leblanc D, Bouchez JL. 1997. Variscan granites of the Pyrenees revisited: their role as syntectonic markers of the orogen. Terra Nova 9: 38–41. [CrossRef] [Google Scholar]
  • Godard G. 1990. Découverte d’éclogites, de péridotites à spinelle et d’amphibolites à anorthite, spinelle et corindon dans le Morvan. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 310: 227–232. [Google Scholar]
  • Godard G. 2001. The Les Essarts eclogite-bearing metamorphic Complex (Vendée, Southern Armorican Massif, France). Géologie Fr. 1: 19–51. [Google Scholar]
  • Gordon SM, Whitney DL, Teyssier C, Fossen H, Kylander-Clark A. 2016. Geochronology and geochemistry of zircon from the northern Western Gneiss Region: Insights into the Caledonian tectonic history of western Norway. Lithos 246–247: 134–148. [CrossRef] [Google Scholar]
  • Gourgaud A. 1973. Granites et migmatites du Forez au sud de Montbrison. Université de Clermont Ferrand. [Google Scholar]
  • Grolier J, Letourneur J. 1968. L’évolution tectonique du grand Sillon Houiller du Massif Central francais, paper presented at XXIII International Geological Congress, Int. Union Geod Geophys, Florence, Italy. [Google Scholar]
  • Guergouz C, Martin L, Vanderhaeghe O, Thébaud N, Fiorentini M. 2018. Zircon and monazite petrochronologic record of prolonged amphibolite to granulite facies metamorphism in the Ivrea-Verbano and Strona-Ceneri Zones, NW Italy. Lithos 308–309: 1–18. [CrossRef] [Google Scholar]
  • Guillot PL, Doubinger J. 1971. Découverte d’Acritarches dans les schistes sériciteux de Génis (Dordogne). C R Hebd Séanc Acad Sci Paris 272: 2763–2764. [Google Scholar]
  • Guillot PL, Lefevre J. 1975. Découvertes de conodontes dans le calcaire à entroques de Génis en Dordogne (série métamorphique du Bas-Limousin). Comptes Rendus L’Académie Sci. Paris D 280: 1529–1530. [Google Scholar]
  • Hamet J, Mattauer M. 1977. Age hercynien, déterminé par la méthode Rb-Sr du granite de l’Aigoual, conséquences structurales. C R Somm Soc Géol Fr 2: 80–84. [Google Scholar]
  • Hamilton MA, Murphy JB. 2004. Tectonic significance of a Llanvirn age for the Dunn Point volcanic rocks, Avalon terrane, Nova Scotia, Canada: implications for the evolution of the Iapetus and Rheic oceans. Tectonophysics 379: 199–209. [CrossRef] [Google Scholar]
  • Hanmer SK. 1977. Age and tectonic implications of the Baie d’Audierne basic-ultrabasic complex. Nature 270: 336. [CrossRef] [Google Scholar]
  • Harrison TM, Watson EB. 1983. Kinetics of zircon dissolution and zirconium diffusion in granitic melts of variable water content. Contrib. Mineral. Petrol. 84: 66–72. [CrossRef] [Google Scholar]
  • Hasalová P, Schulmann K, Lexa O, Štípská P, Hrouda F, Ulrich S, et al. 2008. Origin of migmatites by deformation‐enhanced melt infiltration of orthogneiss: A new model based on quantitative microstructural analysis. J. Metamorph. Geol. 26: 29–53. [CrossRef] [Google Scholar]
  • Hasalova P, Weinberg RF, MacRae C. 2011. Microstructural evidence for magma confluence and reusage of magma pathways: implications for magma hybridization, Karakoram Shear Zone in NW India. J. Metamorph. Geol. 29: 875–900. [CrossRef] [Google Scholar]
  • Heilimo E, Halla J, Hölttä P. 2010. Discrimination and origin of the sanukitoid series: Geochemical constraints from the Neoarchean western Karelian Province (Finland). Lithos 115: 27–39. [CrossRef] [Google Scholar]
  • Henk A. 2000. Foreland-directed lower-crustal flow and its implications for the exhumation of high-pressure-high-temperature rocks. Geol. Soc. Lond. Spec. Publ. 179: 355–368. [CrossRef] [Google Scholar]
  • Henk A, von Blanckenburg F, Finger F, Schaltegger U, Zulauf G. 2000. Syn-convergent high-temperature metamorphism and magmatism in the Variscides: a discussion of potential heat sources. Geol. Soc. Lond. Spec. Publ. 179: 387–399. [Google Scholar]
  • Henry P, Le Pichon X, Goffé B. 1997. Kinematic, thermal and petrological model of the Himalayas: constraints related to metamorphism within the underthrust indian crust and topographic elevation. Tectonophysics 273: 31–56. [CrossRef] [Google Scholar]
  • Holliger P, Cuney M, Friedrich M, Turpin L. 1986. Age carbonifère de l’unité de Brame du complexe granitique peralumineux de Saint-Sylvestre (NO Massif Central) défini par les données isotopiques U-Pb sur zircon et monazite. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 303: 1309–1314. [Google Scholar]
  • Hottin A-M, Calvez JY. 1988. Résultats analytiques K-Ar et Rb-Sr sur quelques minéraux du forage de Sancerre-Couy. BRGM 137. [Google Scholar]
  • Houseman GA, McKenzie DP, Molnar P. 1981. Convective instability of a thickened boundary layer and its relevance for the thermal evolution of continental convergent belts. J. Geophys. Res. Solid Earth 86: 6115–6132. [CrossRef] [Google Scholar]
  • Huerta AD, Royden LH, Hodges KV. 1996. The interdependence of deformational and thermal processes in mountain belts. Science 273: 637–639. [CrossRef] [Google Scholar]
  • Isnard H. 1996. Datation par la méthode U–Pb sur monazites des granites du Mont Lozère et de l’Est de la Margeride (laccolites de Chambon-le- Château et de St-Christophe d’Allier). [Google Scholar]
  • Janousek V, Holub FV. 1997. Two distinct mantle sources of Hercynian magmas intruding the Moldanubian unit, Bohemian Massif, Czech Republic. J. Geosci. 42: 10–0. [Google Scholar]
  • Janoušek V, Braithwaite CJ, Bowes DR, Gerdes A. 2004. Magma-mixing in the genesis of Hercynian calc-alkaline granitoids: an integrated petrographic and geochemical study of the Sázava intrusion, Central Bohemian Pluton, Czech Republic. Lithos 78: 67–99. [CrossRef] [Google Scholar]
  • Joanny V, Lardeaux J-M, Trolliard G, Boudeulle M. 1989. La transition omphacite→ diopside + plagioclase dans les éclogites du Rouergue (Massif Central Français): un exemple de précipitation discontinue. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 309: 1923–1930. [Google Scholar]
  • Joly A. 2007. Relations plutons et discontinuités lithosphériques: approche pluridisciplinaire de la mise en place de plutons granitiques le long du Sillon Houiller (Massif Central Français): apports des études de terrain et des données gravimétriques, magnétiques et ASM pour des modélisations 3D. PhD Thesis, Orléans. [Google Scholar]
  • Joly A, Martelet G, Chen Y, Faure M. 2008. A multidisciplinary study of a syntectonic pluton close to a major lithospheric-scale fault – Relationships between the Montmarault granitic massif and the Sillon Houiller Fault in the Variscan French Massif Central: 2. Gravity, aeromagnetic investigations, and 3-D geologic modelling. J. Geophys. Res. Solid Earth 113. [CrossRef] [Google Scholar]
  • Joly A, Faure M, Martelet G, Chen Y. 2009. Gravity inversion, AMS and geochronological investigations of syntectonic granitic plutons in the southern part of the Variscan French Massif Central. J. Struct. Geol. 31: 421–443. [CrossRef] [Google Scholar]
  • Jung J. 1953. Zoneographie et age des formations cristallophylliennes des massifs hercyniens francais [with discussion]. Bull. Société Géologique Fr. 6: 329–343. [CrossRef] [Google Scholar]
  • Kelsey DE, Clark C, Hand M. 2008. Thermobarometric modelling of zircon and monazite growth in melt-bearing systems: Examples using model metapelitic and metapsammitic granulites. J. Metamorph. Geol. 26: 199–212. [CrossRef] [Google Scholar]
  • Koné M. 1985. Mise en évidence de cisaillements ductiles tangentiel et décrochant dans le sud Livradois (Massif Central Français). Conséquences pour le « métamorphisme Livradois ». Comptes Rendus Acad. Sci. 301: 189–193. [Google Scholar]
  • Kössler P, Tait J, Bachtadse V, Soffel HC, Linnemann U. 1996. Paleomagnetic investigations of Lower Paleozoic rocks of the Thüringer Schiefergebirge. Terra Nostra – Schriften Alfred-Wegener-Stift. 96: 115–116. [Google Scholar]
  • Kossmat F. 1927. Gliederung der varistischen Gebirgsbaues. Abh Sächs Geol Land 1–39. [Google Scholar]
  • Kosztolanyi C. 1971. Géochronologie des gisements uranifères français par la méthode uranium-plomb. Influence du déséquilibre radioactif sur les résultats. Université de Nancy 1. [Google Scholar]
  • Kroner U, Romer RL. 2013. Two plates – Many subduction zones: the Variscan orogeny reconsidered. Gondwana Res. 24, 298–329. [CrossRef] [Google Scholar]
  • Kruckenberg SC, Vanderhaeghe O, Ferré EC, Teyssier C, Whitney DL. 2011. Flow of partially molten crust and the internal dynamics of a migmatite dome, Naxos, Greece: internal dynamics of the naxos dome. Tectonics 30. [Google Scholar]
  • Kusbach V, Ulrich S, Schulmann K. 2012. Ductile deformation and rheology of sub-continental mantle in a hot collisional orogeny: Example from the Bohemian Massif. J. Geodyn. 56: 108–123. [CrossRef] [Google Scholar]
  • Kusbach V, Janoušek V, Hasalová P, Schulmann K, Fanning CM, Erban V, et al. 2015. Importance of crustal relamination in origin of the orogenic mantle peridotite-high-pressure granulite association: example from the Náměšt’ Granulite Massif (Bohemian Massif, Czech Republic). J. Geol. Soc. 172: 479–490. [CrossRef] [Google Scholar]
  • Labrousse L, Prouteau G, Ganzhorn A-C. 2011. Continental exhumation triggered by partial melting at ultrahigh pressure. Geology 39: 1171–1174. [CrossRef] [Google Scholar]
  • Lafon JM. 1986. Géochronologie U-Pb appliquée à deux segments du massif central français, le Rouergue oriental et le Limousin central (PhD Thesis). [Google Scholar]
  • Lafon J-M, Respaut J-P. 1988. Géochronologie U-Pb et leucogranites varisques: cas des massifs de Grandrieu (Lozère) et de la Porcherie (Limousin), Massif Central français. Bull. Minéralogie 111: 225–237. [Google Scholar]
  • Lagarde J-L, Dallain C, Ledru P, Courrioux G. 1994. Strain patterns within the late Variscan granitic dome of Velay, French Massif Central. J. Struct. Geol. 16: 839–852. [CrossRef] [Google Scholar]
  • Lameyre J, Durandeau A, Laurent O, Sagon JP, Vauchelle L, Duzelier D, et al. 1988. Démonstration, par sondage, de la présence du gneiss d’Aubusson sous les granites du batholite de Guéret (Massif Central français) et de la nature tectonique du contact. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 307: 2077–2083. [Google Scholar]
  • Lapierre H, Basile C, Berly T, Canard E. 2008. Potassic late orogenic Stephanian volcanism in the Southwest French Massif central (Decazeville, Figeac, Lacapelle-Marival basins): an example for mantle metasomatism along strike-slip faults? Bull. Société Géologique Fr. 179: 491–502. [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 J-M, Dufour E. 1987. Champs de déformation superposés dans la chaîne varisque. Exemple de la zone nord des Monts du Lyonnais (Massif central français). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 305: 61–64. [Google Scholar]
  • Lardeaux JM, Ledru P, Daniel I, Duchene S. 2001. The Variscan French Massif Central—a new addition to the ultra-high pressure metamorphic “club”: exhumation processes and geodynamic consequences. Tectonophysics 332: 143–167. [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. [Google Scholar]
  • Lasnier B. 1968. Découverte de roches éclogitiques dans le groupe leptyno-amphibolique des monts du Lyonnais (Massif central français). Bull. Société Géologique Fr. 7: 179–185. [CrossRef] [Google Scholar]
  • Lasnier B. 1971. Les peridotites et pyroxenolites a grenat du Bois des Feuilles (Monts du Lyonnais) (France). Contrib. Mineral. Petrol. 34: 29–42. [CrossRef] [Google Scholar]
  • Lasserre M, Tempier P, Philibert J. 1980. Géochronologie Rb/Sr d’une intrusion cambrienne de la région de Saint-Flour (Massif Central français). Comptes Rendus Académie Sci. 737–740. [Google Scholar]
  • Laumonier B, Marignac C, Kister P. 2010. Polymetamorphism and crustal evolution of the eastern Pyrenees during the Late Carboniferous Variscan orogenesis. Bull. Soc. Geol. Fr. 181: 411–428. [Google Scholar]
  • Laurent O, Martin H, Doucelance R, Moyen J-F, Paquette J-L. 2011. Geochemistry and petrogenesis of high-K “sanukitoids” from the Bulai pluton, Central Limpopo Belt, South Africa: Implications for geodynamic changes at the Archaean–Proterozoic boundary. Lithos 123: 73–91. [CrossRef] [Google Scholar]
  • Laurent O, Rapopo M, Stevens G, Moyen JF, Martin H, Doucelance R, et al. 2014. Contrasting petrogenesis of Mg-K and Fe-K granitoids and implications for post-collisional magmatism: Case study from the Late-Archean Matok pluton (Pietersburg block, South Africa). Lithos 196-197: 131–149. [CrossRef] [Google Scholar]
  • Laurent O, Couzinié S, Zeh A, Vanderhaeghe O, Moyen J-F, Villaros A, et al. 2017. Protracted, coeval crust and mantle melting during Variscan late-orogenic evolution: U-Pb dating in the eastern French Massif Central. Int. J. Earth Sci. 106: 421–451. [Google Scholar]
  • Le Breton N, Duthou J-L, Grolier J, Lacour A, Meyer G, Treuil M. 1986. Les diatexites à cordiérite d’Aubusson (Creuse, France) : pétrographie, composition, âge. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 303: 1557–1562. [Google Scholar]
  • Le Corre C, Auvray B, Ballevre M, Robardet M. 1991. Le Massif Armoricain / The Armorican Massif. Sci. Géologiques Bull. 44: 31–103. [CrossRef] [Google Scholar]
  • Ledru P, Autran A. 1987. L’édification de la chaîne varisque dans le Limousin, rôle de la faille d’Argentat à la limite Limousin-Millevaches. Géol Prof Fr. Thème 3: 51–91. [Google Scholar]
  • Ledru P, Marot A, Herrouin Y. 1986. Le synclinorium de Saint-Georges-sur-Loire : une unité ligérienne charriée sur le domaine centre armoricain. Découverte de métabasite à glaucophane sur la bordure sud de cette unité. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 303: 963–968. [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. Société Géologique Fr. 605–618. [Google Scholar]
  • Ledru P, Costa S, Echtler H. 1994. The massif central: structure. Pre-Mesoz. Geol. Fr. Relat. Areas 305–323. [CrossRef] [Google Scholar]
  • Ledru P, Courrioux G, Dallain C, Lardeaux JM, Montel JM, Vanderhaeghe O, et al. 2001. The Velay dome (French Massif Central): melt generation and granite emplacement during orogenic evolution. Tectonophysics 342: 207–237. [CrossRef] [Google Scholar]
  • Legendre C, Briand B, Thierry J, Lebret P, Joly A, Bertin C. 2009. Notice explicative de la carte géologique de Saint-Geniez-d’Olt (861) au 1/50 000. Éditions BRGM Orléans. [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]
  • Lenoir X, Garrido CJ, Bodinier J-L, Dautria J-M. 2000. Contrasting lithospheric mantle domains beneath the Massif Central (France) revealed by geochemistry of peridotite xenoliths. Earth Planet. Sci. Lett. 181: 359–375. [CrossRef] [Google Scholar]
  • Lerouge G, Quenardel J-M. 1988. Les zones de cisaillement carboniferes dans les plutons vendeens et leurs prolongations dans le Nord-ouest du Masif central francais. Bull. Société Géologique Fr. 4: 831–838. [CrossRef] [Google Scholar]
  • Lescuyer J-L, Cocherie A. 1992. Datation sur monozircons des métadacites de Sériès : arguments pour un âge protérozoïque terminal des schistes X de la Montagne Noire (Massif central français). Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 31: 1071–1077. [Google Scholar]
  • Letterrier J. 1978. Aspects chimiques des interactons entre les magmas basiques et leur encaissant pelitique dans le plutonisme. Bull. Société Géologique Fr. 7: 21–28. [CrossRef] [Google Scholar]
  • Lexa O, Schulmann K, Janoušek V, Štípská P, Guy A, Racek M. 2011. Heat sources and trigger mechanisms of exhumation of HP granulites in Variscan orogenic root. J. Metamorph. Geol. 29: 79–102. [Google Scholar]
  • Leyreloup A. 1974. Les enclaves catazonales remontées par les éruptions néogènes de France : nature de la croûte inférieure. Contrib. Mineral. Petrol. 46: 17–27. [CrossRef] [Google Scholar]
  • Linnemann U, Gerdes A, Drost K, Buschmann B. 2007. The continuum between Cadomian orogenesis and opening of the Rheic Ocean: Constraints from LA-ICP-MS U-Pb zircon dating and analysis of plate-tectonic setting (Saxo-Thuringian zone, northeastern Bohemian Massif, Germany). Spec. Pap. – Geol. Soc. Am. 423: 61. [Google Scholar]
  • Lister GS, Baldwin SL. 1993. Plutonism and the origin of metamorphic core complexes. Geology 21: 607–610. [CrossRef] [Google Scholar]
  • Lotout C, Pitra P, Poujol M, Van Den Driessche J. 2017. Ordovician magmatism in the Lévézou massif (French Massif Central): tectonic and geodynamic implications. Int. J. Earth Sci. 106: 501–515. [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]
  • Maillet N, Piboule M, Santallier D, Cabanis B. 1984. Diversité d’origine des ultrabasites dans la série métamorphique du Limousin. BRGM 81: 1–24. [Google Scholar]
  • Malavieille J. 2010. Impact of erosion, sedimentation, and structural heritage on the structure and kinematics of orogenic wedges: Analog models and case studies. GSA Today 4–10. [CrossRef] [Google Scholar]
  • Malavieille J, Guihot P, Costa S, Lardeaux JM, Gardien V. 1990. Collapse of the thickened Variscan crust in the French Massif Central: Mont Pilat extensional shear zone and St. Etienne Late Carboniferous basin. Tectonophysics 177: 139–149. [CrossRef] [Google Scholar]
  • Maluski H, Monié P. 1988. 40Ar–39Ar laser probe multi-dating inside single biotites of a Variscan orthogeneiss (Pinet, Massif Central, France). Chem. Geol. Isot. Geosci. Sect. 73: 245–263. [CrossRef] [Google Scholar]
  • Maluski H, Costa S, Echtler H. 1991. Late variscan tectonic evolution by thinning of earlier thickened crust. An 40Ar/39Ar study of the Montagne Noire, southern Massif Central, France. Lithos 26: 287–304. [CrossRef] [Google Scholar]
  • Marchand J. 1981. Ecaillage d’un « mélange tectonique » profond : le complexe cristallophyllien de Champtoceaux (Bretagne méridionale). CR Acad Sci Paris 293: 223–228. [Google Scholar]
  • Marignac C, Leroy J, Macaudière J, Pichavant M, Weisbrod A. 1980. Evolution tectonométamorphique d’un segment de l’orogène hercynien: les Cévennes médianes, Massif central français. CR Acad Sci Paris 291: 605–608. [Google Scholar]
  • Martin H, Smithies RH, Rapp R, Moyen J-F, Champion D. 2005. An overview of adakite, tonalite–trondhjemite–granodiorite (TTG), and sanukitoid: relationships and some implications for crustal evolution. Lithos 79: 1–24. [CrossRef] [Google Scholar]
  • Martínez Catalán JR, Arenas R, García FD, Cuadra PG, Gómez-Barreiro J, Abati J, et al. 2007. Space and time in the tectonic evolution of the northwestern Iberian Massif: Implications for the Variscan belt. In: Geological Society of America Memoirs. Geological Society of America, pp. 403–423. [Google Scholar]
  • Martínez Catalán JR, Arenas R, Abati J, Martínez SS, García FD, Suárez JF, et al. 2009. A rootless suture and the loss of the roots of a mountain chain: The Variscan belt of NW Iberia. Comptes Rendus Geosci. 341: 114–126. [CrossRef] [Google Scholar]
  • Mattauer M, Etchecopar A. 1976. Arguments en faveur de chevauchements de type himalayen dans la chaîne hercynienne du Massif Central français. Coll Int CNRS Paris 268: 261–267. [Google Scholar]
  • Mattauer M, Laurent P, Matte P. 1996. Plissement hercynien synschisteux post-nappe et étirement subhorizontal dans le versant sud de la Montagne noire (Sud du Massif central, France). Comptes Rendus Académie Sci. Sér. 2 Sci. Terre Planètes 322: 309–315. [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. 1991. Accretionary history and crustal evolution of the Variscan belt in Western Europe. Tectonophysics 196: 309–337. [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, Lancelot J, Mattauer M. 1998. La zone axiale hercynienne de la Montagne Noire n’est pas un “metamorphic core complex” extensif mais un anticlinal post-nappe à coeur anatectique. Geodin. Acta 11: 13–22. [Google Scholar]
  • Maurel O, Monié P, Respaut JP, Leyreloup AF, Maluski H. 2003. Pre-metamorphic 40Ar/39Ar and U-Pb ages in HP metagranitoids from the Hercynian belt (France). Chem. Geol. 193: 195–214. [CrossRef] [Google Scholar]
  • McCarthy A, Chelle-Michou C, Müntener O, Arculus R, Blundy J. 2018. Subduction initiation without magmatism: The case of the missing Alpine magmatic arc. Geology 46: 1059–1062. [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. Société Géologique Fr. 180: 283–292. [CrossRef] [Google Scholar]
  • Melleton J, Cocherie A, Faure M, Rossi P. 2010. Precambrian protoliths and Early Paleozoic magmatism in the French Massif Central: U-Pb data and the North Gondwana connection in the west European Variscan belt. Gondwana Res. 17: 13–25. [Google Scholar]
  • Ménard G, Molnar P. 1988. Collapse of a Hercynian Tibetan plateau into a late Palaeozoic European Basin and Range province. Nature 334: 235. [CrossRef] [Google Scholar]
  • Mercier JCC, Girardeau J, Prinzhofer A, Dubuisson G. 1985. Les complexes ophiolitiques du Limousin : structure, pétrologie et géochimie. Rapp. GPF2 Thème 3: 95–3. [Google Scholar]
  • Mercier L, Lardeaux J-M, Davy P. 1991a. On the tectonic significance of retrograde P-T-t paths in eclogites of the French Massif Central. Tectonics 10: 131–140. [CrossRef] [Google Scholar]
  • Mercier L, Van Roermund HLM, Lardeaux JM. 1991b. Comparison of PTt paths in allochthonous high pressure metamorphic terrains from the Scandinavian Caledonides and the french Massif Central: contrasted thermal structures during uplift. Geol. Rundsch. 80: 333–348. [CrossRef] [Google Scholar]
  • Mialhe J. 1980. Le massif granitique de la Borne (Cévennes). Etude pétrographique, géochimique, géochronologique et structurale (Doctorat). Université de Clermont Ferrand. [Google Scholar]
  • Michon G. 1987. Les vaugnerites de l’Est du Massif central français : apport de l’analyse statistique multivariée à l’étude géochimique des éléments majeurs. Bull. Soc. Geol. Fr. 8: 591–600. [Google Scholar]
  • Milési J-P, Lescuyer J-L. In press. The chessy zn-Cu-ba massive sulphide deposit and the Devonian brevenne volcanosedimentary belt (eastern massif central, France). BRGM. [Google Scholar]
  • Mollier B, Bouchez JL. 1982. Structuration magmatique du complexe granitique de Brâme-St Sylvestre-St Goussaud (Limousin, Massif Central Français). CR Acad Sci Ser 2 294: 1329–1334. [Google Scholar]
  • Molnar P, Lyon-Caen H. 1988. Some simple physical aspects of the support, structure, and evolution of mountain belts. In: Geological Society of America Special Papers. Geological Society of America, pp. 179–208. [CrossRef] [Google Scholar]
  • Molnar P, England P, Martinod J. 1993. Mantle dynamics, uplift of the Tibetan Plateau, and the Indian Monsoon. Rev. Geophys. 31: 357. [CrossRef] [Google Scholar]
  • Monié P, Respaut JP, Brichau S, Bouchot V, Faure M, Roig JY. 2000. 40Ar/39Ar and U-Pb geochronology applied to Au-W-Sb metallogenesis in the Cévennes and Châtaigneraie districts (Southern Massif Central, France). Orog. Gold Depos. Eur. Doc BRGM 297: 77–79. [Google Scholar]
  • Monier G. 1980. Pétrologie des granites du Sud-Millevaches (MCF). Minéralogie, géochimie, géochronologie (Doctorat 3e Cycle). Université de Clermont Ferrand. [Google Scholar]
  • Montel J-M. 1985. Xénolithes peralumineux dans les dolérites du Peyron, en Velay (Massif Central, France). Indications sur l’évolution de la croûte profonde tardihercynienne. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 301: 615–620. [Google Scholar]
  • Montel JM, Abdelghaffar R. 1993. Les granites tardi-migmatitiques du Velay (Massif Central) : principales caractéristiques pétrographiques et géochimiques. Géologie Fr. 1: 15–28. [Google Scholar]
  • Montel JM, Marignac C, Barbey P, Pichavant M. 1992. Thermobarometry and granite genesis: the Hercynian low-P, high-T Velay anatectic dome (French Massif Central). J. Metamorph. Geol. 10: 1–15. [CrossRef] [Google Scholar]
  • Montel J-M, Bouloton J, Veschambre M, Pellier C, Ceret K. 2002. Âges stéphaniens des microgranites du Velay (Massif central français). Géologie Fr. 15–20. [Google Scholar]
  • Mougeot R, Respaut J-P, Ledru P, Marignac C. 1997. U-Pb chronology on accessory minerals of the Velay anatectic dome (French Massif Central). Eur. J. Mineral. 9: 141–156. [CrossRef] [Google Scholar]
  • Moyen J-F, Laurent O, Chelle-Michou C, Couzinié S, Vanderhaeghe O, Zeh A. 2017. Collision vs. subduction-related magmatism: two contrasting ways of granite formation and implications for crustal growth. Lithos 277: 154–177. [CrossRef] [Google Scholar]
  • Murphy JB, Gutierrez-Alonso G, Nance RD, Fernandez-Suarez J, Keppie JD, Quesada C, et al. 2006. Origin of the Rheic Ocean: Rifting along a Neoproterozoic suture? Geology 34: 325–328. [CrossRef] [Google Scholar]
  • Nance RD, Gutiérrez-Alonso G, Keppie JD, Linnemann U, Murphy JB, Quesada C, et al. 2010. Evolution of the Rheic ocean. Gondwana Res. 17: 194–222. [CrossRef] [Google Scholar]
  • Nicolas A, Bouchez JL, Blaise J, Poirier JP. 1977. Geological aspects of deformation in continental shear zones. Tectonophysics 42: 55–73. [CrossRef] [Google Scholar]
  • Nicollet C. 1977. Une nouvelle éclogite à disthène et corindon primaires dans les complexes leptyno-amphiboliques du Massif central français (Lévezou, Rouergue). Bull Soc Fr Miner. Cristal. 100: 334–337. [Google Scholar]
  • Nicollet C. 1978. Pétrologie et tectonique des terrains cristallins anté-permiens du versant sud du dôme du Lévezou (Rouergue, Massif central). Bull BRGM 3: 225–263. [Google Scholar]
  • Nicollet C, Leyreloup A. 1978. Pétrologie des niveaux trondhjémitiques de haute pression associés aux éclogites et amphibolites des complexes leptyno-amphiboliques du Massif Central français. Can. J. Earth Sci. 15: 696–707. [CrossRef] [Google Scholar]
  • O’Brien PJ, Rötzler J. 2003. High-pressure granulites: formation, recovery of peak conditions and implications for tectonics. J. Metamorph. Geol. 21: 3–20. [CrossRef] [Google Scholar]
  • Pagel M, Costa S, Galibert F, Lancelot J, Maluski H, Meyer A, et al. 1992. Forage scientifique de Sancerre-Couy: Géochronologie Sm-Nd, U-Pb, Ar-Ar, Rb-Sr et chronométrie par traces de fission sur le socle. Géologie Fr. 3: 129. [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]
  • Peiffer M-T. 1986. La signification de la ligne tonalitique du Limousin. Son implication dans la structuration varisque du Massif Central français. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 303: 305–310. [Google Scholar]
  • Perchuk LL, Podladchikov YY, Polyakov AN. 1992. Hydrodynamic modelling of some metamorphic processes. J. Metamorph. Geol. 10: 311–319. [CrossRef] [Google Scholar]
  • Percival JA. 1992. Exposed crustal cross sections as windows on the lower crust. Cont. Low. Crust Dev. Geotecton. 317–362. [Google Scholar]
  • Petitpierre E, Duthou JL. 1980. Age westphalien par la méthode Rb/Sr du leucogranite de Crevant, Plateau d’Aigurande (Massif Central français). CR Acad Sci 291: 163–166. [Google Scholar]
  • Piboule M, Briand B. 1985. Geochemistry of eclogites and associated rocks of the southeastern area of the French Massif Central: origin of the protoliths. Chem. Geol. 50: 189–199. [CrossRef] [Google Scholar]
  • Pin C. 1979. Géochronologie U-Pb et microtectonique des séries métamorphiques anté-stéphaniennes de l’Aubrac et de la région de Marvejols (Massif Central). PhD Thesis, Université Montpellier II-Sciences et Techniques du Languedoc. [Google Scholar]
  • Pin C. 1981. Old inherited zircons in two synkinematic variscan granitoids: the’granite du Pinet’and the’Orthogneiss de Marvejols’(Southern French Massif Central). Neues Jahrb. Für Mineral. Abh. 142: 27–48. [Google Scholar]
  • Pin C. 1990. Variscan oceans: ages, origins and geodynamic implications inferred from geochemical and radiometric data. Tectonophysics 177: 215–227. [CrossRef] [Google Scholar]
  • Pin C, Duthou J-L. 1990. Sources of Hercynian granitoids from the French Massif Central: Inferences from Nd isotopes and consequences for crustal evolution. Chem. Geol. 83: 281–296. [CrossRef] [Google Scholar]
  • Pin C, Lancelot JR. 1978. Un exemple de magmatisme cambrien dans le Massif central; les metadiorites quartziques intrusives dans la serie du Lot. Bull. Société Géologique Fr. 7: 203–208. [CrossRef] [Google Scholar]
  • Pin C, Lancelot J. 1982. U-Pb dating of an early Paleozoic bimodal magmatism in the French Massif Central and of its further metamorphic evolution. Contrib. Mineral. Petrol. 79: 1–12. [Google Scholar]
  • Pin C, Marini F. 1993. Early Ordovician continental break-up in Variscan Europe: Nd Sr isotope and trace element evidence from bimodal igneous associations of the Southern Massif Central, France. Lithos 29: 177–196. [CrossRef] [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]
  • Pin C, Vielzeuf D. 1983. Granulites and related rocks in Variscan median Europe: a dualistic interpretation. Tectonophysics 93: 47–74. [CrossRef] [Google Scholar]
  • Pitra P, Poujol M, Van Den Driessche J, Poilvet J-C, Paquette J-L. 2012. Early Permian extensional shearing of an ordovician granite: The saint-eutrope “c/s-like” orthogneiss (montagne noire, French massif central). Comptes Rendus Géoscience 344: 377–384. [CrossRef] [Google Scholar]
  • Poilvet J-C, Poujol M, Pitra P, Van den Driessche J, Paquette J-L. 2011. The Montalet granite, Montagne Noire, France: An Early Permian syn-extensional pluton as evidenced by new U-Th-Pb data on zircon and monazite. Comptes Rendus Géoscience 343: 454–461. [CrossRef] [Google Scholar]
  • Prelević D, Akal C, Foley SF, Romer RL, Stracke A, Van Den Bogaard P. 2012. Ultrapotassic Mafic Rocks as Geochemical Proxies for Post-collisional Dynamics of Orogenic Lithospheric Mantle: the Case of Southwestern Anatolia, Turkey. J. Petrol. 53: 1019–1055. [CrossRef] [Google Scholar]
  • Quenardel J-M, Rolin P. 1984. Palaeozoic evolution of the Plateau d’Aigurande (NW Massif Central, France). Geol. Soc. Lond. Spec. Publ. 14: 63–70. [CrossRef] [Google Scholar]
  • Quenardel J-M, Santallier D, Burg J-P, Bril H, Cathelineau M, Marignac C. 1991. Le Massif Central / The central massif. Sci. Géologiques Bull. 44: 105–206. [CrossRef] [Google Scholar]
  • Rabin M, Trap P, Carry N, Fréville K, Cenki-Tok B, Lobjoie C, et al. 2015. Strain partitioning along the anatectic front in the Variscan Montagne Noire massif (southern French Massif Central). Tectonics 34: 1709–1735. [CrossRef] [Google Scholar]
  • Ramberg H. 1968. Fluid dynamics of layered systems in the field of gravity, a theoretical basis for certain global structures and isostatic adjustment. Phys. Earth Planet. Inter. 1: 63–87. [CrossRef] [Google Scholar]
  • Rapp RP, Norman MD, Laporte D, Yaxley GM, Martin H, Foley SF. 2010. Continent Formation in the Archean and Chemical Evolution of the Cratonic Lithosphere: Melt-Rock Reaction Experiments at 3-4 GPa and Petrogenesis of Archean Mg-Diorites (Sanukitoids). J. Petrol. 51: 1237–1266. [CrossRef] [Google Scholar]
  • Ravier J, Chenevoy M. 1979. Occurence of granulitic formations denoting the presence of a crustal lineament in the Sioule metamorphic series (Massif Central Français). Comptes Rendus Hebd. Séances Académie Sci. Sér. D 288: 1703–1706. [Google Scholar]
  • Respaut J-P. 1984. Géochronologie et géochimie isotopique U-Pb de la minéralisation aurifère de la mine des Pierres Plantées (Lozère) et de son encaissant : le massif granitique de la Margeride. Montpellier. [Google Scholar]
  • Rey P, Vanderhaeghe O, Teyssier C. 2001. Gravitational collapse of the continental crust: definition, regimes and modes. Tectonophysics 342: 435–449. [CrossRef] [Google Scholar]
  • Rey PF, Teyssier C, Whitney DL. 2009. The role of partial melting and extensional strain rates in the development of metamorphic core complexes. Tectonophysics 477: 135–144. [CrossRef] [Google Scholar]
  • Rey PF, Teyssier C, Kruckenberg SC, Whitney DL. 2011. Viscous collision in channel explains double domes in metamorphic core complexes. Geology 39: 387–390. [CrossRef] [Google Scholar]
  • R’Kha Chaham K, Couturie J-P, Duthou J-L, Fernandez A, Vitel G. 1990. L’orthogneiss ø eillé de l’Arc de Fix : un nouveau témoin d’âge cambrien d’un magmatisme hyper alumineux dans le Massif Central français. Comptes Rendus Académie Sci. Sér. 2 Mécanique Phys. Chim. Sci. Univers Sci. Terre 311: 845–850. [Google Scholar]
  • Robardet M. 2003. The Armorica “microplate”: fact or fiction? Critical review of the concept and contradictory palaeobiogeographical data. Palaeogeogr. Palaeoclimatol. Palaeoecol. 195: 125–148. [CrossRef] [Google Scholar]
  • Robardet M, Blaise J, Bouyx E, Gourvennec R, Lardeux H, Le Herisse A, et al. 1993. Paléogéographie de l’Europe occidentale de l’Ordovicien au Dévonien = Palaeogeography of Western Europe from the Ordovician to the Devonian. Bull. Société Géologique Fr. [Google Scholar]
  • Roger F, Respaut J-P, Brunel M, Matte P, Paquette J-L. 2004. Première datation U/Pb des orthogneiss ø eillés de la zone axiale de la Montagne noire (Sud du Massif central) : nouveaux témoins du magmatisme ordovicien dans la chaı̂ne Varisque. Comptes Rendus Geosci. 336: 19–28. [CrossRef] [Google Scholar]
  • Roger F, Teyssier C, Respaut J-P, Rey PF, Jolivet M, Whitney DL, et al. 2015. Timing of formation and exhumation of the Montagne Noire double dome, French Massif Central. Tectonophysics 640: 53–69. [CrossRef] [Google Scholar]
  • Roig J-Y, Faure M. 2000. La tectonique cisaillante polyphasee du Sud Limousin (Massif central francais) et son interpretation dans un modele d’evolution polycyclique de la chaine hercynienne. Bull. Société Géologique Fr. 171: 295–307. [CrossRef] [Google Scholar]
  • Roig J-Y, Faure M, Ledru P. 1996. Polyphase wrench tectonics in the southern french Massif Central: kinematic inferences from pre-and syntectonic granitoids. Geol. Rundsch. 85: 138–153. [Google Scholar]
  • Roig J-Y, Faure M, Maluski H. 2002. Superimposed tectonic and hydrothermal events during the late-orogenic extension in the Western French Massif Central: a structural and 40Ar/39Ar study. Terra Nova 14: 25–32. [CrossRef] [Google Scholar]
  • Rolin P, Duthou JL, Quenardel JM. 1982. Datation Rb/Sr des leucogranites de Crozant et d’Orsennes : Conséquences sur l’âge de la derniere phase de tectonique tangentielle du Plateau d’Aigurande (NW du Massif Central Français). CR Acad Sci Ser II 294: 799–802. [Google Scholar]
  • Rolin P, Marquer D, Colchen M, Cartannaz C, Cocherie A, Thiery V, et al. 2009. Famenno-Carboniferous (370-320 Ma) strike slip tectonics monitored by syn-kinematic plutons in the French Variscan belt (Massif Armoricain and French Massif Central). Bull. Société Géologique Fr. 180: 231–246. [Google Scholar]
  • Rolin P, Marquer D, Cartannaz C, Rossi P. 2014. Carboniferous magmatism related to progressive pull-apart opening in the western French Massif Central. Bull. Société Géologique Fr. 185: 171–189. [Google Scholar]
  • Roques M. 1971. Structure géologique du Massif central. In: Géologie, géomorphologie et structure profonde du Massif Central Français. Symp J Jung, Clermont-Ferrand, pp. 17–32. [Google Scholar]
  • Rosenberg CL. 2001. Deformation of partially molten granite: a review and comparison of experimental and natural case studies. Int. J. Earth Sci. 90: 60–76. [CrossRef] [Google Scholar]
  • Royden L. 1996. Coupling and decoupling of crust and mantle in convergent orogens: Implications for strain partitioning in the crust. J. Geophys. Res. Solid Earth 101: 17679–17705. [CrossRef] [Google Scholar]
  • Royden LH. 1997. Surface deformation and Lower Crustal Flow in Eastern Tibet. Science 276: 788–790. [CrossRef] [Google Scholar]
  • Rubio Pascual FJ, López-Carmona A, Arenas R. 2016. Thickening vs. extension in the Variscan belt: P–T modelling in the Central Iberian autochthon. Tectonophysics 681: 144–158. [CrossRef] [Google Scholar]
  • Sabatier H. 1991. Vaugnerites: special lamprophyre-derived mafic enclaves in some Hercynian granites from Western and Central Europe. Enclaves Granite Petrol. Elsevier Amst. 63–81. [Google Scholar]
  • Saint-Joanis R. 1975. Étude géologique du socle cristallin du Bas- Livradois (Massif central français) dans le périmètre de la feuille d’Issoire. Université Blaise Pascal, Clermont-Ferrand. [Google Scholar]
  • Sandiford M, Powell R. 1990. Some isostatic and thermal consequences of the vertical strain geometry in convergent orogens. Earth Planet. Sci. Lett. 98: 154–165. [CrossRef] [Google Scholar]
  • Santallier D. 1981. Les roches basiques dans la série métamorphique du Bas-Limousin, Massif Central (France) (PhD Thesis). [Google Scholar]
  • Santallier D, Briand B, Menot RP, Piboule M. 1988. Les complexes leptyno-amphiboliques (CLA) : revue critique et suggestions pour un meilleur emploi de ce terme. Bull. Société Géologique Fr. 4: 3–12. [Google Scholar]
  • Sawyer EW. 1994. Melt segregation in the continental crust. Geology 22: 1019–1022. [CrossRef] [Google Scholar]
  • Sawyer EW. 1998. Formation and evolution of granite magmas during crustal reworking: the significance of diatexites. J. Petrol. 39: 1147–1167. [CrossRef] [Google Scholar]
  • Sawyer EW, Cesare B, Brown M. 2011. When the continental crust melts. Elements 7: 229–234. [CrossRef] [Google Scholar]
  • Schaltegger U, Gebauer D. 1999. Pre-Alpine geochronology of the central, western and southern Alps. Schweiz. Mineral. Petrogr. Mitteilungen 79: 79–87. [Google Scholar]
  • Schenk V. 1980. U-Pb and Rb-Sr radiometric dates and their correlation with metamorphic events in the granulite-facies basement of the Serre, southern Calabria (Italy). Contrib. Mineral. Petrol. 73: 23–38. [Google Scholar]
  • Schenk V. 1981. Synchronous uplift of the lower crust of the Ivrea Zone and of southern Calabria and its possible consequences for the Hercynian orogeny in southern Europe. Earth Planet. Sci. Lett. 56: 305–320. [CrossRef] [Google Scholar]
  • Schenk V. 1989. PTt path of the lower crust in the Hercynian fold belt of southern Calabria. Geol. Soc. Lond. Spec. Publ. 43: 337–342. [CrossRef] [Google Scholar]
  • Schulmann K, Lexa O, Štípská P, Racek M, Tajčmanová L, Konopásek J, et al. 2008. Vertical extrusion and horizontal channel flow of orogenic lower crust: key exhumation mechanisms in large hot orogens? J. Metamorph. Geol. 26: 273–297. [CrossRef] [Google Scholar]
  • Schulmann K, Konopásek J, Janoušek V, Lexa O, Lardeaux J-M, Edel J-B, et al. 2009. An Andean type Palaeozoic convergence in the Bohemian massif. Comptes Rendus Geosci. 341: 266–286. [CrossRef] [Google Scholar]
  • Schulmann K, Catalán JRM, Lardeaux JM, Janoušek V, Oggiano G. 2014. The Variscan orogeny: extent, timescale and the formation of the European crust. Geol. Soc. Lond. Spec. Publ. 405: 1–6. [CrossRef] [Google Scholar]
  • Schulz B. 2009. EMP-monazite age controls on PT paths of garnet metapelites in the Variscan inverted metamorphic sequence of La Sioule, French Massif Central. Bull. Société Géologique Fr. 180: 271–282. [Google Scholar]
  • Schulz B. 2014. Early Carboniferous PT path from the Upper Gneiss Unit of Haut-Allier (French Massif Central)-reconstructed by geothermobarometry and EMP-Th-U-Pb monazite dating. J. Geosci. 59: 327–349. [CrossRef] [Google Scholar]
  • Schulz B, Triboulet C, Audren C, Feybesse J-L. 1996. Zoned garnets in metapelites and PT-deformation path interpretation of the Variscan inverted metamorphic sequence of Haut-Allier, French Massif Central. Z. Dtsch. Geol. Ges. 249–273. [Google Scholar]
  • Schulz B, Triboulet C, Audren C, Feybesse J-L. 2001. PT-paths from metapelite garnet zonations, and crustal stacking in the Variscan inverted metamorphic sequence of La Sioule, French Massif Central. Z. Dtsch. Geol. Ges. 152: 1–26. [Google Scholar]
  • Scotese CR, McKerrow WS. 1990. Revised world maps and introduction. Geol. Soc. Lond. Mem. 12: 1–21. [CrossRef] [Google Scholar]
  • Searle MP, Cottle JM, Streule MJ, Waters DJ. 2009. Crustal melt granites and migmatites along the Himalaya: melt source, segregation, transport and granite emplacement mechanisms. Earth Environ. Sci. Trans. R. Soc. Edinb. 100: 219–233. [Google Scholar]
  • Shail RK, Leveridge BE. 2009. The Rhenohercynian passive margin of SW England: Development, inversion and extensional reactivation. Comptes Rendus Geosci. 341: 140–155. [CrossRef] [Google Scholar]
  • Shaw J, Johnston ST. 2016. Terrane wrecks (coupled oroclines) and paleomagnetic inclination anomalies. Earth-Sci. Rev. 154: 191–209. [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]
  • Sintubin M, Debacker TN, Van Baelen H. 2009. Early Palaeozoic orogenic events north of the Rheic suture (Brabant, Ardenne): A review. Comptes Rendus Geosci. 341: 156–173. [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]
  • Skrzypek E, Schulmann K, Tabaud A-S, Edel J-B. 2014. Palaeozoic evolution of the Variscan Vosges mountains. Geol. Soc. Lond. Spec. Publ. 405: 45–75. [CrossRef] [Google Scholar]
  • Solar GS, Pressley RA, Brown M, Tucker RD. 1998. Granite ascent in convergent orogenic belts: testing a model. Geology 26: 711–714.