Thermo-Physical and Geo-Mechanical Characterization of Faulted Carbonate Rock Masses (Valdieri, Italy)

: Water in rock masses is a key factor in geo-mechanics, hydrogeology, mining, geo-thermics, and more. It is relevant in interpreting rock mass behavior (e.g., water-rock interaction or slope stability), as well as in deﬁning heat transfer mechanisms. Pointing out the contribution of secondary porosity in increasing advective heat transfer instead of the conduction phenomenon, this study aims to highlight a different thermal response of sound rocks and faulted zones. Moreover, it provides some methodological suggestions to minimize environment disturbance in data collection and a robust interpretation of the results. An interesting outcrop was identiﬁed in a carbonate quarry near Valdieri (north-west Italian Alps): it was studied coupling a geo-mechanical and a thermo-physical approach. In particular, geo-mechanical and photogrammetric surveys, InfraRed Thermography (IRT), and Thermal Conductivity (TC) measurements were conducted. The rationale of the research is based on the fact that, when a substantial temperature difference between ﬂowing groundwater and rocks was detected, IRT can reveal information about geo-mechanical and hydrogeological properties of the rock masses such as a degree of fracturing and joint interconnection. A comparative ﬁeld and laboratory analysis using different devices enabled a more detailed insight providing values in both dry and wet conditions. A different thermal response was highlighted for the cataclastic zone as well. IRT results showed an evident inverse relationship among the number of joints per meter and the detected surface temperature. This is probably caused by the higher water ﬂow within the cataclastic fault zone. Moreover, low fractured portions of the rock mass presented higher cooling rates and conducted heat far more than those with poor geo-mechanical characteristics (difference up to 40%). A negligible ratio between wet and dried thermal conductivity (about 1%) was also detected in lab measurements, which conﬁrmed that primary porosity is not usually relevant in inﬂuencing thermal properties of the sound rock.


Introduction
Groundwater flow in fractured rock masses has been studied by many earth scientists. It reveals one of the most important factors to be considered in slope stability problems [1][2][3] as well as in the hydrogeological assessment and in heat transfer mechanisms [4,5]. In particular, understanding geo-mechanical features, thermo-physical patterns, and their relationship with groundwater flow, is particularly important when facing non-homogeneous rock masses (flysches, calcschists, or faulted hard rock lithologies). The presence of groundwater flow can significantly improve the heat transfer mechanism [5], which increases the general heat budget. If the groundwater flow is negligible, the subsurface system is heat-conduction dominated and the amount of heat stored and removed depends only on the rate at which heat is conducted through the geologic media. When the groundwater flow is relevant, heat removal via conduction can be negligibly small because heat is rapidly dissipated through convection or advection [6]. A relevant advective heat transfer can be ruled by the presence of fractures acting as secondary porosity and making the rock mass even much more permeable than the matrix.
In order to understand groundwater paths, as well as its discharge and role, a complete characterization of rock masses is essential, especially regarding the definition of orientation, aperture, persistence, and interconnection of discontinuities [7]. This kind of information is usually obtained through geo-mechanical surveys. On the other hand, since water flow with a different temperature than surrounding medium can produce a thermal anomaly, a supplemental remote survey technique as the Infrared Thermography (IRT) coupled with Thermal Conductivity (TC) measures can provide useful data too.
IRT is widely applied in civil engineering and architecture especially for building diagnostics and insulation inspections [8,9]. In Earth Sciences, it is mainly used for the surveillance of volcanic systems [10], to calculate lava discharge rates [11] and fumarole emissions [12,13]. Reference [14] resume an interesting application, such as the study of heat capacity of rocky walls and monitoring of waste disposal sites. Moreover, some authors [15][16][17] used IRT to characterize the thermal response of rocks and joints, as well as landslide mapping [18][19][20][21]. On the other hand, TC measurement techniques are usually applied in engineering to test buildings materials, such as concrete and insulations [22], in the electronics industry [23], in developing new materials (e.g., nanomaterials and polymers, [24]), and for fire prevention [25]. Furthermore, thermal conductivity methods are widely employed in the Earth sciences such as testing heat transfer from rocks and soils [26] to evaluate thermal evolution of oil-bearing basins [27], to assess the thermal potential for geothermal exploration [28][29][30], or aim at the study of the Earth's thermal regime [31,32]. This study aims to recognize different thermal behaviors between sound rocks and faulted zones.
In this context, an interesting test site located in a carbonate quarry in the Argentera Massif (Southern Piedmont, Italian Western Alps) was identified. In this case, a well exposed fault with a clear cataclastic belt is easy to reach and to photograph, so it is a very good site where both the cataclastic fault zone and the sound rock can be compared. In detail, the following was completed.

