Crystallographic and Seismic Anisotropies of Calcite at Different Depths: A Study Using Quantitative Texture Analysis by Neutron Diffraction

: Eight samples of limestones and marbles were studied by neutron diffraction to collect quantitative texture (i.e., crystallographic preferred orientations or CPO) of calcite deforming at different depths in the crust. We studied the different Texture patterns developed in shear zones at different depth and their inﬂuence on seismic anisotropies. Samples were collected in the French and Italian Alps, Apennines, and Paleozoic Sardinian basement. They are characterized by isotropic to highly anisotropic (e.g., mylonite shear zone) fabrics. Mylonite limestones occur as shear zone horizons within the Cenozoic Southern Domain in Alpine thrust-and-fold belts (Italy), the Briançonnais domain of the Western Alps (Italy-France border), the Sardinian Paleozoic back-thrusts, or in the Austroalpine intermediate units. The analyzed marbles were collected in the Carrara Marble, in the Austroalpine Units in the Central (Mortirolo) and Western Alps (Valpelline). The temperature and depth of development of fabrics vary from <100 ◦ C, to 800 ◦ C and depth from <10 km to about 30 km, corresponding from upper to lower crust conditions. Quantitative Texture Analysis shows different types of patterns for calcite: random to strongly textured. Textured types may be further separated in orthorhombic and monoclinic (Types A and B), based on the angle deﬁned with the mesoscopic fabrics. Seismic anisotropies were calculated by homogenizing the single-crystal elastic tensor, using the Orientation Distribution Function calculated by Quantitative Texture Analysis. The resulting P- and S-wave anisotropies show a wide variability due to the textural types, temperature and pressure conditions, and dip of the shear planes.


Introduction
The quantification of the seismic response of rocks is a fundamental task in understanding Earth's structure, from the core to the surface [1]. The large-scale seismic experiments, now with high-resolution arrays, is the most used approach to image the Earth and to resolve the distribution of natural resources in-depth, e.g., water, ore minerals, and oil & gas [2,3]. However, the interpretation of seismic images strongly relies on the knowledge of the seismic response of the anisotropic aggregates of minerals composing the rocks, which in turn is intimately related to the textures of rocks [4], as well as on their extrinsic shape preferred features [5]. In the last two decades, a great effort has been made to quantify the seismic response of natural aggregates, chiefly by using the 2-Dimensional approach of the EBSD (Electron Back Scattered Diffraction) to reconstruct the Orientation Distribution Function (e.g., [6][7][8]). Carbonate rocks occur at the Earth's surface and upper crust levels, as sedimentary cover, but can be found within the intermediate and lower crust, as marbles. Within the upper crust, the tectonic deformation often localizes within carbonate rocks, producing thrust-folds systems [9,10]. Though less diffuse within the metamorphic basements, carbonate rocks are known to be equally crucial at depth because they likely represent weak horizons, often separating tectonometamorphic units [11,12]. For these reasons, the knowledge of the texture and seismic response of carbonate rocks is essential to interpret seismic anisotropy at various depth within the crust [13]. Several works have been dedicated to studying the behavior of calcite single crystal, limestones, and marbles at different temperatures and pressures [14][15][16][17][18][19]. A general scheme of the expected textures has been developed, combining natural and numerical data (see [20] and references therein), but data are still needed. Moreover, virtually no work has been dedicated to combining the neutron diffraction textural data with the prediction of seismic velocities through the homogenization of the elastic tensor [21,22].
With this contribution, we aim at producing new data related to the textures of carbonate rocks. In particular, we aim at relating different pressure, temperature, and shear geometry conditions of deformation with the type of texture developed. We also aim at investigating the control of texture types on the seismic anisotropies. The presented study will hence provide a reference to interpret measured seismic anisotropy in the continental crust [23], where carbonate rocks might be also involved in accommodating shearing along mylonitic horizons.
For this aim, we selected carbonate rocks characterized by fabrics developed at different temperature and pressure conditions and different strain geometries, from pure to simple shear, representing some of the large number of possible combinations in naturally deformed rocks. The samples have been collected mainly in the Alpine system and display a various degree of planar fabrics, developed at different depths within the crust. In this contribution, we will produce new texture data to be included in the general scheme of texture development. Moreover, we will calculate the component of the seismic anisotropies due to the texture by homogenizing the stiffness tensor using the Orientation Distribution Function obtained with the Quantitative Texture Analysis. Finally, we will investigate the influence of the orientation of the foliation plane on the seismic anisotropy [24] at different depth in the crust and with different textures.

