Hugues Raimbourg

Kristijan Rajič

Benjamin Moris-Muttoni

Vincent Famin

Giulia Palazzin

Donald Fisher

Kristin Morell

Saskia Erdmann

Ida Di Carlo

Clément Montmartin

Quartz vein geochemistry records deformation processes in convergent zones

Key points

1) Variations in Al and Li concentration in hydrothermal quartz span more than 2 orders of magnitude and are up to 3000 and 400 ppm 2) Growth rims with large Al and Li concentrations are present only at seismogenic depth and might reflect large fluid pressure drops 3) Crack-seal microstructures formed at the down-dip limit of the seismogenic zone might be associated with smaller fluid pressure variations

Abstract

In several examples of subduction zones, we compared pairs of quartz veins formed either at the lower temperatures of the seismogenic zone (260°C or below), or at the higher temperatures of its downdip limit (~330°C). All the veins analyzed here are mode I cracks that formed contemporaneously with the host-rock main stage of deformation at peak burial conditions. Lower temperature veins show examples of quartz crystals with euhedral shapes and growth rims, while higher temperature veins contain crack-seal microstructures. In the lower temperature realm, quartz growth rims have alternatingly either: 1) high cathodoluminescence (CL), CL-blue color and high concentration in trace elements and fluid inclusions, or 2) low luminescence, CL-brown color and low concentration in trace elements and fluid inclusions. In contrast, the quartz from higher temperature samples is homogeneously low luminescent and CL-brown, except for very restricted domains of the crack-seal microstructures where patches of CL-blue quartz are present. The highly luminescent quartz contains high concentrations of aluminum and lithium, up to 3000 and 400 ppm, respectively. Variations in Al and Li correlate well, so that Li appears as the main charge-compensating cation for Al. We propose that the incorporation of Al and Li reflects the amplitude of the fluid pressure variations, which control crystal growth rates. Quartz geochemistry might therefore unravel the contrast between the seismogenic zone, where large fluid pressure variations are present, and its downdip limit, where fluid pressure variations are much more limited in amplitude.

Plain language summary

We have examined several examples of veins that formed at large depths in collision and subduction zones. These veins formed as a result of the opening of a crack-shaped cavity and were filled by the precipitation of quartz from dissolved silica. These veins thus have characteristics that reflect the water that was present in the pores and cracks of the rock at depth. The samples we have analyzed come from either the seismogenic zone (for temperatures of the order of ~260°C or lower), i.e.

Introduction

Even if present in small quantities, the fluid in the pores of a rocks exerts a large control on its deformation. An example of such an effect is illustrated by the anticorrelation between coseismic coupling and fluid pressure along the Chile subduction zone [START_REF] Moreno | Locking of the Chile subduction zone controlled by fluid pressure before the 2010 earthquake[END_REF]. Around the downdip limit of the seismogenic zone, episodic tremor and slip (ETS) is shown to be correlated with domains of high fluid pressure [START_REF] Audet | Teleseismic constraints on the geological environment of deep episodic slow earthquakes in subduction zone forearcs: A review[END_REF][START_REF] Audet | Hydrologic control of forearc strength and seismicity in the Costa Rican subduction zone[END_REF][START_REF] Gomberg | Cascadia 2007 and Beyond Working Group (2010) Slow-slip phenomena in Cascadia and beyond: a review[END_REF]Cascadia 2007 andBeyond Working Group, 2010).

Conversely, the effect of deformation of the solid matrix on pore fluid is attested at the surface by variations in water levels in wells and flow rates in springs in the vicinity of earthquake rupture zones [START_REF] Elkhoury | Seismic waves increase permeability[END_REF][START_REF] Jonsson | Post-earthquake ground movements correlated to pore-pressure transients[END_REF][START_REF] Muir-Wood | Hydrological signatures of earthquake strain[END_REF]. Nonetheless, below a few kilometers of depth (i.e., in the depth range where earthquakes nucleate or where slow slip occurs), there are few direct observations of how fluid pressure varies in conjunction with the fluid flow response to the deformation of the solid matrix. Although indirect, the most comprehensive dataset on deep, synkinematic fluids is derived from veins that can be shown to form contemporaneously with the deformation of the surrounding rock [START_REF] Fagereng | Incrementally developed slickenfibers -geological record of repeating low stress-drop seismic events?[END_REF][START_REF] Palazzin | Deformation processes at the down-dip limit of the seismogenic zone: The example of Shimanto accretionary complex[END_REF][START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF]. Some veins trap, a record of the fluid that was present during their growth, either in the vein-forming crystals or in fluid inclusions.

A significant constraint on the fluid pressure conditions at depth can be retrieved from veins through microthermometric studies of fluid inclusions. In convergent settings, fluid pressure changes recorded by concomitant trapping of H2O-rich and CO2-rich fluid inclusions with very contrasting P-V-T properties have been observed in gold-bearing quartz veins and interpreted as a record of fluid pressure variations of the order of ~150 to 200 MPa, at temperature conditions of ~350°C [START_REF] Robert | Gold-quartz veins in metamorphic terranes and their bearing on the role of fluids in faulting[END_REF][START_REF] Wilkinson | Pressure fluctuations, phase separation, and gold precipitation during seismic fracture propagation[END_REF]. In subduction zone settings, fluid pressure variations are also suggested by a large range of isochores of CH4-rich and H2O-rich fluid inclusions [START_REF] Mullis | The system methane-water as a geological thermometer and barometer from the external part of the Central Alps[END_REF][START_REF] Mullis | Rapid subsidence and upthrusting in the Northern Apennines, deduced by fluid inclusion studies in quartz crystals from Porretta Terme[END_REF][START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF]. Assuming growth at a constant temperature (220-230°C and 260-290°C for [START_REF] Mullis | Rapid subsidence and upthrusting in the Northern Apennines, deduced by fluid inclusion studies in quartz crystals from Porretta Terme[END_REF] and [START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF], respectively), the difference in the properties of the fluid inclusions suggests fluid pressure variations of the order of ~150 MPa. These estimates should be considered as order of magnitude, as many reequilibration processes, either by diffusion or by cracking, might affect the record of the fluid inclusions after entrapment [START_REF] Hall | Preferential water loss from synthetic fluid inclusions[END_REF][START_REF] Sterner | Synthetic fluid inclusions -VII. Re-equilibration of fluid inclusions in quartz during laboratory-simulated metamorphic burial and uplift[END_REF][START_REF] Sterner | Compositional re-equilibration of fluid inclusions in quartz[END_REF][START_REF] Vityk | Do fluid inclusions in high-grade metamorphic terranes preserve peak metamorphic density during retrograde decompression?[END_REF][START_REF] Vityk | Synthetic fluid inclusions. XV. TEM investigation of plastic flow associated with reequilibration of fluid inclusions in natural quartz[END_REF][START_REF] Vityk | Fluid inclusions as tectonothermobarometers: Relation between pressure-temperature history and reequilibration morphology during crustal thickening[END_REF].

Fluid pressure conditions at depth can also be characterized based on vein-forming minerals. In goldbearing quartz veins, the co-precipitation of gold and quartz might be explained by an extreme reduction in fluid pressure that leads to flash vaporization of the fluid (Weatherley and Wenley, 2013). Furthermore, crack-seal microstructures have been interpreted as reflecting cyclical processes of deformation and fluid pressure variations [START_REF] Cox | The role of fluids in syntectonic mass transport, and the localization of metamorphic vein-type ore deposits[END_REF][START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF]. Even if other microstructures attest to events of microcracking followed by sealing (Fisher et al., 2019a, b), we restrict here the term crack-seal microstructure to the typical microstructure of antitaxial veins with

Accepted Article

This article is protected by copyright. All rights reserved. solid inclusion bands or trails such as described by [START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF]. Such crack-seal microstructures have been described in subduction settings in the Shimanto Belt, Japan [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF] and on Kodiak Island in Alaska, USA [START_REF] Fisher | Strain variations within an accreted slate belt, Kodiak Islands, Alaska[END_REF]Brantley, 1992, 2014;[START_REF] Fisher | Cyclic fluid flow through a regionally extensive fracture network within the Kodiak accretionary prism[END_REF], and have been tentatively interpreted as the expression of ETS at the base of the seismogenic zone [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF]. Nonetheless, as for fluid inclusions, the interpretation of crack-seal microstructures as reflecting discontinuous pulses of fluid and quartz growth is controversial, as alternative models of continuous growth can account for the formation of bands of fluid inclusions (see the review by [START_REF] Hilgers | On the arrangement of solid inclusions in fibrous veins and the role of the crack-seal mechanism[END_REF]).

The idea of large and episodic fluid pressure changes is also supported by conceptual models of combined fluid flow and deformation, such as the fault-valve model; earthquakes result in a large and sudden increase in permeability, which can lead to large fluid transfer across a seal between fluid compartments [START_REF] Schwartz | Slow slip events and seismic tremor at circum-Pacific subduction zones[END_REF]. Nonetheless, such a model, if valid for the upper part of the crust, might eventually break down at depths greater than the brittle-viscous transition. Indeed, the hydrological system in the viscous, lower part of the crust is fundamentally different from the seismogenic, upper region of the crust [START_REF] Ague | Fluid flow in the deep crust[END_REF][START_REF] Cox | Fluid flow in mid-to deep crustal shear systems: Experimental constraints, observations on exhumed high fluid flux shear systems, and implications for seismogenic processes[END_REF][START_REF] Walther | Volatile production and transport in regional metamorphism[END_REF]. The transitional domain is precisely the one where ETS is generated, and the hydrology of this transitional zone is enigmatic. Based on seismological arguments, [START_REF] Saffer | Hydrogeology and mechanics of subduction zone forearcs: Fluid flow and pore pressure[END_REF] proposed an evolution from a moderately overpressured seismogenic zone to an extremely overpressured domain at larger depths. This transitional zone is also complex in terms of rheology, as this is the domain where brittle and viscous processes might coexist [START_REF] Handy | Seismicity, structure and strength of the continental lithosphere[END_REF].

The objective of the present contribution is to describe the evolution with depth and temperature of the fluid pressure across the brittle-viscous transition in subduction zones. To do so, we used the geochemistry of synkinematic quartz veins chosen at different temperatures spanning this transition as a tool to infer the fluid properties during vein growth. For all the samples, which come from several worldwide examples of subduction zones, we have applied cathodoluminescence, LA-ICP-MS, and EPMA analyses on the veins and Raman spectroscopy on the carbonaceous matter (RSCM) sampled from the host rock. The contrast in geochemical features across the temperature range studied is used to unravel the very large variations in fluid pressure and temperature conditions, between the seismogenic zone and the zone of ETS below.

Geological settings and choice of samples

We collected veins from two accretionary prisms: the Kodiak archipelago in Alaska, USA, and the Shimanto Belt on the island of Kyushu, Japan (Figure 1). The sample GPS locations are given in Table 1.