•
Geo-mechanical and a photogrammetric survey of the faults and the surroundings to characterize the rock mass; • time-lapse sequences of IRT thermal images in order to detect thermal transient within both the fault zone and the low fractured mass rock sectors; • thermal conductivity measurements to characterize thermal properties of rocks.
In detail, efforts were focused on comparing IRT and TC data with geo-mechanical parameters acquired on two scanlines following a standard ISRM procedure. In addition, a photogrammetric modelling of the area was conducted, and its derivative products such the orthomosaic and the texturized 3D model were used to double check the geo-mechanical features.

Study Site
The quarry is located close to the village of S. Lorenzo di Valdieri, nearby Cuneo village (North-Western Alps, Piedmont, Italy, Figures 1 and 2). The geological units outcropping in the area belong to the Middle Jurassic-Lower Cretaceous Provençal succession located at the North-Eastern side of the Argentera Massif. The quarry produces pure calcium carbonate (99.9% volume of Calcite, Carbocalcio Cuneese S.p.A., and personal communication) treated and selected for sale as granulated and micronized.
The carbonate rock mass (Lausa Limestone) consists of fine-grained limestone, with abundant decimeter-thick beds of polymictic breccias, generally clast-supported, and with millimeter to decimeter sized clasts of mudstones, coarsely crystalline dolostones and finely crystalline dolostones. Lausa Limestones are followed by grey mudstones and crinoid-rich wackestones, in centimeter to decimeter thick beds, with abundant silicified portions. Fault rocks (Carnieules) and Entracque Marl outcrops toward SW while, in the uppermost part of the quarry, Aptian-upper Cretaceous rocks (Puriac Limestones and Marne Nere fm.) are present ( Figure 2). Carbonates are locally affected by a diffuse hydrothermal dolomitization occurring in the Early Cretaceous period at a very shallow burial depth, which is related to the expulsion of hot fluids (about 200 • C) through faults and fractures during episodes of fault activity. The studied outcrop can be considered as belonging to the Valdieri Marbles according to Reference [33].
The quarry is at an elevation roughly between 900 and 1250 m. The front has an envelope slope of approximatively 53 • and the dip direction is southeast (about 125 • ). It is structured in benches (10 m high, 2-2.5 m wide). The dip of the single terrace face is approximatively 62 • -65 • with local situations of greater acclivity. The front perpendicular to the cultivation has an envelope slope of approximatively 45 • and is set on the schistosity planes. The dip direction of the front envelope is towards the South-West (about 200 • ). The dip of the single step elevation is that of the schistosity and is approximately equal to 50 • but is locally higher.

Introduction
Groundwater flow in fractured rock masses has been studied by many earth scientists. It reveals one of the most important factors to be considered in slope stability problems [1][2][3] as well as in the hydrogeological assessment and in heat transfer mechanisms [4,5]. In particular, understanding geo-

Geo-Mechanical Surveys
The outcrop was investigated by two geo-mechanical scanlines (SL1 and SL2), 73 and 20 m long, respectively (Figure 3a), according to the ISRM procedures [7]. Data were elaborated aiming at highlighting joint density (number of joints per meter, [J/m]) across the fault zone.

Geo-Mechanical Surveys
The outcrop was investigated by two geo-mechanical scanlines (SL1 and SL2), 73 and 20 m long, respectively (Figure 3a), according to the ISRM procedures [7]. Data were elaborated aiming at highlighting joint density (number of joints per meter, [J/m]) across the fault zone. Remote Sens. 2019, 11, x FOR PEER REVIEW 3 of 15 In Figure 3b, the 3D reproduction of the test site is shown. The photogrammetric survey consists of a series of overlapping images acquired with a Canon full frame DSLR coupled with a 50 mm lens. Ground control points for model orientation were dislocated over the quarry front and were acquired by using a rugged GPS tablet (Leica CS25 plus) able to reach a positioning error below 2 cm in the real time RTK mode. MicMac, an open source software, was used to perform tie point extraction, bundle block adjustment, dense image matching, and 3D reconstruction. The obtained cloud consisting of 7.7 million points has been interrogated by means of CloudCompare software for checking structural data using the Compass plugin [34]. A high resolution orthophoto ( Figure 3a) with a Ground Sample Distance (GSD) of 1 cm/pixel was produced in order to measure joint spacing in non-reachable portions of the outcrop.

