Flow of partially molten crust controlling construction, growth and collapse of the Variscan orogenic belt: the geologic record of the French Massif Central

– We present here a tectonic-geodynamic model for the generation and ﬂ ow of partially molten rocks and for magmatism during the Variscan orogenic evolution from the Silurian to the late Carboniferous based on a synthesis of geological data from the French Massif Central. Eclogite facies metamorphism of ma ﬁ c and ultrama ﬁ c rocks records the subduction of the Gondwana hyperextended margin. LP felsic and ma ﬁ c granulites in the lower crust. This model points to the primordial role of partial melting and magmatism in the tectonic-geodynamic evolution of the Variscan orogenic belt. In particular, partial melting and magma transfer (i) triggers mechanical decoupling of subducted units from the downgoing slab and their syn-orogenic exhumation; (ii) the development of an orogenic plateau by lateral ﬂ ow of the low-viscosity partially molten crust; and, (iii) the formation of metamorphic core complexes and domes that accommodate post-orogenic exhumation during gravitational collapse. All these processes contribute to differentiation and stabilisation of the orogenic crust. gravitaire accommodé par l ’ activation de détachements à faible pendage et l ’ exhumation-cristallisation de la racine orogénique partiellement fondue formant des dômes d ’ échelle crustale à c œ ur de migmatites de BP entre c. 305 et 295 Ma, concomitante au ﬂ uage latéral des unités de la zone externe de la chaine. Les derniers magmas extraits de la zone de fusion partielle forment des dykes et des laccolithes mis en place dans des détachements à faible pendage laissant derrière eux une croûte inférieure constituée de l ’ accumulation de granulites réfractaires de composition felsique à ma ﬁ que. Ce modèle met en valeur le rôle primordial de la fusion partielle et du magmatisme sur l ’ évolution tectonique-géodynamique de la ceinture orogénique Varisque. En particulier, la fusion partielle et le transfert de magma (i) déclenchent le découplage mécanique entre le panneau plongeant et les unités subductées, permettant ainsi l ’ exhumation de ces dernières, (ii) favorisent le développement d ’ un plateau orogénique par ﬂ uage latéral de la croûte partiellement fondue de faible viscosité sous l ’ effet de la gravité, (iii) conduisent à la formation de metamorphic core complex et de domes qui accommodent l ’ exhumation post-orogénique au cours de l ’ effondrement gravitaire de la chaine. Tous ces processus contribuent à la différenciation et à la stabilisation de la croûte orogénique. enclaves in high-grade nappes the Moldanubian allochthonous terrane The tholeitic to calc-alkaline signature of the LAC is interpreted as re ﬂ ecting an emplacement of the magmatic protoliths in a back-arc or immature oceanic setting. Alkaline magmas intrusive in Ediacarian 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.