The Kodiak Accretionary Complex consists of several juxtaposed belts, with ages that increase away from the modern trench. The Kodiak Formation forms a large proportion of the complex and is composed of folded and faulted turbidites of Late Cretaceous age (Moore, 1969a, b;[START_REF] Sample | Structural style and kinematics of an underplated slate belt, Kodiak and adjacent islands, Alaska[END_REF]. It is subdivided into three belts (Landward, Central and Seaward Belts), on the basis of regional variations in bedding and cleavage attitudes. The strongly deformed Central Belt is characterized by numerous southeast-verging thrust faults [START_REF] Fisher | Strain variations within an accreted slate belt, Kodiak Islands, Alaska[END_REF][START_REF] Sample | Duplex accretion and underplating in an ancient accretionary complex, Kodiak Islands, Alaska[END_REF][START_REF] Sample | Structural style and kinematics of an underplated slate belt, Kodiak and adjacent islands, Alaska[END_REF] and contains many syn-kinematic quartz veins [START_REF] Fisher | Strain variations within an accreted slate belt, Kodiak Islands, Alaska[END_REF]Brantley, 1992, 2014;[START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]. On the northwestern part of the islands of Kodiak and Afognak, the coherent rocks of the Kodiak Formation are juxtaposed along the Uganik Thrust [START_REF] Rowe | Fluid-rich damage zone of an ancient out-ofsequence thrust, Kodiak Islands, Alaska[END_REF]. The Uganik Thrust places the Kodiak Formation in fault contact with the Uyak Complex, a mélange composed of lenses of basalt and cherts in a sedimentary matrix of mid-Permian to mid-Early Cretaceous biostratigraphic age, which was accreted prior to the Kodiak Formation during the Late Cretaceous [START_REF] Connelly | Uyak Complex, Kodiak Islands, Alaska: A Cretaceous subduction complex[END_REF]. Similar to the Kodiak Belt, the Uyak Mélange was strongly deformed

Accepted Article

This article is protected by copyright. All rights reserved.

during top-to-the-southeast underthrusting [START_REF] Moore | Structural fabric of a mélange, Kodiak Islands, Alaska[END_REF], and contains abundant veins formed contemporaneously with deformation [START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF]. The two samples studied here come from the Uyak Complex and the Kodiak Central Belt along a transect in the northeastern side of Afognak Island [START_REF] Brantley | Segregation veins: evidence for the deformation and dewatering of a low-grade metapelite[END_REF][START_REF] Fisher | Structural evolution of underthrusted sediments, Kodiak Islands, Alaska[END_REF]. The Kodiak Central Belt contains the most deformed rocks of the Kodiak Formation that are exposed in a structural culmination [START_REF] Fisher | Strain variation in an ancient accretionary wedge: implications for forearc evolution[END_REF]. The vein samples from the Kodiak Central Belt are similar to the ones described in [START_REF] Fisher | Strain variations within an accreted slate belt, Kodiak Islands, Alaska[END_REF]Brantley (1992, 2014), which contain crack-seal microstructures.

The Shimanto Belt is the on-land, Cretaceous to Miocene equivalent of the modern Nankai Trough accretionary prism [START_REF] Taira | The Shimanto Belt of southwest Japan and arc-trench sedimentary tectonics[END_REF][START_REF] Taira | The Shimanto Belt in Shikoku, Japan-Evolution of Cretaceous to Miocene accretionary prism[END_REF][START_REF] Taira | The geology of the Shimanto Belt in Kochi prefecture, Shikoku[END_REF][START_REF] Tanaka | Geology and mineral resources of Japan[END_REF]. On Kyushu Island, the Shimanto Belt near the city of Nobeoka is composed of a stack of northwest-dipping thrust sheets. The Nobeoka Tectonic Line (NTL) is a major fault zone separating the Cretaceous from the Cenozoic part of the belt [START_REF] Murata | Nappe structures of the Shimanto terrane of the Mikado-Osuzuyama area in East Kyushu[END_REF](Murata, , 1997[START_REF] Murata | Low-angle nappe structures of the Shimanto terrane in Kyushu and Shikoku, southwest Japan[END_REF]Raimbourg et al., 2014a). Structurally above the NTL, the Morotsuka Group comprises strongly sheared horizons including deep-sea sediments and basalt lenses [START_REF] Mackenzie | Geochemical study of the greenstones of the Cretaceous and Paleogene Shimanto accretionary complex in eastern Kyushu: Implications for origin and mode of emplacement[END_REF][START_REF] Sakai | Stratigraphy of the Shimanto terrane and the tectono-stratigraphic setting of the greenstones of the northern part of Miyazaki prefecture[END_REF][START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF]. A strong metamorphic foliation is developed in the basal part of the Morotsuka Group, along with a N-S or NNW-SSE stretching lineation [START_REF] Fabbri | Phase ductile à vergence nord dans la zone Shimanto de Kyushu (Japon SW)[END_REF]Raimbourg et al., 2014a;[START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF]. The Hyuga Group forms the footwall of the NTL. Its upper horizon, the Hyuga Tectonic Mélange, is formed of basalt and sandstone lenses in a pelitic matrix. This melange is deformed by a dense network of topto-the-SE shear bands that are interpreted to have formed during subduction (Raimbourg et al., 2014a). Below the Hyuga Tectonic Mélange, a coherent formation composed of turbidites (Coherent Hyuga), is less deformed, except for a few crosscutting fault zones with the same kinematics as the overlying mélange. The density of these faults increases towards its upper boundary, a thrust fault that was possibly activated during a collisional event in the Early Miocene [START_REF] Murata | Duplex structures of the Uchinohae Formation in the Shimanto Terrane, Kyushu, Southwest Japan[END_REF][START_REF] Murata | Nappe structures of the Shimanto terrane of the Mikado-Osuzuyama area in East Kyushu[END_REF][START_REF] Murata | Duplexes and low-angle nappe structures of the Shimanto terrane, southwest Japan[END_REF]Raimbourg et al., 2014a). The temperature structure of the Shimanto Belt on Kyushu was constrained by earlier studies using either vitrinite reflectance [START_REF] Kondo | Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan[END_REF], RSCM (Raimbourg et al., 2014a;Raimbourg et al., 2017a), or illite crystallinity [START_REF] Hara | Metamorphic cooling history of the Shimanto accretionary complex, Kyushu, southwest Japan: Implications for the timing of out-of-sequence thrusting[END_REF][START_REF] Mukoyoshi | Style of fluid flow and deformation in and around an ancient out-of-sequence thrust: An example from the Nobeoka Tectonic Line in the Shimanto accretionary complex, southwest Japan[END_REF]. All these studies show an increase in temperature towards the inner wedge, with gaps corresponding to the main faults. In the zone of study, the temperature in the Coherent Hyuga, Hyuga Tectonic Mélange and Morotsuka Group are ~200°C, ~250°C and ~330-350°C respectively [START_REF] Raimbourg | Fluid properties and dynamics along the seismogenic plate interface[END_REF]Raimbourg et al., 2017a). The samples studied here come from the Morotsuka Group, the Hyuga Tectonic Mélange and the Coherent Hyuga, in order to analyze a relatively large temperature range, from the seismogenic zone to its downdip limit. The veins from the Morotsuka Group are of the same type as those described in [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF] and [START_REF] Nishiyama | Fluid overpressure in subduction plate boundary caused by mantle-derived fluids[END_REF], which contain crack-seal microstructures.

Methods

Structural analysis of veins

For each area investigated, we measured vein geometry, along with the orientation and kinematics of coexisting shear zones or faults. We carried out the structural analysis over entire cross-sections of the thrust belts, to make sure that the veins selected are representative of a regional pattern and can be interpreted within the context of a tectonic process-based framework.

Temperature of veining from Raman Spectroscopy of Carbonaceous Material

The vein samples studied here are all hosted in mudstone-rich rocks, which contain abundant, very fine (<micron-scale) carbonaceous material appropriate for RSCM [START_REF] Beyssac | Raman spectra of carbonaceous material in metasediments: a new geothermometer[END_REF][START_REF] Lahfid | Evolution of the Raman spectrum of carbonaceous material in low-grade metasediments of the Glarus Alps (Switzerland)[END_REF]. Raman spectra were obtained using a Renishaw InVIA Reflex microspectrometer (ISTO-BRGM;

Accepted Article

This article is protected by copyright. All rights reserved.

Orléans). A laser (514nm) was focused on the sample by a DM2500 Leica microscope equipped with an x100 objective. Instrument control and Raman measurements were performed with the software package Renishaw Wire 4.0. Acquisition time was 3s and spectra were accumulated 30 times. To avoid defects on the carbonaceous material related to thin-section preparation, analyses were all performed below the surface of the section by focusing the laser beam beneath a transparent crystal (i.e. dominantly quartz and occasionally calcite or albite). The laser beam power at the surface of the sample was set to ~0.5mW.

Focusing on the range mostly below 350°C of the crystalline evolution of carbonaceous matter, we used the procedure by [START_REF] Lahfid | Evolution of the Raman spectrum of carbonaceous material in low-grade metasediments of the Glarus Alps (Switzerland)[END_REF], which is valid up to 320°C, to decompose the spectra. We used the software PeakFit 4.12 to: (i) substract the baseline (by using linear baseline defined by points in the 900-1000 cm -1 and 1900-2000 cm -1 ranges) and (ii) fit the resulting spectra by 5 Lorentzian functions. When converted into temperatures, some of our spectra yielded temperatures slightly above 320°C, but for the sake of using a consistent approach, we used only the paleothermometer described in [START_REF] Lahfid | Evolution of the Raman spectrum of carbonaceous material in low-grade metasediments of the Glarus Alps (Switzerland)[END_REF].

Beyond the method used to derive paleotemperatures, a further difficulty arises in attributing the inferred temperature to the formation of the veins, as the veins might have precipitated before or after the peak recorded in the RSCM. We consider that paleotemperatures as representative of vein formation when the structural evidence allows one to ascribe the veining stage to the peak in burial and temperature. Additionally, where possible, the RSCM temperature estimates were compared with the minimum temperature of vein fluids determined from homogenization temperatures of fluid inclusions.

Cathodoluminescence imaging

Thick sections of veins were placed in a vacuum chamber at an Argon pressure of ~ 60.10 -3 mbar. They were irradiated by a cold cathode source electron gun provided by OPEA (France). Standard voltagecurrent beam conditions were maintained at ~13-15 kV and 90-110 µA, respectively. The angle between the gun axis and horizontal was 18°. The chamber was placed on a simplified optical mount designed to minimize light absorption on the beam pathway. The optical stage was equipped with a Retiga 2000R 1394a cooled camera (captor size 1 inch, resolution 1600 x1200 pixels, square pixel size 7.4 µm). Given the low luminescence of quartz, acquisition time was fixed to ca. 10s. Whenever present, the strong luminescence of carbonate veins prevented the analysis of quartz.

EPMA analysis of composition