Samples Description and Geological Setting
Samples were collected in the French and Italian Alps, Apennines, and Paleozoic basement ( Figure 1). Different locations ( Figure 1) were chosen to sample carbonate rocks from different crustal levels, from upper/intermediate crust limestones to lower crust high-temperature marbles (Table 1).  According to Figure 1, upper-crust samples were collected in the Southalpine domain, Sardinia and Helvetic-Dauphinoise (e.g., Triassic and Jurassic limestones). All samples preserve at macroscopic scale their primary features; intermediate crust samples were collected in the Briançonnais domain and in the Central Austroalpine domain, where Mesozoic sediments were involved in the Alpine tectonics at relatively upper levels. However, they still preserve their stratigraphic relations with the surrounding rocks, as sedimentary strata or laminations. The last group of samples, collected in the Carrara metamorphic complex, and in the Western-Central Austroalpine domain, are metamorphic marbles. In this last group, no sedimentary features are preserved.

•
The BAS sample is a light gray crystalline limestone, fine-grained, homogeneous, and free of any planar fabrics at the macroscopic scale, both primary or secondary, as fractures or joints. It is Jurassic in age and it belongs to the External thrust-belt French Subalpine system [25,26]; it is part of the so-called Calcaires "tithoniques" [27]. It was sampled close to Grenoble (France), La Bastille, the southernmost part of the Chartreuse massive. According to the tectonic reconstruction and the local geologic map [27,28], the thrust system developed a km-scale anticline, associated with brittle faults and fracture systems. The BAS sample has been collected within this anticline but far from the brittle system. The fold and thrust system occurred at <3 km depth, at temperature <100 • C.

