Beaudoin Alexandre
email: alexandre.beaudoin@univ-orleans.fr

Alexandre Beaudoin

Stéphane Scaillet

Nicolas Mora

Laurent Jolivet

Romain Augier

In situ and step-heating 40 Ar/ 39 Ar dating of white mica in low-temperature shear zones (Tenda Massif, alpine Corsica, France)

Keywords: 40 Ar/ 39 Ar dating, in situ and step-heating combined technics, white mica, deformation, shear zone, Tenda massif, Alpine Corsica

Key points -Data are interpretable in terms of age-deformation relationships and capture a broad-scale localization of strain in time and space -At smaller scale, relating age with structures is more difficult because of greater complexity in 40 Ar/ 39 Ar systematics -The shear zone has a duration of ~12 Myr during exhumation and deformation ended in the upper part of the viscous/brittle transition

Introduction

The significance of 40 Ar/ 39 Ar ages in deformed rocks has been the focus of many studies since 1990. A central motivation of such studies is the potential of the 40 Ar/ 39 Ar dating to provide (i) temporal control on deformation structures and crustal shear zones (e.g., [START_REF] Kligfield | Dating of deformation phases using K-Ar and 40Ar/39Ar techniques: results from the Northern Apennines[END_REF][START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF][START_REF] Isik | The 40Ar/39Ar age of extensional ductile deformation and granitoid intrusion in the northern Menderes core complex: implications for the initiation of extensional tectonics in western Turkey[END_REF][START_REF] Turrillot | Late orogenic exhumation of the Variscan high-grade units (South Armorican Domain, western France), combined structural and 40Ar/39Ar constraints[END_REF][START_REF] Rolland | Late Palaeozoic to Mesozoic kinematic history of the Talas-Ferghana strike-slip fault (Kyrgyz West Tianshan) as revealed by 40Ar/39Ar dating of syn-kinematic white mica[END_REF][START_REF] Vincenzo | 40Ar-39Ar laser dating of ductile shear zones from central Corsica (France): Evidence of Alpine (middle to late Eocene) syn-burial shearing in Variscan granitoids[END_REF], (ii) determine and correlate the sequence of events in polydeformed rocks (e.g., [START_REF] Hames | On the loss of 40Ar* from muscovite during polymetamorphism[END_REF][START_REF] Augier | Exhumation, doming and slab retreat in the Betic Cordillera (SE Spain): in situ 40Ar/39Ar ages and P-T-d-t paths for the Nevado-Filabride complex[END_REF][START_REF] Condon | Laser ablation 40Ar/39Ar dating of metamorphic fabrics in the Caledonides of north Ireland[END_REF][START_REF] Castonguay | Dating polyphase deformation across low-grade metamorphic belts: An example based on 40Ar/39Ar muscovite age constraints from the southern Quebec Appalachians, Canada[END_REF][START_REF] Beltrando | Dating microstructures by the 40Ar/39Ar step-heating technique: Deformation-pressuretemperature-time history of the Penninic Units of the Western Alps[END_REF], and (iii) infer the kinematics, evolution and partition of shear in tectonic systems undergoing progressive deformation (e.g., [START_REF] Dunlap | Ages of deformation from K/Ar and 40Ar/39Ar dating of white micas[END_REF][START_REF] West | Progressive localization of deformation during exhumation of a major strike-slip shear zone: Norumbega fault zone, south-central Maine, USA[END_REF][START_REF] Rolland | Constraining deformation stages in brittle-ductile shear zones from combined field mapping and 40Ar/39Ar dating: The structural evolution of the Grimsel Pass area (Aar Massif, Swiss Alps)[END_REF][START_REF] Sanchez | Dating low-temperature deformation by 40Ar/39Ar on white mica, insights from the Argentera-Mercantour Massif (SW Alps)[END_REF][START_REF] Schneider | Dating the longevity of ductile shear zones: Insight from 40Ar/39Ar in situ analyses[END_REF][START_REF] Bellanger | Shortening of the European Dauphinois margin (Oisans Massif, Western Alps): New insights from RSCM maximum temperature estimates and 40Ar/39Ar in situ dating[END_REF]Kellet et al., 2016;[START_REF] Cardello | Initiation and development of the Pennine Basal Thrust (Swiss Alps): a structural and geochronological study of an exhumed megathrust[END_REF]. Practically, the question of dating deformation in such contexts is whether or not the 40 Ar/ 39 Ar targets (typically micas) can actually record deformation as a result of structural-stoichiometric recombination (recrystallization and neoblastesis), and whether these phenomena are kinetically more efficient than temperature for Ar exchange and resetting (e.g., [START_REF] Scaillet | Mg/Fe and [(Mg,Fe)Si-Al2] compositional control on argon behaviour in high-pressure white micas: A 40Ar/39Ar continuous laser-probe study from the Dora-Maira nappe of the internal western Alps, Italy[END_REF][START_REF] Vincenzo | Evidence of a close link between petrology and isotope records: constraints from SEM, EMP, TEM and in situ 40Ar-39Ar laser analyses on multiple generations of white micas (Lanterman Range, Antartica)[END_REF], 2004, 2006). Several works have emphasized that isotopic closure by cooling is prevalent when peak-metamorphic conditions during deformation are hot enough to cause the system to close after the syn-kinematic fabric has developed, resulting in closure ages potentially decoupled from the deformation itself [START_REF] Goodwin | Effects of progressive mylonitization on Ar retention in biotites from the Santa Rosa mylonite zone, California, and thermochronologic implications[END_REF][START_REF] Markley | The relation between grain size and 40Ar/39Ar date for Alpine white mica from the Siviez-Mischabel Nappe, Switzerland[END_REF][START_REF] Mulch | In-situ UV-laser 40Ar/39Ar geochronology of a micaceous mylonite : an example of defect-enhanced argon loss[END_REF][START_REF] Mulch | Recrystallization or cooling ages: in situ UV-laser 40Ar/39Ar geochronology of muscovite in mylonitic rocks[END_REF][START_REF] Scharf | Grain-size effects on the closure temperature of white mica in a crustal-scale extensional shear zone -Implications of in-situ 40Ar/39Ar laser-ablation of white mica for dating shearing and cooling (Tauern Window, Eastern Alps)[END_REF].

The case of white mica is particularly emblematic of such issues. Its 40 Ar retention kinetics have been empirically [START_REF] Hames | An empirical evaluation of the argon diffusion geometry in muscovite[END_REF][START_REF] Kirschner | Staircase 40Ar/39Ar spectra of fine-grained white mica: Timing and duration of deformation and empirical constraints on argon diffusion[END_REF] and experimentally [START_REF] Harrison | Diffusion of 40Ar in muscovite[END_REF] found to straddle the viscous/brittle transition in mid-crustal levels (330-450 °C) across which many fabrics targeted for 40 Ar/ 39 Ar dating develop.

Forming solid solutions with several end-members, white mica is also stable across a wide range of P-T variations in a large variety of lithologies. This potentially induces coupled Arstructural-stoichiometric changes according to the P-T-t-strain path undergone by the host rock and the time spent through the transitional 40 Ar retention zone while the fabric dynamically develops (e.g., [START_REF] Scaillet | Mg/Fe and [(Mg,Fe)Si-Al2] compositional control on argon behaviour in high-pressure white micas: A 40Ar/39Ar continuous laser-probe study from the Dora-Maira nappe of the internal western Alps, Italy[END_REF][START_REF] Vincenzo | Evidence of a close link between petrology and isotope records: constraints from SEM, EMP, TEM and in situ 40Ar-39Ar laser analyses on multiple generations of white micas (Lanterman Range, Antartica)[END_REF][START_REF] Putlitz | Prograde mica 40Ar/39Ar growth ages recorded in high pressure rocks (Syros, Cyclades, Greece)[END_REF]. Such a record is reflected in complex major-elements and 40 Ar/ 39 Ar variations that have defied attempts at extracting useful geochronologic information from such data (e.g., [START_REF] Scaillet | K-Ar (40Ar/39Ar) geochronology of ultrahigh pressure rocks[END_REF][START_REF] Vincenzo | Evidence of a close link between petrology and isotope records: constraints from SEM, EMP, TEM and in situ 40Ar-39Ar laser analyses on multiple generations of white micas (Lanterman Range, Antartica)[END_REF], 2006). The problems fundamentally raised by such systems may be stated as follows (see review in [START_REF] Scaillet | K-Ar (40Ar/39Ar) geochronology of ultrahigh pressure rocks[END_REF][START_REF] Vincenzo | Comparison of 40Ar-39Ar and Rb-Sr Data on Phengites from the UHP Brossasco-Isasca Unit (Dora Maira Massif, Italy): Implications for Dating White Mica[END_REF][START_REF] Warren | When can muscovite 40Ar/39Ar dating constrain the timing of metamorphic exhumation[END_REF][START_REF] Smye | The signature of devolatisation: Extraneous 40Ar systematics in high-pressure metamorphic rocks[END_REF]:

(1) How well are 40 Ar mica diffusion rates constrained and how do these interact with P-Tt induced crystallographic changes, either dynamic or static?

(2) Which are the mechanisms involved at the crystallographic scale and how do they relate with lithology, rheology, and strain partitioning in a strain gradient?

(3) How does kinetic hindrance and transformation rates affect the 40 Ar/ 39 Ar record as a function of starting mineralogy and associated rheology? (4) How does the catalytic action of metamorphic fluids combine with solid-state processes during fluid/solid interactions?

The cocktail of physics implied by such phenomena is common to natural systems, inducing complex feedback relationships between mineralogy, rheology, reactivity, fluid flow, metasomatism, and transient permeability. This has so far impaired the establishment of universal guidelines to interpret 40 Ar/ 39 Ar white mica ages in shear zones (see review in [START_REF] Oriolo | Geochronology of shear zones -A review[END_REF], to a large extent undermining their constraining power for quantitative modeling of thermo-mechanical and rheological processes.

In this study, we tackle the problem the other way round. Barring a satisfactory calibration of the physical parameters collectively involved through [START_REF] Cavazza | Apatite fission-track analysis of Neogene exhumation in northern Corsica (France)[END_REF] and ( 4) above, we simplify the problem by studying a pre-kinematic intrusion that is:

(1) spatially homogeneous in terms of starting protolith with a regionally coherent P-T-t evolution,

(2) sufficiently old to enhance the Ar isotopic contrast with newly-grown syn-kinematic micas, [START_REF] Smye | The signature of devolatisation: Extraneous 40Ar systematics in high-pressure metamorphic rocks[END_REF] spatially heterogeneous in finite-strain pattern to allow the effects of variable deformation to be examined.

Our target is the Tenda massif in Alpine Corsica (France), a regionally-sized coherent unit belonging to the European margin buried and later exhumed during the Apulia-Adria convergence. The massif consists of a ~290 Ma Late-Variscan composite intrusion [START_REF] Rossi | Caractère volcano-plutonique du magmatisme calco-alcalin composite d'âge Stéphanien supérieur-Permien inférieur en Corse[END_REF] that was extensively yet heterogeneously reworked during the Alpine orogeny [START_REF] Mattauer | Transverse lineation and large-scale structures related to Alpine obduction in Corsica[END_REF][START_REF] Jourdan | Balagne orientale et massif du Tenda (Corse septentrionale): étude structurale, interprétation des accidents et des déformations, reconstitutions géodynamiques[END_REF][START_REF] Jolivet | Ductile extension in alpine Corsica[END_REF]1991;[START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF][START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Molli | Orogenic processes and the Corsica/Apennines geodynamic evolution: insights from Taiwan[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF][START_REF] Maggi | Feedback between fluid infiltration and rheology along a regional ductile-to-brittle shear zone: The East Tenda Shear Zone (Alpine Corsica)[END_REF][START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF][START_REF] Gueydan | Sequential extension as a record of Corsica Rotation during Apennines slab roll-back[END_REF]. Large-scale geometry, kinematics and strain gradients have been recently studied [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF] along with the regional P-T evolution [START_REF] Tribuzio | Blueschist facies metamorphism of peralkaline rhyolites from the Tenda crystalline massif (northern Corsica): evidence for involvement in the Alpine subduction event?[END_REF][START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF], while recent Rb-Sr ages can be used to independently constrain spatialtemporal crystallization patterns [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF].

Our strategy combines extensive 40 Ar/ 39 Ar dating with systematic microstructural and petrological characterization of white mica sampled along major and small-scale strain gradients recognized and mapped through the massif [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. By using complementary in situ and step-heating techniques, we document and explore the link between age, grain size, chemical composition, and attempt to relate them to multi-scale strain gradients and structures.

2. Geological setting 2.1. Regional-scale nappe structure of Alpine Corsica

The Tenda massif belongs to Alpine Corsica (Fig. 1), a Cenozoic stack of tectonometamorphic units forming a westward-verging HP-LT subduction wedge overthrust on the autochthonous crystalline basement of Variscan Corsica (e.g., [START_REF] Mattauer | Transverse lineation and large-scale structures related to Alpine obduction in Corsica[END_REF][START_REF] Lahondère | Le métamorphisme éclogitique dans les orthogneiss et les métabasites ophiolitiques de la région de Farinole (Corse)[END_REF][START_REF] Fournier | Alpine Corsica Metamorphic Core Complex[END_REF][START_REF] Caron | Metamorphism and deformation in Alpine Corsica[END_REF]Vitale Brovarone et al., 2013). In this wedge, the Tenda massif corresponds to the distal continental European margin involved in subduction down to blueschist-facies conditions [START_REF] Gibbons | Alpine metamorphim of Hercynian hornblende granodiorite beneath the blueschist facies schistes-lustrés nappe of NE Corsica[END_REF][START_REF] Tribuzio | Blueschist facies metamorphism of peralkaline rhyolites from the Tenda crystalline massif (northern Corsica): evidence for involvement in the Alpine subduction event?[END_REF][START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF]Vitale Brovarone et al., 2013). It is bounded to the east by the East Tenda Shear Zone (ETSZ; Fig. 1), that puts it in contact structurally below the Schistes Lustrés nappe, an internal ocean-derived unit of Ligurian affinity that recorded eclogite-facies conditions [START_REF] Caron | High-pressure-low-temperature metamorphism and polyphase alpine deformation at Sant'Andrea di Cotone (eastern Corsica, France)[END_REF][START_REF] Caron | The transition beetween blue-schist and lawsonite bearing eclogites on the example of corsican metabasalt[END_REF][START_REF] Lahondère | Le métamorphisme éclogitique dans les orthogneiss et les métabasites ophiolitiques de la région de Farinole (Corse)[END_REF][START_REF] Jolivet | Ductile extension in alpine Corsica[END_REF]1991;[START_REF] Fournier | Alpine Corsica Metamorphic Core Complex[END_REF][START_REF] Caron | Metamorphism and deformation in Alpine Corsica[END_REF][START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF][START_REF] Jolivet | Midcrustal shear zones in postorogenic extension: Example from the northern Tyrrhenian Sea[END_REF][START_REF] Ravna | Cold subduction and the formation of lawsonite eclogiteconstraints from prograde evolution of eclogitized pillow lava from Corsica[END_REF]Vitale Brovarone et al., 2011a;[START_REF] Vitale Brovarone | Inherited Ocean-Continent Transition zones in deeply subducted terranes: Insights from Alpine Corsica[END_REF]. To the west, it is separated by the Ostriconi fault (Fig. 1) from the Nappes Supérieures of Balagne that essentially correspond to low-grade Jurassic ophiolites and Cretaceous flysch of Ligurian affinity [START_REF] Beccaluva | K-Ar age determination on some Tethyan ophiolites[END_REF][START_REF] Ohnenstetter | Crystallization and age of zircon from Corsican ophiolitic albitites: consequences for oceanic expansion in Jurassic times[END_REF][START_REF] De Wever | Datations des radiolarites post-ophiolitiques de Corse alpine à l'aide des radiolaires[END_REF].

Tectonometamorphic evolution of the Tenda massif

The Tenda massif consists of a Variscan basement portion overlain by a Meso-Cenozoic sedimentary cover (Fig. 2; [START_REF] Jourdan | Balagne orientale et massif du Tenda (Corse septentrionale): étude structurale, interprétation des accidents et des déformations, reconstitutions géodynamiques[END_REF][START_REF] Rossi | Carte Géologique de France (1/50 000[END_REF]Vitale Brovarone et al., 2013). Its western part is represented by Proterozoic to Paleozoic gneisses and micaschists unconformably overlain by Late-Variscan volcano-sedimentary rocks. The eastern part consists of Late-Carboniferous/Early-Permian granite intrusions [START_REF] Rossi | Caractère volcano-plutonique du magmatisme calco-alcalin composite d'âge Stéphanien supérieur-Permien inférieur en Corse[END_REF]1994a) including, from east to west: the Mt Asto monzogranite, the Casta granodiorite and the Mt Genova leuco-monzogranite. Few sub-volcanic facies such as the Saleccia rhyolite are also present [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF].

The finite structure of the Tenda massif is the result of polyphazed deformation. The first event (D1) is associated with the development of a non-coaxial planar-linear fabric (S1/L1) with a predominant top-to-the-SW sense of shear coeval with burial under the Schistes Lustrés nappe down to blueschists-facies conditions [START_REF] Mattauer | Transverse lineation and large-scale structures related to Alpine obduction in Corsica[END_REF][START_REF] Jourdan | Balagne orientale et massif du Tenda (Corse septentrionale): étude structurale, interprétation des accidents et des déformations, reconstitutions géodynamiques[END_REF][START_REF] Jolivet | Ductile extension in alpine Corsica[END_REF]1991;[START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF][START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Molli | Orogenic processes and the Corsica/Apennines geodynamic evolution: insights from Taiwan[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF][START_REF] Maggi | Feedback between fluid infiltration and rheology along a regional ductile-to-brittle shear zone: The East Tenda Shear Zone (Alpine Corsica)[END_REF][START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF][START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF][START_REF] Gueydan | Sequential extension as a record of Corsica Rotation during Apennines slab roll-back[END_REF]. Three models have been proposed for the post-D1 evolution:

-For [START_REF] Jolivet | Ductile extension in alpine Corsica[END_REF][START_REF] Jolivet | Geometry and kinematics of extension in Alpine Corsica[END_REF], [START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF], [START_REF] Gueydan | Analysis of continental midcrustal strain localization induced by microfracturing and reaction-softening[END_REF] and [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF], a D2 event occurred mainly as a post-orogenic exhumation until semi-brittle conditions. A pervasive deformation reworked former structures. It includes D2-folds and crenulation cleavages associated with a second planar-linear fabric (S2/L2) preserved in more coaxially deformed domains in between top-to-the-NE deformation zones. This deformation shows a large-scale strain gradient toward the core of the ETSZ (Fig. 2; [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. The feldspar-to-phengite reaction is proposed to be the weakening mechanism responsible for strain localization [START_REF] Gueydan | Analysis of continental midcrustal strain localization induced by microfracturing and reaction-softening[END_REF].

-According to [START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF] and [START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF], crenulation cleavages in the ETSZ may be partly interpreted as low-strained domains preserved between zones of localized top-to-the-SW shearing during D1. D2 is associated with the same structures as for the previous model, but is ascribed to syn-convergence exhumation of both the Schistes Lustrés nappe and the Tenda massif. This was accommodated by a basal top-to-the-SW shear zone below and a top-to-the-NE shear zone at the top of the coupled units, with an otherwise internal distributed deformation. The ETSZ was finally reactivated locally by top-to-the-NE shearing in a post-orogenic context and until brittle conditions at the end of exhumation (D3).

- [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF][START_REF] Maggi | Feedback between fluid infiltration and rheology along a regional ductile-to-brittle shear zone: The East Tenda Shear Zone (Alpine Corsica)[END_REF] and [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF] also proposed that exhumation was mostly syn-orogenic but with a still active top-to-the-SW shearing in the ETSZ. Strain gradients from gneisses to phyllonites developed during this phase and are proposed to be due HP-LT conditions are indicated by ubiquitous celadonite-rich phengite mineral assemblages locally including sodic blue amphibole (riebeckite-ferroglaucophane), jadeitebearing aegirine or even clinopyroxene + rutile (Gibbons and Horak, 1994;[START_REF] Tribuzio | Blueschist facies metamorphism of peralkaline rhyolites from the Tenda crystalline massif (northern Corsica): evidence for involvement in the Alpine subduction event?[END_REF][START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF]Rossetti et al., 2015). Corresponding peak-metamorphic conditions fall in the range of 0.8-1.2 GPa for 350-450 °C [START_REF] Molli | Shear zones and metamorphic signature of subducted continental crust as tracers of the evolution of the Corsica/Northern Apennine orogenic system[END_REF][START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF], with maximum temperature rated at 360 ± 50 °C from Raman Spectroscopy on Carbonaceous Material in metasediments of the cover series (Vitale Brovarone et al., 2013). The retrograde path is characterized by an overall temperature decrease during exhumation running parallel to the phengite Si 4+ isopleths and allowing the preservation of highly substituted phengites [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF].

Stratigraphic and radiometric constraints

Stratigraphic data (Fig. 3) constrain the timing of events in the distal zones of the subduction wedge. In northern Corsica, metamorphism affects Bartonian sediments, and is thus post-40 Ma [START_REF] Bézert | Sur l'âge post-bartonien des événements tectono-métamorphiques alpins en bordure orientale de la Corse cristalline (Nord de Corte)[END_REF]. Priabonian (34-37 Ma) sediments are either interpreted as post- [START_REF] Ferrandini | Définition et datation de la Formation de Venaco (Corse) : dépôt d'origine gravitaire d'âge Priabonien[END_REF] or pre-metamorphism [START_REF] Bézert | La déformation progressive de l'Eocène de la région de Corte: nouvelles données pétrostructurales et conséquences pour la tectogenése alpine en Corse[END_REF].

Nappe emplacement occurred after 41 Ma, as shown by reworked Lutetian sediments under the Santa Lucia nappe [START_REF] Rossi | Carte Géologique de France (1/50000), feuille de Corte[END_REF][START_REF] Ferrandini | Définition et datation de la Formation de Venaco (Corse) : dépôt d'origine gravitaire d'âge Priabonien[END_REF] and the Balagne nappe [START_REF] Nardi | Carta geologica della Balagne «sedimentaria» (Corsica) alla scala 1:25 000 e note illustrative[END_REF]. Toward more internal units, metamorphism in the Tenda massif is considered to be post-Paleocene (~53 Ma), as Eocene sedimentary rocks (dated by correlation with similar facies in Alpine Corsica) are affected by this metamorphism [START_REF] Rossi | Carte Géologique de France (1/50 000[END_REF].

Finally, the orogenic wedge structure is unconformably overlain by the Middle Burdigalian (~18 Ma) extensional basins of Francardo and St-Florent, respectively in internal and external position [START_REF] Ferrandini | Chronostratigraphie et paléomagnétisme des dépôts oligo-miocènes de Corse : implications géodynamiques pour l'ouverture du bassin liguro-provençal[END_REF][START_REF] Cavazza | The Miocene Saint-Florent Basin in northern Corsica: stratigraphy, sedimentology, and tectonic implications[END_REF].

Radiometric constraints are quite controversial in Alpine Corsica (Fig. 3). The youngest ages of HP-LT metamorphism in the Schistes Lustrés nappe are constrained between 37 and 34 Ma (U/Pb on zircon, [START_REF] Martin | Late Eocene lawsoniteeclogite facies metasomatism of a granulite sliver associated to ophiolites in Alpine Corsica[END_REF]Lu/Hf on garnet and lawsonite, Vitale Brovarone and Herwartz, 2013). These ages are consistent with a minimum HP-LT age of 34 Ma inferred by [START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF] from 40 Ar/ 39 Ar dating on phengite, with older ages up to ~60 Ma interpreted to reflect excess Ar. Taken at face value, the Upper Cretaceous age from Lahondère and Gerrot (1997; Sm/Nd on whole rock, garnet, glaucophane and clinopyroxene) appears questionable and problematic. Increments of exhumation are inferred from 40 Ar/ 39 Ar [START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF], and Rb-Sr [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF] ages on phengite around 25 Ma and 20 Ma, respectively, and interpreted as crystallization ages during retrogressive shearing.

Accurate age constraints of HP-LT metamorphism in the Tenda massif are still lacking while existing radiometric data are scattered. In the ETSZ, the timing of prograde top-to-the-SW shear is constrained between 54 ± 8 Ma and 48 ± 18 Ma by U/Pb data from syn-kinematic rutile, coexisting acmite-phengite, and coatings of oxides or sulphides [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF]. Late-Eocene/Early-Oligocene 40 Ar/ 39 Ar phengite ages have been reported between 37-34 Ma, [START_REF] Mailhé | Application des méthodes de datation 40Ar/39Ar et traces de fission à l'étude du déroulement de l'orogenèse alpine en Corse[END_REF], and 35-31 Ma [START_REF] Jourdan | Balagne orientale et massif du Tenda (Corse septentrionale): étude structurale, interprétation des accidents et des déformations, reconstitutions géodynamiques[END_REF]. These fall in the same range than in situ 40 Ar/ 39 Ar ages on a single HP-LT phengite grain distributed between 36.5 ± 0.4 and 30.8 ± 0.4 Ma (mean age = 35.1 ± 0.4 Ma) and interpreted to date crystallization during top-to-the-SW shearing in HP-LT conditions (see Fig. 1; [START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF]. Phengite step-heating analyses also yield older ages up to 47 Ma with internally discordant spectra potentially reflecting mixing between different reservoirs possibly including extraneous Ar.

Younger, more concordant 40 Ar/ 39 Ar spectra near ~25 Ma were also obtained on highly strained rocks associated with the retrogressive top-to-the-NE motion along the ETSZ, and considered as the final increments of deformation (see Fig. 1; [START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF]. A Rb-Sr study on phengite (see Fig. 1; [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF] provided ages of 30.8 ± 4.2, 31.8 ± 1.7, 31.5 ± 2.5 Ma (top-to-the-SW shear zone), 27.1 ± 0.3 Ma (top-to-the-SW shear zone partially overprinted by semi-brittle top-to-the-NE shearing; not considered further in the text), and 20.2 ± 2.2 Ma (top-to-the-NE shear zone in its most deformed domains). This temporal sequence is considered as recording the end of the two shearing events. Final exhumation stages are recorded by zircon and apatite fission track data (ZFT and AFT). ZFT ages between 25 and 17 Ma in the Tenda massif and the Schistes Lustrés nappe are interpreted to date cooling with no significant age difference from either sides of the ETSZ [START_REF] Fellin | The thermal evolution of Corsica as recorded by zircon fission-tracks[END_REF]. AFT ages between 25 and 10 Ma [START_REF] Cavazza | Apatite fission-track analysis of Neogene exhumation in northern Corsica (France)[END_REF][START_REF] Zarki-Jakni | Cenozoic denudation of Corsica in response to Ligurian and Tyrrhenian extension: Results from apatite fission track thermochronology[END_REF][START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF] are comparable between the different units, indicating no differential cooling between units during Miocene denudation.

The recent 40 Ar/ 39 Ar study of Di [START_REF] Vincenzo | 40Ar-39Ar laser dating of ductile shear zones from central Corsica (France): Evidence of Alpine (middle to late Eocene) syn-burial shearing in Variscan granitoids[END_REF] has provided ages of Ma in the Variscan basement of central Corsica, interpreted to constrain syn-burial crystallization of phengite near 35 Ma, and syn-exhumation re-equilibration after 33 Ma at a depth of ~18 km. These ages raise again the question of the timing of HP-LT metamorphism in the external domains, diversely interpreted as post-or pre-Priabonian (38-34 Ma) by stratigraphic data.

Sampling strategy

The sampling strategy is based on the recent work of [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF] where strain gradients were mapped. They defined 6 strain grades encompassing the different deformation events described above and locally overlapping as a result of spatial-temporal variations in strain distribution. As the S2/L2 and top-to-the-NE shear bands are considered as the main rocks fabrics in the ETSZ in their model, the strain grades are defined relative to D2 and are based on structural relationships and mineral deformation/recrystallization textures between overlapping fabrics. This model, together with the other two structural/rheological models presented in section 2.2, will be discussed in the light of the new age data. Strain The strain grade map interpolated in [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF] is reproduced in Fig. 2 to show the regional strain pattern and gradients. Note that strain gradients are heterogeneous such that rocks from different strain grades can be locally encountered in a single strain grade zone.

From the initial sampling, 11 samples visually devoid of any sign of fluid circulation/weathering were selected to avoid secondary complications unrelated to deformation. These are mostly located along the main road and the sea shore where the rocks are the freshest. The distribution of the selected samples (Fig. 2; Table 1) encompasses localized shear zones and associated strain gradients (outcrop and hand-sample scales) throughout the ETSZ while ensuring sufficient regional coverage to allow correlation of deformation-age relationships across the ETSZ. Table 1 provides the correspondence between strain grades and sample numbering. Because it is the most prevalently exposed fabric, sampling was mainly concentrated across D2 strain gradients, except sample TE84, sampled within a D1 top-to-the-SW shear zone (with three zones I, II and III, from core to border) previously studied by [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF]. Sampling was focused on the compositionally homogeneous Mt Asto monzogranite protolith, except for two meta-aplites ©2020 American Geophysical Union. All rights reserved.

(TE09 and TE41), a meta-basite (TE84-zone I), and a meta-rhyolite (TE52), this last from the Saleccia shear zone. Microstructural characteristics and chemical composition of white mica aggregates were checked for all dated samples. Structural formula derived from microprobe analysis were calculated with a minimum of ~50 punctual analyses per sample in order to have a statistically representative estimate of the compositional homogeneity of aggregates.

Furthermore, compositional maps were acquired in order to image spatial variations in chemical composition across aggregates potentially containing several mineral generations.

Analytical procedure

Petrographical observations were carried out on polished sections with an optical microscope and a back-scattered electron imaging detector (BSE) from a Merlin Compact ZEISS scanning electron microscope (SEM) operating at 15 kV (ISTO, University of Orléans, France). Compositional maps were acquired with an associated X-ray energy dispersive system (EDS). Electron microprobe analyses (EMPA) were carried out using a CAMECA SX Five electron microprobe equipped with five wavelength-dispersive spectrometers operating at 15 kV and a current of 6nA (ISTO).

For in situ 40 Ar/ 39 Ar dating, rock sections of 1 cm diameter were drilled from XZ hand-sample sections, then polished and metallized for imaging. Sections were then repolished and washed ultrasonically through three successive ultrasonic baths (acetone, ethanol and pure water), each one followed by drying in a 50 °C stove. No glue was used through the entire protocol. For step-heating 40 Ar/ 39 Ar dating, small-sized (< 2 mm) cohesive white mica aggregates were picked up from XZ hand-sample sections. Unlike standard concentrates homogenized over several cm 3 , these provide a bulk (population-like) age with mm-scale textural control (Fig. 4). After very careful inspection under a binocular microscope to remove any adhering phase/contaminant, these were cleaned through three successive ultrasonic baths (acetone, ethanol and pure water), each one followed by drying in a 50 °C stove, before final optical inspection and weighing.

For irradiation, in situ sections were stacked in a customized Al can with 1.5 mm thick Al disk holders in between each sample and loaded with Fish Canyon sanidine as the irradiation monitor (FCS = 28.02 ± 0.28 Ma, [START_REF] Renne | Intercalibration of standards, absolute ages abd uncertainties in 40Ar/39Ar dating[END_REF]. Each monitor holder contains up to 13 loading positions (one in the center and 6 distributed radially in two concentric circles) to provide 3D neutron flux control. Alignment between polished sections and standard levels was checked before and after irradiation. Aggregates used for the bulkwise step-heating approach were individually wrapped in Al foils and stacked in a dedicated, 5 mm OD diameter, Al can with monitor packets (each loaded with 2 Fish Canyon sanidine grains) interspersed every 4 samples. Simultaneous irradiation of both cans lasted for 5 hours in the Cd-lined CLICIT port of the OSU reactor facility (Corvalis, USA).

In situ analyses were performed with a LSX 213-G2 UV deep-UV CETAC ® (213 nm) pulsed laser operating at 20 Hz and 3.5 mJ/pulse. A minimum spot size of 50 µm (max. 100 µm) was applied to ensure a favorable sample/blank ratio while allowing spatially-resolved inter-grain targeting (mean crystal size was ~50-100 µm, precluding intra-grain variations to be resolved).

Step-heating analyses were conducted with a continuous-wave 25 W Synrad ® CO2 laser [START_REF] Scaillet | Mg/Fe and [(Mg,Fe)Si-Al2] compositional control on argon behaviour in high-pressure white micas: A 40Ar/39Ar continuous laser-probe study from the Dora-Maira nappe of the internal western Alps, Italy[END_REF] defocused to twice as much as the size of the sample lying on a stainless steel holder underneath a differentially pumped ZnS viewport, each step consisting in 30 sec exposure to the beam (at sequentially higher laser-power steps until complete melting).

The extracted gas was purified for 6 minutes using a cold trap held at -127 °C, two air-cooled GP50 SAES getters at 250 °C, and two hot Ta filaments held at 1500 °C, before admission into one of the three static noble gas Helix-SFT Thermo Scientific ® mass-spectrometers operated at ISTO with a mass resolution in excess of 750. Ar and Cl isotopes were sequentially measured in 20 cycles with an electron multiplier ( 35 Cl, 36 Ar ± 1 H 35 Cl, 37 Ar ± 37 Cl, 38 Ar ± 1 H 37 Cl, 39 Ar, 40 Ar isotopes) and a Faraday cup ( 40 Ar only) by peak-switching.

Gas handling and mass-spectrometer operations are fully automated using in-house, standalone LabView-based software interfacing gas extraction and data collection. Procedural blanks were analyzed every sample gas admission (in situ experiments) or every third heating step (step-heating experiments) in the same conditions as the sample analysis. Typical values were 0.1-0.05 fA and ~0.001 fA for m/e = 40 and m/e = 39, 38, 37, 36, respectively. Data regression and age calculations/corrections were made following [START_REF] Scaillet | Numerical error analysis in 40Ar/39Ar dating[END_REF]. Regressed 40 Ar/ 39 Ar isotopic data are tabulated either as 40 Ar/ 39 Ar step-heating analyses and in situ 40 Ar/ 39 Ar analyses in the supporting information. Individual age errors include propagation of all instrumental and procedural uncertainties. Total sample (bulk) ages are reported as integrated (inverse-variance weighted mean) and total-gas (individually summing the Ar isotopes of all steps) ages quoted at ± 1σ.

Microstructural analysis

In all samples, the main metamorphic fabric is composed of alternating greenish micalayers and quartz + feldspar layers (magmatic K-feldspar relics and ubiquitous metamorphic albite) ± chlorite, accessory epidote, Fe-Ti oxides, apatite, monazite, and zircon. Metabasite TE84-zone I stands apart with white mica (> 80 %) and clinopyroxene as main minerals. Sample characteristics, including strain grade, structural domains, type of protolith, mineralogical assemblage, structures, microstructures, grain size with maximum/minimum length along basal section (001), and average referred as fine (<100 µm), medium (100-150 µm) and coarse (>150 µm) and chemical zoning are listed in Table 1.

Grade 0W

Meta-granite TE61 is characterized by S1 and C1 shear bands localized either in white mica or quartz-feldspar aggregates (Fig. 4). Microscopically, S1 is lined by white mica aggregates (Fig. 5a). Crystals are on average coarse-grained, even if fine-grained zones exist.

They show evidence of intracrystalline deformation such as undulose extinction and gentle (001) bending with grain segmentation linked to viscous/brittle behavior (Fig. 5b and5c).

Large grains can be segmented via kink-bands or cracks sub-normal to (001), that sometimes appear serrated, separating slightly misoriented domains, resulting in the reduction of the average grain size (Fig. 5b and5c).

Meta-basite TE84 is composed of three different zones (Fig. 4). Zone I (core of shear zone) is characterized by regular C1/S1 structures in white mica-clinopyroxene aggregates (Fig. 5d). Grains are totally transposed into S1 with an average grain size smaller than 50 µm and show rare evidence of chlorite precipitation along cleavage or along boundaries (see white arrow; Fig. 5d). Zone III (outside the shear zone) displays structures comparable to sample TE61 described above. Zone II is located in between and resembles zone III except that it is characterized by a higher content of white mica and the lack of relic K-feldspar.

Incipient crenulation is seen to affect the D1 fabric as a result of D2 reworking (see Fig. 8 for the outcrop description).

Grade 1

Grade 1 (meta-granites TE38, TE47 and TE62) is characterized by folded S1 microlithons between S2 crenulation cleavage (Figs. 4 and5e). S2 fabric is planar and defined either by a sharp change in S1 orientation or by crystallization of new white mica crystals (Fig. 5f). Rare evidence of chlorite precipitation along cleavage or grain boundaries are found in TE38 and TE62. Bending of S1 into S2 reflects internal deformation (undulose extinction) but is revealed by SEM examination to involve small grains progressively misoriented from each other with no apparent internal deformation. White mica is fine-to medium-grained and grain size is on average smaller and thinner than in Grade 0W for the same protolith (TE61).

Grains marking the S2 are smaller than the grains forming S1 (Fig. 5f).

Grade 3

Grade 3 (meta-granite TE43) is mostly characterized by low-strained zones that preserve crenulation structures between sheared zones displaying C2/S2 structures (Fig. 4).

Crenulated zones have the same microstructural characteristics as Grade 1. In more deformed zones, C2 shear bands lie at shallow angle to S2 (< 10°; Fig. 4), rendering their recognition difficult (Fig. 5g and5h). Here, the grain size of white mica is comparable to previous strain grades but crystals appear almost totally transposed.

Grade 4

The meta-granites (TE77, TE86) and meta-rhyolite (TE52) are characterized by homogeneous C2/S2 structures (Fig. 4). Microscopically, white mica aggregates are either transposed in C2 shear bands or define sigmoidal microlithons in between (Fig. 5i and5j).

Grains within shear bands have a coarser grain size (until >400 µm) but a thinner thickness than in microlithons (Fig. 5i and5j). Grains oriented perpendicular to shear bands can be fractured (Fig. 5j).

Meta-aplites (TE09 and TE41) are characterized macro-(Fig. 4) and microscopically (Fig. 5k and5l) by concordant white mica and quartz-feldspar layers lacking clear shearing sense due to decreasing angle between C2 and S2. White mica grains are totally transposed in the aggregates (Fig. 5l).

Average grain size in Grade 4 is similar to the other strain grades, showing that neither systematic nor consistent grain size reduction accompanied the increase in deformation intensity.

White mica composition

Compositional data show that white micas are phengites that broadly plot along the Tschermak substitution line in the Fe + Mg vs Si plot (Fig. 6a). They show highly variable celadonite content, with Si and Fe + Mg ranging between 3.35-3.75 and 0.35-0.95 a.p.f.u., respectively. Analyses are clustered independent from strain grade, but primarily as function of lithology (bulk-rock control). Phengites with the lowest celadonite content (i.e., Si = 3.35-3.45 a.p.f.u.) belong to the relatively Fe-Mg poor meta-aplite TE09, while phengites with the highest celadonite content (i.e., Si = 3.65-3.75 a.p.f.u.) belong to the Fe-Mg rich meta-basite TE84-zone I. All remaining granite-derived samples plot between 3.45 and 3.65 a.p.f.u. On average, the data slightly plot off the ideal Tschermak substitution line (Fig. 6a) with an amount of pyrophyllite component between 0 and 15 % (Fig. 6b) revealed by K contents comprised between 0.70 and 1.05 a.p.f.u. (see table "chemical composition of white mica" in supporting information). Two samples show a distinctly higher pyrophyllite component, between 10 and 40 % (meta-rhyolite TE52) and between 20 and 35 % (meta-aplite TE41).

The XMg is comprised between 0.45 and 0.55 in most samples, but four are out of this range (Fig. 6c). Zone III of meta-basite TE84 contains phengites with relatively low mean XMg at 0.36 ± 0.04, while zones I and II display relatively high mean XMg at 0.65 ± 0.06 and 0.60 ± 0.04, respectively. Meta-granite TE86 also has a high mean XMg at 0.59 ± 0.05. X-ray compositional maps show different types of chemical zoning in Fe, Al and K (Fig. 7), listed in Table 1 for all samples. Half of the samples are devoid of significant chemical zoning (Fig. 7a, 7b, 7c and Table 1). Other samples show variable zoning patterns marked by anti-correlated Al-rich (dark grey) zones and Fe-rich (light grey) zones in SEM pictures (Fig. 7). Phengites from Grade 0W can show clear core-rim zonations (Fig. 7d,7e,7f and Table 1). A general trend can be seen between Al-rich cores and Fe-or Mg-(sometimes Si-) rich rims, with a possible decorrelation between Fe and Mg. Crenulated zones in Grades 1-3 display correlation between phengite composition and textural habitus with S2-micas distinctly showing an apparent depletion in K, Fe, and an enrichment in Al (Fig. 7g,7h,7i and Table 1). Chemical zoning of phengite is not systematic in these samples (see Table 1).

When present, it is irregular in the more strained ones (Grades 3-4) with no clear core-rim trend (Fig. 7j). In these cases, it is systematically marked by Al-rich (dark grey) zones and Fe-rich (light grey) zones, such as in Grade 0W (Fig. 7j and7k). Samples displaying a marked pyrophyllite component show the same Al-Fe zoning along with an interfering K zoning, also irregular and not systematically correlated with other elements (Fig. 7l, 7m, 7n).

40 Ar/ 39 Ar results

Results are presented in tables " 40 Ar/ 39 Ar step-heating analyses" and " 40 Ar/ 39 Ar in situ analyses" in supporting information along with 38 Ar/ 39 Ar (Cl/K) and 37 Ar/ 39 Ar (Ca/K) ratios not displayed in the figures. These were all close or within error of zero, indicating no significant contribution from Ca-or Cl-rich impurities or contaminating phases. This cross section, along the main road to Saint-Florent, cuts across the entire strain grade evolution (Fig. 8a), from Grade 0W (lower parts and low-strain domains), through Grade 1 (middle parts), and Grade 4 (topmost parts). Data are presented in two figures, one for the (western) basal parts of the cross section (Grades 0W and 1; Fig. 8) and one for the (eastern) top parts (Grade 4; Fig. 9). Grade 0W was not sampled in the lowest part of the section due to extensive weathering but was selected into low-strain boudins within Grade 1 zone (Fig. 8a). Two are from the same outcrop, in the core (meta-basite) and the border (meta-granite) of the top-tothe-SW shear zone partially overprinted during D2 (zone I and III from sample TE84, respectively; Fig. 8b). The third lies about 50 m further east (TE61; meta-granite) and resembles TE84-zone III in terms of structures. Phengites from TE84-zone I (TE84A; Figs. 4 and8d) yielded an oscillating spectrum between 34 and 38 Ma. The total gas age is 36.8 ± 0.2 Ma. Phengites from TE84-zone III (TE84Y; Figs. 4 and8d) also yielded an oscillating spectrum between 28 and 38 Ma with a younger total gas age of 33.9 ± 0.2 Ma. Phengites from TE61, sampled along S1 (TE61B; Fig. 4) were analyzed twice (Fig. 8d). The first run yielded a relatively concordant spectrum with a nearly flat segment between 31 and 32 Ma over 90 % of the total 39 Ar released, followed by increasing ages climbing to ~36 Ma with a total gas age of 32.2 ± 0.1 Ma. The second run yielded a similar shape with older ages and a nearly flat segment between 34 and 35 Ma over 92 % of the total 39 Ar released, followed by an increase to ~37 Ma until the end of the spectrum. The total gas age is 34.6 ± 0.2 Ma.

Two meta-granites from Grade 1, TE38 and TE62, were collected 3 m apart in the same outcrop (Fig. 8c). Phengites from TE38, sampled in folded S1 (TE38C; Figs. 4 and8d) yielded an upward convex spectrum with minimum apparent age at ~32 and 35 Ma, respectively at the beginning and the end of the analysis, and a maximum with a nearly flat segment at ~39 Ma between 25 and 85 % of the total 39 Ar released. The total gas age is 37.8 ± 0.2 Ma. Phengites from TE62 sampled in folded S1 (TE62B; Figs. 4 and8d), also yielded an upward convex age spectrum with minimum apparent ages at ~34 Ma and 37 Ma, respectively at the beginning and the end of the analysis, with a maximum apparent age at ~40 Ma between 58 and 82 % of the total 39 Ar released. The total gas age is 38.6 ± 0.2 Ma.

Toward the northeast, deformation evolves to Grade 4 (Fig. 9a). Two samples were selected (meta-aplite TE41 and meta-granite TE86), separated by ~1 km from each other by a late sinistral brittle fault (see Figs. 2 and9a).

TE41 is an aplite-derived ultramylonite (Fig. 9b). Phengites sampled along C2 planes (TE41D; Fig. 4) were analyzed twice (Fig. 9d). The first run yielded an asymmetric upward convex age spectrum with a minimum apparent age at ~34 Ma at the beginning, increasing after to a maximum age of ~44 Ma between 15 and 35 % of 39 Ar release, then falling back to ~36 Ma at the end. The total gas age is 40.8 ± 0.2 Ma. The second run has a similar shape but rising to ages about 2 Myr younger. The total gas age is nearly the same at 40.3 ± 0.2 Ma. These ages are significantly older than those from samples from the less deformed parts of the same section (indicated by the blue shaded domain in Fig. 9d).

TE86 is a highly deformed meta-granite mylonite (Fig. 9c). Phengites sampled in C2 shear band or transposed S2 (TE86X; Figs. 4 and9d) yielded a relatively concordant age spectrum with a nearly flat segment between ~28 and 30 Ma over the first 95 % of the total 39 Ar released. The total gas age is 29.3 ± 0.1 Ma, younger than Grade 0W and Grade 1 samples of the same cross section (blue shaded area in Fig. 9d).

Cross section 2: Punta di Cepo

This cross section (Fig. 10a) cuts across a large-scale strain gradient including the bottom (meta-granite TE47, Grade 1) and the top of the ETSZ (meta-granite TE43, Grade 3).

TE47 is characterized by recumbent microfolds of S1 with a larger amplitude (up to 3 cm) than other Grade 1 samples (generally < 1 cm), resulting in the obliteration of the top-tothe-SW shearing criteria (Fig. 10b). Phengites sampled in folded S1 (TE47B; Fig. 10b) yielded a 40 Ar loss-like spectrum between ~31 and 36 Ma from 3 to 24 % of the total 39 Ar released, followed by a flat segment at ~36 Ma until 98 % of the total 39 Ar released (Fig. 10d). The total gas age is 35.3 ± 0.2 Ma.

TE43 is characterized by crenulated domains traversed by top-to-the-NE shear zones (Fig. 10c). Four phengite aggregates were sampled, two in crenulated domains (TE43A and TE43B), and two in shear zones (TE43B and TE43C). They all provide ages on average 8 Myr younger than aggregate TE47B (Fig. 10d). TE43A was analyzed twice. The first run yielded a 40 Ar loss-like spectrum starting at ~26 Ma and stabilizing around ~27-28 Ma between 10 % and 96 % of the 39 Ar released. The total gas age is 27.8 ± 0.1 Ma. The second run also yielded a 40 Ar loss-like spectrum starting at ~23 Ma stabilizing between ~26 and 27 between 20 % and 98 % of the 39 Ar released. The total gas age is 26.2 ± 0.2 Ma. TE43B yielded a relatively concordant age spectrum between 27 and 28 Ma from 12 to 94 % of the total 39 Ar released. The total gas age is 27.6 ± 0.1 Ma. TE43C provided an upward convex age spectrum with minimum apparent age at ~25 and 23 Ma, and a maximum apparent age at ~28 Ma near 90 % of the total 39 Ar. The total gas age is 26.9 ± 0.1 Ma. TE43D was analyzed twice. The first run yielded an upward convex age spectrum with minimum apparent age at ~24-24.5 Ma (< 20 % and > 96% of total 39 Ar), and a maximum apparent age at ~25. 5-26.0 Ma in between. The total gas age is 25.3 ± 0.1 Ma. The second run yielded a 40 Ar loss-like pattern between ~25 and 29 Ma. The total gas age is 26.7 ± 0.1 Ma.

Cross section 3: Saleccia and Loto beaches

This cross section (Fig. 11a) exposes two main D2 shear zones where strain is strongly localized, including the Saleccia shear zone (meta-rhyolite TE52; Fig. 11b) and the ETSZ (meta-aplite ultramylonite TE09 and meta-granitic mylonite TE77, Fig. 11c), the last providing a test on the possible influence of variable lithology in a single outcrop (Grade 4).

Phengites from meta-rhyolite TE52 (Saleccia shear zone) was sampled along S2 (TE52C; Fig. 4), and analyzed twice (Fig. 11d). The first run yielded an age spectrum including a first flat segment at ~27.5 Ma followed by variable ages between ~28.5 and 29.5

Ma. The total gas age is 28.6 ± 0.1 Ma. The second run displays a 40 Ar loss-like pattern between ~27 and 30 Ma over the whole gas release. The total gas age is 28.2 ± 0.1 Ma.

Phengites from meta-granite TE77 (ETSZ), sampled along C2 shear band or transposed S2 (TE77C; Fig. 4) were analyzed twice and globally yielded older ages than TE52C (Fig. 11d). The first run yielded a 40 Ar loss-like pattern from ~28 to 31 Ma between 7 and 90 % of the total 39 Ar released, the end of the spectrum being highly discordant between ©2020 American Geophysical Union. All rights reserved.

34 and 39 Ma. The total gas age is 30.8 ± 0.1 Ma. The second run yielded a similar 40 Ar losslike spectrum from ~26 to 31 Ma with a total gas age of 29.3 ± 0.1 Ma.

Two phengite aggregates from the meta-aplite TE09 (ETSZ), sampled along a C2 plane (TE09A and TE09Y; Fig. 4) provided younger ages than TE77C and TE52C (Fig. 11d).

The first one yielded a very discordant age spectrum with a 40 Ar loss-like profile from ~23 and 29 Ma with an overall complex shape. The total gas age is 25.7 ± 0.1 Ma. The second one was more concordant, yet with still a complex shape between 25 and 28 Ma. The total gas age is 26.8 ± 0.1 Ma.

In situ data

In situ 40 Ar/ 39 Ar dating was applied to statistically quantify ages vs. microstructure relationships according to strain grade, and investigate the relationships with the age variations revealed by the corresponding bulk 40 Ar/ 39 Ar step-heating data. Three polished sections from Grades 0W, 1 and 4 (TE84, TE62, TE77, respectively) were analyzed. Data are presented below either in the form of in situ age maps or sorted in situ ages spectra. 7.2.1. Grade 0W (meta-granite TE84) Section TE84 (Fig. 12) was drilled out across zones II and III in a top-to-the-SW localized shear zone (see Fig. 8b). S1 is marked by green phengite aggregates alternating with white quartz-feldspar layers (Fig. 12a). Phengite is more abundant in zone II which lacks relic K-feldspar clasts as found in zone III. 68 analyses were gathered on phengite (43 in zone II, 25 in zone III) and 61 on K-feldspar (zone III). Phengite ages range between 27.2 ± 0.3 and 36.9 ± 2.9 Ma (Fig. 12c). Statistically, older ages are prevalently found in zone II, hosting 84% of the ages older than 32 Ma, while zone III displays 89% of the ages younger than 32 Ma. The spatial age distribution is otherwise random in each sub-zone (Fig. 12a). The first part of the spectrum of sorted ages, dominated by zone III ages, shows discordant ages between ~29 to 33 Ma between 2 and 34 % of cumulative increase (Fig. 12c). The remainder, mostly corresponding to zone II, is more concordant with a weak increase from ~33 to 35 Ma between 34 and 99 % of cumulative increase. In situ ages in K-feldspar define a much broader, almost continuous, age range between 34.0 ± 0.5 and 141.9 ± 0.5 Ma. The spatial age distribution occurs without apparent microstructural control either between and within the clasts, but with older ages (> 77 Ma) clearly associated with the larger crystals (Fig. 12a).

Noteworthy, the younger limit of K-feldspar age distribution nearly coincides with the oldest phengite component found at 35 Ma (Fig. 12b).

Grade 1 (meta-granite TE62)

Section TE62 (Fig. 13) exposes S1 formed by alternating phengite and quartz-feldspar layers, crenulated between S2 cleavage oriented at 80° clockwise from the horizontal of the picture (Fig. 13a). S2 is thicker (20-50 µm) in the bottom left corner of the analyzed zone, and displays recrystallized phengites (see BSE image, Fig. 13b) enriched in Al compared to S1phengite (see Al composition map in Fig. 13c). 100 analyses were acquired on phengite, 50 along S1 and 50 along S2. In situ ages range from 33.4 ± 1.0 to 42.6 ± 0.5 Ma (Fig. 13a). A slight correlation can be observed between age and microstructures with the oldest ages Ma) are evenly distributed among S1 and S2. The spectrum of sorted ages is discordant and can be divided in two parts (Fig. 13d). The first one increases regularly from ~35.5 to 38.5

Ma between 1 and 82 % of cumulative increase. The spectrum then increases nonlinearly from ~38.5 to 43 Ma between 82 and 100 % of cumulative increase. Past 80 %, the spectrum is mostly dominated by S2-tied phengite ages.

Grade 4 (meta-granite TE77)

Section TE77 (Fig. 14) is made of finely alternating phengite and quartz-feldspar layers that experienced significant simple shear. Two orientations can be observed, with the dominant S2 fabric defining the compositional layering and C2 shear bands developed at an angle of 20° (Fig. 14a). Crystals are almost completely transposed in those structures but some of them appear oblique (see BSE image in Fig. 14b and14d). These have a similar chemical composition to layer-parallel phengite grains (see Al composition map in Fig. 14c). 74 analyses were acquired on phengite, 23 in layer-parallel crystals carried by S2, 46 in C2parallel crystals, and 5 in oblique crystals. Ages span from 23.7 ± 1.6 to 33.6 ± 0.5 Ma with no apparent correlation with microstructures (i.e., S2 versus C2). Note however that ages in oblique crystals belong to the oldest interval (30.5-34.0 Ma; Fig. 14e). Spectrum of sorted ages increases regularly from ~27 Ma to 31.5 Ma between 2 and 90 % of cumulative increase (Fig. 14e). Beyond, the apparent ages rise more sharply to ~34 Ma.

Discussion

Potential Ar reservoirs

Before establishing the relationships between deformation, texture and 40 Ar/ 39 Ar ages, the coexistence of possibly different Ar reservoirs must be addressed, especially in the case of the step-heating data from the phengite aggregates. While the principal Ar reservoir is obviously phengite, other components may also contribute to the degassing, including in particular the K-feldspars represented by large igneous porphyroclasts (Fig. 12) and fine (<50 µm) microcline rims around porphyroclasts of metamorphic albite in Grade 0W samples (see [START_REF] Maggi | Feedback between fluid infiltration and rheology along a regional ductile-to-brittle shear zone: The East Tenda Shear Zone (Alpine Corsica)[END_REF]. The incorporation of igneous K-feldspar in phengite aggregates used for step-heating can be excluded due to its large size. Likewise, the presence of metamorphic microcline appears unlikely because of its scarcity (in any case, its contribution should be very minor considering the volume of the analyzed aggregates). This is even more so in the case of the higher strain grades, where the igneous K-feldspar is totally consumed and metamorphic microcline occurs as localized traces. Other potential contaminants are Cl-or Ca-rich minerals including chlorite or epidote (or fluids trapped as inclusions) as found in some samples. These occur mostly as very small, rare and localized grains rather than disseminated crystal spatially associated with mica. Experimental 38 Ar/ 39 Ar (Cl/K) and 37 Ar/ 39 Ar (Ca/K) ratios are close or within error of zero however, ruling out significant contribution from such phases (see supporting information). We infer that the age spectra are largely dominated by the degassing from phengite and that the contribution from subsidiary reservoirs such as K-feldspar or impurities is negligible.

Bracketing constraints on age-deformation relationships

In the light of potential contamination with excess Ar and/or inheritance, establishing whether the observed ages variations are directly interpretable in terms of age-deformation relationships first requires to bracket the temporal range allowed by independent dating constraints in the Tenda massif (see Fig. 3).

A maximum age for alpine metamorphism can be estimated near the end of the Paleocene (~53 Ma), as Eocene sedimentary rocks (dated by correlation with similar facies in Alpine Corsica) are metamorphosed [START_REF] Rossi | Carte Géologique de France (1/50 000[END_REF]. This is consistent with the oldest radiometric ages obtained in the Tenda massif such as the U/Pb-rutile age at 48 ± 18 Ma, together with the aegyrine-phengite U/Pb systematics at 53.7 ± 7.7 Ma of [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF].

Besides, Rb-Sr ages on phengites from top-to-the-SW shear zones constrain the end of D1 at 30.8 ± 4.2, 31.8 ± 1.7, and 31.5 ± 2.5 Ma (see Fig. 1; [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF]. All these data bracket the permissible age span for D1 between 50 and 27 Ma. Concerning D2, [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF] have provided a 20.2 ± 2.2 Ma age from a greenschist-facies, highly sheared top-tothe-NE shear zone, constraining the end of the top-to-the-NE deformation stage (see Fig. 1). This is consistent with the cooling tempo provided by low-temperature thermochronology.

ZFT and AFT yielded scattered ages between 25 and 20 Ma and between 23 and 14 Ma, respectively [START_REF] Cavazza | Apatite fission-track analysis of Neogene exhumation in northern Corsica (France)[END_REF][START_REF] Zarki-Jakni | Cenozoic denudation of Corsica in response to Ligurian and Tyrrhenian extension: Results from apatite fission track thermochronology[END_REF][START_REF] Fellin | The thermal evolution of Corsica as recorded by zircon fission-tracks[END_REF][START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF], indicating a fast cooling event during final exhumation through the closure-T for the two systems (310 and 110°C, respectively; [START_REF] Stockli | Application of Low-Temperature Thermochronometry to Extensional Tectonic Settings[END_REF]. Considering the paucity of brittle structures and the lack of brittle increment of deformation over the main shear zones [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF], these data bracket the end of D2-deformation in the ETSZ between 25 and 20 Ma at the upper limit of the viscous-brittle transition near 300 °C. fracturing, rather than by kink band development or high amplitude folding [START_REF] Wilson | Deformation of biotite and muscovite: optical microstructure[END_REF][START_REF] Bell | Deformation of biotite and muscovite: TEM microstructure and deformation model[END_REF][START_REF] Behrmann | A study of white mica microstructure and microchemistry in a low grade mylonite[END_REF][START_REF] Bell | Kinks in mica: role of dislocations and (001) cleavage[END_REF][START_REF] Mares | Experimental deformation of muscovite[END_REF]. The lattice preferred orientation is stronger in more deformed sample (compare fig.

5g/5h with fig. 5a). This cannot be achieved by dislocation creep because dislocation climb is limited in mica due to the limited number of slip systems. Moreover, slight bending or kinking as observed in Grade 0W is no longer observed in stronger strain grades, indicating that dislocation glide is not an important mechanism. As an example, the folded S1 in Grade 1 is not represented by folded/kinked grains but by a sequence of strain-free grains that are progressively misoriented from each other toward S2 (Fig. 5f). This observation is reported by [START_REF] Bell | Deformation of biotite and muscovite: TEM microstructure and deformation model[END_REF] in biotite where fracturing along and across ( 001) is responsible for grain segmentation without dislocations rearrangement, achieving the same result as recovery (i.e., segmentation of several grains misoriented from each other). The fact that grain thickness (but not grain size, as measured along (001)) is reduced implies a strong component of segmentation along (001). In parallel, aggregates record crystallization of new grains and overgrowth, as shown by grains oriented along S2 (Fig. 5f) or thin but large grains (up to 400 µm) aligned along shear bands (Fig. 5i). This calls for an important participation of dissolution/precipitation processes. These are evidenced by serrated boundaries and the nucleation/growth of mica in boundaries and triple junction of albite and quartz (Fig. 5h,5j, 5l and 7a), attesting to solute transfer by advective transport of components for mica nucleation. Consequently, in Grades 1 to 4, competing grain size growth and reduction mechanisms led to a continuous alignment of (001) along structures (i.e., folded S1, S2 and C2) but without significant reduction in overall grain size distribution (Table1).

Tracking sequential changes in the fabric development through the deformation sequence (and Grades) is rendered difficult due to the isochemical behavior of white mica across during deformation and the lack of other diagnostic phases to delineate discrete changes through the P-T-strain history. The different phengite generations appear compositionally similar with inter-sample variations chemically buffered by the bulk-rock composition rather than by temporal P-T-strain changes (the more ferromagnesian the protolith, the more ferromagnesian the phengite; see Fig. 6). The lack of extra ferromagnesian phase for cation exchange in an overall mineralogically stable assemblage partly explains such a bulk-compositional effect (Table 1). However, mass-balance calculations indicate pervasive fluid-rock exchange in the study area [START_REF] Maggi | Feedback between fluid infiltration and rheology along a regional ductile-to-brittle shear zone: The East Tenda Shear Zone (Alpine Corsica)[END_REF], the phengite composition primarily reflecting the mineralogical transformations attending D1 under the (local) control of protolith variations, P-T conditions, and external fluid availability. These data indicate that the bulk mineralogy, the phengite composition, and the phengite amount were acquired early during D1 after the total consumption of magmatic ferromagnesian minerals. Once consumed, the lack of such reactants precluded further mineralogical exchange, explaining why samples deformed during D2 dominantly consist of aggregates of pre-existing D1 grains. Relic core compositions in core-rim zonations (Fig. 7; Table 1) may be linked to early local equilibrium reflecting early (pre-to syn-D1) mineral transformations.

These microtextures are clear, homogeneous and concentric in D1 samples (Fig. 7d,e,f), indicating that early crystal growth is the main fabric-forming process of mica aggregates at this stage. Conversely, chemical zonation appears patchy, heterogeneous (asymmetric and discontinuous) in D2 samples (Fig. 7j), and is interpreted as truncated relic of early (syn-D1) growth zonation. Accordingly, we infer that recrystallization processes during D2 deformation occurred without major within-grain chemical equilibration but rather via dissolution-precipitation assisted by some grain-boundary fluid. Accordingly, nucleation of new grains does not necessitate a change in chemical composition as a driving force [START_REF] Behrmann | A study of white mica microstructure and microchemistry in a low grade mylonite[END_REF].

Summarizing, textural relationships across the deformation sequence indicate that after a first phase of crystal growth and coarsening during D1, phengite aggregates recorded continuous deformation until almost complete transposition of former fabrics during D2 without significant crystal-scale re-equilibration, nor bulk grain size reduction past Grade 0W. How these mineral deformation and recombination processes interacted with the 40 Ar/ 39 Ar isotopic record is addressed in the next sub-section.

Mechanisms of 40 Ar isotopic record during progressive deformation

Four main mechanisms can be responsible for the age younging with deformation, as evidenced by our data: (1) Ar transfer via fast-pathway pipe diffusion or cross-lattice defects activated at lower temperature than volume diffusion (e.g., [START_REF] Lee | Multipath diffusion in geochronology[END_REF][START_REF] Lo | Argon release mechanisms of biotite in vacuo and the role of short-circuit diffusion and recoil[END_REF][START_REF] Kramar | Heterogeneous 40Ar* distributions in naturally deformed muscovite: in situ UV-laser ablation evidence for microstructurally controlled intragrain diffusion[END_REF]2003;[START_REF] Hames | Single-crystal 40Ar/39Ar age variation in muscovite of the Gassetts Schist and associated gneiss, Vermont Appalachians[END_REF][START_REF] Cosca | 40Ar* loss in experimentally deformed muscovite and biotite with implications for 40Ar/39Ar geochronology of naturally deformed rocks[END_REF]. ( 2) Dynamic sub-grain size refinement reducing the effective diffusion length scale and closure-T (e.g., [START_REF] Goodwin | Effects of progressive mylonitization on Ar retention in biotites from the Santa Rosa mylonite zone, California, and thermochronologic implications[END_REF][START_REF] Dunlap | Argon loss during deformation of micas: constraints from laboratory deformation experiments[END_REF][START_REF] Mulch | In-situ UV-laser 40Ar/39Ar geochronology of a micaceous mylonite : an example of defect-enhanced argon loss[END_REF][START_REF] Cosca | 40Ar* loss in experimentally deformed muscovite and biotite with implications for 40Ar/39Ar geochronology of naturally deformed rocks[END_REF]. [START_REF] Smye | The signature of devolatisation: Extraneous 40Ar systematics in high-pressure metamorphic rocks[END_REF] Partial to full recrystallization inducing variable resetting and mixing with incompletely reset precursory grains (e.g., [START_REF] Scaillet | 40Ar/39Ar laser-probe dating by step heating and spot fusion of phengites from the Dora Maira nappe of the western Alps, Italy[END_REF][START_REF] West | Dating mylonitic deformation by the 40Ar-39Ar method: An example from the Norumbega Fault Zone, Maine[END_REF][START_REF] Kirschner | Staircase 40Ar/39Ar spectra of fine-grained white mica: Timing and duration of deformation and empirical constraints on argon diffusion[END_REF][START_REF] Dunlap | Neocrystallization or cooling? 40Ar/39Ar ages of white micas from low-grade mylonites[END_REF][START_REF] Dunlap | Argon loss during deformation of micas: constraints from laboratory deformation experiments[END_REF][START_REF] Augier | Exhumation, doming and slab retreat in the Betic Cordillera (SE Spain): in situ 40Ar/39Ar ages and P-T-d-t paths for the Nevado-Filabride complex[END_REF]. (4) Crystal-scale dissolution/precipitation replacive changes assisted by an intergranular free fluid phase (fluidassisted diffusional creep, see [START_REF] Scaillet | Excess 40Ar transport scale and mechanism in high-pressure phengites: A case study from an eclogitized metabasite of the Dora-Maira nappe, western Alps[END_REF][START_REF] Scaillet | K-Ar (40Ar/39Ar) geochronology of ultrahigh pressure rocks[END_REF].

As discussed in section 8.3, there is no evidence for higher density of defects at increasingly higher strain grade. Basal parting along (001) dominates over grain segmentation across (001), and cannot bring about a reduction of effective diffusion length (and closure-T) if volume diffusion occurs primarily along the mica interlayer. Thus, deformation-enhanced Ar loss by mechanisms ( 1) and ( 2) cannot explain alone the age variations between the different strain grades, calling for mechanisms ( 3) and ( 4) as the dominant processes to explain the patchy and irregular age distribution and trends observed at all scales. This is consistent with an important participation of dissolution-precipitation in deformation of mica aggregates, as already evidenced in section 8.3 and discussed in recent studies [START_REF] Airaghi | Microstructural vs compositional preservation and pseudomorphic replacement of muscovite in deformed metapelites from the Longmen Shan (Sichuan, China)[END_REF][START_REF] Berger | Microstructures, mineral chemistry and geochronology of micas along a retrograde evolution: An example from the Aar massif (Central Alps, Switzerland)[END_REF]. This may allow compositional re-equilibration at temperatures below which static (solid-state) diffusion of major elements becomes ineffective (e.g., [START_REF] Dempster | Zoning and recrystallization of phengitic micas: implications for metamorphic equilibration[END_REF]. A related form of fluid-phase diffusional creep (or fluid-assisted dissolution-precipitation creep) operating in the ductile regime has already been argued as a major driver for Ar re-mobilization with partial or total resetting of 40 Ar/ 39 Ar ages below the nominal closure-T for Ar loss by diffusion in micas [START_REF] Scaillet | Excess 40Ar transport scale and mechanism in high-pressure phengites: A case study from an eclogitized metabasite of the Dora-Maira nappe, western Alps[END_REF][START_REF] Scaillet | K-Ar (40Ar/39Ar) geochronology of ultrahigh pressure rocks[END_REF].

Multi-factorial analysis of age vs. deformation relationships

In order identify possible cause-and-effect relationships between field observations, finite-strain fabrics, and apparent ages, we devise below a multi-parametric analysis relating all observables with the mapped age trends and local-scale 40 Ar/ 39 Ar systematics.

Probability density distribution (PDD) analysis

Data were plotted as a probability density distribution (PDD) summing the individual Gaussian density of every datum (Fig. 15), providing a continuous kernel density spectrum from which characteristic peaks, or frequency groups, can be tentatively assigned to one (or more) strain grade(s). Because the system studied is temporally dynamic and spatially heterogeneous, we devise below a global trend analysis based on such empirical PDD patterns to correlate the age variations with the observables (structural position, strain grade, grain size and geochemistry) collected via our multi-scale, multi-parametric approach.

In attempting to integrate the observations into a coherent model of age vs. deformation relationships, we must account for the following features.

(1) The majority of the age spectra is variably and erratically discordant with no coherent trend or systematic increase during degassing (Figs. 8 to 11). The lack of systematic correlation between internal discordance and 40 Ar/ 39 Ar spectrum shape prohibits assigning a single mechanism, like partial loss/resetting by diffusion, to the trends. We interpret such spectra as mixed-reservoirs patterns reflecting the differential degassing of structurally-thermally distinct domains rather than diffusion-driven experimental kinetics.

(2) There is an overall agreement in terms of internal discordance between the in situ data and step-heating systematics (compare them in Figs. 12 to 14), indicating that the latter do preserve a faithful record of internal 40 Ar/ 39 Ar gradients.

(3) As noted before, thermal-diffusional effects during the HP-LT cycle were unable to completely erase the pre-Alpine signature of magmatic K-feldspar while the phengite ages are distinctly Alpine in origin in the same lithology; the stark age contrast between these two phases directly matches the overprinting textural relationships, indicating that 40 Ar/ 39 Ar systematics in the phengites are kinetically controlled by recrystallization in the matrix phase rather than by diffusion alone.

Motivated by observation (2), we assume in what follows that any mixing effect of withingrain components (reservoirs) degassed during step-heating is subordinate and that the internal age gradients are preserved such that each individual step-age indirectly relates to a true discrete event (e.g., closure/crystallization or resetting age). Minimal experimental mixing (i.e., averaging) is also suggested by the mutually discordant replicates from the same samples (e.g., TE61B in Fig. 8) indicating that the step-heating data do retain a representative record of internal 40 Ar/ 39 Ar discordance. Accordingly, the 40 Ar* gradients manifested by the discordant age spectra are approximated to reflect distinct (discrete) phengite generations or subdomains with a characteristic closure or (re)crystallization age. Specifically, we evaluate a model whereby the experimental 40 Ar/ 39 Ar spectra reflect the time-integrated, incremental growth/recystallization stages experienced by the samples and that the pooled analysis of the individual age distributions across the sample suite provides a way to extract regionally meaningful trends beyond sample-specific or local effects. The relevance and implications of such assumptions are discussed later.

The full PDD spectrum relative to step-heating ages (Fig. 15, left column) pooling all strain grades displays two peaks culminating near 36 Ma (Mode I) and 27 Ma (Mode II), with second-order satellite peaks erratically distributed throughout. Grade 0W is clustered in the mid of the range (30-40 Ma), with a multi-modal pattern showing no dominant component (Fig 15a). Grade 1 displays a spectrum coinciding with mode I of the pooled spectrum (Fig. 15b). Grade 3 is concentrated at the younger end of the record with a dominantly unimodal and symmetric peak contributing exclusively to Mode II (Fig. 15c). In contrast, Grade 4 is dominantly bimodal and skew-asymmetric with a tail of apparent ages fading away from both modes toward younger ages (Fig. 15d). While the younger maximum coincides and contributes significantly to the pooled Mode II, the second peak at ~41 Ma is distinctly older than Mode I and is unique to this category. This older component is carried by one single sample, TE41, displaying the oldest ages of the suite (Fig. 9).

Applied to the in situ data, the spectral PDD analysis reveals a less noisy pattern over a duration spanning the same length (20-44 Ma, Fig. 15, right column). Grade 0W shows a dominant mode peaking at ~33-34 Ma centrally located relative to the full distribution, with a second order peak near 29-32 Ma (Fig. 15e). Grade 1 plots at the older end with a dominantly unimodal peak at 37-38 Ma with a neutral polarity trend (Fig. 15f). Grade 4 defines the youngest component of the full record sharply peaking at ~28 Ma (Fig. 15g) with a second order peak near 32-34 Ma.

Mixed-coefficient parameter (MCP) analysis

To inquire further into possible age vs. deformation relationships, the PDD analysis can be coupled with mixed-coefficient parameter spectra to illustrate how much a given sample characteristic or category (strain grade, geochemistry, sample location and grain size) is associated to a given pooled age -or temporal bandwidth -across the PDD. MCP spectra (Fig. 16) are computed by weighing pointwise a parameter or category attached to each individual age by the corresponding age-density kernel across the PDD, then by summing over the full age distribution and normalizing by the local PDD. As such, they provide a density-weighted estimate of a given parameter contribution across the PDD spectra using a continuous ordinal scale. While intermediate values across this scale do not strictly represent true (continuous) transitional variations between potentially discrete or categorical values (e.g., like the ordinal strain grade), such a construct potentially provides a semi-quantitative estimate of the mixed-magnitude contribution of a given parameter provided the age depends continuously on that parameter. For instance, strain grades are defined by discrete, ordered, values (categories). Implicit in the transition from Grade 0W through Grade 4 is the concept of cumulative finite strain due to the incremental (continuous) development of the shear fabric. Grades are ordered through a relative scale scoring the magnitude of the associated deformation through such a continuum and the continuous MCP plot smoothes out the transition between these categories according to the age probability attached to a given grade.

For this as for any physically continuous parameter, this display affords a practical way to reveal coupled age-parameter variations (or lack thereof) that would be difficult to correlate across the temporal series otherwise.

Overall, the MCP compositional plots for XMg (Fig. 16a) and K20 (Fig. 16b) indicate no correlation with age variation. One local exception is suggested by the in situ ages in section TE84 displaying older high-XMg (0.65 ± 0.06, ~ 34 Ma) phengites in zone II and younger low XMg (0.36 ± 0.04, 29-34 Ma) phengites in zone III (Fig. 12). Mg-rich micas are empirically considered more retentive than Fe-rich endmembers [START_REF] Harrison | Diffusion of 40Ar in biotite: Temperature, pressure and compositional effects[END_REF][START_REF] Scaillet | Mg/Fe and [(Mg,Fe)Si-Al2] compositional control on argon behaviour in high-pressure white micas: A 40Ar/39Ar continuous laser-probe study from the Dora-Maira nappe of the internal western Alps, Italy[END_REF][START_REF] Dahl | The crystal chemical basis for Ar retention in micas: inferences from interlayer partitioning and implications for geochronology[END_REF] and the older ages at ~34 Ma of the Mg-rich specimens of zone II could potentially point to lower-retention closure ages. However, this effect has been noted for this sample only and cannot explain the regional trends. Such a lack of control by XMg, along with the lack of correlation between grain size and ages (Fig. 16c), strengthen our assumption that Ar lattice (volume) diffusion did not dominate the phengite 40 Ar/ 39 Ar record, as stated earlier (Section 8.4 and 8.5.1). Likewise, a K-loss (pyrophyllite) effect on 40 Ar/ 39 Ar ages could be suggested by sample TE62 (Fig. 13) which displays syn-S2, pyrophyllite-rich, phengites (see depletion in K, Fig. 7g) older than earlier syn-S1 phengites. Such reverse discordance could be attributed to recoil loss of 39 Ar from the K-poor pyrophyllite intermixed within phengite. In the same way as for chloritized biotite, this could bias the apparent age via a combination of 39 Ar recoil into K-poor domains [START_REF] Lo | 39Ar recoil artifacts in chloritized biotite[END_REF][START_REF] Vincenzo | The effect of chlorite interlayering on 40Ar-39Ar biotite dating: an 40Ar-39Ar laser-probe and TEM investigations of variably chloritised biotites[END_REF] with 39 Ar loss from the low-retentive low-K domains during irradiation (e.g., [START_REF] Hess | Kinetics of Ar isotopes during neutron irradiation: 39Ar loss from minerals as a source of error in 40Ar/39Ar dating[END_REF]. We note however that such a correlation is not systematic. For instance, while phengites TE41 and TE52 are characterized by a high pyrophyllite component (Fig. 6b), only the first is associated with apparently older ages relative to the rest of the dataset (see sub-section 8.5. [START_REF] Cavazza | Apatite fission-track analysis of Neogene exhumation in northern Corsica (France)[END_REF]and Figs. 9 and 11). Thus, a systematic and direct correlation between K content and age can also be dismissed.

In contrast, a control of deformation is clearly highlighted by the temporal strain distribution revealed by the MCP analysis. High strain grades (Grades 3 and 4) and low strain grades (Grades 0W and 1) rocks, respectively localized toward the (northeast) upper parts and the (southwest) lower parts of the ETSZ, are respectively associated with the youngest (< 30 Ma) and oldest (> 30 Ma) ages (Figs. 8 to 11). This is clearly expressed by the MCP spectrum of strain grades (Fig. 16d) with the first half (youngest ages) dominated by oscillations between Grades 3 and 4, the second one (oldest ages) dominated by Grades 0W and 1. The MCP plot of the structurally-projected distance (Fig. 16e) of each sample onto a virtual cross-sectional plane oriented SW-NE (with an origin arbitrarily set to the deepest sample) reveals a clear polarity trend with an oscillatory younging toward the northeast.

While the present-day relative position in the structural stack may be locally offset by late faults or may not be representative of a purely cylindrical deformation field, such a pattern is strongly suggestive of the progressive localization of deformation in that direction through the 40-20 Ma timespan. Interestingly, a bulk antipathetic trend is manifested by the Siphengite MCP plot showing an oscillatory-like decrease in Si-content in this age range (Fig. 16f), consistent with progressively shallower lithostatic pressures in the corresponding time interval, although the Si-content may be partly buffered by the host rock composition.

Regional implications

The PDD analysis indicates a deformation continuum between 40 and 20 Ma with two major episodes culminating near ~35-36 Ma (Mode I) and ~27-28 Ma (Mode II) and secondary events scattered throughout (Fig. 15). Compositional and grain size effects, as well as structural inheritance from earlier pre-Alpine (magmatic) fabric, appear quite marginal in controlling the bulk phengite age distribution. The latter appears dominated by the spatiotemporal development of the deformation sequence recorded across the whole Tenda unit.

Mode I and Mode II are respectively represented by Grades 0W-1 and Grades 3-4. These are respectively characterized by S1/L1/top-to-the-SW (and minor S2) and S2/L2/top-to-the-NE structures. These relationships bear on the attribution of D1 and D2 structures respectively to burial and exhumation, and the timing of their transition. Mode I and Mode II are separated by a probability deep at 33-34 Ma in the full PDD spectra (Fig. 15). Noteworthy, part of the internal Schistes Lustrés nappe was in the burial stage close to the peak-pressure at this stage [START_REF] Martin | Late Eocene lawsoniteeclogite facies metasomatism of a granulite sliver associated to ophiolites in Alpine Corsica[END_REF][START_REF] Vitale Brovarone | Timing of HP metamorphism in the Schistes Lustres of Alpine Corsica: New Lu-Hf garnet and lawsonite ages[END_REF]. This was the case of the Tenda massif as well assuming that top-to-the-SW shearing in the ETSZ accommodated burial (rather than exhumation and cooling) contemporaneously in the crustal wedge. That this date is possibly very close to the age of maximum burial is suggested by the minimum 40 Ar/ 39 Ar ages near 34-35 Ma in K-feldspar of sample TE84 (Fig. 12). These are equal to the maximum age of crystallization of the coexisting Alpine phengite and thus correspond to the most severely (if not totally) reset K-feldspar ages during peak-T metamorphic conditions prior to subsequent closure and accumulation of 40 Ar due to the combined effects of temperature decrease during incipient exhumation and cessation of local deformation. On the other hand, the start of exhumation in the most external alpine deformation zones (i.e., Variscan basement of central Corsica) is constrained at ~33 Ma, marking the beginning of exhumation of the whole nappe stack, including the Tenda massif [START_REF] Vincenzo | 40Ar-39Ar laser dating of ductile shear zones from central Corsica (France): Evidence of Alpine (middle to late Eocene) syn-burial shearing in Variscan granitoids[END_REF]. Following the interpretation of [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF] and in light of these new data, D1 structures (S1/L1/top-to-the-SW) and D2 structures (S2/L2/top-to-the-NE) can be respectively attributed to burial and exhumation culminating at ~35-36 Ma (Mode I) and ~27-28 Ma (Mode II), with a switch at ~34 Ma. This is consistent within error bars with the top-to-the-SW/top-to-the-NE switch independently constrained by Rb-Sr (Rossetti et al., 2017). This tectonic inversion/reactivation occurred when Corsica passed from an east-dipping Alpine subduction geometry to a back-arc domain setting with the development of the west-dipping Apennine subduction system near 35-30 Ma [START_REF] Molli | Orogenic processes and the Corsica/Apennines geodynamic evolution: insights from Taiwan[END_REF][START_REF] Malusà | Contrasting styles of (U)HP rock exhumation along the Cenozoic Adria-Europe plate boundary (Western Alps, Calabria, Corsica)[END_REF].The coincidence between the D1/D2 switch and the low probability of the full PDD spectrum at 33-34 Ma (Fig. 15) strongly suggests a low amount of strain and recrystallization at this time consistent with more distributed deformation during the tectonic switch [START_REF] Molli | Deformation and metamorphism at the eastern border of the Tenda Massif (NE Corsica): a record of subduction and exhumation of continental crust[END_REF][START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. The preservation of such an age sequence also suggests minimal mixing of the different age components during incremental 40 Ar/ 39 Ar degassing, as previously argued.

While age mixing is believed to be very limited, we recognize that it can affect to a variable extent the data, primarily in the form of grain-scale disequilibrium brought about by incomplete recrystallization in the absence of efficient thermal drive for Ar diffusion and resetting. The crux observations in this connection are twofold. [START_REF] Cavazza | Apatite fission-track analysis of Neogene exhumation in northern Corsica (France)[END_REF] There is no correlation between extent of deformation (grade) and internal 40 Ar/ 39 Ar scatter. (2) The latter is preserved both as aggregate-scale 40 Ar/ 39 Ar gradients related to aggregate-scale textures (see the relative old oblique S1 crystals in TE77 in Fig. 14) and as grain-scale 40 Ar/ 39 Ar gradients unrelated to grain-scale textures due to dissolution-precipitation replacive changes.

Accordingly, we interpret the progressively reduced maximum and minimum ages with decreasing total-gas age (Fig. 17a) as the effect of a juvenile fabric being progressively formed with a decreasing contribution from early-formed micas. Considering that 40 Ar/ 39 Ar systematics in the phengites ( 1) are kinetically dominated by recrystallization (rather than by diffusion) in the matrix phase and (2) that they reflect mixing of distinct phengite generations, the youngest steps give a maximum age for the end of deformation in a given sample while the oldest ones gives a minimum age for the start of deformation (e.g., [START_REF] Wijbrans | 40Ar/39Ar dating of white micas from an Alpine highpressure metamorphic belt on Naxos (Greece): the resetting of the argon isotopic system[END_REF][START_REF] Beltrando | Dating microstructures by the 40Ar/39Ar step-heating technique: Deformation-pressuretemperature-time history of the Penninic Units of the Western Alps[END_REF].

Accordingly, the ages of the youngest steps recorded herein at ~22 Ma indicate the maximum age for the end of deformation in the ETSZ, which again is slightly older (but still consistent within error bars), with the Rb-Sr ages and interpretation of [START_REF] Rossetti | Pressure-temperature-deformation-time of the ductile Alpine shearing in Corsica: From orogenic construction to collapse[END_REF].

Considering ( 1) that such a young age is consistent with low-temperature thermochronologic data between 300 and 100 °C, and (2) the paucity of brittle structures, these data indicate that the deformation in the ETSZ ended in the uppermost parts of the viscous/brittle transition.

Taken together, the D1/D2 switch at ~34 Ma and the end of D2 at ~22 Ma return a ~12 Myr timelife for the ETSZ as a top-to-the-NE extensional shear zone. This is a minimal duration for regional strain localization. This timelife falls in the range of reasonably known durations of shear zone activity constraints (e.g., [START_REF] Schärer | Duration of strike-slip movements in large shear zones: The Red River belt, China[END_REF][START_REF] Brichau | Constraining the long-term evolution of the slip rate for a major extensional fault system in the central Aegean, Greece, using thermochronology[END_REF][START_REF] Duchêne | Pressure-Temperature-Time Evolution of Metamorphic Rocks from Naxos (Cyclades, Greece): Constraints from Thermobarometry and Rb/Sr dating[END_REF][START_REF] Turrillot | Late orogenic exhumation of the Variscan high-grade units (South Armorican Domain, western France), combined structural and 40Ar/39Ar constraints[END_REF][START_REF] Charles | Timing, duration and role of magmatism in wide rift systems: Insights from the Jiaodong Peninsula (China, East Asia)[END_REF][START_REF] Schneider | Dating the longevity of ductile shear zones: Insight from 40Ar/39Ar in situ analyses[END_REF][START_REF] Beaudoin | The Ikaria hightemperature Metamorphic Core Complex (Cyclades, Greece): Geometry, kinematics and thermal structure[END_REF]. In keeping with previous models [START_REF] Jolivet | Ductile extension in alpine Corsica[END_REF][START_REF] Jolivet | Geometry and kinematics of extension in Alpine Corsica[END_REF][START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF][START_REF] Gueydan | Analysis of continental midcrustal strain localization induced by microfracturing and reaction-softening[END_REF][START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]; see section 2.2), this duration, together with the D2 structures pervading the ETSZ (Fig. 2; [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF], denotes this structure as a major top-to-the-NE post-orogenic extensional shear zone rather than as a local short-lived reactivation.

In terms of strain localization, the spatial distribution of the youngest steps of Grade 3-4 samples (Fig. 17b) suggests that the latest strain increments during D2 migrated progressively toward the northeast upper levels. This regional trend follows the large scale take place decoupled from the local solid-state finite strain record, leaving a texture of undeformed (strain-free) sub-grains or crystal domains formed in an otherwise dynamically forming fabric.

The meaning of the age discordance characterizing such systems may be elusive in settings lacking extensive (in situ and regional) 40 Ar/ 39 Ar sampling coverage as deployed here. In particular: do such 40 Ar/ 39 Ar systematics represent a continuum of deformation or a blurred record of discrete events overlapping in time that are successively overprinted to leave just a vestige of mixed ages? The answer probably lies somewhere in between. The fact that the most deformed samples here are not completely homogenized clearly argues in favor of such 40 [START_REF] Fusseis | Networking of shear zones at the brittle-to-viscous transition (Cap de Creus, NE Spain)[END_REF][START_REF] Fusseis | Micromechanisms of shear zone propagation at the brittleviscous transition[END_REF][START_REF] Pec | Semi-brittle deformation of granitoid gouges in shear experiments at elevated pressures and temperatures[END_REF][START_REF] Marti | Experimental investigation of the brittle-viscous transition in mafic rocks -Interplay between fracturing, reaction, and viscous deformation[END_REF].

However, as it is difficult at small scale to relate ages retrieved in some (micro)structures with the true age of these (micro)structures, along with the mixed ages that might produce intermediate values with no direct geologic meaning, we keep this exploration qualitative.

Conclusion

In the light of our 40 Ar/ 39 Ar multi-scale analysis of the East Tenda shear zone, the following conclusions can be drawn:

1) Intrinsic parameters involved in Ar diffusion including mica composition and grain size do not control the bulk 40 Ar/ 39 Ar ages, indicating a minimal participation of diffusion. Along with the lack of inherited Ar and the general agreement of ages with independent age constraints, we uncover a direct link between deformation/recrystallization processes and the 40 Ar/ 39 Ar record at the regional scale.

2) 40 Ar/ 39 Ar spectra are interpreted as mixed-reservoirs patterns reflecting the differential degassing of structurally-thermally distinct domains rather than diffusion-driven experimental kinetics. The agreement of internal discordance between in situ and step-heating ages suggests that the later preserve a faithful record of internal 40 Ar/ 39 Ar gradients and that individual step ages may be indirectly associated to discrete (re-crystallization) events, strengthening the use of the step-heating technique for tracking the deformation tempo recorded by the mica fabric.

3) The data uncover an isotopic record of deformation lasting for over 20 Myr. Key dates round ~34 Ma and ~22 Ma can be assigned to the end of burial and cessation of deformation during subsequent exhumation, respectively; these allow to infer a ~12 Myr timelife of extensional motion along ETSZ until the Tenda Massif reached the uppermost parts of the viscous/brittle transition. Regionally, younging trends reveal spatio-temporal variations in strain distribution with progressive localization toward the upper contact of the unit and smaller-scale localization in weak lithologies, notably meta-aplites. Thus, methodically, 40 Ar/ 39 Ar dating as deployed here offers interesting perspectives if coupled with rheological purposes.

4) The lack of systematic correlation between i) ages and microstructures and ii) ages and strain grades indicates that the deformation mechanisms controlling 40 Ar/ 39 Ar systematics can be decoupled from local solid-state finite strain record. Along with the effects of mixed ages and partial Ar resetting, this invites utmost caution in relating ages to microstructures and inferring rates of strain localization. [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. The geological outlines in the background are from Rossi et al. (1994a;2001) and [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. Cross sections presented in figures 8, 9, 10 and 11 (step-heating results) are indicated. The asterisks indicate samples also dated by the in situ approach.

Figure 3: Biostratigraphic and radiometric constraints on tectonic and metamorphic events in Alpine Corsica. See text for further explanations. 1) [START_REF] Bézert | Sur l'âge post-bartonien des événements tectono-métamorphiques alpins en bordure orientale de la Corse cristalline (Nord de Corte)[END_REF]; 2) [START_REF] Bézert | La déformation progressive de l'Eocène de la région de Corte: nouvelles données pétrostructurales et conséquences pour la tectogenése alpine en Corse[END_REF][START_REF] Ferrandini | Définition et datation de la Formation de Venaco (Corse) : dépôt d'origine gravitaire d'âge Priabonien[END_REF]3) Nardi et al., 1978;[START_REF] Rossi | Carte Géologique de France (1/50000), feuille de Corte[END_REF][START_REF] Ferrandini | Définition et datation de la Formation de Venaco (Corse) : dépôt d'origine gravitaire d'âge Priabonien[END_REF]4) [START_REF] Rossi | Carte Géologique de France (1/50 000[END_REF]5) [START_REF] Maggi | Clinopyroxene-rutile phyllonites from the East Tenda Shear Zone (Alpine Corsica, France): pressure-temperature-time constraints to the Alpine reworking of Variscan Corsica[END_REF]6) Mailhé, 1982;7) Jourdan, 1988;8) Brunet et al., 2000;9) Rossetti et al., 2015;10) Vitale Brovarone and Herwartz, 2013 ; 11) [START_REF] Martin | Late Eocene lawsoniteeclogite facies metasomatism of a granulite sliver associated to ophiolites in Alpine Corsica[END_REF]12) Di Vincenzo et al., 2016;13) Fellin et al., 2006;14) Cavazza et al., 2001;Zarkni-Jakni et al., 2004;[START_REF] Daniel | Crustal-scale strain partitioning: footwall deformation below the Alpine Oligo-Miocene detachment of Corsica[END_REF]15) [START_REF] Ferrandini | Chronostratigraphie et paléomagnétisme des dépôts oligo-miocènes de Corse : implications géodynamiques pour l'ouverture du bassin liguro-provençal[END_REF][START_REF] Cavazza | The Miocene Saint-Florent Basin in northern Corsica: stratigraphy, sedimentology, and tectonic implications[END_REF] Figure 4: Pictures of samples and in situ location of white mica population sampled for stepheating experiments. Each analyzed aggregate is identified by a bold letter. The scale and the orientation is the same for all pictures. Some dated aggregates were not sampled directly from the rock sections as displayed but a few mm apart (as indicated by the symbol outside the section limits).

  to transient rheology with softening due to feldspar-to-micas reaction and hardening by neoblastesis of alkali feldspar (albite porphyroclast rimmed by microcline) during fluid metasomatism. Only the final stages of exhumation are recorded during post-orogenic extension, where top-to-the-SW phyllonites are the only structures reactivated semi-brittlely by top-to-the-NE shearing.

  grades are summarized below: Grade 0s.s.: undeformed to slightly deformed rocks (white mica crystallization only incipient; no 40 Ar/ 39 Ar sampling), Grade 0W: D1 deformation with S1, L1 and top-to-the-SW shear bands (C1 shear bands). Grade 1: D2 deformation with folded S1 between S2 crenulation cleavage, associated with L2 with a large amount of pure shear component, Grade 2: local top-to-the-NE shear bands (C2 shear bands) in an otherwise crenulationdominated fabric (no 40 Ar/ 39 Ar sampling), Grade 3: top-to-the-NE mylonitic fabric completely transposing S1 with local preservation of D2 folds, Grade 4: higher-strain top-to-the-NE mylonites displaying more tightly spaced (<1 cm), lowangle shear bands (< 30° to S2) relative to Grade 3.

7. 1 .

 1 Step-heating data40 Ar/ 39 Ar age spectra are presented along three W-E to SW-NE cross sections over which structural position of the samples is projected, showing the transition from deep to shallower levels of the Tenda Massif toward the ETSZ (see Fig. 2 for location). Reported synthetic cross sections display the structural-spatial relationships along with outcrop pictures and a synthetic spectrum plot collecting the relevant step-heating results. The same color code is used to identify samples across the different displays. 7.1.1. Cross section 1: southern area

(

  40.5-43.0 Ma) spatially associated with S2-type phengite and the youngest ages (33.0-35.5 Ma) to S1-dominated microlithons. Intermediate ages follow this trend. 70 % of the data fall between 38.0-40.5 Ma and are related to S2 structures. Younger intermediate ages (35.5-38.0

  Ar/ 39 Ar systematics mostly reflecting intermediate ages of mixed origin. The question thus refers to how well one can relate ages with texture/structure. Our data amply demonstrate that such connection is blurred beyond what standard optical, SEM, and UVlaser probe techniques may presently resolve. The level of crystal-isotope interactions in such systems is rooted into the mica interlayer seen as the preferential pathway for chemical reactivity. Combined with grain-scale gradients in fluid activity and transient changes in permeability/porosity, this greatly contributes to make such systems particularly challenging to date. In 300-400°C fluid-present systems deforming near the viscous/brittle transition, the concept of dynamic closure supersedes the classic view that a K-Ar system may be simply defined as closed or open with respect to a static, nominally defined, closure temperature. While necessarily empirical, such a conceptual extension opens useful prospects in terms of 40 Ar/ 39 Ar dating of deformation by incorporating common phenomena such as the catalytic action of fluids, fluid-mineral equilibria, and fluid-rock-deformation feedback relationships that are impossible to rationalize otherwise. In the Tenda case, our integrated 40 Ar/ 39 Ar approach appears to uncover an isotopic record of deformation lasting for as long as 20 Myr in the temporal range recorded by white mica. It permits to put important constraints in the temporal range recorded by white mica, most notably on the ~12 Myr timelife inferred for extensional motion along the ETSZ. Younging trends have permitted to capture and explore spatio-temporal variations in strain localization at different scales, that are required in rheological works, considering that deformation mechanisms can be strain and timedependent

Figure 1 :

 1 Figure 1: Geological context of the studied area. (a) Tectonic map of the western Mediterranean region. (b) Geological map of Alpine Corsica (modified after Vitale Brovarone et al., 2013) associated with a representative cross-section of the Tenda massif (modified after Beaudoin et al., 2017). Previous 40 Ar/ 39 Ar and Rb-Sr ages are projected on the cross section. Ce: Centuri continental unit; Fa: Farinole continental unit; In: Inzecca unit; Ma: Macinaggio unit; Ne: Nebbio unit; Sa: Sampolo unit; SDP: Serra di Pigno continental unit; SL: Santa Lucia; Pi: Pineto ophiolitic unit; Po: Popolasca unit.

Figure 2 :

 2 Figure 2: Deformation intensity map and samples location (modified after[START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. The geological outlines in the background are fromRossi et al. (1994a; 2001) and[START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF]. Cross sections presented in figures 8, 9, 10 and 11 (step-heating results) are indicated. The asterisks indicate samples also dated by the in situ approach.

Figure 5 :

 5 Figure 5: Microstructural observations of white mica aggregates. Strain grade is indicated for each pictures. (a) Lattice preferred orientation with (001) parallel to S1 in Grade 0W. (b) Bending and kinking of (001) associated (top arrow) or not (bottom arrow) with cracking normal to (001). (c) Grain segmentation. (d) C1/S1 structures in clinopyroxene/white mica aggregate in TE84-zone I. (e) Folded S1 between S2 cleavage (arrows) in Grade 1. (f) Recrystallized grains in S2. Note that folding in S1 is achieved through different strain free grains that are progressively misoriented from each other.

Figure 5 :

 5 Figure 5: Microstructural observations of white mica aggregates (continued). (g and h) Transposition of grains in C2 shear zones in Grade 3. (i and j) Coarse grains aligned in C2 shear bands between microlithons domains with oblique grains (white arrow) in Grade 4. (k and l) Almost total transposition of grains in discrete layers in meta-aplites of Grade 4.

Figure 6 :

 6 Figure 6: Chemical composition of white mica. Data are differentiated by colors and marker shapes, indicating both samples and strain grade (see legend). (a) Fe + Mg vs Si plot. Data are broadly aligned along the Tschermak substitution but are clustered by sample, reflecting bulk-rock compositional control. Data plotting off the Tschermak substitution line indicate the existence of other components (substitutions). (b) Muscovite-celadonite-pyrophyllite ternary diagram showing a pyrophyllite substitution, in particular for samples TE41 and TE52 where the pyrophyllite component reaches 20-40 %. (c) Average XMg of phengite. Brackets represent standard deviation.

Figure 7 :

 7 Figure 7: SEM and XRD compositional maps in phengite aggregates. (a), (b) and (c) illustrate homogeneous Al and Fe concentration in some samples (XY sections). (d), (e) and (f) Core (Al-rich)-rim (Fe-rich) zoning in a grain from Grade 0W. (g), (h) and (i) Relative enrichment in Al and depletion in K and Fe in S2-phengite (arrows) compared to S1-phengite in Grade 1. (j) and (k) Al, Mg and Fe zoning that is not associated with core-rim textures anymore (XY section). (l) Heterogeneous distribution of K interfering with other zonations involving Al (m) and Fe (n) in samples with pyrophyllite-rich phengite.

Figure 8 :

 8 Figure 8: Step-heating analyses from Grades 0W and 1 from cross section 1 (see Fig. 2). (a) Cross section showing deformation features and samples location. (b) Outcrop picture of sample TE84 (Grade 0W). Top-to-the-NE indicators show a weak reworking of this thick shear band localized in a meta-basite, but the thin section shows pervasive top-to-the-SW shear bands (see Fig. 5d). (c) Outcrop picture of sample TE38 (Grade 1) showing crenulated structures. (d) Step-heating spectra.

Figure 9 :

 9 Figure 9: Step-heating analyses from Grade 4 from cross section 1 (see Fig. 2). (a) Cross section showing deformation features and samples location. (b) and (c) Outcrop pictures of meta-aplite TE41 and meta-granite TE86 showing parallelization of S2 and C2 shear bands. Few top-to-the-NE indicators remain. (d) Step-heating spectra.

Figure 10 :

 10 Figure 10: Step-heating analyses from cross section 2 (see Fig. 2). (a) Cross section showing deformation features and samples location. (b) Hand-sample picture of sample TE47 (Grade 1) showing a strong crenulation structure. (c) Hand-sample picture of sample TE43 (Grade 3) showing crenulated structures preserved between highly sheared zones where C2 shear bands and S2 tend to parallelize. In situ location of analyzed aggregates are indicated. (d) Stepheating spectra.

Figure 11 :

 11 Figure 11: Step-heating analyses from cross section 3 (see Fig. 2). (a) Cross section showing deformation features and samples location. (b) Outcrop picture of sample TE52 (Grade 4), deformed rhyolite from the Saleccia shear zone. (c) Outcrop picture of sample TE09 (ultramylonitic aplite) and an equivalent of sample TE77 (meta-granite mylonite; Grade 4), sampled few meters away, showing S2/C2 structures. (d) Step-heating spectra.

Figure 12 :

 12 Figure 12: in situ 40 Ar/ 39 Ar dating of phengite and K-feldspar on section TE84, Grade 0W. (a) Distribution of in situ ages, color scaled for phengite and grey scaled for K-feldspar (grains are contoured in red). In situ ages in phengite are older (~32-37 Ma) in zone II, characterized by the highest mean XMg than in zone III (~27-32 Ma), characterized by the lowest mean XMg. Besides, they are randomly distributed in each zone. In situ ages in Kfeldspar are randomly distributed but older (~35-150 Ma) than ages in phengite. (b) Spectra of sorted in situ ages for K-feldspar and phengite. (c) Close-up of the phengite spectrum. In red, the step-heating spectrum of TE84Y (zone III) for comparison.

Figure 13 :

 13 Figure 13: in situ 40 Ar/ 39 Ar dating of phengite on section TE62, Grade 1. (a) Distribution of in situ ages shows statistically older ages (~38-43 Ma) in S2 than in S1 (~33-38 Ma), particularly in the framed bottom left zone (b), the two structures being characterized by slightly different chemical composition exemplified by Al (c), more concentrated in S2. (d) Spectrum of sorted in situ ages. In red, the step-heating spectrum of TE62B for comparison.

Figure 14 :

 14 Figure 14: in situ 40 Ar/ 39 Ar dating of phengite on section TE77, Grade 4. (a) Distribution of in situ ages that are globally randomly distributed in transposed grains in either S2 or C2 shear bands (~25-35 Ma). Oblique grains in un-transposed zone surrounded in red (b and d) are however associated with the oldest ages (~30-35 Ma) without significant variation in chemical composition (c). (e) Spectrum of sorted in situ ages. In red, the step-heating spectrum of TE77C for comparison.

Figure 15 :

 15 Figure 15: Results of the PDD analysis per strain grade for step-heating (left column) and in situ (right column) dating. See text for further explanations.

Acknowledgments and data availability

Federico Rossetti and two anonymous reviewers are thanked for their constructive comments and suggestions that improved the manuscript. This work has received funding from the European Research Council (ERC) under the Seventh Framework Programme of the ©2020 American Geophysical Union. All rights reserved. European Union (ERC Advanced Grant, grant agreement No 290864, RHEOLITH) and the Institut Universitaire de France. It is also contribution of the Labex VOLTAIRE. The 40 Ar/ 39 Ar facility at ISTO was funded and is supported by the ERC Advanced grant RHEOLITH, the LABEX project VOLTAIRE (ANR-10-LABX-100-01), the Région Centre project ARGON, and the project EQUIPEX PLANEX (ANR-11-EQPX-0036).

Datasets for this research are available at: https://zenodo.org/record/4014973

Altogether, these bracketing constraints indicate that the ~20 Myr age range spanned by our (45-22 Ma) and previous 40 Ar/ 39 Ar data (37-34 Ma;[START_REF] Mailhé | Application des méthodes de datation 40Ar/39Ar et traces de fission à l'étude du déroulement de l'orogenèse alpine en Corse[END_REF]35-31 Ma;[START_REF] Jourdan | Balagne orientale et massif du Tenda (Corse septentrionale): étude structurale, interprétation des accidents et des déformations, reconstitutions géodynamiques[END_REF]47-25 Ma;[START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF] potentially represents genuine alpine ages unscathed from excess or inherited 40 Ar. [START_REF] Brunet | Migration of compression and extension in the Tyrrhenian Sea, insights from 40Ar/39Ar ages on micas along a transect from Corsica to Tuscany[END_REF] discussed the possible occurrence of excess argon due to local (in situ) 40 Ar inheritance commonly affecting HP-LT reworked pre-Alpine intrusives (e.g., [START_REF] Monié | Preservation of Hercynian 40 Ar/ 39 Ar ages through high-pressure low-temperature Alpine metamorphism in the Western Alps[END_REF][START_REF] Arnaud | Evidence for excess argon during high pressure metamorphism in the Dora Maira Massif (western Alps, Italy), using an ultra-violet laser ablation microprobe 40Ar/39Ar technique[END_REF][START_REF] Scaillet | K-Ar (40Ar/39Ar) geochronology of ultrahigh pressure rocks[END_REF].

However, this does not seem to apply to our data as shown by the stark age contrast between the matrix phengites and the relic K-feldspars in meta-granite TE84 (Fig. 12). These do not overlap at all, with a neatly defined knickpoint point at ~35 Ma separating the two in the spectrum of sorted ages (Fig. 12). This strongly suggests that the Variscan 40 Ar* component released from magmatic K-bearing minerals was efficiently driven off during partial recrystallization at HP-LT conditions.

Based on these premises, we discuss in the two next sub-sections the mechanisms of mica deformation and their possible interactions on the phengite 40 Ar/ 39 Ar record.

Mechanisms of mica deformation

Despite the intense deformation recorded by the Tenda massif, phengite crystals appear poorly internally deformed (Fig. 5). They dominantly reflect syn-tectonic growth producing newly-formed aggregates with a crystallographic preferred orientation [START_REF] Etheridge | Experimentally produced preferred orientation in synthetic mica aggragates[END_REF][START_REF] Vernon | Microfabric of mica aggregates in partly recrystallized biotite[END_REF], in keeping with studies showing that fabric-forming micas in mylonites are generally strain-free [START_REF] Etheridge | Experimentally produced preferred orientation in synthetic mica aggragates[END_REF][START_REF] Wilson | Deformation of biotite and muscovite: optical microstructure[END_REF][START_REF] Behrmann | A study of white mica microstructure and microchemistry in a low grade mylonite[END_REF]. Indeed, in this configuration, (001) is sub-perpendicular to the maximum stress axis promoting internal deformation by dislocation glide on (001) associated with gently bended/folded (001), kink boundaries between misoriented (001) domains (Fig. 5b) or strain gradient [START_REF] Beaudoin | Deformation behavior of continental crust during subduction and exhumation: Strain distribution over the Tenda massif (Alpine Corsica, France)[END_REF], but does not reproduce it exactly. This is attested by Grade-3 (but not Grade-4) samples populating the youngest peak in the full PDD spectrum (compare Fig. 15c and15d). This is also true at the sample scale. For example, resetting down to 22-24 Ma in sample TE43 is not restricted to the more deformed zones but also occurs in early low-strained crenulated zones (Fig. 10). Thus, it is worth to note that without age constraints, the use of regional strain gradients as a proxy for strain localization in time may be valid at large scale but it may be locally invalid within individual shear zones. Besides the regional strain localization in space and time highlighted by combining structural analysis and 40 Ar/ 39 Ar systematics, the data also indicate that strain localization follows local effects controlled notably by the type of protolith. This is exemplified by the younger ages in meta-aplite TE09 compared to the nearby meta-granite TE77 (Fig. 17b).

Thus, 40 Ar/ 39 Ar dating appears strongly useful to explore strain localization patterns in shear zones characterized by a heterogeneous nature of protolith, as it is difficult in such context to recognize strain gradients, and then infer on strain localization in space. 8.7. Dynamic closure and 40 Ar/ 39 Ar dating of deformation across the viscous/brittle transition The analysis above inherently carries with itself the concept of dynamic closure whereby discrete micas generations may locally record the timing of growth/replacement by dissolution-precipitation, resulting in a mozaic of texturally and isotopically complex crystals bearing no systematic apparent relationship to mineral composition, microstructure, or overprinting textures. Particularly relevant here is the observation that the inferred grain or domain-scale age distribution cannot be directly tied to a specific deformation stage or texture recorded at the sample scale via the in situ approach unless the regional record is properly understood. This is because the processes involved (dissolution-precipitation) can Table 1: Synthesis of structural and petrological characteristics of samples. Grain size is characterized by the minimal and maximal apparent diameter of crystal flakes along with the average grain size parallel to the (001) cleavage in white mica aggregates. "Coarse" designates grain size greater than ~150 µm, "medium" between 100-150 µm, and "fine" smaller than 100 µm and are average grain sizes. Grain size were checked on transversal planes from XZ thin section. It was also checked on XY thin sections for 8 samples in order the get true diameters on basal planes. Since structures are undulated, those thin sections display only few basal sections, that globally give the same diameter as XZ sections, as the apparent aspect ratio on basal plane is near 1:1.

Sample