Aluminum concentration in quartz samples was analyzed by compositional mapping with electron probe micro analysis (EPMA). The maps were acquired on Cameca SX Five housed at ISTO-BRGM joint facilities in Orleans. The voltage was set at 15 kV and the current at 200 nA for all maps (Figure 10, Figure 12, Figure 13), but for one where it was set at 250 nA (Figure 11). Step size and dwell time were set at 1 µm and 150 ms for all maps, but for the map in Figure 11 where they were set at 2 µm and 270 ms.

LA-ICP-MS analysis

We carried out our analyses on double-polished thick sections (~150 to 200 µm) of vein material, which were, prior to ablation, removed from the glass slide by immersion for more than 24 h in acetone and then carefully cleaned with acetone for remains of glue. The margins of the thick sections were subsequently fixed to a new glass plate with cyanoacrylate glue. On these sections, we restricted our

Accepted Article

This article is protected by copyright. All rights reserved.

analyses to the domains of clear quartz, i.e. free or with a very low density of fluid inclusions. These domains are rare in the veins, which limits the amount of available material.

The LA-ICP-MS analyses were conducted at the CNRS, Institut de Sciences de la Terre d'Orléans (ISTO), using a RESOlution-SE 193 nm ArF excimer laser with a S155 ablation cell coupled to an Agilent 8900 QQQ mass spectrometer. The ablation chamber of the laser ablation system was flushed with He (350 ml/min) and Ar (15 l/min), to which further Ar (~0.9 l/min) was added on the way to the ICP-MS. The ICP-MS was tuned on NIST612 to low production rates of ThO and doubly-charged 42 Ca of <0.15% and <0.30%, respectively, and to a U/Th ratio of 100±1%.

For the quartz analyses, the laser was operated at a frequency of 10 or 5 Hz, set to deliver a fluence of 8 J/cm 2 on the sample surface, using an analytical spot size 50 µm in diameter. These conditions (especially the lower ablation frequency of 5 Hz) were found to provide good ablation results. Background and ablation signals were acquired for 30 and 40 s, respectively. For each spot, the concentrations of 7 Li, 23 Na, 27 Al, 28 Si and 39 K were determined using dwell times of 50 ms for each isotope. Before and after each set of ~10-12 sample analyses, 3 analyses of NIST610 and 2 analyses of NIST612 have been performed. NIST610 was used as an external standard for calibration of all acquired data. SiO2 was used as an internal standard following the assumption that each point contains close to 100 wt%. NIST 612 was analyzed for quality control, which showed that accuracy and precision were ≤2-5 %.

For data reduction, we used the software GLITTER TM (Van Achterbergh et al., 2001), where we have rigorously inspected and manually selected background and ablation signals, applying a linear fit to all standard data. Despite performing ablation at optimized conditions, initial ablation on quartz commonly caused catastrophic failure (spallation of larger quartz fragments), typically for the first <5 s. Such catastrophic ablation resulted in artificially heightened concentration peaks for Na, K and possibly affecting other elements. In our data reduction, we excluded such concentration peaks, and accordingly excluded the first 5 to 10 s of ablation signal in the integration.

Results

Vein structures/microstructures 1) Structures :

To study comparable structures throughout the different field examples, we focused here on mode I cracks perpendicular to the macroscopic stretching direction. Typical mode I quartz vein structures are shown in Figure 2 and Figure 3 and summarized as a sketch in Figure 4, while the corresponding structural data are presented in Figure 5.

Deformation occurs by simple shear in the Coherent Hyuga and in the Hyuga Tectonic Mélange, with a dense network of shear zones with top-to-the-SE sense of shear (Figure 5). In the Morotsuka Group, the deformation seems mainly coaxial, with a foliation close to horizontal and quartz vein fibers orientated NNW-SSE. Evidence for a contribution of simple shear is nonetheless present in the Morotsuka Group, mostly in the form of en échelon sets of veins indicating top-to-the-S sense of shear (Figure 2B).

In the Kodiak accretionary prism, the Uyak Complex has a spatially variable pattern of deformation. Some outcrops have conjugate sets of shear bands symmetrical about the foliation, indicating NW-SE extension. Other outcrops show a dominant simple shear, with top-to-the-SE shear bands (Figure 5). In the Kodiak Central Belt, the area investigated shows also either top-to-the-NW or top-to-the-SE shear bands filled with quartz and associated with mode I veins. In addition, en-échelon veins are widespread over the area, and indicate a top-to-the-SE sense of shear (Figure 2D). In summary, the strain ellipsoid in all these examples is either simple shear or a combination of pure shear and simple shear (Figure 4). The contribution of pure shear is larger in the higher temperature examples in the Shimanto Belt and the Kodiak accretionary prism. In all the examples, the mode I veins considered are perpendicular, or at a high angle to the foliation, and perpendicular or at high-angle to the X finite strain axis.

In homogeneous, clay-rich lithologies, Mode I veins cut across the metamorphic foliation at high angle, sometimes with an en-échelon geometry, and are also folded and shortened perpendicular to the foliation (Figure 2). In more heterogeneous lithologies, Mode I veins are principally developed in the boudinaged portions of rocks more rigid than the shale matrix, e.g. sandstone layers or lenses embedded in shales, at high angle to the foliation (Figure 3). In many instances, Mode I veins are mode abundant in the vicinity of quartz-filled shear bands (Figure 2C and Figure 3A andB). The shear bands have a low to moderate angle to the foliation. In general, they have a fibrous structure, with fibers parallel to the stretching lineation present on the foliation plane. In some instances, the shear bands are the result of the evolution of an en-échelon set of mode I veins (Figure 2D).

2) Microstructures :

For the present sudy, we subdivided the microstructures of quartz into two categories: growth rims and crack-seals. Growth rim microstructures (Figure 6) were observed in the Coherent Hyuga, in the Hyuga Tectonic Mélange in Japan and in the Uyak Complex in Alaska (Figure 3). The quartz crystals forming the veins occasionally exhibit crystal facets pointing towards the inner axis of the crack.

Under the optical microscope, the growth rims are in most instances underlined by sharp contrasts in the density of fluid inclusions (Figure 6 left). In CL, the growth rims appear as alternations of CL-blue, high luminescence and CL-brown, low-luminescence domains within the crystal (Figure 6 right). In many instances, the CL-blue rims correlate with domains of quartz with a high density of fluid inclusions (Figure 6A, C, D, Figure 7), but other examples show no correlation (Figure 6B). The CL-blue luminescence is evanescent, i.e. it disappears after a few tens of seconds of beam exposure and the CL-color then changes to brown.

Crack-seal microstructures were observed in quartz veins from the Kodiak Central Belt (following [START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]) and from the Morotsuka Group in the Shimanto Belt (following Ujiie et al. ( 2018)) (Figure 8). The crack-seal microstructures that we analyzed are located in mode I veins at high-angle to the foliation S1 (Figure 2).

The crack-seal microstructures consist of bands of solid inclusions between bands of quartz (Figure 8 left). The spacing of these solid inclusion bands is often ~10-20 µm but relatively irregular. The crackseal microstructures have a limited lateral extent and most often are restricted to the interior of a single, elongate quartz grain (Figure 8B,C). In most instances, the quartz grains hosting the crack-seals are elongate perpendicular to the crack-seal bands (Figure 8C).

The quartz crystal surrounding the crack-seal has a homogeneous, CL-brown (or purple) color of luminescence, and has a low intensity of luminescence. The quartz within the crack-seal microstructure has the same luminescence (Figure 8C), except for restricted portions showing a much higher luminescence with CL-blue colors (Figure 8A and B, see also Figure 12 and Figure 13). The solid inclusions within the quartz are non-luminescent.

Temperature of deformation and veining

The paleotemperature data derived from RSCM and associated with the veins are detailed in Table 1, while typical spectra are shown in Figure 9. In the Kodiak Accretionary Complex, additional

Accepted Article

This article is protected by copyright. All rights reserved.

temperature data are provided for the NE transect along NE Afognak Island, near Big waterfall Bay (Table 2).

There is a significant temperature difference, illustrated in the shape of Raman spectra, between the lower and the higher temperature units of the Shimanto Belt and the Kodiak Accretionary Complex. In these two examples, the crack-seal structures with homogeneous brown CL are found in the higher temperature units (T~300-350°C), while the growth rim structures, with contrasts in luminescence, are found in the lower temperature units (T~200-260°C).

EPMA analysis of quartz composition

Lower-temperature microstructures (growth rims):

Quartz vein samples precipitated at temperature below or or on the order of 250°C show growth rims, evidenced by sharp variations in fluid inclusion density and in luminescence (Figure 6). EPMA maps of these growth rims show that the CL-blue, high-luminescence correlates with high aluminum content, so that the growth rims visible in CL are also apparent in Al concentration maps (Figure 10 and Figure 11).

Two growth rim examples from the Shimanto Belt (Coherent Hyuga) and the Kodiak accretionary prism (Uyak Complex) that showed the presence of CL zoning were analyzed by EPMA (see locations in Table 1). These samples show zonings of Al, while all the other elements, such as K or Na, do not show any zoning. As a general rule, there is a much larger contrast in Al content in Shimanto Belt samples (Figure 10) than in Kodiak samples (Figure 11). Despite these differences, in both cases, Al zonings correlate with cathodoluminescence zonings, so that bright blue CL zones correspond to high Al concentration zones and dark brown CL zones to low Al concentration zones (Figure 10 and Figure 11).

The Al-rich rims have a relatively complex structure, with large heterogeneties in Al content. A network of low-luminescence channels isolate regions with a large Al content (Figure 10). In addition, spots of very high Al and K contents are present in the network of low-luminescence channels. These K-and Alrich spots correlate with individual fluid inclusions and might correspond to small clay particles wrapping the cavities.

Higher-temperature microstructures (crack-seal):

Two crack-seal examples were analyzed by EPMA from the Shimanto Belt and the Kodiak accretionary prism, both of which show the presence of patches of high-luminescence, CL-blue quartz, in addition to the dominant, low-luminescence, CL-brown quartz. Both samples display variations in Al concentration that correspond to crack-seal microstructures, whereas other elements like K do not correlate to these microstructures (Figure 12 and Figure 13).

In the two examples, the CL-blue quartz within the crack-seal microstructures has a higher Al content than the surrounding, CL-brown quartz. The contrast in luminescence/Al concentration occurs within a single quartz grain (Figure 12).

In one example, the quartz between the solid inclusion bands has a homogeneous Al concentration (Figure 12). In another example, the luminescence and Al content within the quartz domain of the crack-seal microstructure have a more complex distribution (Figure 13). The zones of blue luminescence and high Al content constitute local patches but not all of the crack-seal quartz. At the higher resolution given by EPMA, contrasts in quartz Al content, in the form of series of bands of variable content, occur within a single quartz band bordered by two successive inclusion bands (Figure 13). The solid inclusions forming the crack-seal inclusion bands are enriched in K and Al and depleted

Accepted Article

This article is protected by copyright. All rights reserved.

in Si compared to the host quartz of the crack-seal bands and correspond to small clay particles (Figure 12 and Figure 13).

LA-ICP-MS analysis of trace elements in quartz