•
The SOD sample was collected in the Southalpine domain. The sample is part of the Triassic cover, which crops along the entire Southalpine domain. The general interpretation infers the thrust and fold system to convergent tectonics active during the formation of the Alpine system (Figures 1 and 2). It was collected within the meter-scale mylonitic horizons associated with a meter to tens of meters asymmetric folds [29]. According to these authors, the thrust horizons occurred at relatively upper depth, not exceeding 5-7 km, corresponding to a maximum of 200 • C. The studied sample shows an abrupt change in fabric from undeformed to highly strained domains (Figure 2), often associated with grain-size reduction, pressure-solution, and plastic deformation in carbonate grains, testified by mechanical twinning. These mesoand microstructures support temperatures as high as 200 • C [30,31]. The kinematics of these thrust horizons have been reconstructed based on stratigraphic and structural constraints, being characterized by strong simple shear components also testified by classical shear indicators, as cmto mm-scale drag folds and porphyroclasts geometries [32]. • SA1 sample was collected in the Sardinia basement. It is part of the Ceroidi Limestone (Gonnesa Fm.), pre-Ordovician protolith of the External Zone. This metasedimentary cover has been deeply involved in the Paleozoic-Variscan tectonics [33]. During the Variscan tectonic activity, green-schists metamorphic conditions were attained [34], locally associated with high strain shear horizons, from meters to tens of meters thick [35,36]. The SA1 sample was collected in one of these high strain horizons, made by meta-carbonate limestone. The Variscan deformation developed as pervasive folding systems. However, the SA1 sample was collected far from folds and within the domain where the mylonitic foliation is well visible at the macroscopic scale ( Figure 2). As qualitative estimates for pressure and temperature, the intermediate limit can be put following the metamorphic evolution of the inner part of the belt, as described by Elter et al. [37]. Consequently, a limit of 350 • C and 4-5 kbar, corresponding to roughly 12-15 km depth can be used to constrain the development of the fabric on this horizon. Moreover, mesoscopically no asymmetric shear sense indicators have been recognized, supporting the microstructural and textural observations [35], which suggested a pure shear strain. • STE sample was collected in the Central Alps, within the intermediate Austroalpine domain [38].
Here the Mesozoic sedimentary cover was involved in the Alpine tectonics, occurring as thick-skinned thrusting of cover and its pre-alpine basement [39]. Thrusting localized along high strain horizons within cover and basement [40]. The temperature conditions of this event were estimated at T = 300-350 • C at depth = 15-20 km [41,42]. We collected a sample within the mylonite horizons in the Mesozoic cover, namely the Fraele Fm. It is characterized by a strongly developed mylonite fabric, associated with a meter-to tens of meter folds. A strong transposition of the lithostratigraphic features has been recorded along this horizon [43]. A simple shear component has been described for this deformation based on meso-and microstructural analysis [41]. Sample scale features, as mm-sized drag faults (Figure 2), support the simple shear component. • CC sample was collected in the Western Alps, in the sub-Briançonnais domain. It consists of Jurassic limestone strongly deformed during Alpine thrusting of the external part of the chain [44]. According to [45,46] the External Briançonnais units equilibrated at a pressure not exceeding 5 kbar, corresponding to depths <15 km, and T = 300-350 • C. Moreover, overall simple shear kinematics is constrained by map-to micro-scale indicators, extensively described in the area [47][48][49]; however, the studied sample does not show clear microstructures that could be uniquely related to a simple shear geometry ( Figure 2). • P1 was collected in the well-known Carrara marbles quarries. They consist of the Mesozoic cover of the Autochthon of the Tuscan Units metamorphosed during the Alpine evolution [50,51]. Km-to meter-scale folding systems have been mapped, and several kinematics reconstructions have been proposed [52] suggesting large scale shear, related to the thrusting of the deep units, resulted in meso-and micro-scale simple shear fabrics (e.g., [20,53]). However, large scale strain partitioning also occurred, allowing a large volume of marbles to escape the development of planar or linear fabrics [54]. In fact, the studied sample is characterized by isotropic equigranular texture of recrystallized calcite (Figure 2), suggesting a static recrystallization process associated with grain-boundary migration both leading to grain size increase [32]. Within the Alpi Apuane Metamorphic Complex the estimated conditions of metamorphism point to T = 400-500 • C at 8-10 kbar [50,55,56], corresponding to 25-30 km depth. • VP3 and MA1 were both collected in the Austroalpine domain of the Alps [44]. The Austroalpine domain is interpreted as fragments of the pre-alpine continental crust involved in the Alpine subduction-collision system. The collected samples are marbles whose metamorphism has produced during high-temperature Permian-Triassic extension [57]. VP3 was collected in the Valpelline Series, Austroalpine of the Western Alps [58,59], while MA1 was collected in the Languard-Campo, Austroalpine domain of the Central Alps [60]. They are characterized by 0.5-1.5 mm (VP3) to <0.5 mm (MA1) grain-size and strong planar fabric, marked by the shape preferred orientation of calcite and local diopside and quartz ( Figure 2). They display diffuse mechanical twinning and undulose extinction, both suggesting grain-scale plastic deformation [32]. Shear indicators are not univocal, MA1 often shows asymmetrical bending of calcite tails around porphyroblast suggesting a component of simple shear, but pure shear cannot be ruled out. Similarly, VP3 displays a less tendency to develop asymmetric microstructures, point to a greater contribution of pure shear. Metamorphic conditions have been quantitatively determined for both samples, VP3 developed at T = 700-800 • C and P > 6 kbar, corresponding to >18 km depth, while MA1 fabric was estimated to form at T = 600-750 • C and P = 6.5-7.5 kbar, corresponding to 20-24 km depth [58,60].