Principles
Thermography is based on the physical principle stating that, any object having temperature higher than absolute zero, produces electromagnetic radiation. The electromagnetic spectrum is divided arbitrarily into several wavelength regions, called bands, which are distinguished by the methods used to produce and detect the radiation. Thermal radiation occurs in the portion of the spectrum between 0.1 and 100 μm, including parts of the ultraviolet along with visible and infrared regions. However, most of the thermal radiation falls within the IR band [15,35,36]. Thermal cameras for a remote sensing purpose are usually designed and calibrated for the Long-Wavelength IR band (8-14 μm) due to the peak energy emissions for ambient Earth surface temperatures and to the In Figure 3b, the 3D reproduction of the test site is shown. The photogrammetric survey consists of a series of overlapping images acquired with a Canon full frame DSLR coupled with a 50 mm lens. Ground control points for model orientation were dislocated over the quarry front and were acquired by using a rugged GPS tablet (Leica CS25 plus) able to reach a positioning error below 2 cm in the real time RTK mode. MicMac, an open source software, was used to perform tie point extraction, bundle block adjustment, dense image matching, and 3D reconstruction. The obtained cloud consisting of 7.7 million points has been interrogated by means of CloudCompare software for checking structural data using the Compass plugin [34]. A high resolution orthophoto ( Figure 3a) with a Ground Sample Distance (GSD) of 1 cm/pixel was produced in order to measure joint spacing in non-reachable portions of the outcrop.

Principles
Thermography is based on the physical principle stating that, any object having temperature higher than absolute zero, produces electromagnetic radiation. The electromagnetic spectrum is divided arbitrarily into several wavelength regions, called bands, which are distinguished by the methods used to produce and detect the radiation. Thermal radiation occurs in the portion of the spectrum between 0.1 and 100 µm, including parts of the ultraviolet along with visible and infrared regions. However, most of the thermal radiation falls within the IR band [15,35,36]. Thermal cameras for a remote sensing purpose are usually designed and calibrated for the Long-Wavelength IR band Remote Sens. 2019, 11, 179 6 of 23 (8-14 µm) due to the peak energy emissions for ambient Earth surface temperatures and to the presence of atmospheric windows ( [20]). The relation between the total emissive power and an object's temperature is stated by the Stefan-Bolzmann's Law [37,38].
where W, is the total radiant emittance from the body, expressed in Watt m −2 , ε is the object emissivity, σ, is the Stephan Boltzmann constant = 5.6697 × 10 −8 Wm −2 K −4 , and T is the absolute temperature (expressed in K) of the emitting material. When a surface is sensed by a thermal camera, the thermal radiation pattern acquired depends not only on the surface temperature but also on the emissivity of the material, on the atmospheric absorption, as well as on the reflection observed on the radiation surface, due to the surrounding bodies [10,39]. These variables must be taken into account at the time of measurement, and their contribution must be known in order to estimate the real temperature of the object.

Surveys Methodology
Two measurement sessions were carried out in 2018 including the first one on January 23 (winter session) and the second one on June 18 (summer session). The aim was to evaluate the behavior of the rock mass in different seasons, insulation conditions, and water content within the fault zone. In every session, a series of thermal images were acquired at an hourly rate. The instrument employed was a FLIR E8 thermal camera (FLIR Systems Inc, [40]) equipped with an uncooled FPA micro bolometer sensor, with a resolution of 320 × 240 pixel, IFOV = 2.6 mrad. It can perform measurements in the thermal infrared band (between 7.5-13 µm). The device was placed on a stable tripod 80 m away from the target, giving a GSD, Ground Sample Distance, of approximatively 0,18 m/pix. In addition to the thermal image, this device produces an optical image, which helps identify the measured area and interpret more easily the obtained thermal data. Thermal data and frames were imported and calibrated by using Flir Tools software (FLIR Systems Inc, version 5.13.17214.2001). After that, the data matrix and optical images were exported to be aligned and to calculate differentials by using the ArcGis Raster Calculator (ESRI ArcMap 10.2.2).
Images were labelled using the following convention: the letter F means the Frame, followed by two numbers identifying the hour when the image was taken, omitting minutes (e.g., F16 refers to an image taken at 16.00 p.m.). During both measuring sessions, the air temperature and humidity were recorded. The first session, conducted from the late morning, allowed for detecting both the heating and the onset of the cooling phase of the rock mass. The air temperature varied between −1 and 11 • C. The second one recorded only the cooling phase of the rock wall, avoiding the contribution of the solar radiation. The measurement session lasted from 19:00 to 23:00, while the air temperature ranged between 24 and 14 • C. In this session, an appreciable run-off flow coming from the top of the quarry front was detected on the Bench 1, steaming from East to West (Left to right in Figures 1 and 3). The temperature of the water flowing was 16.0 • C. The environmental conditions are shown in Table 1. The emissivity of the rock, the reflected temperature, the shooting length, and the relative humidity have been set appropriately for each shot (except for the winter session, when humidity and reflected temperature were not recorded). In particular, emissivity has been detected by means of lab tests on a rock sample heated in the hoven up to 40 • C. A piece of black electrical tape with a known emissivity of 0.95 was used as a reference, following the "reference emissivity material method" as given in Reference [41]. By assuming emissivity equal to one and by placing a crumpled and re-flattened piece of aluminum foil in front of the camera, the reflected temperature was measured by using the reflector method [42,43]. In particular, the reflector was positioned trying to reproduce the same slope angle and filling the frame. Although outdoor conditions could affect this measure, we consider this approach the one optimizing time saving and accuracy.