The concentrations of the main trace elements of quartz vein samples from the Shimanto Belt and the Kodiak Accretionary Complex are shown in Table 3 andplotted in Figures 14 and15 As expected from EPMA maps, Al is the main trace element in quartz in terms of concentration, with values as high as ~3000 ppm (all data in wt ppm). The other main trace elements are Li, Na and K, in decreasing order, in much lower concentration than Al. These monovalent cations can be incorporated in quartz as charge-compensating cations for Si-Al substitution. The plot of the monovalent cations as a function of Al shows a good correlation for Li in each of the three accretionary prisms, while Na and K do not seem correlated with Al (Figure 14).

Focusing on Li and Al, there is a strong correlation between these two elements in quartz from both accretionary prisms (Figure 15A). This correlation holds despite significant differences in the Li/Al ratio between each accretionary prism (0.037 and 0.109 in Kodiak and the Shimanto Belt, respectively, Figure 15B).

The total Al or Li content incorporated in quartz is also largely different from one accretionary prism to another (Figure 15). The vein quartz from the Kodiak Accretionary Complex and the Shimanto Belt incorporate Al (Li) up to 445ppm (17ppm) and 2985 ppm (421ppm), respectively (Table 3). The vein quartz from the Shimanto is therefore more enriched in trace elements than the Kodiak Accretionary Complex.

As already seen with EPMA for Al, the luminescence is strongly correlated with the Al and Li content. In each of the three geological examples, the CL-brown, low-luminescence quartz has a much lower Al and Li concentration than the CL-blue, high-luminescence quartz (Figure 15). The comparison of the two examples shows a more complex pattern, as the low-luminescence quartz from the Shimanto Belt has a larger Al and Li concentration than the high-luminescence quartz from Kodiak.

Discussion

A) Tectonic settings of veining in subduction zones

A-1) Temperature of veining An underlying assumption of this work is that the temperature of veining, which is unknown, can be inferred from paleotemperatures determined in the host rock (Table 1). This hypothesis is relevant only if the thermal history is relatively simple, with no more than a single peak in P-T conditions, and if the geometry and kinematics of the veins allow one to connect them, using metamorphic and structural arguments, to the peak P-T conditions of the tectonic evolution.

Higher temperature samples: The veins from the Morotsuka Group in the Shimanto Belt and from the Kodiak Central Belt in Alaska show similar features: mode I veins, sometimes with en-échelon geometry, cutting across relatively homogeneous and fine-grained pelites, and associated with shear bands at low-angle to the metamorphic foliation S1. The veins cut across S1, at the same time they are folded and shortened parallel to the Z-kinematic axis of S1 (Figure 2). In the Kodiak Central Belt, the progressive shearing of the veins parallel to S1 has also been described [START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]. In both the Japanese and Alaskan cases, the veins are therefore cutting across S1 and deformed by it, so the veins and the foliation must be contemporaneous. In addition, in the Morotsuka Group of the This article is protected by copyright. All rights reserved.

Shimanto Belt, S1 is formed of the highest grade metamorphic minerals (actinolite and white mica [START_REF] Toriumi | Tectono-metamorphism of the Shimanto Belt[END_REF]), and therefore is associated with the peak temperature recorded by RSCM. In the case of Kodiak, the Central Belt corresponds both to the lowest structural level and to the highest temperature (Table 2). It is also the zone of highest strain, as apparent with the rotation of the foliation towards horizontal as a result of shear [START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]. It is therefore likely, following previous authors, that the deformation structures observed in the Central Belt of Kodiak (including the metamorphic foliation and the contemporaneous veins), were formed during peak temperature conditions when the unit was at the depths of the accretionary wedge.

Lower temperature samples: In the Hyuga Tectonic Mélange, veining is contemporaneous with the formation of a penetrative network of top-to-the-SE shear zones, kinematically consistent with deformation during subduction (Figure 5) [START_REF] Palazzin | Deformation processes at the down-dip limit of the seismogenic zone: The example of Shimanto accretionary complex[END_REF][START_REF] Raimbourg | Distributed deformation along the subduction plate interface: The role of tectonic mélanges[END_REF]. Furthermore, these shear zones contain chlorite, which is the sole metamorphic mineral. The association between deformation microstructures (including shear zones and mode I veins) and the precipitation of metamorphic minerals suggests that deformation occurred for peak temperature conditions. Another argument stems from the homogenization temperatures of fluid inclusion contained in the veins, which provide a lower bound to the temperature of formation. In the Uyak Complex of the Kodiak accretionary prism and the Hyuga Tectonic Mélange, homogenization temperatures are of the order of 290 +/-20°C [START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF] and in the 200-280°C (Raimbourg et al., 2014b) -175-305°C [START_REF] Kondo | Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan[END_REF] range, respectively. These fluid inclusion data further show that the veins formed for temperature conditions on the order of the maximum temperatures recorded by RSCM on the host rock.

Finally, the strong top-to-the-SE shear deformation and associated pervasive veining observed in the rocks from the coherent Hyuga unit (Figure 5) is restricted to the upper portion of this unit, in the footwall of a fault separating it from the overlying Hyuga Tectonic Mélange (Raimbourg et al., 2014a). This fault contact was activated during a collisional stage in the Early Miocene (Raimbourg et al., 2017a). This stage is the main tectonic event during the Cenozoic and was responsible for thickening of the accretionary prism along several major faults. From this tectonic scenario, we propose therefore that the veining in the Coherent Hyuga unit also occurred at peak conditions of temperature.

In summary, the temperatures estimated by RSCM on host rock in the Shimanto Belt and Kodiak examples are a reliable estimate of vein formation temperatures.

A-2) Relationship between veining and subduction zone deformation

Beyond the temperature of veining, further discussion is needed to associate veining with the different domains of a subduction zone where deformation occurs. As can be inferred from structural data (Figure 5) and from their close association in space (Figures 2 and3), the mode I veins considered here are contemporaneous with the formation of the shear zones (and the foliation in the higher temperatures cases).

The strain pattern is principally simple shear in the lower temperature examples of Hyuga Tectonic mélange and Coherent Hyuga formation, with top-to-the-trench sense of shear (Figure 5). As discussed above, the deformation/veining in the Coherent Hyuga is most likely associated with the motion along the out-of-sequence thrust faults that were activated in Early Miocene (Figure 18) (Fisher et al., 2019c;Raimbourg et al., 2017). A similar interpretation was proposed to account for deformation within the the Hyuga Tectonic mélange, which is located in the vicinity of a major fault, the Nobeoka Tectonic Line, activated during the Early Miocene collisional event [START_REF] Kondo | Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan[END_REF]. Based on mapping arguments that show that the top-to-the-trench deformation within the Hyuga Tectonic mélange predates the Early Miocene deformation stage [START_REF] Murata | Duplexes and low-angle nappe structures of the Shimanto terrane, southwest Japan[END_REF][START_REF] Raimbourg | Distributed deformation along the subduction plate interface: The role of tectonic mélanges[END_REF], an alternative interpretation is that the mélange internal deformation occurred during underthrusting, along the

Accepted Article

This article is protected by copyright. All rights reserved.

subduction interface, similarly to other tectonic mélanges from the Shimanto Belt [START_REF] Kitamura | Dynamic role of tectonic mélange during interseismic process of plate boundary mega earthquakes[END_REF][START_REF] Onishi | Deformation history of tectonic melange and its relationship to the underplating process and relative plate motion: An example from the deeply buried Shimanto Belt, SW Japan[END_REF][START_REF] Raimbourg | Distributed deformation along the subduction plate interface: The role of tectonic mélanges[END_REF]. Similarly, the top-to-the-trench shear in Uyak Complex was interpreted as resulting from deformation during underthrusting [START_REF] Byrne | Evidence for a weak and overpressured decollement beneath sediment-dominated accretionary prisms[END_REF][START_REF] Moore | Structural fabric of a mélange, Kodiak Islands, Alaska[END_REF]. In addition to simple shear, the additional contribution of layerparallel extension to bulk strain in Uyak Complex (Figure 5) has been already reported in the same geological unit [START_REF] Byrne | Evidence for a weak and overpressured decollement beneath sediment-dominated accretionary prisms[END_REF], as well as in other examples [START_REF] Onishi | Deformation history of tectonic melange and its relationship to the underplating process and relative plate motion: An example from the deeply buried Shimanto Belt, SW Japan[END_REF]. It is interpreted as reflecting the effect of a very weak décollement on top of the underthrusted sediments.

In higher temperature examples, the strain pattern includes a larger contribution of coaxial stretching (Figure 4 and4). Nonetheless, top-to-the-trench simple shear is present, principally in the form of enéchelon set of veins, which locally evolve into shear bands (Figure 2). The respective contribution of pure and simple shear is difficult to evaluate.

In the Morotsuka formation, top-to-the-trench deformation is reported along the coast [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF], while more coaxial fabrics are described in the equivalent formation further west [START_REF] Charvet | Vue générale sur l'orogenèse Shimanto et l'évolution tertiaire du Japon sud-ouest[END_REF][START_REF] Fabbri | Phase ductile à vergence nord dans la zone Shimanto de Kyushu (Japon SW)[END_REF]Raimbourg et al., 2014). There, flattening is also apparent in isoclinal folds whose limbs are parallelized to the flat-lying, metamorphic foliation (see the structures of deformation in (Raimbourg et al., 2014)). Strain (including in particular the development of a metamorphic foliation) and metamorphic temperature are highest in the structurally lowest portion of the Morotsuka Group (Raimbourg et al., 2014;[START_REF] Toriumi | Tectono-metamorphism of the Shimanto Belt[END_REF]. We interpreted the peak of deformation and metamorphism as contemporaneous with deformation within the upper plate, near the plate interface. A similar interpretation, i.e. deformation within the upper plate near the décollement, was proposed to account for the veins and associated shear bands in the Kodiak Central Belt [START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]. Compared to their interpretation, our results show, though over a more restricted area, a higher contribution of coaxial extension.

All the examples presented show domains of intense and localized strain, compared to surrounding terranes. The respective locations, within an idealized subduction zone system, of these different examples, is summarized in Figure 18. Uyak Complex and Hyuga Tectonic mélange deformation is associated with underthrusting along the plate interface, Coherent Hyuga deformation with damaging in the footwall of out-of-of-sequence fault zones, Makimine and Kodiak Central Belt deformation with shear and flattening, at higher temperature, at the base of the upper wedge. We propose therefore that veining in these examples is associated with deformation along the main tectonic contacts, across a large depth range including the seismogenic zone and its downdip limit.

Beyond the spatial association between veins and zones of concentrated strain, quartz precipitation has been interpreted as directly resulting from earthquake in tectonic mélanges [START_REF] Kondo | Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan[END_REF][START_REF] Yamaguchi | Dynamic changes in fluid redox state associated with episodic fault rupture along a megasplay fault in a subduction zone[END_REF]Yamaguchi et al., 2012). Similarly, at higher temperature, the two examples of veins showing crack-seal microstructures in the Shimanto Belt and the Kodiak Complex have been proposed as possibly resulting from ETS [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF][START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF]. The general process to connect vein growth to pulses of deformation (for example in the form of earthquakes or ETS) is, as suggested in [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF] and [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF], that deformation of the solid results in dilatant opening of veins and fluid pressure drop, oversaturating the solution in Si and leading to quartz precipitation. In certain cases, the fluid pressure drop resulting from slip on fault is possibly so large as to result in vaporization of part of the fluid (Weatherley and Wenley, 2013).

