Fluid-Rock Interactions in a Paleo-Geothermal Reservoir (Noble Hills Granite, California, USA). Part 2: The Inﬂuence of Fracturing on Granite Alteration Processes and Fluid Circulation at Low to Moderate Regional Strain

: Fracture connectivity within fractured granitic basement geothermal reservoirs is an important factor controlling their permeability. This study aims to improve the understanding of ﬂuid–rock interaction processes at low to moderate regional strain. The Noble Hills range (Death Valley, CA, USA) was chosen as a naturally exhumed paleo geothermal reservoir. A series of petrographic, petrophysical, and geochemical investigations, combined with a fracture distribution analysis, were carried out on samples collected across fracture zones. Our results indicate that several generations of ﬂuids have percolated through the reservoir. An increase of (1) the alteration degree; (2) the porosity values; and (3) the calcite content was observed when approaching fracture zones. No correlation was identiﬁed among the alteration degree, the porosity, or the calcite content. At a local scale, samples showed that the degree of alteration does not necessarily depend on the fracture density or on the amount of the strain. It is concluded that the combined inﬂuence of strain and coeval ﬂuid–rock processes drastically inﬂuence the petrophysical properties of fracture zones, which in turn impact geothermal production potential.


Introduction
Long-term exploitation of geothermal resources is closely linked to reservoir rock petrophysical properties and regional geological settings [1]. Due to the low matrix porosity and permeability of granitic rocks, underground granitic units are considered as a reservoir only when fractures are present. These fractures provide the essential reservoir permeability and porosity for fluid flow [2,3] and are therefore of primary importance regarding geothermal exploitation [4][5][6][7][8]. These structures control the deep flow at the reservoir scale, in case of good connectivity [2,3,[9][10][11], which is one of the most important controls on the permeability [12]. Several granitic reservoirs, as Soultz-sous-Forêts in the Upper Rhine Graben, France [6,13] or the Cooper Basin in Australia [14], give invaluable experience in terms of exploration and exploitation feedback. In addition to these datasets gained from data limited to boreholes and indirect geophysical methods, studies of surface reservoir analogues are common [15][16][17]. The MEET H2020 project (Multidisciplinary and multi-context demonstration of EGS exploration and Exploitation Techniques and potentials) [18] aims to develop enhanced geothermal systems throughout Europe. Within this project, the Noble Hills (NH) range, located in the southern termination of the Death Valley (DV, CA, USA), has been chosen as an analogue of fractured granitic basements in a context of transtensional deformation [19]. Klee et al. (2021) [19] highlighted numerous evidences of hydrothermal alterations: (1) a propylitic alteration affecting pervasively a large volume of rock during the cooling of the pluton and (2) an argillic alteration, also called vein alteration [19][20][21][22][23], inducing changes of the bulk-rock chemical and mineralogical compositions and of physical properties [7,24].
A vast amount of literature [9,11,25] has proven that increasing amounts of strain within fault zones drastically change their petrophysical properties. The present paper focuses on the influence of fractures on the fluid circulation and alteration processes at low to moderate regional strain.
The present study aims to: 1.
Characterize the relations among the varying amounts of strain, fracture densities, and alteration degrees at the NH scale, as well as the sample scales through the case studies.

3.
Identify the different fluid circulation episodes through the granite body.
Macroscopic and microscopic petrographic studies, XRD mineralogical characterization of whole rock and clay minerals, bulk rock chemical analyses by inductively coupled plasma (ICP), mass spectrometry (MS), atomic emission spectrometry (AES), scanning electron microscopy coupled with energy dispersive spectrometry (SEM-EDS) for structural observation and local chemical analyses, mass balance calculations, fracture density calculations, calcimetry, and porosimetry, were performed on samples collected in the vicinity of fracture zones within areas of the NH affected by low to moderate regional transcurrent strains. The results will be discussed and compared with the protolith studied by [19].

Geological Setting
The Death Valley (DV, Figure 1a) is located in the core of a Cenozoic distributed system of dextral strike-slips, about 700 km long, comprising the Walker Lane Belt and the Eastern California Shear Zone (ECSZ/WLB) [26][27][28][29]. This narrow northwest-trending system, lying between the Basin and Range region to the east and the Sierra Nevada batholith to the west [28], today accommodate~25% of the Pacific-North America relative motion [26,30]. DV is a structural depression, NNW-SSE oriented, bounded by the Black and Funeral Mountains to the west and by the Panamint Range to the east [31]. It has been formed by a right-lateral movement giving a pull-apart structure [32].
The area of interest for this study is the Noble Hills range (NH, Figure 1b). It is located in the southernmost part of the DV region and trends parallel to the NW-striking SDVFZ at its southern end. Geological markers along the SDVFZ trace [33] suggest that the NH correspond to a transported fragment of the frontal part of the Owlshead Mountains (OM), a Cretaceous (~95 Ma, [34]) granitic pluton at a 40-41 km distance to the SE. Several attempts have been made and discussed in the literature to give insight to the geological setting of the NH and structural relationships between SDVFZ and GFZ [35][36][37]. Particular emphasis has been given to a detailed description of sedimentary sequence deposits on each side of the NH Crystalline Bedrock Slice (CBS) [38,39]. The CBS is composed of Proterozoic sediments with upward younging direction the Crystal Spring (CS) quartzite, CS dolomite, detrital flysh, and carbonate sequences possibly part of the Pahrump Group, intruded by 1.1 Ga diabase sills, the whole intruded by Mesozoic granite [40]. However, a detailed structural analysis of the CBS itself is missing. Thus, Section 4.1 of this study will give new elements, improving the observations made by [38,41] concerning the NH structure. A precise fracture pattern characterization of the NH was performed by [42], through a wide-ranging analysis scale from the microscopic scale to the regional scale. These authors showed that the NH fracture network geometry has been controlled by the SDVFZ and the GFZ.  [19,38,41]. The quaternary is not displayed, but available on the map provided by [38]. Sample locations are represented by the white dots. The asterisk shows the samples used as case studies in this work. NHF-Noble Hills Formation.
The NH granite was poorly studied in terms of microstructure, petrography, and geochemistry. A recent study conducted by Klee et al. (2021) [19] characterized the granite as a leucocratic equigranular monzogranite (S-type), ranging in the high-K calc-alkaline and peraluminous domains. The NH granite is composed of primary plagioclase, quartz, K-feldspar, and biotite. Klee et al. (2021) [19] have also shown that the NH granite underwent two alteration stages, forming secondary minerals that recorded the chemical and paleo-thermal conditions of the system: (1) a pervasive, propylitic alteration linked to the cooling of the pluton; and (2) a more local alteration corresponding to the argillic alteration (illitization and sometimes kaolinitization), which overprints the propylitic alteration. Among the primary minerals, optical observations and geochemical data show that only plagioclase and biotite are affected by those alteration processes.

Material and Sample Selection
Digital field mapping techniques have been used in the field; these include the use of portable rugged tablet laptop with internal GPS. The QGIS software has been used for sample location, geological digitization, and structural data acquisition. Georeferenced topographic maps and ortho-imagery were the initial input of the database.
A detailed petrographic, mineralogical, geochemical, and petrophysical characterization was conducted on 25 samples from the NH granite, based on fieldwork and laboratory analyses.
Samples were selected by targeting fracture zones located in CBS areas that have been affected by low to moderate strain, aimed at characterizing the influence of fracturing on fluid circulations and the associated argillic alteration, mentioned only by [19] until now. Thus, the collected samples consisted of altered granite presenting open fractures, veins, reactivated veins, and breccias. All samples were georeferenced for precise locations as well as database supply (Figure 1b). The selection took into account ranges of fracture density (F d ) defined by [19] from scanlines realized on thin sections, based on [43]: Samples selected for this study had a fracture density higher than 2694 fracs/m, meaning they fell into categories from Fd2 to Fd4. Fd4 corresponds to highly strained zones in which granite becomes a breccia. In order to correlate the fracture density with the amount of alteration, three hand specimens and seven thin sections, located in the granite body, were used for fracture extraction and analysis. These fracture characterizations are based on the scanline method, described in Section 3.1. Data from Chabani et al. (this issue, under review) [44] will be used to complete the characterization.

Fracture Network Parameters
Fracture networks can be characterized by their spatial arrangements [45]. The fracture network geometry is used to predict fluid circulation [45] and evidence of structural growth processes [46]. Spacing measurements through the 1D-scanline method are widely used to characterize arrangements [47][48][49][50]. They consist of the digitization of fractures along those lines, in order to calculate the linear fracture density P 10 characterized by the number of fractures per length calculated along the scanline [47,48]. The Terzaghi correction has been applied [51]. Two parallel scanlines were performed on the hand specimens, perpendicular to a major fracture. For each oriented thin section mosaic, two scanlines were realized perpendicular to the main structures. For each scanline, stick plots and cumulative frequency diagrams were realized to describe the fractures spatial distribution [50]. A coefficient of variation C v was computed for each scanline in order to quantify the fracture distribution [48]: C v < 1 indicates a regular fracture spacing, C v~1 indicates a random distribution, and C v > 1 indicates a clustered distribution. Fracture density data obtained by [42] on selected samples will be used in this study.

Petrographic Characterization
Twenty-four samples were selected to prepare thirty-one covered and polished thin sections, for petrographic observations. Optical microscopy was used to study the mineralogical assemblage, the alteration paragenesis, the microfabric, and the degree of microfracturing. The thin sections were observed under a Leica DM4500-P optical microscope, equipped with a Leica DFC450C camera at Institut Polytechnique UniLaSalle (ULS, Beauvais, France). Large field area imaging under polarized-analyzed and polarized-nonanalyzed light was conducted using Leica automatized stage facility and Leica Application Suite (LAS) v4.11.0 software [19].
Based on optical microscopical observations, microsites were selected on 9 thin sections for energy dispersive X-ray spectrometry (EDS) analyses [52] performed at ULS. These were conducted with a Hitachi S-3400N scanning electron microscope (SEM) equipped with a Thermo Ultradry EDS [53] and associated with NSS thermal scientific software. These analyses aimed to obtain qualitative and semi-quantitative chemical compositions and to characterize the microstructure of different selected phases. The analytical conditions consisted of a 50 µA beam current, an accelerating voltage of 20 kV, and an acquisition time of 30 s. A NORAN-type correction procedure was used.
In addition, cathodoluminescence (CL) analyses [54] was performed on 8 thin sections from 4 samples. CL imaging was performed (1) at the University of Göttingen using a "hot-cathode" cathodoluminescence microscope HC3-LM apparatus and (2) at CY Cergy Paris University in the Geosciences and Environment Cergy (GEC) laboratory, using a cathodoluminescence Olympus BX50. CL of minerals is predominantly a "defect luminescence". This allowed the detection of distribution patterns of certain trace elements, such as iron (Fe), which is the most efficient quencher element and manganese (Mn), which is the most important activator element. Many minerals show visible CL colors like quartz (blue-purple), K-feldspar (blue when fresh and brownish when weathered), plagioclase (green and blue-purple when highly affected by hydrothermal alteration), calcite (yellow or yellow-orange when resulting from hydrothermal alteration), dolomite (orange-red), magnesite (red), apatite (yellow), kaolinite (dark blue) [54,55]. Illite shows no luminescence [54,55].
Most of the mineral abbreviations used in this paper refer to [56] mineral symbols and only a few others were defined by the authors.

X-ray Diffraction (XRD)
XRD analyses were carried out to identify and determine the mineral phases (semiquantitative, around 3-5%). The analyses were performed at (1) ULS using a D8-Advance Bruker-AXS (Siemens) diffractometer with a Ni-filtered CuKα radiation at 40 kV and 40 mA, a primary Soller slit of 2.5 • , divergence slit of 0.6 mm, and a secondary Soller slit of 2.5 • , with a detector slit of 0.1 mm and an anti-scattering slit of 0.6 mm; and (2) the University of Göttingen, using a PHILIPS PW1800 diffractometer with a Cu-anode and an automatic divergence slit. Whole rock powders and oriented clay fractions (2-6 µm and <2 µm) analyses were performed on 12 samples. Quantitative phase analysis based on reference intensity ratio values were performed on randomly oriented whole-rock powders with a step length of 0.5 • and a scan speed of 0.014 • /s over the range 3 • -70 • 2θ for whole rocks composition. The uncertainty is estimated to be ±5%. Clay mineral separation was conducted using a technique described by [19], based on [57,58], according to the standard techniques suggested by [59]. For the XRD analyses of both fractions, oriented specimens were measured at air-dried (AD), solvated with ethylene glycol (G), and after heating at 550 • C (H) conditions. These three analytical steps were routinely used to better determine clay minerals content and swelling properties [7]. The clay minerals identification, which is based on d-values and the relative intensity of their 00l reflections, was undertaken, referencing [60,61]. These measurements were performed with a step length of 0.5 • and a speed of 0.01 • /s per step over the range 3 • -35 • 2θ. The interpretations of the data obtained at ULS were performed using the DIFFRAC EVA v4.2 (©Bruker AXS) software. Two fractions were collected in order to (1) separate the last produced or neoformed clay minerals (fraction <2 µm), which could be assimilated to the result of the last fluid circulation event, from the fraction 2-6 µm assimilated to old grain recrystallization or to a possible mix between detrital and neoformed clay mineral; (2) obtain the respective clay composition; and (3) obtain the temperature conditions.

Kübler Index and Kaolinite Crystallinity Index determination
The Kübler Index (KI) was used to define the limits of metamorphic zones (diagenetic zone, anchizone, and epizone) [19,62], following the recommendations for Kübler Index calibration of [59] and the CIS-KI transformation formalism of [63], as well as the temperature of illite formation [64,65]. The KI was calculated from the illite crystallinity (IC), which is defined as the full width at half maximum (FWHM) of the (001) 10Å peak of illite, on the AD oriented clay fractions. A recent study [66] has confirmed that IC provides a useful method for characterizing regional grades of diagenesis and low-grade metamorphism. Temperature is thought to be the main factor controlling IC, but other parameters, such as lithology, also have important effects [67][68][69]. Working at constant lithology allows this effect to be avoided. In some conditions, an altered granite shows similarities with diagenetic reactions present in feldspathic sandstones [20]. Thus, it is possible to characterize temperature ranges by using the KI values obtained for granitic rocks. The values obtained by IC were thus standardized using the crystallinity index-standard (CIS) samples provided by [70] in order to calculate the KI. The KI values of raw data expressed in ∆ • 2θ, were measured into three slots, corresponding to different campaigns, which induced three standardizations given by [19]. The Kaolinite Crystallinity Index (KCI) was calculated and put in parallel with the KI in order to identify any correlation between both, and to determine temperature evolution as an indicator for samples free of illite. The KCI is defined as the FWHM of the (001) 7Å peak of kaolinite on the AD oriented clay fraction.

Manocalcimetry
Manocalcimetry is used as an indicator of the total calcite content in a rock sample and to assess its possible influence on permeability and consequently better understand the hydrothermal sealing of a reservoir [71]. Calcite content was determined using an OFITE 152-95 manocalcimeter composed of a glass flask and a high precision manometer. The analyses were performed on 15 samples at CY-GEC. High quality results were obtained with an accuracy of ±0.5 wt.%. Prior to measurement, sample preparation was conducted according to [71]. Two replicates were performed for each sample, following the protocol described by [19], in order to check the reproducibility of the results. This procedure was considered to be good when the difference between the two results was lower than 0.5 wt.%, corresponding to the precision mentioned above. The CaCO 3 percentage was calculated according to [19]. The average calcite content of a fresh granite is 0.252 wt.% and does not exceed 1.8 wt.% [72]. As a consequence, measurements above this last value can be regarded as a calcite anomaly [71], due to hydrothermal alteration.

Ethanol Saturation Porosimetry
Porosity measurement quantify the available volume for fluid storage [73]. The connected porosity was measured on three samples by the triple weighing method [19,74] defined by the RILEM standard (test n • I. 1,1978). The measurements were carried out at CY Cergy Paris University in the Mechanics and Material for Civil Engineering laboratory (L2MGC). They consist in the saturation of the samples after vacuum degassing. The analyses and calculations were conducted using the technique described by [19]. Ethanol has been chosen instead of water in classical methods, in order to avoid possible clay swelling [75], which could lead to the destruction of the sample and bias to porosity values. Even though ethanol (0.469 nm) is a larger molecule than water (0.343 nm), the pore volume is not estimated to be under-evaluated, because the pore size is likely much larger than that of ethanol molecules.

Bulk Rock Geochemical Analyzes
Chemical analyzes of the major element oxides, rare earth elements, and trace elements were conducted on five samples, selected as case studies in this paper. Analyzes were performed at the Bureau Veritas Minerals (Vancouver, BC, Canada) using inductively coupled plasma emission spectrometry or mass spectrometry (ICP-ES and ICP-MS, respectively). Samples were crushed and mixed with LiBO 2 /LiB 4 O 7 flux. Crucibles were fused in a furnace at 980 • C. The obtained cooled bead was dissolved in ACS grade nitric acid and analyzed. Loss on ignition (LOI) was determined by igniting the samples split then measuring the weight loss. These chemical analyses were performed in order to determine the element transfers and the degree of alteration.

Mass-Balance Using Gresens' Method
Bulk-rock mass-balance were used to establish element transfers during hydrothermal alteration by applying Gresens' mass-balance procedure [76], consisting of the quantification of gains and losses of major elements by comparing unaltered and altered samples. NH_2 is considered as the reference because it is the freshest sample collected in the area by [19] and 3 altered samples were selected from this study (NH_20, NH_28, NH_32_3, and NH_37). The following equation defined by [76] relates the composition and volume of altered rocks to unaltered/fresh rocks: With X n corresponding to the gain or loss of a considered element n (absolute mobility (%)), F v the volume factor, d A and d F the density (in g/cm 3 ) of the altered (A) and the fresh (F) rocks, respectively, and C n A and C n B , the oxide percentage of the element n for the altered (A) and the fresh (F) rocks, respectively, which were given by the bulk geochemical analyses, recalculated without the LOI. The densities were calculated from the sample masses obtained during the porosity measurements. The F v was calculated for each altered sample by considering X Al2O3 = 0 following the analysis of the values. Xn was calculated for each element. Then, the obtained Fv values being very close to 1, the Xn values were recalculated by considering F v = 1 (constant volume). Hence, X n values > 0 (positive values) represent the gains and X n values < 0 (negative values) represent the losses for the considered elements. This calculation was applied for the major elements in each sample.

Noble Hills Structural Overview
We investigated the geometry of outcropping Proterozoic sequences between Pipeline and Cave Spring washes ( Figure 2a) and confirm at first the plutonic nature of the contact at the base of the Proterozoic sequence with the underlying Mesozoic granite by metamorphic halos [40]. The whole architecture of the NH presents signs of deformation affecting both Proterozoic and granitic units. Proterozoic units present much evidence of deformation with (1) several tectonically induced duplications of the stratigraphic sequence ( Figure 2b) and (2) a lateral bending, stretching, and boudinage of this Proterozoic cover nappe stack [41]. Boudinage was identified in the field and in map view where progressive dismantling of stratigraphic markers is observed. (Figure 2b,c). Development of synthetic, oblique dextral shears (Figure 2c), offsetting along the CBS strike and accommodating the deformation, have also been identified. The age and tectonic significance of the nappe stack described in (1) is still unclear at present and possibly predates Mesozoic magmatic intrusion, as it does not appear intensively deformed along the basal contacts of the nappes. Lateral bending and stretching mentioned in (2) and depicted in Figure 2b,c are related to SDVFZ activity, since all kinematic indicators are consistent with progressive southeast oriented dextral shearing. This feature is ubiquitous within the NH as described by [33]. Given these new findings dealing with internal CBS structural organization, it is expected that areas along the rear southern limit of the CBS have recorded less transcurrent shear compared to areas situated along its northern front. Thus, the CBS gives an opportunity to study the effects of progressive transcurrent shearing within a granitic geothermal analogue. The above considerations have important implications regarding strain distribution within the NH and, in turn, concerning the sampling strategy. Given the above-mentioned findings, we consider the CBS as an exceptional example of large-scale cataclastic flow in which rock masses of the northern domain are dragged as a continuous body into a transcurrent deformation corridor (Figure 2a). No signs of a discrete, large-scale deformation structure appear on the map presented in Figure 2, where yellow dashed lines show instead a continuous flow accommodated by brittle deformation processes.

Petrographic Description
As described by [19], the primary assemblage of the NH granite is made of plagioclase (mainly oligoclase in composition), quartz, K-feldspar (perthitic orthoclase and sometimes microcline), and biotite. The granitic samples collected for this study show strong evidence of intense alteration, especially on plagioclase and biotite. Plagioclase transformed into illite, kaolinite and/or calcite and biotite into illite and oxides. In rare samples, K-feldspar perthite could be partially altered into illite. However, the magmatic texture is almost preserved, except in the case of breccias, which are not necessarily highly altered. For all samples, the alteration product is always the same, but the relative proportions of the different secondary phases might vary. The mineralogical composition and the degree of fracturing of each sample is given in Table 1. The clay mineral composition of some samples has been precisely determined for the <2 µm and 2-6 µm (Table 1). It reveals the presence of illite/smectite (I/S) mixed-layers in three samples (NH_17, NH_24, and NH31_1_2) in both fractions when it was measurable. The I/S identified is always illite-rich (R3), with more than 90% of illite [60], where R represents the Reichweite parameter [60]. A small amount of corrensite (trioctahedral variety of regular 50:50 chlorite/smectite mixed-layer (C/S) [77]) was also identified in the samples NH_16, NH_27, and NH_32. Both mixed-layers (I/S and C/S) were already identified by [19] in the sampled granite away from the fracture zones. Two new patterns were identified in this study ( Figure 3):

•
The first pattern ( Figure 3a) shows illite characterized by peaks at 10.1 Å and 4.99 Å, which support the previous results [19]. Illite is present in almost all the samples except NH_28, NH_36, NH_38, and NH_39. Kaolinite is also present and could be associated to illite or not.

•
The second pattern (Figure 3b) shows the presence of dioctahedral vermiculite. It is characterized by a peak at 14.32Å in air-dried condition, which slightly swells up to 14.59Å after glycol solvation and slightly collapses to 12.26Å after heating. It was identified only in samples NH_24, NH_26, and NH_33, where it is associated to illite and kaolinite. Table 1. List of samples collected in the field with their mineralogical composition determined after microscopic observations (optical-microscope and SEM-EDS) and XRD analysis (primary minerals and secondary minerals), their degree of microfracturing after [43], and their clay mineral composition for the fractions <2 µm and 2-6 µm. Abbreviations after [56], except for the following: Olg-oligoclase, Ox-oxide, I/S-illite/smectite, C/S-corrensite, n.m.-not measured, n.a.-not analyzed.  Figure 3. Two examples of XRD patterns obtained for the clay fraction <2 µm in air-dried (AD-dark), after glycol solvation (G-blue) and heated (H-red) conditions for the NH granite, which were not identified in the fresh granite described by [19], and showing different clay compositions. The high proportions of illite and kaolinite are characteristic of the argillic alteration [19]. Argillic alteration is prevalent in the vicinity of fracture zones (Table 1). Vermiculite can be interpreted as either hydrothermal alteration or weathering. This will be discussed in Section 5.1.
Microscopic observations and SEM-EDS analyses have shown that, in almost all the samples, carbonates (dolomite, siderite, ankerite, and calcite) and oxides were present as veins and or in replacement of plagioclase. Moreover, XRD analyses revealed the presence of gypsum and/or halite in minor amounts in few samples (Table 1). Halite only occurred in sample NH_36 and gypsum was found in samples NH_17, NH31_3, NH_34, and NH_36.

Fracturing and Fluid Circulation
The NH granite presents ubiquitous signs of both fracturing and alteration, due to its location along the SDVFZ major fault corridor. The area lies along a deformation corridor governed by transcurrent deformation along the SDVFZ. In addition to the present contribution, Klee et al. [41] argue in favor of a deformation gradient affecting the CBS. At low to moderate strains, the system presents fracture zones composed of a fault core in which fractures are branching and where most of the displacement is concentrated (Figure 4a). These fracture zones affect the surrounding rock, creating an important fracture network (the damage zone, DMZ). Open fractures and veins are observed in this DMZ ( Figure 4b).
Fault zones may act as a channel when connected and open fractures are present or as a barrier when fracture are sealed by mineralization. The distinction is made between (1) zones of opening related structures (i.e., where no signs of displacement are recognized between the two borders of the vein, Figure 4b); and (2) zones of lateral displacement (i.e., indicative of a shear displacement is observable) through three samples used as case studies (NH_31 (Figure 5a), NH_36, and NH_23 (Figure 5b)). Note that these samples are spaced at 2 m and are almost perpendicular to each other.

Opening Related Structures with Minimum Shear Displacement
A significant number of veins, veinlets, and microfractures were identified in the different samples and outcrops investigated. Veinlets of various mineralization natures (≤100 µm wide) were also identified in the samples. This variety of fracture infills shows different fluid circulation episodes. Rare veinlets of quartz were observed and intersected by illite veinlets (Figure 6a), which can be contemporaneous to kaolinite veinlets. Illite development is dependent on host mineral properties in terms of mechanical resistance and chemical stability. In that sense, quartz and K-feldspar, which remains unaltered, presents sharp open fractures in which illite precipitates (Figure 6b). Plagioclase, which is altered, presents blurred vein borders and start to spread out pervasively, due to the secondary porosity created by the alteration process (Figure 6b). These veinlets are generally intersected by carbonate veins (Figure 6c

Infill of Fractures Developed with a Confirmed Lateral Shear Displacement
The previous section shows that all of the samples in the area present many fractures where fluid has circulated, precipitating or not various secondary minerals. However, the area is constrained by a shear component. Three samples were selected in order to describe the relationship between the amount of the strain and fluid circulation.

Sample NH_31
The outcrop of the sample NH_31 ( Figure 5a) presents a vein around 10 cm thick crosscutting the granite and with a pinkish halo developed on each border. This halo represents the hydrothermalized zone of the granite. The sample NH_31 (Figure 7a) is composed of (1) a brecciated vein presenting brecciated quartz remnants, which were overprinted by a carbonate matrix and granitic clasts transported by the carbonate matrix; (2) a hydrothermalized zone, which corresponds to the pinkish halo; and (3) the porphyritic granite crosscut by veinlets filled by carbonates coming from the brecciated vein. Quartz veins are scarce in the NH granite. Fractures are dominantly rich in carbonates. It can be observed through this sample that carbonates crosscut the quartz.
• Spatial fracture distribution: two scanlines were realized in the granitic part of the sample (Figure 7a), from the brecciated vein towards the host rock, in order to evaluate the evolution of the fracture density along those two profiles. The spatial analysis is summarized in Table 2.
The cumulative frequency plotted against the fracture-projected position is presented in Figure 7b. The fracture frequency regularly and slowly increases over the first half part of the diagram curve, then increases regularly but more significantly. The Cv of 0.67 and 0.60 for each scanline shows a regular to random fracture arrangement, which is confirmed by the stick plots in Figure 7c. Both scanlines present an equivalent fracture density of 829 and 882 fracs/m. Moreover, three thin sections were gained from the sample and one scanline was realized for each thin section (Figure 8a,b). As for the sample scanlines, the thin section scanlines were realized perpendicular to the vein and towards the host rock ( Figure 8c). SL_NH_31_1 shows a fracture frequency that increases following a random arrangement (Cv = 1.03) (Figure 8c,d).
The SL_NH_31_2, parallel to the previous scanline, also presents an irregular distribution. The frequency increases slowly at the beginning, followed by a high fracture density zone around 0.7-1.05 (fracture cluster), then a more significant increase comprising a new fracture cluster around 2-2.2 ( Figure 8c). The slightly higher Cv compared to the previous scanline (Cv = 1.07) and the fracture distribution ( Figure 8c) indicate a random to clustered arrangement of the fractures.
The last scanline SL_NH_31_3, which is the furthest from the vein, shows a greater overall increase of the fracture frequency compared to the two previous scanlines (Figure 8c). Only one fracture cluster was identified around 1.7-1.9 cm even if the Cv is higher (Cv = 1.11). The fracture arrangement is again considered as random to clustered.
Cv values, as well as fracture density values, increase from the vein towards the host rock, from 1.03 to 1.11 and from 2997 to 5084 fracs/m, respectively (Figure 8d).
By comparing fracture analyses between sample scale and thin section scale, a clear difference is observed among P 10 , the mean spacing, and Cv values ( Table 2). At thin section scale, the fracture density is significantly higher (×4.7) than at sample scale and the mean spacing is six times lower. At sample scale, the Cv is lower than 1 indicating an almost regular spacing between the fractures, even if there is a change of the slope (Figure 7c). Whereas, at the thin section scale, the Cv is around 1 or slightly higher, indicating a global random distribution of the fractures along the scanlines with the appearance of a few clusters (Figure 7c). For both, differences in the fracture distribution can be linked to a subjective bias during the data collection [78], but also to the image resolution, which prevents from seeing all the fractures and so induced a bias. Thin section fracture analyses are, thus, more precise.  • Petrographic and petrophysical characterization: the brecciated vein presents a large variety of mineralogical phases. As shown in Figure 7a, this zone shows a carbonate matrix containing clasts of granite, which overprints a brecciated quartz vein remnant showing an undulatory extinction but no evidence of dynamic recrystallization. The carbonate matrix is complex, composed of several phases (Figure 9a-g). Phase 2 is composed of a matrix rich in carbonates (calcite, dolomite, ankerite) transported clasts of quartz, K-feldspars and few biotite, but also presents barite precipitation patches (Figure 9e). This phase crosscuts the quartz and the calcite veins.

3.
Phase 3 consists again of calcite veins. However, they crosscut the whole rock, and appear as a dark single phase (Cal C) under CL (Figure 9b).
Evidence of hydraulic fracturing, as presented in Figure 6e, have also been observed in this part of the sample.
By focusing in the carbonate matrix, at the contact with the granite, a new phase is observed, composed of dolomite embedding small clasts (mainly quartz). This phase is intersected by the phase 2, which is intersected by the phase 3 (Figure 9e,f). Examining the hydrothermalized granite along the brecciated vein (NH_31_1 and NH_31_2), the granite shows a preserved magmatic texture with plagioclase fully altered into illite, ± kaolinite, and calcite, completely altered biotite into illite and oxides, and unaltered perthitic K-feldspar (Figure 9h,i). Hematite and numerous dolomite/ankerite and calcite veinlets are present. Moving away from the brecciated vein (NH_31_3), the granite appears less affected by the alteration. Plagioclase and biotite are only partially altered and K-feldspar and quartz are unaltered (Figure 9j,k). Unlike the hydrothermalized zone, the granite presents only a few veinlets of calcite.

Sample NH_36
This sample (Figures 1b and 10a) shows a fractured granite with a preserved magmatic texture intersected by brecciated zones. As for sample NH_31, the fracture distribution analyses, as well as a petrographic description, were realized.
• Spatial fracture distribution: the spatial fracture distribution analysis was realized based on data from [42]. One scanline was realized through the sample in order to evaluate the evolution of the fracture distribution and density (Table 2, Figure 10a). The cumulative frequency against the fractures projected position presented in Figure 10b shows a fracture frequency slowly increasing, as well as a fracture cluster. The Cv of 1.03 indicates a random arrangement of the fractures along the scanline (Figure 10c). A fracture density of 342 fracs/m was compiled (Figure 10c).
Two thin sections were made from the sample (Figure 10a,d). Two perpendicular scanlines were realized on each thin section (Figure 10d).  Regarding the thin section NH_36_1, the scanline SLH_NH_36_1, perpendicular to the main fractures, shows a fracture frequency increases slowly the first 30%, then increases more strongly by presenting five fracture clusters (Figure 10e). This fracture distribution is highlighted by a Cv = 1.49, confirming a clustered arrangement of the fractures (Figure 10f). The fracture clusters are respectively comprised between, 0.9-1.2 cm, 1.4-1.5 cm, 1.55-1.6 cm, 1.7-1.85 cm, and 2.05-2.1 cm intervals (Figure 10e,f). The scanline SLV_NH_36_1, parallel to the main fractures, shows a different trend. The fracture frequency increases slowly and irregularly (Figure 10e). This distribution is highlighted by a Cv = 1.01, indicating a random arrangement of the fractures (Figure 10f). Considering both scanlines, the average fracture density is of 5474 fracs/m.
Regarding thin section NH_36_2, the scanline SLH_NH_36_2, parallel to the main fractures shows a fracture frequency trending similarly to SLV_NH_36_1. Four fracture clusters were identified (Figure 10e), also visible in the stick plot between 0.95-1.3 cm, 1.85-2.05 cm, 2.5-2.7 cm, and 3.1-3.2 cm (Figure 10f). The Cv = 1.11 confirmed a random to clustered arrangement of the fractures. The scanline SLV_NH_36_2, perpendicular to the main fractures, crosscut the magmatic preserved texture of the granite and a breccia. A change of the curve inclination is visible around 65% of the scanline, where the fracture frequency strongly increases (Figure 10e). It corresponds to the limit between both zones (Figure 10d). The Cv = 0.92 indicates a global regular to random fracture arrangement along this scanline (Figure 10f). Considering both scanlines, the average fracture density is of 4581 fracs/m.
Because fractures are difficult to recognize in the breccia of the sample NH_36_2, the fracture density calculated for both thin sections may be underestimated.
As for sample NH_31, an important difference between the P 10 values of the sample and the thin sections is observed. However, a global random arrangement of the fractures at both scales is observed. Only SLH_NH_36_1 presents a clustered distribution.

• Petrographic and Petrophysical Characterization
The two thin sections (Figure 11a-c) allow petrographic characterization of the sample. The first thin section, NH_36_1 (Figure 11a,b), shows a breccia of the granite with significant variations of the grain size. No crystallographic preferential orientation (CPO) is observed. Some zones show brecciated minerals with small precipitation of carbonates as calcite and dolomite in the microfractures (Figure 11d). Primary minerals are only slightly altered but intensively deformed. Quartz shows a pronounced undulatory extinction and starts to dynamically recrystallize by means of sub-grain rotation processes (Figure 11d). Other zones show preserved minerals, but also a significant amount of carbonate, mainly calcite, between the clasts (Figure 11e). The second thin section, NH_36_2 (Figure 11c), shows a granite with a preserved magmatic texture, but fractured, and a breccia composed of a carbonate matrix and transported clasts coming from the granite with no CPO. The preserved granite part shows a low alteration degree. Quartz, K-feldspars (orthoclase and microcline), and biotite are fractured and unaltered. Plagioclase is fractured, slightly altered, and replaced by kaolinite (Figure 11f). The breccia shows a matrix composed of dolomite and calcite flowing through transported clasts of unaltered quartz and feldspar and lenses of kaolinite (Figure 11g).
Bulk analyses show a low LOI (2.6%), a calcite content of 5.95%, and a porosity of 11.5%. XRD bulk and clay fraction analyses also reveal the presence of gypsum and halite.  • Spatial fracture distribution: two scanlines were realized in the granitic part of the sample NH_23 to evaluate the fracture arrangement from the vein towards the host rock (Figure 12a). The spatial analysis is summarized in Table 2. Fracture distributions presented in Figure 12b,c for each scanline show an evolution in two steps of the fracture cumulative frequency (Figure 12b).
Along the SL1_NH_23 scanline, the frequency slowly increases and shows a regular fracture distribution. At 40% of the scanline, the frequency slightly increases, with a regular arrangement of the fractures. The stick plot Figure 12c highlights this distribution. The SL2_NH_23 scanline presents a similar evolution as SL1_NH_23 scanline, with an increase in the frequency at 3 /4 of the scanline. This slope change is expressed at the end by a fracture cluster around 8.7-10 cm, clearly visible on the stick plot (Figure 12c). The Cv of 0.68 and 0.67 for each scanline shows a regular to random arrangement of the fractures along both scanlines (Figure 12c, Table 2). Both scanlines fracture densities are equivalent, with respectively 1382 and 1322 fracs/m (Figure 12c, Table 2). In order to complete the dataset, the method was repeated on two thin sections from the granite. One is at the border of the vein (NH_23_1) and the second away from it (NH_23_2) (Figure 13a). A scanline was realized in both thin sections (Figure 13b), perpendicular to the major vein. The spatial analysis is presented in Table 2.
Along the SL_NH_23_1 scanline (Figure 13c), a low fracture density zone is present at the beginning of the scanline, then there is an increase of frequency between 0.9-1.4 cm (fracture cluster, Figure 13d), followed by an irregular fracture distribution until the 0.4 last cm where a fracture cluster was identified. Three fracture clusters are present in the SL_NH_23_2 scanline around 1.5 cm, between 1.7-2.1 cm, and around 2.3 cm. From the first fracture cluster, the frequency increases significantly and irregularly. The scanline SL-NH_23_1, with the Cv = 0.77, shows a regular to random arrangement of the fractures, while the Cv = 1.19 of the scanline SL_NH_23_2 indicate a random to clustered arrangement (Figure 13c). The fracture distribution varies with the position of the scanline and shows a P 10 slightly higher close to the vein (4798 fracs/m) than away from it (3879 fracs/m) (Figure 13d).

•
Petrographic and petrophysical characterization: the vein of the sample NH_23 consists of two parts. The major part is mainly composed of quartz, few feldspars, and oxide minerals having no CPO. Ankerite veinlets crosscut it, and are also intersected by calcite veinlets (Figure 14a,b). This vein includes clasts of carbonates having the same characteristics as the carbonate vein at the contact with the granite (Figure 12b). This thin carbonate vein is composed of ankerite layers and granite lenses (Figure 14c). By zooming in on the ankerite layer, some part of it appears as lenses with automorphic crystals in the swell, which are crushed and stretched through the pinches (Figure 14d). Granite is pinched between these ankerite layers. Quartz clasts present in the granite lenses show an undulatory extinction, as well as some evidences of subgrain rotation phenomena. Carbonates can deform plastically and accumulate large amount of strain at relatively low P-T conditions, from a minimum temperature of 180 • C [78][79][80][81]. Quartz crystal plastic deformation is dominant from 600 • C [82]. Hence, in the present situation, carbonates accommodate large amounts of deformation by means of crystal plasticity, while granite is being deformed by cataclastic flow, i.e., a process accommodated by brittle processes. A localized and significant amount of deformation occurred in this vein. The major vein of this sample is thus composed of (1) a highly deformed zone made of ankerite and granite lenses and (2) a siliceous vein, which contains clasts of similar composition to the carbonate vein. We suggest that the carbonate vein was the first event, which was intensively deformed, and then intersected by the siliceous vein, tearing off pieces of the carbonate vein.
Calcite and ankerite veins identified in Figure 14a,b intersect the carbonate vein and the siliceous vein. As presented above, the magmatic texture of the granitic part of the sample is preserved. However, a color change of the granite is observed close to the vein. At a sample scale, the granitic part along the vein (NH_23_1) shows a pronounced pinkish color (Figure 13a). In the thin section, the granite appears almost completely altered. Plagioclase and biotite are no longer recognizable and are replaced by illite and oxides. Perthitic K-feldspar are not completely affected by the alteration, but are partially replaced by illite (Figure 14e). Away from the vein (NH_23_2), the granite appears whitish at a sample scale (Figure 13a). Compared to NH_23_1, plagioclase are completely replaced by illite and kaolinite, but biotite are only partially altered and K-feldspar are unaltered (Figure 14f). This indicates that the granite is more affected by the alteration along the vein, in the hydrothermally altered zone. Bulk analyses performed on the granitic part of the sample show a LOI of 5.2%. In all the granite, calcite veinlets are observed coming from the major vein.

Geochemical Analyses
The geochemistry of major elements, presented in Table 3, allows establishing diagrams representing selected elements versus LOI ( Figure 15) in order to study the alteration degree of the selected samples.

Calcite Content
Calcite is present in the majority of the samples. Samples homogeneously distributed in the studied area were selected for calcite content analyses. The calcite content is represented on the map in Figure 16 by the white dots and figures. The values range between 1.2% and 8.7%. Fracture zones, as well as the surrounding granite, present variable calcite content. However, the calcite content appears higher in the southeastern part of the range than in the northwest. All of the samples, except NH_17, NH_26, and NH_31_3, have a calcite content higher than 1.8%, meaning that they were affected by hydrothermal alteration, according to White et al. (2005) [72].

Porosity and Mass-Balance
Porosimetry measurement was only performed on samples NH_27, NH_31_3, and NH_36 (Table 4). NH_36 is the closest to a fracture zone and has the highest porosity (11.5%). This suggests that the porosity increases with proximity to fracture zones. Porosity measurements were used to better decipher the petrophysical properties of the rock for geothermal purpose and to realize mass-balance calculations. The chemical compositions recalculated after suppression of the LOI, the densities, and volume factors (F v ), assuming X Al2O3 = 0 of the rock samples used for mass balance calculations, are given in Table 4. The F v shows values very close to 1, indicating a small change in volume between the freshest sample and the altered ones.
The results of mass-balance calculations, corresponding to the absolute mobility values of each element for each sample, are given in Table 5, by considering a constant volume, F v = 1, as exposed by [76]. Positive values represent the gains and negative values the losses of the considered elements. Absolute mobility values indicate almost immobility of Al 2 O 3 , TiO 2 , and MnO. They also show losses of SiO 2 and Na 2 O, which are three times more important in NH_27 than in NH_31-3, and twice that of NH_36. The losses of Fe 2 O 3 and MgO are more important in NH_31_3 and less important in NH_36. CaO losses are observed in NH_31_3, whereas CaO gains are observed in NH_36 and a high gain is observed in NH_27. Finally, the results show a gain of K 2 O in samples NH_27 and NH_31_3 and a loss in NH_37.  Table 4. Porosity values, recalculated bulk-rock major element weight percentage, altered rock density/fresh rock density ratio, volume factor (F v ) assuming the immobility of Al 2 O 3 (X Al2O3 = 0).

Temperature Conditions
The clay composition of the fractions <2 µm and 2-6 µm was given in Table 1. For some of the samples containing enough illite, Kübler Index (KI) was measured for both fractions ( Table 6). The fraction <2 µm likely is supposed to represent neoformed clay minerals, thus corresponding to the youngest event. By contrast, the fraction 2-6 µm might contain either more developed illite crystals by inherited grains recrystallization, or a mix of detrital and neoformed minerals. The KI gives information concerning the degree of metamorphism for each fraction. In some fraction, it was not possible to measure the KI due to the too low amount of illite. KI values range from 1.32∆ • 2θ to 0.29∆ • 2θ for the <2 µm fraction and from 1.27∆ • 2θ to 0.21∆ • 2θ for the 2-6 µm fraction ( Table 6). The spatial distribution of the samples and their associated KI values are shown in Figure 17 by the dots. In both fractions, most of the samples present KI values corresponding to the diagenetic zone, except for two samples (NH_32 and NH_34). These samples are located in the northwestern part of the main granitic body (Figure 7a,b), at the contact with the Proterozoic basement. They present lower KI values of low anchizone and epizone. No significant variations in KI or trends are observed between both fractions. Some samples contain only kaolinite and no illite. Positive trends have been identified between KI and KCI for both clay fractions <2 µm and 2-6 µm ( Figure 18). Samples NH_6, NH_8, and NH_31_1_2, with out of range values, due to the presence of I/S and/or corrensite interfering with the 7Å peak, were excluded from the correlation. The trends observed in the two different fractions are similar and give more confidence into the concomitant evolution of KI and KCI. The correlation provides temperature indications for illite-free samples. Thus, samples NH_26, NH_33, and NH_36 show values equivalent to diagenetic zone and anchizone conditions.

Argillic Alteration Dominance
The characterization of the NH granite alteration processes in the vicinity of fracture zones reveals the presence of various secondary minerals, such as clay minerals (kaolinite, illite, I/S mixed layer, corrensite, and vermiculite), carbonates (calcite, dolomite, ankerite, siderite), and oxides. Thus, two successive types of hydrothermal alteration events, characterized by [19] in the protolith, were observed:

1.
The propylitic alteration, which is an early stage of alteration affecting pervasively the granite during the cooling of the pluton [83]. It involves mainly the formation of corrensite and epidote considered as the major indicators of the propylitic alteration [19,22,84]. Only plagioclase and biotite are partially affected.
However, compared to [19], the samples are more altered and corrensite is not as evident as in the samples from this study. Kaolinite and illite are the prevailing clay minerals. This confirms a significant fluid-rock interaction near fractures and a predominant argillic alteration, which has overprinted almost completely the propylitic alteration. Kaolinite, being the dominant clay mineral, indicates either a more important leaching of the rock or the circulation of a more acidic fluid [90].
The samples NH_24, NH_26, and NH_33 contain a small amount of vermiculite, which can be associated to hydrothermal alteration at low temperature or to weathering. It is commonly thought that most of the vermiculite is formed under supergene conditions [91,92]. In these samples, the amount of vermiculite is low and the weathering contribution already discarded by [19].

Thermal Evolution toward Fracture Zones
Studied samples present KI values characteristic of the diagenetic zone, meaning temperatures lower than 200 • C [93]. However, KI values of samples NH_32 and NH_34 have lower KI values with anchizonal to epizonal conditions, meaning temperatures around 300 • C [93]. These two samples present the highest fracture density (Fd4) suggesting that shear heating could contribute to a local increase of temperature [94]. Moreover, KI values are similar between both clay fractions indicating the predominance of neoformed illite crystallization close to fracture zones. It is highlighted by the presence of illite veins observed microscopically. By comparing KI values from this study and from [19], the KI tends to decrease approaching fracture zones. Thus, the greater the fracturing, the higher was the temperature, except for NH_31_1_2 and NH_24. These two samples present higher KI values, meaning lower temperatures. It can be explained by the presence of I/S and vermiculite, forming at lower temperature than illite and kaolinite [91,92,95]

A Multi-Stage Paleo-Fluid Circulation
Hydrothermal alteration in crystalline basement rocks induces the precipitation of secondary minerals that can seal the fractures [2]. The argillic facies, described above, is also characterized by fractures filled by various kinds of secondary minerals, as a result of different fluid generations. Veins are formed from fluids that had reacted with granite [96] and transported various chemical elements. Through the petrographic analyses of the whole samples, a relative chronology between the different veins can be determined. Six vein generations, following each other, were identified (Figures 6 and 9): 1.
Quartz veins resulting in the precipitation of secondary quartz due to primary silicate partial dissolution.

2.
Illite veinlets, which have different behaviors, according to the mineral crossed. Indeed, as presented by [19,97], plagioclase is the main pathway for fluid flow due to their abundance in the rock and to dissolution pit porosity, allowing the interconnection between the pores. In quartz and K-feldspar; however, only microfractures create the porosity allowing the fluid to circulate. This explains why veinlets look like straight lines in quartz and K-feldspar and are twisted and blurred lines in altered plagioclase. Illite veinlets can be contemporary to kaolinite veinlets.

3.
A dolomitic brecciated vein embedding essentially quartz clasts, such as in sample NH_31.

4.
Carbonate veins with different compositions: contemporary ankerite/dolomite veins, intersected by siderite veins and the whole intersected by later calcite veins. They are preferentially oriented NW-SE, according to the direction of the NH. It is suggested that the fluids, having precipitated these carbonate veins, have circulated through the fracture network formed by the activity of the SDVFZ.

5.
A carbonate brecciated vein that is likely to have transported clasts of quartz, K-feldspar, and biotite, and is presenting precipitations of barite. This phase, but also barite itself, is only present in sample NH_31. This phase was probably due to a later event that reactivated the main fracture composing the sample and let a new fluid circulate. 6.
Calcite veins, as shown in the sample NH_31, which are of a different compositions from generation 4.
Two types of calcite veins and the carbonate vein of phase 5 were, thus, identified through the sample NH_31, allowing the definition of two generations of calcite. However, they do not allow a decision on if the calcite veins identified in the other samples correspond to the generation 4 or 6. CL analyses are ongoing to identify the type of calcite generation. Sample NH_31 and the multiple vein generations show how complex the fluid circulation history was in this area.

Fluid Circulation and Argillic Alteration Effects on Petrographic and Petrophysical Behavior
The effects of the alteration on the element transfers are described by the diagrams plotting major element oxides against LOI ( Figure 15) and by mass-balance calculations (Table 5) using the Gresens method [76] to quantify the losses and gains of elements during hydrothermal alteration. To this end, data from this study and from [19] were used. LOI has been defined as an indicator of the alteration degree [19,98,99] supporting optical observations. In this study, LOI values range from 2.5% to 10.9% (Table 3), whereas, Klee et al. (2021) [19] obtained values from 1.1% to 2.8% for the protolith. These values confirm a more pronounced alteration of the granite in the vicinity of fracture zones.
Concerning element transfers, only minor variations of Al 2 O 3 were observed compared to the LOI, also confirmed through mass-balance calculations. Al is thus considered as immobile. This study has shown a significant SiO 2 negative correlation with LOI. This loss is confirmed by mass-balance calculations. It indicates the partial silicate alteration, explaining the presence of quartz veins and veinlets through the granite. A negative correlation with LOI was observed for Na 2 O, confirmed by the loss of Na 2 O obtained by mass-balance calculation. Indeed, altered plagioclase are depleted in Na 2 O. No Na-bearing newly formed minerals (except a little I/S and C/S) was encountered. It is likely that Na was exported out of the alteration zone. When the alteration of plagioclase is less pronounced, as in samples NH_31_3 and NH_36 (Table 3), the Na 2 O depletion is lower. However, a global enrichment of K 2 O is observed linked to illite crystallization. A positive correlation is observed between the amount of alteration and illite formation. K-feldspar being rarely altered, chemical modelling would be necessary to determine whether the amount of K release by biotite alteration would be sufficient to allow the rather important formation of illite. MgO and Fe 2 O 3 present no clear correlation with LOI. However, low depletion of Fe 2 O 3 and MgO compared to the protolith were identified by mass-balance calculations. It is suggested that these depletions are linked to the alteration of biotite but they are compensated by their precipitation as oxides or in carbonates. Thus, it is likely that no Fe or Mg is exported out of the alteration zones. CaO and calcite content present no correlation with the LOI. Samples NH_19, NH_36 (this study), and NH_3 (in [19]) present a similar calcite content, respectively, of 6.3%, 6.0%, and 6.5%, whereas they present different LOI values (4.1%, 2.6%, and 2.8%, respectively). The case of a high calcite content associated to a low LOI can be explained by (1) physicochemical fluid composition bringing Ca, and not allowing a complete alteration of plagioclase and biotite into clay minerals [100]; (2) the residence time of the fluid not being sufficient [101]; or (3) a non-sufficient amount of fluid to allow the transformation of plagioclase and biotite into clay minerals. CaO losses are attributed to the alteration of plagioclase. However, the amount of Ca released by the alteration of plagioclase is certainly not sufficient to explain the crystallization of large amounts of carbonate into the granite veins. We infer that Ca is related to an external source as it is gained by the system (open system). As carbonate veins are related to the SDVFZ activity, we suggest that the source of the Ca is the Proterozoic series, which contain dolomite and carbonate sequences.
Fracture zones created porosity through microfractures and spaces between secondary minerals (Figure 6f) precipitated during primary mineral alteration. When fracturing increases, alteration increases too, as well as porosity [10]. However, the subsequent mineralization and chemical alteration can either decrease or increase the porosity of the rock [2]. For a similar LOI, porosity values vary significantly (e.g., NH_3 [19], NH_23, and NH_31_3). No correlation is visible between porosity and alteration or calcite content.

Relation among Fluid Circulation, Alteration, and Fracture Density
It is expected that areas along the rear southern limit of the CBS have recorded less transcurrent shear than areas situated along its northern front (Section 4.1). Recent work by Chabani et al. (2021) [42] has shown a structural compartmentalization occurring at the NH scale with a varying intensity that may influence fluid flow through the fracture network. Likewise, they identified a complex network of joints at outcrop scale, also encountered in this study at the sample and thin section scales, playing a key role in the fracture connectivity. The authors identified several fracture sets, among which, the NW/SE and the E/W oriented sets are predominant. The NW/SE trend, following the SDVFZ direction, controls the geometry in the whole range, whereas the E/W set, characterized by several short fractures, is mostly present in the central and southeastern part of the NH. The geometrical analysis of both fracture sets showed that the connectivity is ruled by fractures of different sizes. Regarding the fracture distribution, the highest fracture density is recorded within the rear southern part with values five times higher than in the northern part. This can be explained by the complex tectonic setting and the gradient of deformation observed in the entire NH range, with evidence of extreme shearing in the southern part as already presented. Moreover, through scanline analyses along drone photogrammetric outcrop profiles and ground outcrop profiles, Chabani et al. (this issue, under review) [44] show a clear correlation between fracture density and distance from a fault zone. The closer the fault zone, the greater the fracture density. There is no straightforward correlation between amounts of deformation and fracture density. Fracture density depends on a series of parameters notably rock competence, anisotropy of units and possible presence of detachment zones concentrating deformation in narrow areas.
This study, together with that of Klee et al. (2020) [19], confirm the increase of fracture density towards fault zones through several samples collected at different distances from the faults. Moreover, at the NH scale, the LOI increases as the fracture zones are approached (Section 5.4). The more intense the fracturing, the higher the alteration. Likewise, porosity values obtained for samples from Klee et al. (2021) [19] range from 2.2% to 5.2%, whereas, porosity values for this study samples range from 10.1% to 14.0%. Close to fracture zones, samples have more than twice the porosity than the samples away from them. Calcite content ( Figure 16) also show higher values near to fracture zones, highlighting the importance of the presence of fractures allowing fluid circulations and thus carbonate precipitation, which mainly filled and sealed the fractures [71,102]. Carbonate veins are present in the whole range, preferentially orientated following the NH strike, probably related to the SDVFZ activity. However, calcite content is generally higher in the southeastern part of the range. This part is characterized by a high fracture density due to complex tectonic setting, related to the GFZ activity [42]. The connectivity within this part is ruled by the small and large fractures, and an additional complex joint network leading the fluid supply toward the fault zone.
The comparison of the data from Klee et al. (2021) [19] with those in this study at the NH scale shows that the increase in fracture density is correlated with the increases in alteration degree, porosity, and calcite content, respectively. However, no correlation is observed between alteration, porosity, and calcite content all together. The sample NH_36 ( Figure 11) has shown that calcite content is not related to the intensity of alteration (Section 5.4). This highly fractured sample has been characterized by a specific fracture arrangement, composed mainly of fracture clusters, following the SDVFZ direction. It allows the fluid to circulate following these clusters, and then creating a flow anisotropy. However, only a low alteration degree is observed, which can either be due to (1) specific physicochemical fluid composition [100]; (2) an equilibrium between the fluid and the surrounding rock; (3) an insufficient residence time of the fluid [101]; or (4) a non-sufficient amount of fluid. Fracture density increase does not necessarily induce the increase of alteration (Figure 19a,b), whereas at the NH scale, a general correlation between fracture density and alteration was shown. Sample NH_36 shows that fluid can circulate through a highly fractured granite without necessarily producing alteration (Figure 19b).
Numerous fractures and veins, showing a high fluid-rock ratio, crosscut the granite altering it consequently (Figure 20a,b). Faulkner et al. (2010) [11] have shown that within crystalline rocks, the flow can be ruled by a small number of fractures within the rock surrounding a fracture zone. The fracture interconnections constitute the main parameter, which can enhance the fluid circulation. Samples NH_31 and NH_23 showed the relation between fracturing and alteration at local scale (Section 4.3.2). Thin section fracture analyses showed that the fracture density average is similar in both samples. NH_31 consists of a brecciated vein developing an alteration gradient but no fracturing gradient (Figure 20c), while NH_23, perpendicular to NH_31, consists of a mylonite developing an alteration and a fracturing gradient (Figure 20d). Regarding the stress axis orientation, it is suggested that the main stress axis σ3 is at a very high angle to the NH_31 initial quartz vein. The close to normal angle between the initial quartz vein and the stress axis σ3 was favorable to reactivation, inducing a brecciated vein without developing a significant additional fracture gradient in the surrounding rock (Figure 20c). It has been shown that remnants of a quartz vein composed the NH_31 brecciated vein. This quartz vein, developed during a past event of unknown age and origin, possibly magmatic (last crystallization fluids). It is suggested that this vein was reactivated, letting a new fluid rich in Ca circulates. Then carbonates crystallized giving a brecciated vein as a final product. The alteration gradient observed in the wall rock is especially pronounced at the borders of the brecciated vein. The proximity with the main fluid pathway induced an important leaching by the fluid. Within NH_23 (Figure 20d), it is supposed that the quartz-feldspathic unit was initially a dyke crosscutting the granitic pluton. Such dykes have been observed in the surroundings. During SDVFZ activity, reactivation of this zone of weakness permitted percolation of a Ca-rich fluid. An intensive shearing induced dynamic recrystallization developing a mylonite in the borders and its associated fracturing gradient in the wall rock. In that case, the main stress axis σ1 would be characterized by an intermediate angle to the shear plane. A fracturing gradient is observed in the surrounding rock. The fracture density induced in the surrounding rock does not further increase. Asymmetric clasts and shear bands are observed within this mylonitic layer. Fluid circulation and induced fracturing have promoted fluid-rock interactions and the granite alteration, which is more pronounced in NH_23 than in NH_31 (Figures 9h-k and 14e,f). Those two samples showed that at sample and outcrop scales, fracturing and alteration are heterogeneous, but also that the alteration degree is not always related to the fracture density (Figure 20c,d). These sample analyses allowed us to study a fault zone development depending on strain, fluid-rock ratio, and material rheology. It is shown that at low to moderate regional strain, incipient local high strain concentration is present ( Figure 20). It has been observed that:

1.
The fluid interacts with the surrounding rock, altering it pervasively, and fractures are sealed by secondary minerals (Figure 20a,b). Those veins can be reactivated by shearing [103]. Thus, it creates a brecciated vein (e.g., NH_31) or even a mylonite (e.g., NH_23), developing an alteration gradient. Sample NH_23 shows that mylonitic deformation appears at low to moderate regional strain when carbonates are involved. Fault reactivation tends to cause a mineral fill breakage and reopens the fracture [2,104].

2.
Open fractures let fluid circulate without interaction with the surrounding rock. Shearing creates a breccia without further alteration within the granite, such as in NH_36 (Figure 19).
At any scales, fracture density promotes fluid flow through the fracture network provided a fluid is present [5,11]. However, the alteration degree is not correlated to the fracture density.

Implication of Alteration in Terms of Geothermal Reservoir Properties
Fresh granite has a very low initial matrix porosity (<1%), which "does not allow" the fluid circulation, hence inducing a low permeability [9]. A granitic geothermal reservoir is considered as exploitable when it presents a connected fracture network increasing the permeability and in which a sufficient amount of hot fluid circulates [2,3,9,11]. By flowing through fractures, the fluid interacts with the rock, increasing the matrix porosity and permeability promoted by the dissolution of the primary minerals and crystallization of newly-formed phases [13]. However, newly-formed minerals can also crystallize in the fractures and seal them, reducing the porosity and the permeability [105,106]. Thus, sealed fractures are transformed from conduits into barriers to fluid flow [83,107,108] and impact the geothermal production by decreasing the connectivity between the fractures and hence the permeability [2,109].
The NH altered granite shows an important fracture network and a high matrix porosity, which allowed fluids to circulate and to interact with the surrounding rock. Indeed, numerous veins with various infills crosscut the granite indicating several successive fluid circulation events. Carbonates occupy a prominent part of the fracture system, which can be easily dissolved thanks to acid injections in order to connect the boreholes to major conductive fractures for geothermal exploitation. It is suggested that, when NH were an active reservoir at depth, fractures were only partially filled, and the fluid could flow through the fracture network. It has been seen that the stress field has a major impact on the fluid circulation. Indeed, drilling into fractures at a very high angle from the main stress axis σ1, such as exposed by sample NH_31, is favorable to allow fluid to circulate. In the event of geothermal exploitation, only minor chemical stimulation, if any, would have been necessary [13,71,110]. Moreover, the alteration amount is not necessarily related to the degree of fracturing and, therefore, the fracture network influence on fluid circulation is hardly predictable. Thus, due to the conduit-barrier role, the deformation gradient, the degree of fracturing, the alteration processes, and the relationship between them, have a large impact on geothermal production [109].

Conclusions
Noble Hills (NH) is a newly studied area, in terms of an exhumed granitic geothermal reservoir. It provides an excellent opportunity to give fundamental scientific input in 3D, allowing for better understanding of granitic reservoir behavior in a trans-tensional context. This analogue shows how complex a granitic reservoir can be, in terms of structures, fluid circulation, and fluid-rock interactions. This study proposes a geometric, petrographic, petrophysical, and geochemical description in order to characterize the influence of fracturing on fluid circulation and alteration processes.
Approaching fracture zones at large scale, the NH granite shows signs of several generations of fluid circulations resulting in successive veins of various mineralization. Fluid circulation being more important in the vicinity of fracture zones, a stronger fluidrock interaction is observed. Thus, argillic alteration prevails compared to in the protolith. It is highlighted by the increase of the LOI correlated to a Na depletion due to plagioclase alteration and a K enrichment associated to illite precipitation. Likewise, the porosity, the calcite content, and the temperature increase nearer to fracture zones. However, no correlation exists among LOI, porosity, and calcite content altogether. Moreover, a high fracture density does not necessarily imply a strong alteration (e.g., sample NH_36).
This relation among fluid circulation, alteration, and fracturing is also visible at a sample scale. The higher the fracture density, the more pronounced the alteration. However, samples NH_31 and NH_23 showed how complex this relation could be: • NH_31 shows a reactivated vein giving a brecciated vein (quartz and granitic clasts in a carbonate matrix), which induced no fracture gradient in the surrounding porphyritic granite, but developed an alteration gradient. • NH_23 shows a carbonate mylonite creeping around deformed granite lenses, which induced a fracture gradient, as well as an alteration gradient in the surrounding porphyritic granite.
These two samples, showing different deformation features as a result of a different orientation within the stress field, present a similar fracture density. Strain was accumulated within a carbonate mylonite within sample NH_23. The alteration gradient is more pronounced in this sample, but the fracture density does not change. Sample NH_36, consisting of a granite composed of a cohesive breccia, whose matrix is made of dolomite, shows a low alteration. Even locally, alteration does not always depend on the deformation gradient or on the fracture density.
This study shows that the deformation gradient observed at a large scale is also visible locally in the context of a low to moderate strain. It also highlights the importance of mineral crystallization in a geothermal reservoir. Fluid flow depends on a connected network of open and permeable fractures. In reservoirs at depth, fractures can be sealed and act as a barrier. Stimulation techniques are needed in order to reopen them so that a drain allowing a new fluid circulation is created. The influence of fracture zones on fluid flow and alteration is difficult to predict and yet impacts the production from the reservoir.
In order to complete this work, a future study will focus on the influence of a fault zone on fluid circulation and alteration processes, in the NH, at high strain conditions. Moreover, laboratory investigations could be performed to better understand the controlling parameters of the ongoing processes and fracturing of rocks.