Sample Reference System and Methods
The samples were cut in cubes of ≈1 cm edge ( Figure 3); the three axes of the cube correspond to the orthonormal reference of the sample. In this orthogonal sample reference system, the plane XY corresponds to the stratification or shear plane, while X direction is taken parallel to the flow direction or lineation if it exists ( Figure 3). The references in Figure 3b are used in all the Pole figures (Figures 4-6). Pole figures display pole densities in multiples of a uniform distribution (mud).    We used the non-disruptive method of the Quantitative Texture Analysis (QTA) by neutron diffraction. This method allows the evaluation of the Texture (also called Crystallographic Preferred Orientation, CPO) of the samples. It has been successfully applied to various types of rocks and synthetic materials: monomineralic quartzites [61], marbles [20], limestones [62], dunite [63], or glaucophanite [64], in poly-phasic rocks, as amphibolite [65,66], quartz-feldspatic mylonitic orthogneisses and gabbros [67,68], subaqueous lavas [69] or sandstones [70]. This procedure uses the high penetration and high flux of neutrons available at the nuclear reactor at the Institute Laue-Langevin (Grenoble, France) allowing the measurement of samples with volumes of approximately 1 cm 3 in one to eight hours [71]. Here, we present the results from the D1B, D19, and D20 diffractometers (http://www.ill.eu). The raw data have been analyzed using the software package MAUD (Materials Analysis Using Diffraction [72]): diffracted intensities are used to calculate the Orientation Distribution Function (ODF [73]), from which the representative lattice planes are extracted and represented as Pole Figures (PF). PFs were represented showing the reference axes from the Sample Reference System (X, Y, Z in Figure 3). The coverage of the three-dimensional sample space has been acquired differently at the three beamlines (D1B, D19, and D19), due to their technical specifications. Briefly, D1B works at lambda = 2.52 Å, equipped with a position-sensitive detector, covering a 2θ range of 128 • ; D19 is equipped with a very large (120 • × 30 • ) 2-dimensional position-sensitive detector, and works best in the wavelength range 0.8 to 2.4 Å; D20 offers wavelength range 1.4 to 2.1 Å and the position sensitive detector covers 153.6 • . The respective angular settings are reported in Table 2. The ODF was calculated using the E-WIMV approach [73]. The reliability of the refinement is shown by the refinement parameters reported in Table 3. Calcite lattice parameters used for the refinements are a = 4.9849 Å c = 17.0479 Å. 2Theta offset, sample position, and background parameters were also refined.  The seismic properties were obtained using petrofabric data, as extensively explained by [7]. Here we calculated the sample seismic properties by averaging the single-crystal elastic tensor [77] weighted by the ODF, as obtained by quantitative texture analysis by neutron diffraction. In practice, the ODF calculated using the Maud procedure was used to homogenize the single-crystal tensor in Maud. Results for ODF refinements are reported in Table 3. Although various averaging procedures can be adopted [78], we used the arithmetic mean of Hill [79]. The resulting homogenized tensors (Table 4) were used in MTEX Matlab Toolbox [80] to compute and plot seismic properties [81] shown in Figure 6. The same averaged tensors have been then used to calculate the seismic anisotropies at different dip angles of the foliation plane, simply rotating the tensors and changing the plane of projection for the pole Figure in  Figure 7 are plotted in a way that only the Y-axis remains constant and the observer is looking down from the Earth's surface having the macroscopic foliation plane rotating with respect to the Y-axis.  As a general consideration, this approach to computing seismic velocities generally produces higher values when compared to natural samples. The most likely reason for this is that other important factors, intrinsic and extrinsic, controlling seismic velocities, and anisotropies are not included in the calculation. These factors are fluid pressure, crack density and their orientation, porosity, grain boundary geometry, and their shape preferred orientation [5,53,82,83].
A thermodynamic approach [66,84] was also used to calculate the seismic properties of a modeled carbonate rock at different pressure and temperature conditions, roughly corresponding to the estimate conditions of texture development of the studied samples. We used the PerpleX software package (version 6.8.6 [85]) to compute the seismic parameters as dependent on the thermodynamic function G (molar Gibbs free energy), which is minimized to establish phases, amounts, and compositions stable as a function of pressure and temperature. The software requires a chemical composition of the system as input. Software and data used for the calculation of phase relations and seismic velocities are available at http://www.perplex.ethz.ch. The graphical outputs were obtained using pyWerami source [86]. The calculated seismic properties correspond to those of a polycrystal rock with no preferred orientations. An averaged composition of carbonate rocks has been used as the model chemical system [87].  [88] and Perplex datafile repository. Holland and Powell (1998) [88] and successive update [89] (datafile hp62ver at http://www.perplex.ethz.ch/datafiles/) was used as a thermodynamic database. A typical output ( Figure 8A) consists of pressure and temperature grid (in this calculation, P = 1-15 kbar and T = 100-1000 • C) where the stable phases and assemblages are shown (i.e., pseudosection). Using the application Werami, which is part of the PerpleX package [85], specific properties were investigated in the P-T space and plotted in Figure 8B,C. Namely, we investigated the variations of modes, density, Vp and Vs of all phases within the investigated P-T ranges. Input data for PerpleX calculation may be found as Supplementary Materials.  Figure 4 reports the calculated PFs for a single crystal-like of calcite, an aggregate of equigranular randomly-disposed crystals, and an aggregate of equigranular grains, by a statistical Fisher-Bingham distribution.

Textures
They show the contributions to Texture expected by the crystallographic symmetry only or by textures: (i) single crystal distribution is characterized by a strong c-axis maximum and a three-fold axes distribution in the plane normal to the c-axis; (ii) oppositely, a random or uniform distribution of equigranular aggregate will reduce or obliterate the crystal symmetry anisotropy, producing a random statistical distribution, in any directions of the 3d space; (iii) lastly, an imposed (e.g., tectonic stress) distribution, here represented as the Fisher-Bingham distribution [91], will produce strong maxima that follow neither the crystal symmetry and nor a random distribution.  (100)) and m-planes (110) since these orientations better describe the deformation mechanisms active during plastic deformation at lattice scale and are commonly reported in the literature. Table 3 reports the R factors and F 2 factor, commonly used in texture analysis to assess the quality of ODF refinement. In general, the samples show different texture distributions and densities but they can be divided into 2 principal types: a random distribution (01BAS, 06P1); c-axis textured distribution (02SOD, 03SA1, 04STE, 05CC2, 08VP, 09MA). The c-axis textured distribution may be further divided into orthorhombic (02SOD, 05CC2, 09MA) and monoclinic symmetry (03SA1, 04STE, 05VP) types. The orthorhombic type (A) is characterized by the parallelism between the (001) poles and the Z direction. On the other hand, the monoclinic type (B) shows an angle between the (001) poles cluster and the Z-direction. Symmetrically, the (100) poles, which generally display a girdle distribution, may lay orthogonal to the Z direction or with an angle, generally 30 to 45 • . F 2 values (Table 3) may be used as the overall index to assess the randomness of the texture [92,93]. As expected for poorly textured samples, 01BAS, 06P1 have F 2 close to 1 mrd 2 , even though 01BAS displays a minimum amount of preferred orientation. 02SOD, 04STE, 05CC2, and 09MA are characterized by well-developed textures, though the F 2 index is relatively low, being between 1.10 and 1.26 mrd 2 . Higher F 2 values are shown by 03SA1 and 09MA. Figure 6 reports the calculated seismic velocities. Namely, P-wave Vp (km/s), S-wave anisotropy and Vs1 and V2 polarization. Their distributions closely resemble the main texture feature ( Figure 5). In fact, since the c-axis direction is the slowest Vp in the single crystal setting, and the a-axes directions are close to fastest (Figure 3), the general feature that shows up is that the pole to (001) distribution correspond to the slowest directions in the Vp plots. So, the prediction of the slowest direction in carbonate rocks is relatively straightforward since it is influenced by the c-axis distribution. In contrast, S-wave anisotropies are more complex to predict since they result from a combination of textures ( Figure 5). Specifically, 02SOD, 06P1, 08VP, and 09MA display a girdle distribution of the poles to (100) and (110), which correspond to a pronounced direction of fastest propagation. Conversely, 01BAS, 03SA1, and 05CC2 are characterized by single or multiple clusters for S-wave anisotropies, giving a more complex pattern of polarization. 04STE displays a mixed distribution, where a weak fastest S-wave girdle links the two clusters.