In summary, the veins considered here were formed within the domains of concentrated deformation in subduction zones, either plate interface or out-of-sequence faults, and the precipitation of quartz within the veins is the result of the deformation events within the surrounding sedimentary rocks.

Accepted Article

This article is protected by copyright. All rights reserved.

B) Evolution in luminescence and microstructure as a function of temperature

In all the geological examples of mode I quartz veins described here, there is large evolution in luminescence with temperature: the luminescent, CL-blue quartz that can be observed in growth rims of euhedral quartz from Mode I veins below 300°C (Figure 6) is almost absent in veins formed at temperature in the range ~300-350°C (Figure 8). In higher temperature samples, the precipitation of luminescent, CL-blue quartz is restricted to a small proportion of crack-seal microstructures (Figure 12 and Figure 13), which otherwise show the same low luminescence and CL-brown colors as the surrounding vein quartz (Figure 8).

A further difference is the fact that crack-seal microstructures (in the sense of antitaxial veins with bands of solid inclusions [START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF]) are observed only in the higher temperature examples, the Makimine Group of the Shimanto Belt [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF] and the Kodiak Central Belt [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF][START_REF] Fisher | Cyclic fluid flow through a regionally extensive fracture network within the Kodiak accretionary prism[END_REF][START_REF] Fisher | The character and distribution of mineralized fractures in the Kodiak Formation, Alaska: Implications for fluid flow in an underthrust sequence[END_REF]. A systematic correspondence between crack-seal microstructures and temperature of veining is difficult to establish, because of the lack of temperature constraints. In Variscan syn-tectonic quartz veins in Aachen that contain crack-seal microstructures, the fluid inclusion homogenization temperatures range from 280 to 390°C [START_REF] Becker | Crack-seal microstructure evolution in bimineralic quartz-chlorite veins in shales and siltstones from the RWTH-1 well, Aachen, Germany[END_REF]. In the quartz-gold veins from Val d'Or, Abitibi, the conditions of metamorphism are lower to middle greenschist-facies conditions. From the analysis of fluid inclusion records, a temperature of formation of ~350°C is inferred [START_REF] Robert | Gold-quartz veins in metamorphic terranes and their bearing on the role of fluids in faulting[END_REF]. One of the classical crack-seal microstructures described in [START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF] comes from the Morcles nappe, but its precise location is not given, and the paleotemperature field is not known in the corresponding area. RSCM data in nearby areas show temperatures between 250 and 350°C [START_REF] Boutoux | Thermal and structural evolution of the external Western Alps: Insights from (U-Th-Sm)/He thermochronology and RSCM thermometry in the Aiguilles Rouges/Mont Blanc massifs[END_REF][START_REF] Girault | The 3-D Thermal Structure of the Helvetic Nappes of the European Alps: Implications for Collisional Processes[END_REF]. The strong thermal gradient observed from the internal to the external part of the mountain belt [START_REF] Girault | The 3-D Thermal Structure of the Helvetic Nappes of the European Alps: Implications for Collisional Processes[END_REF] prevents any precise determination of the temperature of growth of the sample described in [START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF]. Another crack-seal example in [START_REF] Ramsay | The crack-seal mechanism of rock deformation[END_REF] comes from Jurassic limestones from the Windgällen area forming the autochtonous cover of the Aar massif. RSCM temperature estimates in the area yielded temperatures between 300 and 400°C [START_REF] Nibourel | Large vertical displacements of a crystalline massif recorded by Raman thermometry[END_REF].

In spite of these uncertainties, in all examples where a temperature can be associated with confidence to the vein formation, the typical crack-seal structures are formed for temperatures of the order of 300°C or higher.

C) Luminescence and trace element incorporation in quartz

C-1) Origin of the evancescent, blue luminescence in lower-temperature veins

The brown or purple luminescence is commonly observed in metamorphic quartz [START_REF] Götze | Origin, spectral characteristics and practical applications of the cathodoluminescence (CL) of quartz -a review[END_REF]. This is the CL color we observed in all samples, at least for part of the quartz (Figure 6, Figure 8, Figure 10, Figure 11). In addition, a bright, evanescent, CL blue color is specifically observed in growth rims of quartz crystals that grew at temperature of ~250°C or lower (Figure 6, Figure 10, Figure 11), and in very restricted domains of a few crack-seal microstructures observed in quartz that grew at temperature higher than ~300°C (Figure 8, Figure 12, Figure 13).

The short-lived blue luminescence is related to two wavelengths, around 390 and 500 nm [START_REF] Götze | A comparison of optical microscope-and scanning electron microscope-based cathodoluminescence (CL) imaging and spectroscopy applied to geosciences[END_REF][START_REF] Götze | Origin, spectral characteristics and practical applications of the cathodoluminescence (CL) of quartz -a review[END_REF]. This luminescence is typically observed in low-temperature hydrothermal quartz with a large Al content [START_REF] Lehmann | Growth related zonations in authigenic and hydrothermal quartz characterized by SIMS-, EPMA-, SEM-CL-and SEM-CC-imaging[END_REF][START_REF] Perny | Microdistribution of Al, Li, and Na in α quartz. Possible causes and correlation with short-lived cathodoluminescence[END_REF][START_REF] Rusk | Trace elements in hydrothermal quartz: Relationships to cathodoluminescent textures and insights into vein formation[END_REF]. Accordingly, the EPMA maps in our lower temperature samples show a consistent match between CLblue luminescence and high Al content (Figure 10 and Figure 11), confirmed by LA-ICP-MS analyses

Accepted Article

This article is protected by copyright. All rights reserved.

(Figure 15). The alkali-(or hydrogen-) compensated [AlO4/M + ] centre (where M + is a chargecompensating, monovalent cation) is likely responsible for this luminescence.

The correlation between the aluminum concentration and luminescence is nevertheless not very robust. For example, the band at 390nm was detected in ultrapure quartz with <0.1 ppm Al [START_REF] Gorton | Experimental analysis of the composite blue cathodoluminescence emission in quartz[END_REF]. The analysis by [START_REF] Ramseyer | Factors influencing short-lived blue cathodoluminescencoef αquartz[END_REF] of a large range of vein quartz from the Alps shows that only quartz crystals containing more than 50 ppm Al display a short-lived blue cathodoluminescence. Above this threshold, there is no correlation between the intensity of the luminescence and Al concentration. Similarly, [START_REF] Lehmann | Growth related zonations in authigenic and hydrothermal quartz characterized by SIMS-, EPMA-, SEM-CL-and SEM-CC-imaging[END_REF] distinguish a low-luminescent, Alpoor quartz from a high-luminescent, Al-rich one, with a threshold at 500 ppm. In the high Al concentration range, the density of luminescent intrinsic defects (non bridging O-, O-and Si-vacancies) is decorrelated with the Al concentration.

The respective role of the different charge-compensating cations (Na, K, Li, H) is also difficult to determine. On the basis of comparison of different microstructures of quartz, [START_REF] Ramseyer | Factors influencing short-lived blue cathodoluminescencoef αquartz[END_REF] discarded Li+ as possibly causing the luminescence. On the other hand, Na+ and K+ do not correlate with Al concentration in the hydrothermal quartz considered here (Figure 14) or in the quartz crystals from volcanic rocks [START_REF] Watt | Cathodoluminescence and trace element zoning in quartz phenocrysts and xenocrysts[END_REF]. In vein quartz from the Swiss Alps, Li + and H + are the main compensating cations for Al [START_REF] Perny | Microdistribution of Al, Li, and Na in α quartz. Possible causes and correlation with short-lived cathodoluminescence[END_REF]. In our samples, there is a clear correlation between Al and Li (Figure 14) and the luminescence (Figure 15). Compiling the literature data from hydrothermal veins confirms this conclusion (Figure 16). In most samples, high luminescence is correlated with high level of Al 3+ + (Li + /H + ) substitution in quartz. Furthermore, similar to our samples, the high luminescence is restricted to the lower-temperature range of the samples investigated.

C-2) Trace element incorporation in quartz:

Because of charge compensation, the concentrations of the main impurity cations in quartz follow the relationship Li + + H + + K + + P 5+ = Al 3+ + B 3+ [START_REF] Baron | OH-point defects in quartz in Band Li-bearing systems and their application to pegmatites[END_REF][START_REF] Müller | Hydrogen speciation and trace element contents of igneous, hydrothermal and metamorphic quartz from Norway[END_REF]. In high-T high-P experiments by [START_REF] Frigo | OH defects in quartz in granitic systems doped with spodumene, tourmaline and/or apatite: experimental investigations at 5-20 kbar[END_REF], aluminum is the major cation, with a concentration that correlates to some extent with Li, especially at high Al content, with a Li/Al atomic ratio of ~0.6. The large compilation of metamorphic, magmatic and hydrothermal quartz from Norway shows a large scatter in Li vs. Al plots [START_REF] Müller | Hydrogen speciation and trace element contents of igneous, hydrothermal and metamorphic quartz from Norway[END_REF]. Restricting the analysis to hydrothermal quartz, the most abundant elements are Al and Li [START_REF] Müller | Trace elements in quartz -a combined electron microprobe, secondary ion mass spectrometry, laser-ablation ICP-MS, and cathodoluminescence study[END_REF]. Similarly, Al and Li are the main impurities in hydrothermal vein quartz from the Alps [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF]. Our measurements confirm that Al and Li are the cations preferentially incorporated in quartz during low-temperature hydrothermal growth, with concentrations of Al that can reach 3000 ppm (Figure 14).

Additionally, in hydrothermal quartz from the Alps, Al and Li concentrations show a relatively good correlation, with Li/Al atomic ratios between 0.5 and 0.7 or 0.8 [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF]. In our study, the three cases investigated exhibit a very good correlation between Al and Li concentration, but the atomic ratio is quite variable from one case to another (Figure 15). The atomic ratios are between 0.14 to 0.54 (a factor of 3.86 applies to convert wt ratio to atomic ratio), i.e. significantly lower than in [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF]. The plot of all literature data shows that Li/Al atomic ratios are quite variable from one example to another (Figure 16), but in any case lower than one.

Li is not able to compensate for the charge deficit of Si-Al 3+ substitution, while Na and K are in small concentrations and uncorrelated to Al (Figure 14). Therefore, according to the equation above, the rest of the Al 3+ -related charge compensation must result from H + incorporation, as reported in other studies [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF][START_REF] Perny | Microdistribution of Al, Li, and Na in α quartz. Possible causes and correlation with short-lived cathodoluminescence[END_REF]. Nonetheless, the respective incorporation of H + and/or Li + to compensate for the Si-Al charge deficit is dependent on the Al concentration (compare Figure 17A andB). The samples with high Al concentration show the highest Li/Al atomic ratio (Figure respect to H is observed in our samples (compare Figure 15A andB).

In summary, Al is the main impurity incorporated in hydrothermal quartz. The charge imbalance associated with this incorporation is compensated by the uptake of Li and H, with a variable ratio.

D) Interpretation of the luminescence in terms of growth conditions