Thermal Conductivity Analysis
Thermal conductivity defines the ability of a material to transfer heat corresponding to the amount of heat transferred from a body per unit area, as it is expressed by the Fourier's Law of thermal conduction defined using the formula below.
where q is the heat flow in the direction of X, and T is the temperature. The coefficient of proportionality λ is the thermal conductivity [43]. A sequence of field and laboratory measures adopting two different methods such as Transient Line Source (TLS, [44]) and Transient Divide Bar (TDB, [43]), was, therefore, carried out on a meaningful number of samples. This is important in order to understand the influence and the accuracy of some experimental conditions, as well as to compare obtained results with known mixing models for measuring bulk thermal conductivity, as proposed in the literature [45].
Aiming at obtaining precise and reliable values of thermal conductivity, density, and porosity parameters were also measured in intact rock specimens. Because of the irregular shapes of the collected samples, bulk density was obtained by water saturation and by weighing all the collected rocks, which refers to Archimede's Principles [23]. Each specimen was then oven-dried at 70 • C for one day in order to obtain the dry mass and to infer porosity [46].

Transient Line Source (TLS)
K2DPro, which is a thermal property analyzer developed by Decagon Devices, compliant to standards ASTM D5334 and IEEE 442 [47,48] was used. It relies on the Transient Line Source (TLS) method, assuming a heat production at a constant rate by a heated wire and predicting an infinitely long and small source inside an infinite and homogeneous medium.
The obtained temperature, can be found in (Equation (3)).
where q (W m −1 ), is the specific rate at which the heat is generated, T 1 and T 2 ( • C) are the recorded temperatures at time steps t 1 and t 2 (s), respectively, and λ is the thermal conductivity. The device uses a specific algorithm based on References [49,50] line heat source analysis. This permits measurements during a heating and a cooling interval, and the separation of the effects of the heat pulse from ambient temperature changes. The algorithm fits time and temperature data using an exponential integral function through a non-linear least squares method. Based on this theory, the K2Dpro device allows for measures with an accuracy of ±10% [47]. It consists of four specific sensors designed for well-defined samples. In this study, the most proper single needle probe (RK-1 sensor, 0 cm long, 3.9 mm diameter) was directly inserted into the investigated materials. With regard to field measures, since thermal conductivity gives a punctual value, two "end members" were tested: sound rock and cataclastic zone. The needle probe was easily introduced in the abundant fine-grained matrix of the cataclastic zone, ensuring a good coupling. Regarding the sound rock, a 4 mm hole was drilled on the surface, and a low quantity of thermal grease was introduced to improve thermal contact between the needle and the rock. Since thermal conductivity is strictly linked to physical, mechanical, and mineralogical characteristics of the material [51,52], laboratory measures both in dry and in wet conditions were conducted too. Due to the difficulty to collect undisturbed samples of the cataclastic fault zone, only sound rock was tested in the laboratory. Each final value achieved is an average over 10 measurements.

Transient Divided Bar (TDB)
The Transient Divided Bar (TDB) is an implemented apparatus designed by Reference [44]. It consists of a stack of two copper blocks acting as a heat sink (the upper one) and a heat source (the lower one), respectively, and, in between the cylindrical rock specimen, is interposed. It is based on Reference [53], adapted for defining thermal conductivity of semiconductors. It predicts that thermal conductivity is obtained through temperature changes between the heat sink and the heat source during the test as well. Temperature changes are recorded and defined by the amount of heat removed from the upper copper block C u ∆T u in a specific period of time ∆T, according to the Fourier's postulate. This is given by (Equation (4)).
is the cross-sectional area of the rock sample respectively, T is the temperature, and u and l are suffixes that refer to the upper and lower block. Details of the TDB device as well as the procedure followed to conduct tests, are described in Reference [43]. In particular, TDB is designed for laboratory analysis on sound rocks. For this study, it was applied as a comparative methodology with the TLS apparatus, in order to obtain more reliable thermal conductivity data.
Cylindrical specimens 25 ± 0.1 mm in diameter and 20 ± 0.5 mm in thickness, and with parallel, flat (within 0.1 mm) and smooth (within 0.03 mm) surfaces, were prepared. A film of silicone paste (0.1 mm) was applied to improve thermal contact between the rock sample and copper blocks of the used device. Similar to TLS, each final value obtained is averaged over 10 measurements.