Texture Types Evolution with Temperature, Pressure, Shear Geometry, and Strain
Texture data reported in Figure 5 fall into the types of calcite texture most commonly described for natural rocks [94] and reported in experimental results [16,17,95,96]. In particular, we separated two types: random, where no or a little preferred orientation is found, and textured, characterized by a strong c-axis maximum and a-axes girdle distribution. Textured types have been further separated in Type A and Type B, to emphasize the relations between Texture and shear plane, as described for pure shear (i.e., orthorhombic geometry, Type A) and simple shear (i.e., monoclinic geometry, Type B) [20,94,97]. Type A and B have been used to infer the deformation regime and the sense of shear (e.g., [98][99][100][101]). However,   [102] demonstrated the risk of simplifying those textural relations with deformational regimes. In fact, the most common calcite deformation mechanisms may allow, even at low temperature, a quick transition from monoclinic to orthorhombic symmetries, even under a simple shear regime and at low strain (>2). Similar observations have been made on naturally deformed marbles [103]. Accordingly, the studied samples also show that discriminating between high temperature (HT) and low temperature (LT) textures is not an easy task. For example, 02SOD, 05CC2, and 09MA1 samples, though their textures developed respectively at ≈200 • C, 300-350 • C, and 600-750 • C (Table 1), they all display a similar orthorhombic symmetry; specularly, 04STE, 03SA1, 08VP share the same texture but developed at different temperatures, 300-350 • C, ≈350 • C, and 700-800 • C (Table 1). Literature data, between 100 and 400 • C at increasing shear strain, suggest a general evolution from Type A texture, generally far from the shear zone core, to Type B textures and, more often, the coexistence of Type A and B [82]. Type A seems to become dominant at very high shear strain values and very high pressures [16,104]. At low temperature (<300 • C) and low pressure (<3 kbar), as in carbonate fault gouges or thrust systems, Type B is generally favorite [105][106][107] while A and B types start to compete where temperature, pressure, and shear strain all overcome certain thresholds. These threshold values are hard to quantify; most likely pressure should exceed 3-4 kbar and shear strain 1-2 [13,16,82,97,104,107]. Similar competing textures were also described for dolomite deforming between 240-300 • C under simple shear geometry [108]. Besides, the recent experimental work of Schuster et al. (2019) [16] shows that very high shear strains (>80) tent to stabilize Type A texture that replaces early A or B types. This experimental work also strongly supports the findings by Ebert et al. (2007) [104] that constrained similar evolution from field observations. In particular, Ebert et al. (2007) [104] showed evolving random to Type A and B textures as temperature and shear strain increase. Ebert et al. (2007) [104] also showed textures similar to those of 04STE and 03SA1 samples where the angle between the c-axis maximum and the Z-axis is much higher than typically described.
Similar considerations may be done on metamorphic marbles produced at intermediate temperature (≈400-600 • C) and pressures related to subduction-collisional tectonics (≈10-20 kbar). In particular, the large literature available on the Carrara Marble allows us to constrain textural types with shear strain geometry and grain size. Large (1-2 mm) to intermediate (0.1-0.5 mm) grain size marbles, characterized by granoblastic textures may show random texture [53,109] as well as Type A [52]. The first most likely due to static annealing or larger scale strain-partitioning and the second associated with grain boundaries migration and subgrain rotation dynamic recrystallization processes [32]. In marbles showing shape preferred orientations and microscopic foliation, the Type A is better developed [20] suggesting a high shear strain condition. Experimental work on the Carrara marbles under torsion [110], ranging between 500 and 700 • C and shear strain from 0 up to 50, support the observations [111] that Type A easily stabilizes at high strain conditions replacing Type B, no matter if coaxial and non-coaxial. Again, this process seems to be more likely occurring at high temperatures (>500 • C) [111]. Barnhoorn et al. (2004) [110] also showed that at low temperature (≈500 • C) random texture may replace Type B textures increasing shear strain. Other intermediate temperature marbles (400-500 • C), from elsewhere (e.g., [11,103,112,113]) support this shift from B to A types texture. However, no shear strain constraints are available to include these observations in the above scheme. Higher temperature (>600 • C) marbles, similar to samples 08VP3 and 09MA1, naturally [22] or experimentally [14,110] deformed do not show a strong shift from B to A type texture. These observations may be taken to support the idea that shear strain most likely controls the textural type, more than temperature. However, more observations are needed to build a valid general scheme. More observations are also needed to produce a solid scheme to relate deformation mechanisms with the observed textures. The texture observed in these samples does not uniquely relate to specific deformation mechanisms. Type A and B textures are present in samples where intra-crystalline slip mechanisms are relevant, as 08VP3 and 09MA1. On the other hand, where twinning is present, together with grain size reduction deformation mechanism (02SOD, 04STE, 05CC2), both Type A and Type B may develop [101]. Moreover, the combination of sub-grain rotation mechanisms, known to be important at similar conditions [16], and grain-boundary migration may have contributed to the observed random texture in P1.