The coexistence of two types of luminescence in alternating growth rims poses the questions of which physico-chemical parameters control the luminescence and how do these parameters vary during crystal growth.

Neither the pressure, the temperature, or the nature of the solvent in the fluid seem to influence the intensity of the short-lived blue luminescence [START_REF] Ramseyer | Factors influencing short-lived blue cathodoluminescencoef αquartz[END_REF][START_REF] Rusk | Trace elements in hydrothermal quartz: Relationships to cathodoluminescent textures and insights into vein formation[END_REF]. This statement is further justified by the same homogenization temperatures of fluid inclusions contained in CL-blue and CL-brown growth rims [START_REF] Raimbourg | Fluid circulation in the depths of accretionary prisms: an example of the Shimanto Belt, Kyushu, Japan[END_REF].

Other studies have pointed out the influence of the composition of the fluid. Al variable incorporation in quartz might be controlled by changes in the pH of the solution, as it strongly influences Al complexation and its concentration in the solution [START_REF] Mookherjee | Aluminum speciation in aqueous fluids at deep crustal pressure and temperature[END_REF][START_REF] Perny | Microdistribution of Al, Li, and Na in α quartz. Possible causes and correlation with short-lived cathodoluminescence[END_REF][START_REF] Rusk | Trace elements in hydrothermal quartz: Relationships to cathodoluminescent textures and insights into vein formation[END_REF]. But the large variations in pH advocated in [START_REF] Rusk | Trace elements in hydrothermal quartz: Relationships to cathodoluminescent textures and insights into vein formation[END_REF] are relevant in the context of magma-derived hydrothermal alteration and associated ore genesis. In the low grade domains of subduction zones and accretionary prisms, such large changes in pH are less likely, so that the pH does not seem to be the controlling factor for the alternating precipitation of CL-blue and CL-brown quartz.

Alternatively, the growth rate and the growth mode (i.e., stationary or cyclic) have a large influence on the luminescence [START_REF] Götze | Origin, spectral characteristics and practical applications of the cathodoluminescence (CL) of quartz -a review[END_REF][START_REF] Perny | Microdistribution of Al, Li, and Na in α quartz. Possible causes and correlation with short-lived cathodoluminescence[END_REF][START_REF] Ramseyer | Factors influencing short-lived blue cathodoluminescencoef αquartz[END_REF]. [START_REF] Ihinger | Determination of relative growth rates of natural quartz crystals[END_REF] have shown that faster-growing crystal facets incorporate more Al than slower-growing ones and concluded that there is a relationship between Al incorporation and crystal growth rate. This effect of growth rate is further supported, in the samples considered here, by the very large variations in the density of fluid inclusions that correlate with the luminescence: the CL-blue growth rims are much richer in fluid inclusions than the CL-brown ones (Figure 6, Figure 10, see also [START_REF] Raimbourg | Fluid properties and dynamics along the seismogenic plate interface[END_REF][START_REF] Raimbourg | Fluid circulation in the depths of accretionary prisms: an example of the Shimanto Belt, Kyushu, Japan[END_REF]). This contrast in fluid inclusion density points to significantly higher growth rates of the CL-blue rims, which prevents defect-free growth of the crystal.

Additionally, further insights can be gained from the comparison of growth zoning patterns in hydrothermal quartz with those that can be observed in olivine crystallized from a melt. Patterns of phosphate zoning in olivine grown from a melt have been interpreted as the result of rapid growth of the crystal (i.e., surface-controlled growth), overcoming the diffusion of P away from the growing crystal [START_REF] Welsch | Dendritic crystallization: a single process for all the textures of olivine in basalts?[END_REF][START_REF] Welsch | Phosphorus zoning reveals dendritic architecture of olivine[END_REF]. In contrast, low P concentrations are observed in the case of diffusion-controlled growth of olivine, when the growth of the crystal is slow enough for the efficient diffusion of the P away from the crystal. Similarly, framboidal pyrite with a large concentration of impurities are grown from a fluid phase in surface-controlled régime, while euhedral pyrite with a lower impurity content have formed in diffusion-controlled régime [START_REF] Merinero | Nucleation, growth and oxidation of framboidal pyrite associated with hydrocarbon-derived submarine chimneys: lessons learned from the Gulf of Cadiz[END_REF].

The same concept might apply here: the CL-blue quartz grew in surface-controlled régime, with high growth rates and abundant impurities (fluid inclusions, Al+compensating cation substitution for Si), while the CL-brown quartz grew at a much lower rate, in a diffusion-controlled regime, where efficient Al and Li diffusion away from the growing crystal occurs.

Accepted Article

This article is protected by copyright. All rights reserved.

The reason for the differences in growth rate of quartz might be related to variations in super saturation levels, following the models of olivine or pyrite crystals [START_REF] Merinero | Nucleation, growth and oxidation of framboidal pyrite associated with hydrocarbon-derived submarine chimneys: lessons learned from the Gulf of Cadiz[END_REF][START_REF] Welsch | Phosphorus zoning reveals dendritic architecture of olivine[END_REF]. Similarly, [START_REF] Mullis | Fluid inclusion studies during very low-grade metamorphism[END_REF] postulated that the rapid quartz growth was caused by a sudden drop in fluid pressure, causing phase separation in the fluid (the single-phased mixed fluid separates into methane-rich and water-rich fluids) and precipitation of quartz. The equilibrium concentration of dissolved silica decreases with pressure decrease at constant temperature [START_REF] Anderson | The solubility of quartz in supercritical water[END_REF][START_REF] Fournier | An equation correlating the solubility of quartz in water from 25" to 900°C at pressures up to 10,000 bars[END_REF][START_REF] Morey | The solubility of quartz in water in the temperature interval from 25° to 300°C[END_REF]. Therefore, a sudden pressure drop leads to supersaturation of silica and enhanced growth of quartz from the solution.

In conclusion, the formation of Al-and Li-rich, CL-blue quartz is probably related to enhanced growth rates, which involves a large increase in the magnitude of the driving force for precipitation. In turn, the most likely scenario to account for such stages of rapid growth is rapid pressure drops, leading to supersaturation of silica and increased rates of quartz precipitation. E) Model of pore fluid pressure variations in subduction zones E-1) Evidence for fluid pressure variations:

The precipitation of highly luminescent, CL-blue, and Al-Li-rich quartz, is best explained, as seen above, as the result of rapid fluid pressure drops. For cyclical precipitation of the CL-blue growth rims to occur, they must alternate with periods of recovery of the fluid pressure, where CL-brown growth rims would precipitate, at a much slower pace.

Further arguments for fluid pressure variations during precipitation of hydrothermal quartz are provided by trapped fluid inclusions. Several studies have argued for large (up to 150-200 MPa) variations in the paleopressure conditions of the fluid that is contained as inclusions within hydrothermal quartz grown in veins, in low-grade terranes from collision or subduction settings [START_REF] Herrington | Colloidal gold and silica in mesothermal vein systems[END_REF][START_REF] Mullis | Rapid subsidence and upthrusting in the Northern Apennines, deduced by fluid inclusion studies in quartz crystals from Porretta Terme[END_REF][START_REF] Robert | Gold-quartz veins in metamorphic terranes and their bearing on the role of fluids in faulting[END_REF][START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF]. Fluid inclusion data might be questionable, because post-entrapment reequilibration processes can affect the P-V-T properties of the fluid [START_REF] Raimbourg | Fluid properties and dynamics along the seismogenic plate interface[END_REF][START_REF] Sterner | Synthetic fluid inclusions -VII. Re-equilibration of fluid inclusions in quartz during laboratory-simulated metamorphic burial and uplift[END_REF][START_REF] Vityk | Do fluid inclusions in high-grade metamorphic terranes preserve peak metamorphic density during retrograde decompression?[END_REF][START_REF] Vityk | Fluid inclusions as tectonothermobarometers: Relation between pressure-temperature history and reequilibration morphology during crustal thickening[END_REF]. Methane-rich fluid inclusions are much less prone to reequilibration than water-rich fluid inclusions [START_REF] Vityk | Fluid inclusions as tectonothermobarometers: Relation between pressure-temperature history and reequilibration morphology during crustal thickening[END_REF]. When restricting the study to methane-rich fluid inclusions, large (~150 MPa) variations in fluid pressure are still recorded in the inclusions (Raimbourg et al., 2017b;[START_REF] Vrojlik | Tectonically-driven fluid flow in the Kodiak accretionary complex, Alaska[END_REF]. The geochemical scenario developed above, where quartz precipitation in veins results from fluid pressure variations, is therefore supported by fluid inclusion data.

E-2) Depth changes in the amplitude of fluid pressure variations

The samples considered here span a temperature range from 200 to 350°C. Such range corresponds to the transition, in terms of mechanical behavior, from the seismogenic portion to the deeper, aseismic portion of the plate interface, i.e. the downdip limit of the seismogenic zone [START_REF] Hyndman | The seismogenic zone of subduction thrust faults[END_REF][START_REF] Oleskevich | The updip and downdip limits to great subduction earthquakes: Thermal and structural models of Cascadia, south Alaska, SW Japan, and Chile[END_REF]. The base of the seismogenic zone is also the locus of ETS generation [START_REF] Audet | Teleseismic constraints on the geological environment of deep episodic slow earthquakes in subduction zone forearcs: A review[END_REF]. As discussed above, all the veins considered were formed in zones of localized strain, such as the plate interface or out-of-sequence thrust, and quartz precipitation is probably the result of the deformation of the surrounding rock matrix. In particular, the higher temperature samples studied here contain crack-seal microstructures. They are similar and were collected in the same area as [START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF] and [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF], in the Morotsuka Group of the Shimanto Belt and in the Kodiak Central Belt, respectively. Crack-seal microstructures were interpreted in these studies as a possible record of ETS based on kinematics, temperature of formation, and recurrence interval.

Accepted Article

This article is protected by copyright. All rights reserved.

The comparative study of lower (~200-250°C) and higher (~300-350°C) temperature quartz veins shows, from the contrast in luminescence, large differences in the incorporation of Al and Li in hydrothermal quartz as a function of temperature. Veins formed at ~250°C show two types of luminescence, including a bright, blue luminescence that corresponds to high Al and Li contents [START_REF] Raimbourg | Fluid properties and dynamics along the seismogenic plate interface[END_REF]. Such blue luminescence is the most visible when mode I veins contain euhedral crystals with growth rims, such as described in the present work, but more generally most of quartz mode I veins show such a luminescence [START_REF] Raimbourg | Fluid circulation in the depths of accretionary prisms: an example of the Shimanto Belt, Kyushu, Japan[END_REF]. In contrast, veins formed at ~350°C are composed principally of a homogeneous, CL brown quartz (Figure 6, Figure 8 [START_REF] Raimbourg | Fluid properties and dynamics along the seismogenic plate interface[END_REF]). In a similar study, periodic patterns of growth are observed in quartz from samples precipitated at temperature conditions of ~300°C, while in the sample precipitated at 400°C, only a single type of crystal is formed, with a low luminescence and a low concentration in impurity [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF].