Geo-Mechanical Surveys
Bench 1 and 2 are characterized by four main intersecting discontinuity systems (Figure 3a) wich are the schistosity (188/63) and three persistent joint systems, K1, K2, and K3 (respectively, 170/50, 208/79 and 022/81). Fractures are generally closed or with aperture <2 mm while the filling, present as an abundant fraction only in the fault zone, is mostly soft and made of sandy/silty weathered material. No water flows were reported in the winter survey session (probably due to freezing), while a certain humidity degree was present in the summer session. Surface alteration increases approaching to the fault zone. The main geo-structural parameters are summarized in Table 2 and Figure 4. The number of fractures per linear meter (J/m) provides a good estimate for the potential water circulating through the rock mass. This indicator is consistent with more standardized indexes, such as RQD [54], GSI [55], and RMR [56], as can be seen in Figure 5.
With respect to the Hoek & Brown classification, the GSI value was inferred according to rock mass interlocking and surface conditions, while UCS and m i values were chosen based on an empiric field observation, according to authors' suggestions [55].
Considering arbitrary joint density thresholds, four lithotechnical units are identified and labelled RT1, RT2, RT3, and RT4, respectively. Quality (latu sensu) of the rock mass decreases from RT1 to RT4 (from Blocky, to Very blocky, to Disturbed to Disintegrated) since the number of joints per meter increases, and RQD, GSI, and RMR decrease. Figure 5 shows the location of the four litho-technical units on the two surveyed benches. In Figure 6, some geo-mechanical parameters [55] are calculated for every lithotechnical unit, by using the RocData software (Version 5.009, RocScience Inc.) (Figure 7). fault zone. In every session, a series of thermal images were acquired at an hourly rate. The instrument employed was a FLIR E8 thermal camera (FLIR Systems Inc, [40]) equipped with an uncooled FPA micro bolometer sensor, with a resolution of 320 × 240 pixel, IFOV = 2.6 mrad. It can perform measurements in the thermal infrared band (between 7.5-13 μm). The device was placed on a stable tripod 80 m away from the target, giving a GSD, Ground Sample Distance, of approximatively 0,18 m/pix. In addition to the thermal image, this device produces an optical image, which helps identify the measured area and interpret more easily the obtained thermal data. Thermal data and frames were imported and calibrated by using Flir Tools software (FLIR Systems Inc, version 5.13.17214.2001). After that, the data matrix and optical images were exported to be aligned and to calculate differentials by using the ArcGis Raster Calculator (ESRI ArcMap 10.2.2).
Images were labelled using the following convention: the letter F means the Frame, followed by two numbers identifying the hour when the image was taken, omitting minutes (e.g., F16 refers to an image taken at 16.00 p.m.). During both measuring sessions, the air temperature and humidity were recorded. The first session, conducted from the late morning, allowed for detecting both the heating and the onset of the cooling phase of the rock mass. The air temperature varied between −1 and 11 °C. The second one recorded only the cooling phase of the rock wall, avoiding the contribution of the solar radiation. The measurement session lasted from 19:00 to 23:00, while the air temperature ranged between 24 and 14 °C. In this session, an appreciable run-off flow coming from the top of the quarry front was 3.3. Thermal

IRT Winter Field Acquisition
The winter measurement session was conducted from the late morning to the late afternoon on Bench 1. Figure 8 shows raw images and differences between every image and the first one. The solar radiation affects heavily the temperature pattern over the rock outcrop: a quite homogeneous increasing in temperature is shown, with some difference due to the orientation of the surfaces with respect to the position of the sun. For example, the darker loose deposit on the benches (C) has the higher surface temperatures at the start of the survey, and the lower ones at the end. They show a

IRT Winter Field Acquisition
The winter measurement session was conducted from the late morning to the late afternoon on Bench 1. Figure 8 shows raw images and differences between every image and the first one. The solar radiation affects heavily the temperature pattern over the rock outcrop: a quite homogeneous increasing in temperature is shown, with some difference due to the orientation of the surfaces with respect to the position of the sun. For example, the darker loose deposit on the benches (C) has the higher surface temperatures at the start of the survey, and the lower ones at the end. They show a lower thermal inertia than the rocks outcropping in the quarry front. Moreover, it is evident that shaded sectors tend to heat less than the surroundings (cfr. Point S in Figure 8). This is consistent with literature data, which highlights the massive disturbance induced by sun heating [15,17,39]. Only the cooling down phase seems to highlight some peculiar pattern of the fault zone. However, images acquired after the sunset are not enough for providing data necessary for retrieving a cooling rate.