Seismic Anisotropy
Seismic velocities and anisotropy of carbonate rocks in the studied samples change with texture types (Figure 6). Vp velocities vary from a minimum of 6.1 km/s to a maximum of 6.8 km/s; AVp anisotropy varies between 0.3 and 10.4 ( Figure 6). Vs varies from 3.3 and 3.6 km/s, while Vs1 anisotropy varies between 0.3 and 6.8. The direction of maximum Vp varies from parallel to X direction of the fabrics (05CC2, 02SOD, 09MA), corresponding to Type A texture, to more than 45 • in the XZ plane (04STE, 03SA1, 06P1, 08VP3), for Type B texture. Vp displays a general girdle distribution that roughly describes a plane of max Vp. Max Vp plane changes its orientation even within the same group of Vp max direction. Vp max in Type B samples has an angle with X direction that varies between 20 to 50 • . More interestingly, the Vp maximum frequently dips about 20-30 • away from the circumference (XZ plane). Similarly, the Vp max girdle plane dips between 20 to 50 • from the plane of shear. Besides, Vs velocities frequently show a double to triple maxima distributions. A large literature is available for sedimentary rocks, though the large majority is dedicated to shales and sandstones (see complete review [7]). Similar Vp and Vs values have been measured and calculated for calcite [13,22,82,114,115] at different temperature and depth conditions. Vp velocities most frequently show a girdle distribution of maxima variably dipping to the plane of shear, recalling Type A and Type B distributions. In general, the Vp and Vs velocities overlap with those calculated for our samples, a part for those samples from very shallow depth (<1000 m) [114,115] where only after reducing porosity and increasing cementation the Vp values go from low (<5 km/s) to higher values, more similar to those here measured and calculated (>6 km/s). These observations support the general idea that carbonate rocks anisotropy also contributes to the overall anisotropy of the crust [7,22,82], at any depths. To further investigate the nature of these anisotropies we calculate the seismic isotropic component by the thermodynamic modeling described in chapter 3. The thermodynamic modeled seismic velocities, Vp and Vs (Figure 8), describe a range of Vp wave velocities from ≈ 6.0 km/s to 6.5 km/s. Vs velocities vary from 2.6 to 3.0 km/s. The P-T plots in Figure 8B,C are characterized by three main parts where Vp and Vs vary progressively as a function of the changing mineral assemblages and mode proportion, as described in the pseudosection of Figure 8A. The three sectors are characterized by (i) the presence of aragonite instead of calcite a high pressure and low temperature; (ii) the association of calcite and dolomite at an intermediate temperature and pressure from low to high; (iii) the stabilization of clinopyroxene at high temperature, for low to high pressures. In general, Vp and Vs wave velocities decrease moving toward higher temperatures and increase moving toward higher pressures. According to the pressure and temperature constraints for the studied samples (Table 1) we are below the aragonite phase transition and the Vp should vary from 6.0 to 6.5 km/s, and Vs from 2.8 to 3.0 km/s. The calculated values from texture analysis [68,84] as well as those reported in the literature, differ by about 10% from these calculated values. These values are generally closer to measured seismic velocities taken orthogonal to the maximum shear direction (i.e., lineation) [53] or in the more isotropic samples or very shallow sample, dominated by low porosity (<10%) [115]. Consequently, these Vp and Vs can be used as a reference for the isotropic component of carbonate rocks, to investigate the amount of anisotropy to be expected when observing real seismic data for similar rock composition. In general, the analyzed samples show an increase of seismic velocities due only to textural anisotropies that can be about 0.2 to 0.8 km/s for Vp and 0.3 to 0.6 Vs. These results show that seismic properties (e.g., velocities and anisotropies) even for a single rock type (e.g., carbonate rocks) may strongly change in the crust, as a function of temperature, pressure, and shear strain, as a result of texture types that strongly control seismic propagation. Similar considerations have been proposed for other rock compositions. In these rock types, specific mineral phases showed preferred orientations that produced relevant deviations from the isotropic component of the seismic properties (e.g., [7,13,23,84]). This evidence suggests that at depth there is important overlapping between different rock types in term of their seismic properties, such that unambiguous determination of geological information from natural seismic data imply a combined use of many parameters [116]. In addition, it is clear that also the orientation of the macroscopic fabric, shear foliation, lineation, or any fabric associated with texture in rock-forming minerals play an important role in defining seismic anisotropies [24,84].
To further investigate the relevance of the orientation of the shear plane at different depth in the crust, we used the approach of Ko and Jung (2015) [24] by rotating the shear flow plane around the Y-axis. In Figure 7 the results are shown by changing the plane of view, from XZ, used in Figures 3-6, to XY. Figure 7 shows the progressive variations of the seismic properties by changing the dip angle from 0 to 90 • (0, 30, 45, 70, and 90). From the resulting rotated tensors, the components of Vp, Vs1, and Vs2 were read parallel to the horizontal component of the X direction (X' in Figure 7), as well as parallel to the normal to the Earth's surface (T in Figure 7). Random textures are characterized by a constant Vp and Vs2 maxima parallel to the Y direction and the lowest direction being within the X'T plane ( Figure  9). The lowest velocity is always recorded parallel to the vertical T direction, at any dip. A relevant deviation is found for the Vs1 where the highest velocities parallel to T are recorded at dip > 45 • , while parallel X' at dip < 45 • (Figure 9).  Figure 9).  Figure 10 represents the Vp velocities at different depths, at varying dip. From 0 to 15 • dipping of the shear planes, as those generally developed along the flat part of thrust systems [9], the Type A textures produce higher Vp velocities propagating horizontally, parallel to X' component ( Figure 7) and slower parallel to the normal to the topographic surface (T component in Figure 7). This setting does not change with depth, while the fastest Vp velocities are recorded at depth (i.e., lower crust). This pattern is reversed at different dip angles: upper crust at 50 • , intermediate crust between 20 and 30 • , and the lower crust at 80 • . At intermediate crust conditions, this flip may correspond to the evolution of the thrust system from flats to ramps. The change of dip corresponds to a substantial increase of Vp velocities parallel to X', from 6.4 to 6.7 km/s. At the same depths, it corresponds to a decrease parallel to T for rocks having a Type B texture pattern. Oppositely, at depth, both Type A and B produce an overall decrease of horizontally propagating velocities, from 6.8 to 6.3 km/s and from 6.6 to 6.2 km/s, respectively. Type B reaches its minimum at 50 • and then increases slowly approaching vertical dips. In general, steep dipping shear zones, as those in transtensional or transpressional systems (e.g., positive and negative flower structures), will produce slower or at maximum equal Vp velocities propagating horizontally (X') than shear zone within the same system but with lower dipping angles ( Figure 10). Figure 11 represents the evolution of Vp seismic velocities (vertical component T) with depth at changing dip angles. Figure 11 also reports the predicted P-wave velocities in the continental crust as modeled by Lloyd et al. (2011) [23] for "rock recipes". Vertical red boxes correspond to Vp ranges for middle to lower crust proposed by Rudnick and Fountain (1995) [117]. Two aspects arise from this image. Firstly, carbonate rock may overlap with average rock types or "recipes" at various depth in the continental crust. Large discrepancies to these overlapping may arise by the dipping of the shear planes. At similar depths, Vp velocities may differ of about 0.3-0.5 km/s, only due to different dipping. Secondly, a specific pattern is reproduced by Type A texture, at any depth. The type A texture shows that it is characterized by an increase of Vp velocities as the dipping angle increases. This pattern is reproduced at the three different depths our samples represent. Moreover, Types B displays more variability in the way seismic velocities change with the dipping of the shear plane.   [23]. Vertical red boxes correspond to Vp ranges for middle to lower crust Rudnick & Fountain (1995) [117]. Random, A and B texture Types described in the text.

Conclusions
• Eight samples of limestones and marbles were studied by neutron diffraction. We collected Texture of naturally deformed calcite at different depths in the crust, from upper to lower crust. • Different Texture patterns were recognized, from random to strongly textured (Type B); the latter, with orthorhombic (Type A) or monoclinic (Type B) symmetries. Seismic anisotropies were calculated using the Orientation Distribution Function to homogenize the elastic tensor of calcite. • A wide variability of seismic anisotropies arise from the various textural types, crustal positions, and dipping of the shear planes.

•
These results may be valuably used to guess geometry of shear zones at upper to depth for carbonate rocks or at deeper crustal levels where thick carbonatic horizons are involved localizing deformation ( Figure 11). • This work suggests that at depth, even for carbonate rocks, there is an important overlapping between different rock types in term of their seismic properties ( Figure 11); unambiguous determination of geological information from natural seismic must consider the orientation of the macroscopic fabric, shear foliation, lineation, or any fabric associated with texture playing a role in developing seismic anisotropies [24,84].

•
More work still needs to be done to include other parameters influencing seismic anisotropies, intrinsic and extrinsic, to build a database of natural rocks properties [5][6][7]13,23].
• Lastly, texture types may produce unique seismic velocities pattern. For example, Type A ( Figure 11) produces a unique pattern at any depth, characterized by an increase of Vp velocities (e.g., 6.3 to 6.8 at 20-25 km depth) as the dipping angle increases. Future studies may show if these singular patterns might be used to better constrain the interpretation of seismic profiles.