orogenic belt. In particular, partial melting and magma transfer (i) triggers mechanical decoupling of subducted units from the downgoing slab and their syn-orogenic exhumation; (ii) the development of an orogenic plateau by lateral flow of the low-viscosity partially molten crust; and, (iii) the formation of metamorphic core complexes and domes that accommodate post-orogenic exhumation during gravitational collapse. All these processes contribute to differentiation and stabilisation of the orogenic crust.

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;Sawyeret 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;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 Weinberg and Mark, 2008) leading to the formation of orogenic plateaux 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 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., 2008Schulmann et al., , 2014Lardeaux 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 Couzinié et al., 2017Couzinié et al., , 2019, (ii) U-Pb dating of zircon and monazite by LA-ICP-MS on the Carboniferous plutonic rocks Laurent et al., 2017) and (iii) a detailed petrogenetic model for these granitoids . 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(Matte, , 1991Franke and Engel, 1986;Franke, 1989Franke, , 2000Kroner 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., 2006Berger et al., , 2010aLardeaux, 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).

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(Chantraine et al., , 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 lowgrade 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 Faure, 1995;Gardien et al., 1997;Bellot, 2007).

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 crosscutting 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., 2017Couzinié et al., , 2019Lotout 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.

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., 1999Barbey et al., , 2015Villaros 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 calcalkaline 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.
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 LREErich 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., 2014Couzinié et al., , 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., 2011Laurent et al., , 2014Prelević 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 sek 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 "mantlederived" granitoids.

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(Engel et al., , 1981Engel, 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 lowangle 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.

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  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). 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, 1991Matte, , 2001Lardeaux, 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., 1990Gardien et al., , 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, strikeslip 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. 8-10).

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 40 Ar/ 39 Ar 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  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(Faure et al., , 2009b.  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. 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(Burg et al., , 1989bBallèvre and Balé, 1992). Amphibolites yield an Rb-Sr isochron of 387 ± 2 Ma interpreted as dating high-grade metamorphism. 40 Ar/ 39 Ar 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., 2010aBerger et al., , 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 . The relationship between this nappe pile and the one exposed . 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.
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 40 Ar/ 39 Ar 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 highpressure 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., 1990Cuney et al., , 1993. These rocks first yielded TIMS U-Pb zircon ages of 379 ± 19 Ma and 355 ± 2 Ma  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 40 Ar/ 39 Ar 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., 2007Gébelin et al., , 2009Rolin et al., 2009Rolin et al., , 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 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., 2006Gébelin et al., , 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(Gébelin et al., , 2009. The footwall of these detachment zones represented sites of strong (meteoric) fluid-rockdeformation 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., 2008Joly et al., , 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 highgrade fabric parallel to lithological-tectonic contacts identified at the surface. These reflectors are only crosscut and offset by highangle faults and low-angle detachments. For example, highangle 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 brittleductile 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.

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. 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. 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. The top of the nappe pile is made of cordieritebearing 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. 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. 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 syntectonic 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., 1991aMercier et al., , 1991bGardien 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(Duthou et al., , 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., 1990Gardien et al., , 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 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. 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 40 Ar/ 39 Ar on recrystallized muscovite , 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  (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) 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 highpressure 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  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 . The underlying paragneisses display low-pressure granulite facies metamorphism associated with widespread partial melting dated at ca. 308 Ma . 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-garnetwhole 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., 1999Barbey et al., , 2015Fig. 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 (Laurent et al., 2014). This diagram separates granitoids related to different sources , 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). 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. 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 . 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 midamphibolite 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 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 Laurent et al., 2017) point to a protracted history of zircon growth and/or dissolutionreprecipitation 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  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., 1978Engel et al., , 1981Demange, 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).

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.

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(Matte, , 1991(Matte, , 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 (1986Matte ( , 1991Matte ( , 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 (1986Matte ( , 1991Matte ( , 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. Faure et al. (1997) refined the geodynamic model elaborated by Matte (1986Matte ( , 1991Matte ( , 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 highpressure 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 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.

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., 2007Gébelin et al., , 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 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 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 earlymid 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.

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(Matte, , 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. Highpressure 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 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., , 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 mantlederived 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 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 ribbonshaped 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 tectonicgeodynamic evolution of the Variscan belt.
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 hyperextended margin occurred as a consequence of convergence between Laurussia and Gondwana from the Silurian to the Carboniferous (Figs. 17C-E).

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(Lasnier, , 1971Pin and Vielzeuf, 1983;Bouchardon et al., 1989;Gardien et al., 1990). The presence of UHP mineral assemblages in the LAC of the Lyonnais  and Limousin (Berger et al., 2010a(Berger et al., , 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., 1989aBurg et al., , 1989bGardien, 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. 18A-C) (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., 1987Ballèvre et al., , 2002Godard, 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. 40 Ar/ 39 Ar 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(Bosse et al., , 2005Paquette 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, garnetlherzolite 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(Duthou et al., , 1994Pin 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 premetamorphic 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  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.

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., 1989aBurg et al., , 1989bGardien et al., 1990;Mercier et al., 1991aMercier et al., , 1991bLardeaux 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., 1989aBurg et al., , 1989bSchulz et al., 2001) suggests 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 tholeitic 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 Ediacarian 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. 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(Faure et al., , 2009aBurg 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(Faure et al., , 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 steeplydipping 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  and the Bohemian Massif .
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., 1997Lardeaux 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 highpressure (10-20 kbars) of paragneisses, orthogneisses and amphibolites c. 20 to 30 Ma after subduction (Nicollet and Leyreloup, 1978;Duthou et al., 1981Duthou et al., , 1994Pin 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 northvergent 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.

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 sek et al., 2004;Schulmann et al., 2009). The positive e 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 crustderived 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 lowtemperature 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(Bosse et al., , 2005Paquette 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 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., 1997Faure et al., , 2009aLardeaux, 2014;Lardeaux et al., 2014). The cooling ages obtained on the UGU of the Aigurande plateau indicates that at least some of the highgrade 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. 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., 2007Gébelin et al., , 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., 2009Rolin et al., , 2014. As presented above, these strikeslip 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., 1993Faure et al., , 2009bRoig 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., 1996Barbey 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 SiO 2 , 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., 2014Couzinié 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 crustand mantle melting .
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 . 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 sek 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  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.

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 Van Den Driessche and Brun, 1992;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 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 crustalscale 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 Vanderhaeghe, 2009). Such subdomes have been described in Naxos Aegean domain and interpreted as reflecting crustal-scale convection 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., 1999Barbey et al., , 2015Ledru 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 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., 2008Joly et al., , 2009Lapierre 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(Bruguier et al., , 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 topto 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(Schenk, , 1981(Schenk, , 1989).

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 geodynamictectonic 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).

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 . 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 calcalkaline 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).