IRT Summer Field Acquisition
The summer measurement session was conducted entirely after the sun had set, with a decreasing air temperature. An appreciable flow of water (T ~ 16 °C) was present on the top the bench (not visible from the IR camera position due to the presence of the loose deposit), and a certain degree of humidity was visible in the filling of the fault zone fractures.
Looking at the raw images ( Figure 9) from a to e (from 19:00 to 23:00), it can be noted that RT1 and RT2 portions are systematic with a higher surface temperature than RT3, RT4, and C portions (loose deposits with vegetation). The temperature trend over time shows an evident decrease, in every sector of the frame. Making some spot measurements (Sp1 to Sp6), it is possible to quantify the cooling rate of the single measurement points along a meaningful profile. The measured points are chosen along a line and include portions of sound rock, wet sectors, and fault rocks. Although a punctual measure may be affected by errors due to tripod instability and camera vibrations, this does not undermine measurement representativeness. Actually, images were aligned before measuring. Furthermore, variations in terms of the temperature of the adjacent pixel are not so significant. The

IRT Summer Field Acquisition
The summer measurement session was conducted entirely after the sun had set, with a decreasing air temperature. An appreciable flow of water (T~16 • C) was present on the top the bench (not visible from the IR camera position due to the presence of the loose deposit), and a certain degree of humidity was visible in the filling of the fault zone fractures.
Looking at the raw images ( Figure 9) from a to e (from 19:00 to 23:00), it can be noted that RT1 and RT2 portions are systematic with a higher surface temperature than RT3, RT4, and C portions (loose deposits with vegetation). The temperature trend over time shows an evident decrease, in every sector of the frame. Making some spot measurements (Sp1 to Sp6), it is possible to quantify the cooling rate of the single measurement points along a meaningful profile. The measured points are chosen along a line and include portions of sound rock, wet sectors, and fault rocks. Although a punctual measure may be affected by errors due to tripod instability and camera vibrations, this does not undermine measurement representativeness. Actually, images were aligned before measuring. Furthermore, variations in terms of the temperature of the adjacent pixel are not so significant. The error committed could be regarded as negligible.   Full lines represent the temperature measured at every time step, while the dashed line correspond to the measured air temperature over time. Wet portions (Sp3 and Sp6, detected during the field survey) are characterized by a lower surface temperature in a, and show a lower average cooling rate, with respect to drier ones (Sp1, Sp2, Sp4, and Sp5). Moreover, all the curves tend to have a lower steepness in the right part of the graph, showing a decrease in the cooling rate over time, as illustrated in Figures 9 and 10.  Figure 11 shows thermal differentials between consecutive images ΔTH. In general, hourly variations correspond to ΔTH between −2.5 °C and +0.1 °C. The hourly cooling rate is not constant over time, but generally decreases with a fluctuating trend. This fact could be due to equipment sensitivity, environmental disturbance, and position of the spot measures. In fact, the environmental temperature constantly decreases over time. Wet and dry areas do not show peculiar behavior with respect to this parameter.  Figure 10 shows thermal differentials of every image compared to the first one, ∆T C . Portions of the bench with different thermal behavior are emphasized. In general, ∆T C between 0 • C and −8 • C are shown over the whole frame. Among the six measured spots, the wet portions (Sp3 and Sp6) have a lower final thermal differential (in e → ∆T C Sp3~− 4.0 • C and ∆T C Sp6~− 3.5 • C) while the dry rock areas show more consistent thermal losses (e → ∆T C Sp2~− 6 • C). Figure 11 shows thermal differentials between consecutive images ∆T H . In general, hourly variations correspond to ∆T H between −2.5 • C and +0.1 • C. The hourly cooling rate is not constant over time, but generally decreases with a fluctuating trend. This fact could be due to equipment sensitivity, environmental disturbance, and position of the spot measures. In fact, the environmental temperature constantly decreases over time. Wet and dry areas do not show peculiar behavior with respect to this parameter.

Thermal Properties
The main results of thermal conductivity, porosity, density, and temperature analysis are summarized in Table 3. The field five measurements on the main representative portion of the sound rock mass (EM1) as well as on the cataclastic zone (EM2) were conducted ( Figure 12). Laboratory tests in dry and wet conditions (Table 3) were also performed. In both cases, each test represents an average of five values over 10 measures.