On the basis of (i) the inferred relationships between quartz growth and fluid pressure variations and (ii) the tectonic interpretation of the veins in the framework of subduction zone, we propose in Figure 18 a model of fluid pressure variations across depth. In the lower temperature régime, the large variations in Al incorporation in veins are the expression of large fluid pressure variations, possibly as a result of earthquakes (Figure 18). In contrast, the quasi-absence of CL-blue quartz in higher temperature veins is somehow puzzling, as the crack-seal microstructures they contain [START_REF] Fisher | The role of silica redistribution in the evolution of slip instabilities along subduction interfaces: Constraints from the Kodiak accretionary complex, Alaska[END_REF][START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF], as well as their fluid inclusion record [START_REF] Nishiyama | Fluid overpressure in subduction plate boundary caused by mantle-derived fluids[END_REF][START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF], suggest cyclical changes in fluid pressure. Compared to lower temperature samples, the absence of CL-blue quartz could be related to several factors, such as the higher temperature, a depletion in Al in the fluid, or smaller fluid pressure variations.

The difference in temperature between lower and higher T samples is relatively small, ~75°C, so that its effect on the diffusion coefficient of solutes in the fluid is limited. For example, the diffusion coefficient of silica in water is only increased by a factor of 2 or less [START_REF] Nakashima | Diffusivity of ions in pore water as a quantitative basis for rock deformation rate estimates[END_REF], i.e. much lower than the increase of the trace element content by up to 3 orders of magnitude between CL-brown and CL-blue quartz. Another explanation could be a large change in the Al-content of the fluid between lower and higher temperature samples. There is nevertheless no large mineralogical change with increasing temperature that would imply a significant depletion in Al in the fluid in higher temperature samples [START_REF] Fagereng | Fluid-related deformation processes at the up-and downdip limits of the subduction thrust seismogenic zone: What do the rocks tell us?[END_REF]. For example, the reaction described in [START_REF] Kameda | Low-grade metamorphism around the down-dip limit of seismogenic subduction zones: Example from an ancient accretionary complex in the Shimanto Belt, Japan[END_REF], in similar rocks and for similar temperatures, does not induce any change in Al content of the fluid nor any change in the phyllosilicates forming the parageneses (chlorite and mica). It follows that the absence of CL-blue quartz in the higher temperature samples might result from a smaller amplitude of fluid pressure variations than in lower temperature samples (Figure 18). Even within the crack-seal microstructures, the CL-blue patches (the triangles in Figure 8 A andB) are very restricted, so that the fluid pressure cycles accompanying cracking then sealing cycles might have limited amplitude compared to lower temperature samples. Nonetheless, crack-seal microstructures are the only domains where the CL-blue luminescence is locally present, which could correspond to events of high strain during ETS, whereas the rest of the vein would grow in steady-state conditions.

The conclusions about the limited amplitude of fluid pressure variations in the higher temperature samples are to some extent contradictory to the large range of vapor/liquid ratio or homogenisation temperatures in fluid inclusions from the Morotsuka unit [START_REF] Nishiyama | Fluid overpressure in subduction plate boundary caused by mantle-derived fluids[END_REF][START_REF] Ujiie | An explanation of episodic tremor and slow slip constrained by crack-seal veins and viscous shear in subduction mélange[END_REF]. But, as discussed above, water-rich fluid inclusions are possibly affected by various reequilibration processes, even more in deformed zones at temperature where plastic deformation is activated. Postentrapment reequilibration is therefore a possible interpretation for the large variations in fluid inclusion data, which would deviate then from the record of past pressure.

Accepted Article

This article is protected by copyright. All rights reserved.

The model proposed in Figure 18 summarizes these interpretations about the fluid pressure variations in zones of concentrated deformation of a subduction system, across the depth range encompassing the seismogenic domain and its downdip limit. This model, which compares different depth domains, is only qualitative. What is shown are fluid pressure variations in time and not fluid overpressure amplitude. The amplitude of fluid overpressure may be smaller at seismogenic depths than in the deeper domain of ETS generation [START_REF] Saffer | Hydrogeology and mechanics of subduction zone forearcs: Fluid flow and pore pressure[END_REF].

The smaller fluid pressure fluctuations near the downdip limit of the seismogenic zone might reflect the onset of quartz plastic deformation around ~300°C [START_REF] Behr | A naturally constrained stress profile through the middle crust in an extensional terrane[END_REF][START_REF] Stipp | The eastern Tonale fault zone: a "natural" laboratory for crystal plastic deformation of quartz over a temperature range from 250°C to 700[END_REF], which prevents large differential stress (as might result from a large drop in pressure within rock cavities) to be sustained. In other words, a plastically weak rock matrix would deform easily enough not to let the fluid within its cavity drop to a pressure that deviates significantly from the ambient solid pressure. Although the precise range of the brittle to viscous evolution of quartz is not well constrained, our own analysis of Shimanto Belt rheological behavior across the observed temperature gradient shows that samples deformed at T~250°C do not deform plastically, while the quartz crystals from domains at T~350°C show microstructures of recrystallization and the development of a crystallographic preferred orientation [START_REF] Palazzin | Brittle-ductile transition in subduction zones: the role of quartz[END_REF]. Alternately, smaller fluid pressure variations at ETS might also results from the lower strain amplitude of deformation events, or to much lower variations in stresses [START_REF] Ide | A scaling law for slow earthquakes[END_REF][START_REF] Schwartz | Slow slip events and seismic tremor at circum-Pacific subduction zones[END_REF], compared to normal earthquakes in the seismogenic zone. Conclusions 1-Quartz veins from a range of geological examples of low-grade terranes in subduction settings have geometries and kinematics that allow most of them to be associated with peak burial or with the deep evolution in the upper wedge. In such cases, host rock maximum temperature derived from RSCM can therefore be assigned with confidence to vein formation conditions. 2-Quartz veins from these various examples show geochemical characteristics that vary systematically as a function of temperature of formation. Lower-temperature (~200-260°C) quartz showing growth microstructures, exhibits alternations of CL-blue, high luminescence, Al-and Li-rich, fluid-inclusionsrich rims and CL-brown, low luminescence, trace element-poor, fluid-inclusions-poor rims. Higher temperature (~300-350°C) quartz is almost exclusively CL-brown, with a low luminescence, and is poor in trace elements and fluid inclusions.

3-Crack-seal textures in the higher temperature quartz are the only microstructures where the luminescent and trace element-rich quartz are present. Yet this type of quartz is not observed in all crack-seal microstructures, and when present, its spatial extent is very limited. 4-The trace element analysis of quartz shows a first-order clear correlation between luminescence patterns and Al and Li concentration variations. Al and Li concentration are correlated, with a Li/Al ratio that increases when trace element concentration increases. Li concentration is not sufficient to compensate for all the charge deficit of the substitution of Si by Al; we suggest that H + compensates for the difference. 5-The geochemical features of quartz are interpreted as reflecting the conditions of precipitation, in particular the fluid pressure conditions. The fluid inclusions-rich and trace elements-rich quartz is interpreted as reflecting out of equilibrium, rapid growth, following a large pressure drop, while the fluid inclusions-poor and trace elements-poor quartz has grown under steady-state conditions. 6-Following this interpretation, the geochemical contrast between quartz formed at lower and higher temperature reflects a large contrast in the hydrological system along depths in subduction zones, in the domains of where strain is concentrated (i.e. the subduction interface and out-of-sequence thrusts). At seismogenic depths, the fluid pressure undergoes large variations, possibly as a result of earthquakes, resulting in the precipitation of the two types of quartz. In the transition zone at the downdip limit of the seismogenic zone, the fluid pressure is relatively constant. The only exceptions are crack-seal microstructures, which may reflect fluid pressure changes associated with ETS. The amplitude of such fluid pressure changes in the ETS domain is probably much more limited compared to the seismogenic domain. 1) and (1') Structures of mode I quartz veins studied here. Mode I veins are orientated perpendicular to bedding or the metamorphic foliation plane and perpendicular to the X-kinematic axis. They form preferentially in the necked domains of sandstone lenses (a), sometimes in correspondence with shear zones (b). Veins also form in the vicinity of quartzinfiltrated shear bands, branching from the main shear plane (c), as sets of en-échelon veins (d) or as isolated veins (e). Quartz precipitates within shear bands, in which case it shows a fibrous structure parallel to the direction of motion (f). In (1), the bulk strain is simple shear, while in (1') strain is pure shear. At the microscopic scale, growth rims form in euhedral quartz in the core of Mode I veins filling boudin necks (2), while crack-seal textures are observed in veins cutting across the shale matrix (3). (top) and Kodiak island (bottom), for both lower temperature (left) and higher temperature (right) units. All the samples studied are mode I veins at a large angle to the foliation plane. In the Kodiak Central Belt, the sample is a small tension crack perpendicular to a macroscopic quartz-filled shear zone indicated in the plot. For the Uyak Complex site 2019KO24, the left plot shows data in the present geometry, while the right plot show data after rotating the foliation plane to horizontal. The number on each plot refers to the number of veins. All stereoplots in lower hemisphere, with arrows pointing to the motion of the hanging wall.

Figures captions

Accepted Article

This article is protected by copyright. All rights reserved.

to domains with zero or low density of fluid inclusions, in order to not mix the fluid and host quartz trace element signals.

Figure 16 : Compilation of Li vs. Al concentrations (in wt ppm) of hydrothermal quartz measured by LA-ICP-MS, from literature data and this study. Similar to the examples studied here, the low-luminescence quartz has a lower concentration in Al and Li than the high luminescence quartz. Additionally, the hydrothermal quartz formed at relatively higher temperature (i.e. 400°C and above) has a low concentration in Al and Li. (A) All data points and (B) data points for Al < 500ppm. References: Alpine samples (Gigerwald, Grosstal and Vättis): [START_REF] Jourdan | Evidence of growth and sector zoning in hydrothermal quartz from Alpine veins[END_REF]; Norway quartz: [START_REF] Müller | Hydrogen speciation and trace element contents of igneous, hydrothermal and metamorphic quartz from Norway[END_REF]; Hakos mountains (Greece): [START_REF] Müller | Trace elements in quartz -a combined electron microprobe, secondary ion mass spectrometry, laser-ablation ICP-MS, and cathodoluminescence study[END_REF]; accretionary prisms : this study. Figure 18 : Model of vein quartz geochemical variations as a function of temperature, in subduction zones. Vein quartz at temperatures below 300°C show alternation of Al-rich and Al-poor growth rims. Vein quartz formed at temperatures above 300°C show an homogeneous, Al-poor quartz. The crack-seal microstructures, which were observed only in the higher temperature samples, show also a low concentration in Al, except in restricted domains. We interpreted these variations in terms of fluid pressure variations. At lower temperature, corresponding to seismogenic zone, large fluid pressure variations occurs, possibly as a result of earthquakes. At the downdip limit of the seismogenic zone, variations in fluid pressure are of much smaller amplitude, even if cyclical microstructures such as crack-seal microstructures, are observed. These cyclical microstructures might be associated to cyclical deformation events, such as ETS.