Thermal Properties
The main results of thermal conductivity, porosity, density, and temperature analysis are summarized in Table 3. The field five measurements on the main representative portion of the sound rock mass (EM1) as well as on the cataclastic zone (EM2) were conducted ( Figure 12). Laboratory tests in dry and wet conditions (Table 3) were also performed. In both cases, each test represents an average of five values over 10 measures. Table 3. Results from field and laboratory activities. Keys: TLS (Transient Line Source methods) and TDB (Transient Divided Bar technique), EM2 and EM1, fault and rock respectively. The symbol percentage indicates the percentage between two measurements of thermal conductivity. T is the temperature ( • C), ϕ indicates the porosity (% vol), and defines density (g/cm 3 ). λb, is the bulk of thermal conductivity (W m −1 K −1 ). St.dev., the standard deviation of λb.      Figure 13).
Taking into account only the sound rock, if we compare averages values of thermal conductivity for field and laboratory data, the underestimation of TLS with respect to TDB decreases if the TLS technique is used in the field. Moreover, a clear difference of the EM2 with respect to EM1, can be further highlighted. It behaves differently with both TDB and TLS rock data. Moreover, very low standard deviations (0.021 to 0.110) can be observed (Table 3).
Looking at the more representative portions of the sound rock mass, relationships of wet thermal conductivity with porosity (η, vol %) and density (ρ, g/cm 3 ) are summarized in Figure 14. Figure 13. Field data relative to thermal conductivity (above) and temperature (below) both for cataclastic fault zone (Orange) and low fractured rock mass (Blue). Estimated errors in the measurement are represented by vertical error bars.  (Table 3, Reference [31]).
Taking into account only the sound rock, if we compare averages values of thermal conductivity for field and laboratory data, the underestimation of TLS with respect to TDB decreases if the TLS technique is used in the field. Moreover, a clear difference of the EM2 with respect to EM1, can be further highlighted. It behaves differently with both TDB and TLS rock data. Moreover, very low standard deviations (0.021 to 0.110) can be observed (Table 3).
Looking at the more representative portions of the sound rock mass, relationships of wet thermal conductivity with porosity (η, vol %) and density ( , g/cm 3 ) are summarized in Figure 14.

Discussion
As a general point of view, good agreement between geo-mechanical characteristics of rock masses and their thermophysical behavior is highlighted. Thermal imaging results obtained during rock mass heating and cooling phases, differ in significance depending on several reasons. First of all, the heating trend depends more on the sun position than on the rock mass characteristics. In fact, the relatively uniform and constant slope aspect and surface roughness do not affect substantially thermal measurements.
Then, a reflected and scattered solar radiance play a major role during the day, and considerably affect the rock mass emission term. Thermal contrast between rock mass and air temperature is higher after the sunset. This highlights the rock mass behavior (Figure 8). This is consistent with literature Density ranges from 2.54 to 2.74 g/cm 3 ( Figure 14A, Table 3), which are quite high values for a carbonatic rock. Since the negligible porosity (0.30 to 0.48 vol %, Figure 14B, Table 3) and a lack of interconnection between factures, water content does not affect measurements.