Accepted Article

This article is protected by copyright. All rights reserved.

Tables Table 1 : Location of the sample analyzed in this study and associated paleatemperatures determined by RSCM, following the calibration by [START_REF] Lahfid | Evolution of the Raman spectrum of carbonaceous material in low-grade metasediments of the Glarus Alps (Switzerland)[END_REF].

Table 2 : Temperatures derived from RSCM analysis along the Big Waterfall Bay transect, NE of Afognak Island, Kodiak Accretionary Complex, Alaska.

Table 3 : Trace element concentration, in wt ppm, of Li, Al, Na and K, determined by LA-ICP-MS, of vein quartz from the Helminthoides Flysch, from the Shimanto Belt and from Kodiak Island. The analysis given here were restricted to volumes of quartz with zero or low density of fluid inclusions. The sign "-" refers to an analysis that is either below detection limit, or that was excluded as the ablation profile showed irregular spikes.

  , while comparison with other examples of veins from the literature is shown in Figures 16 and 17.

Figure 1 :

 1 Figure 1 : Geological map of the Shimanto Belt on Kyushu, Japan (left) and of the Kodiak Accretionary Complex, Alaska, USA (right), with the location of the quartz vein samples analyzed in this study.

Figure 2 :

 2 Figure 2 : Mode I quartz vein structures from Morotsuka Group in the Shimanto Belt (A and B) and from the Kodiak Central Belt in Alaska (C and D). In all examples, veins cut the metamorphic foliation S1 at a high-angle. In (A) and (C), veins are folded as a result of shortening perpendicular to S1. In (B) and (D), veins form en échelon pattern, with rotation of the internal part of each vein, consistent with top-to-the-S or SE sense of shear. In (C), veins form in the vicinity of shear-zones cutting across S1. The shear zones are also filled with quartz. The blue arrows point to veins of the same type as the ones studied here. All photographs in the kinematic X-Z plane. Outcrops : (A)-18NOB23B, (B)-18NOB13, (C)-2019KO32A and (D)-2019KO30.

Figure 3 :

 3 Figure 3: Mode I quartz vein structures from the Coherent Hyuga (A), the Hyuga Tectonic Mélange (B), in Japan, and Uyak Formation (C), in Alaska (shown by the blue arrows). (A) Quartz veins in sandstone lenses deformed by a top-to-the-SE shear zone (B) Quartz veins in boudinaged sandstone lenses, in the vicinity of a top-to-the-SE shear band. (C) Quartz veins in the necked parts of boudinaged chert layers. All photographs in the kinematic X-Z plane.

Figure

  Figure 4 : (1) and (1') Structures of mode I quartz veins studied here. Mode I veins are orientated perpendicular to bedding or the metamorphic foliation plane and perpendicular to the X-kinematic axis. They form preferentially in the necked domains of sandstone lenses (a), sometimes in correspondence with shear zones (b). Veins also form in the vicinity of quartzinfiltrated shear bands, branching from the main shear plane (c), as sets of en-échelon veins (d) or as isolated veins (e). Quartz precipitates within shear bands, in which case it shows a fibrous structure parallel to the direction of motion (f). In (1), the bulk strain is simple shear, while in (1') strain is pure shear. At the microscopic scale, growth rims form in euhedral quartz in the core of Mode I veins filling boudin necks (2), while crack-seal textures are observed in veins cutting across the shale matrix (3).

Figure 5 :

 5 Figure 5: Structural data of veins and associated deformation features in the ShimantoBelt (top) and Kodiak island (bottom), for both lower temperature (left) and higher temperature (right) units. All the samples studied are mode I veins at a large angle to the foliation plane. In the Kodiak Central Belt, the sample is a small tension crack perpendicular to a macroscopic quartz-filled shear zone indicated in the plot. For the Uyak Complex site 2019KO24, the left plot shows data in the present geometry, while the right plot show data after rotating the foliation plane to horizontal. The number on each plot refers to the number of veins. All stereoplots in lower hemisphere, with arrows pointing to the motion of the hanging wall.

Figure 6 :

 6 Figure 6 : (Left) Optical microscope and (right) CL images of quartz veins showing growth rims with alternating CL-blue, high luminescence and CL-brown, low luminescence (indicated respectively by blue and brown triangles). Zones with an extremely high density of fluid inclusions appear much darker in optical microscope images. In many instances (A, C), the CLblue rims coincide with domains of quartz extremely rich in fluid inclusions. (A) Coherent Hyuga from the Shimanto Belt (2018NOB02E), (B) Hyuga Tectonic Mélange from the Shimanto Belt (HN51b), (B) Uyak Complex from the Kodiak Island (2019K032).

Figure 7 :

 7 Figure 7 : Close-up view of the contrast in fluid inclusion density between the growth rims of the sample in Figure 6A. Optical microscope.

Figure 8 :

 8 Figure 8 : (Left) Optical microscope and (right) CL images of crack-seal microstructures in quartz veins. The quartz in these veins has predominantly a CL-brown, low luminescence, except for restricted domains within the quartz bands of crack-seal veins (e.g. the white triangles in (A) and (B)). Note in (C) the presence of a later-stage, CL-blue, crosscutting veinlet. (A) and (B) : Kodiak Central Belt (2019KO-32), (C): Morotsuka Group from the Shimanto Belt (2018NOB23B).

Figure 9 :

 9 Figure9: Typical Raman spectra of carbonaceous matter in lower and higher temperature units of the Shimanto Belt and Kodiak accretionary prism. The measured spectrum (black line) is modelled (red line) as the sum of five Lorentzian functions(gray lines, labelled D4, G, D3, G and D2), following[START_REF] Lahfid | Evolution of the Raman spectrum of carbonaceous material in low-grade metasediments of the Glarus Alps (Switzerland)[END_REF].

Figure 10 :

 10 Figure 10 : (A) Optical microscope and (B) cathodoluminescence image of growth rims in quartz grains from the Shimanto Belt sample 2018NOB-2E. (C) EPMA maps of aluminum (left) and potassium (right) of the rectangle shown in (A) and (B). The low-luminescence rims correspondto a lower aluminum concentration. The fluid inclusions are the most abundant in the CL-blue and brightest rims, which corresponds also to the highest Al concentrations. Note also the brightest spots in the Al map, which correlate with large K concentrations, and also with fluid inclusions.

Figure 11 :

 11 Figure 11 : (A) Optical microscope and (B) cathodoluminescence image of growth rims in quartz grains from the Uyak Complex from Kodiak Island (2019KO-24A).(C) EPMA maps of aluminum (left) and potassium (right) of the rectangle shown in (A) and (B). The low-luminescence rims correspond to a lower Al concentration. A thin growth rim, enriched in Al, is shown by blue arrows in (C). Note the brightest spots in Al map, which correlate with large K concentrations, and also with fluid inclusions. The CL contrast, as well as the Al concentration differences, are much smaller than on the Shimanto Belt sample (Figure10). On (C) left, the lower bound of the greyscale has been adjusted to highlight the variations within the quartz.

Figure 12 :

 12 Figure 12 : Crack-seal microstructure in quartz from the Shimanto Belt (sample 18NOB23B). (A) Optical microscope and (B) CL image. Some of the grains, which are all elongated in the N-S direction of the figure, contain successive bands of solid inclusions in their core, forming crack-seal microstructures. The blue luminescence is restricted to the quartz within the crackseal structures.(C) EPMA map of Al concentration. There is a sharp contrast in Al concentrationbetween the Al-richer quartz within the crack-seal microstructure and the Al-poorer quartz outside, within the same quartz grain. Note that the black color on the table on the right corresponds to a non-zero concentration, in order to highlight the Al concentration contrast within quartz.

Figure 13 :

 13 Figure 13 : Crack-seal microstructure in quartz from the Kodiak Formation (sample 2019KO32A). (A) Optical microscope and (B) CL image and (C) EPMA Al map from the two zones shown in (A). In (C) the red domains (highest Al content) correspond to clay solid inclusion bands, while the quartz shades are between blue and turquoise. The crack-seal structures highlighted by the dark, solid inclusion bands (A) host quartz that shows, for a portion of it, a CL-blue luminescence absent elsewhere in the vein (B). The CL-blue luminescence corresponds to a higher Al content in quartz (C). The CL-blue and Al-rich domains in quartz bands of crackseals show a patchy and complex spatial distribution, with local alternation of rims ofAl-rich and Al-poor quartz (C-bottom).

Figure 14 :

 14 Figure 14: In situ LA-ICP-MS concentration of monovalent cations Na, K and Li as a function of Al, in quartz from veins collected in the Shimanto Belt and Kodiak Island. Li is the monovalent cation with the largest concentration, and the only one whose concentration shows a correlation with Al concentration. All data are listed in Table3.

Figure 15 :

 15 Figure 15: Compilation of Li vs. Al concentration (in wt ppm) of hydrothermal quartz measured by LA-ICP-MS in the Shimanto Belt and Kodiak accretionary prism. (A) all data points, (B) data points for Al < 800 ppm. In all three geological examples, the CL-brown, low luminescence quartz, has lower Al and Li concentration than the CL-blue, high-luminescence quartz. In the two geological examples, the concentrations align along relatively well-defined correlation lines, which vary significantly from one case to another. Note that the analyses were restricted

Figure 17 :

 17 Figure 17 : Li/Al ratio in all hydrothermal quartz measured by LA-ICP-MS from literature data and this study. (A) All data points, (B) data points for Al < 200ppm. Li/Al ratio is higher for larger Al concentration (left) than for low ones (right).

Acknowledgements

The autors acknowledge support from both LabEx VOLTAIRE (ANR-10-LABX-100-01) and EquipEx PLANEX (ANR-11-EQPX-0036) projects. The dataset of sample locations, host-rock temperatures and LA-ICP-MS analysis is available in repository EarthChem (doi:10.26022/IEDA/111566). This is IPGP contribution XXXX .

Accepted Article

This article is protected by copyright. All rights reserved.

Yamaguchi, A., Ujiie, K., Nakai, S. and Kimura, G. (2012) Sources and physicochemical characteristics of fluids along a subduction-zone megathrust: A geochemical approach using syn-tectonic mineral veins in the Mugi mélange, Shimanto accretionary complex. Geochem. Geophys. Geosyst. 13, 1-22.

Accepted Article

This article is protected by copyright. All rights reserved. This article is protected by copyright. All rights reserved.

Accepted Article

This article is protected by copyright. All rights reserved.

Accepted Article

This article is protected by copyright. All rights reserved.