Discussion
As a general point of view, good agreement between geo-mechanical characteristics of rock masses and their thermophysical behavior is highlighted. Thermal imaging results obtained during rock mass heating and cooling phases, differ in significance depending on several reasons. First of all, the heating trend depends more on the sun position than on the rock mass characteristics. In fact, the relatively uniform and constant slope aspect and surface roughness do not affect substantially thermal measurements.
Then, a reflected and scattered solar radiance play a major role during the day, and considerably affect the rock mass emission term. Thermal contrast between rock mass and air temperature is higher after the sunset. This highlights the rock mass behavior (Figure 8). This is consistent with literature data [15,17]. For those reasons, the results obtained under daylight conditions are affected by a not insignificant source of error.
The summer measurement session, aimed at recording the cooling phase, highlights the different thermal response of the fault zone (RT3 and RT4) with respect to the low fractured sectors (RT1 and RT2). In fact, raw images show that the temperature of the fault zone, characterized by a higher degree of fracturing and low GSI, RQD, and RMR, is sensibly lower in every time step. This outcome is dissimilar from what is obtained in Reference [15]. In fact, authors report higher surface temperature for a high fractured portion of the rock mass. The difference can be ascribed to the hydraulic conditions of the discontinuities. In the case examined herein, there is a consistent water circulation, while, in the reference, the rock mass is in dry conditions. Figure 15 shows the comparison between the number of joint per meter and the measured temperature for the F19 frame on the SL1 scanline ( Figure 9a). It is evident that an inverse relationship among these two parameters exist. This is due to the higher water content of RT3 and RT4 units. In fact, while the higher temperature is registered for the sound rock (~32 • C), the water that circulates in fractures at 16 • C cool down the fault zone up to 4 to 5 degrees (F19, Figure 9a). The inverse solution is also possible to apply toward detecting cooler sectors. Thus, the information on the degree of fracturing is possible to retrieve (i.e., SP6, Figure 9). data [15,17]. For those reasons, the results obtained under daylight conditions are affected by a not insignificant source of error. The summer measurement session, aimed at recording the cooling phase, highlights the different thermal response of the fault zone (RT3 and RT4) with respect to the low fractured sectors (RT1 and RT2). In fact, raw images show that the temperature of the fault zone, characterized by a higher degree of fracturing and low GSI, RQD, and RMR, is sensibly lower in every time step. This outcome is dissimilar from what is obtained in Reference [15]. In fact, authors report higher surface temperature for a high fractured portion of the rock mass. The difference can be ascribed to the hydraulic conditions of the discontinuities. In the case examined herein, there is a consistent water circulation, while, in the reference, the rock mass is in dry conditions. Figure 15 shows the comparison between the number of joint per meter and the measured temperature for the F19 frame on the SL1 scanline ( Figure 9a). It is evident that an inverse relationship among these two parameters exist. This is due to the higher water content of RT3 and RT4 units. In fact, while the higher temperature is registered for the sound rock (~32 °C), the water that circulates in fractures at 16 °C cool down the fault zone up to 4 to 5 degrees (F19, Figure 9a). The inverse solution is also possible to apply toward detecting cooler sectors. Thus, the information on the degree of fracturing is possible to retrieve (i.e., SP6, Figure 9). Hot and dry areas-and, thus, less fractured one-are characterized by a higher cooling rate, while the cold and wet ones (with poor geotechnical characteristics) shows lower temperature differentials. This is in agreement with the Newton's law of Cooling, stating that the rate of heat loss of a body is directly proportional to the difference in the temperatures between the body and its surroundings.
The sound rock seems to conduct far more than cataclastic fault (in dry conditions), due to the better thermal conductivity values both from field and laboratory surveys (Table 3). This could be due to many factors ranging from different physical, mechanical, and mineralogical characteristics. In particular, lower thermal conductivities of the cataclastic fault zone could be due to a higher content of clay minerals in its wet thermal conductivity can be observed. However, this underestimation abruptly decreases when comparing field data on  3  6  9  12  15  18  21  24  27  30  33  36  39  42  45  48  51  54  57  60  63  66  69  Hot and dry areas-and, thus, less fractured one-are characterized by a higher cooling rate, while the cold and wet ones (with poor geotechnical characteristics) shows lower temperature differentials. This is in agreement with the Newton's law of Cooling, stating that the rate of heat loss of a body is directly proportional to the difference in the temperatures between the body and its surroundings.
The sound rock seems to conduct far more than cataclastic fault (in dry conditions), due to the better thermal conductivity values both from field and laboratory surveys (Table 3). This could be due to many factors ranging from different physical, mechanical, and mineralogical characteristics. In particular, lower thermal conductivities of the cataclastic fault zone could be due to a higher content of clay minerals in its matrix or to a higher porosity in dry conditions. Laboratory experiments were focused on the characterization of the sound rock including comparative and multiple measurements, which allow the obtainment of more reliable values of thermal conductivity through different methodologies and working conditions. In this case, a constant underestimation of TLS data than TDB (Transient Divided Bar technique) for both dry and wet thermal conductivity can be observed. However, this underestimation abruptly decreases when comparing field data on the same rock type. In this case, thermal conductivity using TLS in the field results in higher values, approaching those from the TDB device during laboratory activity. A negligible λ w /λ d ratio (ratio between wet and dry thermal conductivity) retrieved in comparative analysis of bulk thermal conductivity of sound rock confirm that porosity does not play an important role in the sampled outcrop. As reported in the literature [57][58][59], fractures can highly increase the λ w /λ d ratio, but only when fractures are interconnected and affected by groundwater circulation. On the other hand, some different factors (e.g., air or re-crystallized calcite aggregates or some other low conductive materials fulfilled fractures [60]) can contribute to lowering the λ w /λ d ratio even if fractures are abundant in the rock mass. In spite of the very low fracture porosity, the negligible λ w /λ d ratio can be ascribed to the mono-mineralic content (around 99.9% volume of Calcite, Carbocalcio Cuneese S.p.A., personal communication), as well. The laboratory tests demonstrate, once more, as the thermal conductivity is strictly connected to the porous medium, the structure and the mineralogy of the rock [61]. Reliability of the obtained results is supported thanks to low standard deviations (Table 3) as well.

Conclusions
Both mesoscopic and microscopic characteristics influence thermal properties of the rock masses. The thermal response of an outcrop in heating and in cooling phases, which was confirmed in this study, is caused by different and simultaneous heat transfer processes. At a sample scale, thermal behavior is led by the microscopic features (mineralogical composition, microstructure, interstitial fluid characteristic, etc.) and by the heat conduction process. Thermal conductivity values are systematically higher in the sound rock with respect to the fault zone. A negligible λ w /λ d ratio shows that, in his case, fracture porosity is not relevant in influencing thermal properties of the sound rock. Heat conduction is less efficient at the mesoscopic and macroscopic scale, where convection and advection processes due to water circulation into the fracture network are predominant. Thus, the thermal imaging, which highlights a portion of the outcrop characterized by a higher water content, can be a useful tool for detecting a densely fractured litho-technical unit. This can be used when a direct investigation of rock masses (or part of them) is not possible and can be successful in defining portions characterized by similar mechanical and hydro-geological behavior.