Relationship between Rock Porosity and Infrared Cooling Rate in Non-Standard Specimens of Tuffs Used in the Hungarian Cultural Heritage

: This paper is focused on the application of Infrared Thermography to non-standard rock specimens, in terms of size and deterioration conditions, of Hungarian tuff to monitor their cooling process and to look for a relationship between the rock Cooling Rate Index and the porosity. Literature data agree on the potential of Infrared Thermography for the indirect estimation of rock porosity in fresh specimens through the IRTest, but this technique has never been tested on non-standard specimens. To this purpose, tests on three varieties of Hungarian tuffs were carried out. These materials were selected for their cultural importance linked to their usage as building stones and in other historical applications in Northern Hungary. Tuff specimens underwent a ﬁxed number of salt crystallization cycles. The Cooling Rate Index (CRI) for each specimen was calculated according to the literature experience and correlated to their porosity estimated by water, helium, and mercury intrusion. The results show that the rock cooling process is related to porosity since more porous rocks are characterized by faster cooling. Positive linear trends were achieved for weathered specimens considering 20 min monitoring (CRI 20 ), which is double the time suitable for untreated rocks. The reason should be searched in salt crystallization’s effects on the rock texture, paving the way to further studies on this pioneering branch of technological application.


Introduction
Porosity is a physical property conditioning rock mechanical behavior under stress in terms of strength and deformability [1][2][3][4][5]. The presence of a void network within the rock represents a weakening feature, reducing the contact between grains, lowering rock resistance against axial stresses, and enhancing the rock attitude towards deformation. The solid literature experience in this subject points out that Uniaxial Compressive Strength (UCS) and Young's modulus are inversely proportional to porosity, with statistical trends scattered depending also on the type of voids in the rock, i.e., pores or cracks or both, e.g., [6][7][8]. Porosity is recognized to be also related to rock durability from the point of view of weathering. In fact, the growth of salt and ice crystals in the pore system and their related stresses can be acknowledged among the major deterioration mechanisms of porous materials [9][10][11][12][13][14][15]. This physical property can be referred to as "effective" and "total" depending on whether only the volume of interconnected voids or the total void volume, respectively, is considered. Porosity estimation is therefore relevant for predicting the behavior of rocks used as construction or building materials and specific engineering geological aspects. The International Society for Rock Mechanics (ISRM) and the European Committee for Standardization Comité Européen de Normalisation (CEN) defined Demjén tuff is the most porous and soft, with compressive and tensile strengths of 10 and 1 MPa, respectively. Bogács and Sirok tuffs have better physical and mechanical properties, with compressive strength and tensile strength in the ranges of 26-28 MPa and [3][4] MPa. For all three tuff varieties, water saturation may cause an extreme deterioration of the mechanical performance [23].

Methods
The methodological approach set for this study can be summarized by 4 main points: (1) salt crystallization tests; (2) IRTest for the calculation of the Cooling Rate Indexes; (3) rock porosity estimation under water, helium, and mercury intrusion; (4) statistical data processing.

Salt Crystallization Tests
Salt crystallization tests were performed to reproduce the rock degradation effects naturally arising from the weathering phenomena involving in-pore water circulation. Tuff specimens had a non-standard shape, being cylinders with a 50 mm diameter and a 30 mm height. A total of 72 specimens were prepared. These were then grouped into 6 sets, each undergoing a different number of crystallization cycles (0, 2, 4, 6, 10, 23), with 1 set serving as an untreated reference (Cycle 0). The different steps were chosen based on how quickly the observed state of conservation changed, so as to isolate representative samples for different decay stages for the following IRTests. According to the standard test method EN 12370 (1999), each salt crystallization cycle consisted of the following phases: 2 h total immersion in a 14% Na2SO4 aqueous solution, minimum 16 h drying at 105 °C, and 2 h cooling. After each cycle, the sample mass and macroscopic changes were recorded.

Infrared Cooling Monitoring and CRI Calculation
The IRT technique relies on the principle that every material at a temperature above absolute zero emits thermal radiation. This is characterized by a wavelength mainly falling within the infrared band of the electromagnetic spectrum [26], and it is proportional to the temperature of the emitting body, according to the Stefan-Boltzmann law. Therefore, the surface temperature of a body can be easily estimated by a thermal camera, which

Methods
The methodological approach set for this study can be summarized by 4 main points: (1) salt crystallization tests; (2) IRTest for the calculation of the Cooling Rate Indexes; (3) rock porosity estimation under water, helium, and mercury intrusion; (4) statistical data processing.

Salt Crystallization Tests
Salt crystallization tests were performed to reproduce the rock degradation effects naturally arising from the weathering phenomena involving in-pore water circulation. Tuff specimens had a non-standard shape, being cylinders with a 50 mm diameter and a 30 mm height. A total of 72 specimens were prepared. These were then grouped into 6 sets, each undergoing a different number of crystallization cycles (0, 2, 4, 6, 10, 23), with 1 set serving as an untreated reference (Cycle 0). The different steps were chosen based on how quickly the observed state of conservation changed, so as to isolate representative samples for different decay stages for the following IRTests. According to the standard test method EN 12370 (1999), each salt crystallization cycle consisted of the following phases: 2 h total immersion in a 14% Na 2 SO 4 aqueous solution, minimum 16 h drying at 105 • C, and 2 h cooling. After each cycle, the sample mass and macroscopic changes were recorded.

Infrared Cooling Monitoring and CRI Calculation
The IRT technique relies on the principle that every material at a temperature above absolute zero emits thermal radiation. This is characterized by a wavelength mainly falling within the infrared band of the electromagnetic spectrum [26], and it is proportional to the temperature of the emitting body, according to the Stefan-Boltzmann law. Therefore, the surface temperature of a body can be easily estimated by a thermal camera, which is a device operating within the range of wavelengths as long as 13-14 µm. The use of IRT for rock characterization has increasingly become the subject of numerous scientific studies focused either on rock deterioration and fracturing state [11,[27][28][29][30] or on rocky outcrops [31][32][33][34][35][36].
Recent literature studies shed light on the positive correlation between rock cooling, monitored by IRT, and porosity [19][20][21][22]. Starting from the hypothesis that the presence of voids within the rock texture could affect its thermal behavior, a series of tests on numerous intact rock specimens were performed. These proved that the greater the porosity of the rock, the faster its cooling is.
The IRTest on intact rock specimens is carried out according to the following four main stages: (1) Oven-heating of the rock specimens for 24 h at 70 ± 5 • C.
(2) Monitoring of the rock cooling by acquiring regular IRT images of the upper face of the specimens, which are placed on a smooth laminate worktop in a temperature-controlled room (25 ± 1 • C), for a defined time interval. The initial IRT frame is acquired when the cooling begins, thus representing the reference for the following stage.
(3) Post-processing of the IRT images in order to define the average surface temperature of the framed specimen face; this is the average temperature value calculated within a specific area of the IR image, set to coincide with the specimen's upper face.
(4) According to Equation (1), the Cooling Rate Index (CRI) of the rock is calculated for a defined monitoring time window. CRI is expressed by the ratio between the difference of temperature between two measurements (∆T) within a specific time window (∆t).

CRI = ∆T/∆t
(1) Literature outcomes agree that the best relationship between rock porosity (estimated by buoyancy in water) and rock cooling is achieved by calculating the CRI within a 10 min time window starting from when the cooling begins (CRI 10 ). A satisfactory correlation was found for rocks belonging to different geological formations and for different shapes and sizes of specimens [19][20][21][22].
In this paper, the IRTest was experimentally carried out, on the same specimens prepared for salt crystallization (see Section 3.1), by using a high-sensitivity infrared thermal imaging camera with a 320 × 240 pixels infrared resolution and operating in a range of temperatures from −20 to 650 • C (with ±2 • C accuracy) and a thermal sensitivity < 0.04 • C. The rock emissivity was set to 0.975, suitable for the rock type tested according to the literature [37].
Rock IR monitoring had a 5 min cadence, thus allowing the final computation of CRI within the first 5, 10, 15, and 20 min from the cooling beginning.
Thermal images were post-processed through software capable of handling radiometric images and extracting the average surface temperature values corresponding to the framed rock specimen face ( Figure 2).

Rock Porosity Estimation
The porosity of the reference group of tuff specimens (Cycle 0) was estimated according to EN1936:2006 by the long-established Archimedes buoyancy method on 16 specimens. This method relies on the measure of three parameters of a single specimen: the weight of the dry specimen (w d ); the weight in air of the specimen fully saturated with water (w s ); and the weight of the saturated specimen suspended and fully immersed in water (the Archimedes weight w h ) [38]. From these weights, and after determining the real rock density (ρ r ) through the pycnometer method, total (n) and effective (n' H 2 O ) porosities can be calculated according to Equations (2) and (3), where ρ b is the bulk density. Some specimens started showing a loose structure from the second cycle, sometimes crumbly even to the touch for the most advanced cycles. This made them unsuitable for saturation with water. Therefore, to avoid errors in the porosity estimation due to both possible salt crystal dissolution and partial loss of specimen mass, porosity was computed by a multianalytical approach using structured-light 3D scanning, helium pycnometry, and mercury intrusion porosimetry (MIP) on 10 representative weathered samples from Cycles 2, 4, 6, 10, and 23.
By structured-light 3D scanning, a Scan-in-a-box FX scanner was combined with the software IDEA 1.1 SR 8 FX for building 3D digital models of the samples and calculating their volume, averaging from three independent measurements. The weighing of the samples provided their bulk density.
By helium pycnometry, a Quantachrome Pentapyc 5200e PPY-30T pycnometer was operated at 25 • C constant temperature and 131 kPa maximum pressure for determining the matrix volume and, indirectly, the matrix density of the samples, averaged among three measurements. Total porosity and pore volume were also calculated.

Statistical Data Processing
CRIs within 5, 10, 15, and 20 min of cooling were plotted against the n and n' H 2 O values for the untreated specimens, measured according to the EN1936:2006 specifications. A simple regression analysis was carried out, and the best-fit curves and the coefficients of determination (R 2 ) were calculated by the "least squares" method to determine if the correlation between porosity and CRI could be ensured. The achieved statistical trends were then compared to the literature trends [21], referring to specimens similar in shape and size to those tested herein, to prove the reliability of the testing principle. For the other tuff groups (from Cycle 2), CRIs were correlated to n He and n' Hg porosity to find out if significant correlations can be achieved as well.

Results of Salt Crystallization Tests
A summary of the results of the salt crystallization tests is shown in Figure 3. Demjén tuff has the lowest durability to salt crystallization, crumbling already during the first aging cycles. Differential erosion is the main weathering pattern, highlighting the faster disintegration of the groundmass, compared to the slower detachment of the coarser-size phenocrysts.
The resistance of Bogács and Sirok tuffs is distinctively higher instead. Most specimens show negligible macroscopic changes and mass variations (excluding the weight increases due to the continuous salt loading) throughout the experimental program. Bogács tuffs record just the formation of a few voids on the surface and a slight, occasional peeling, which becomes more evident during the very last cycles. On the other hand, several specimens of Sirok tuff start deteriorating after the 5th cycle, with traces of differential erosion, disintegration, and microcracking, leading to mass decreases up to~10%. These become even clearer during the following cycles, reaching drops of~25%-30% for some specimens at the end of the tests. The strength variability observed in Sirok tuffs, dependent on changing textural characteristics (e.g., different amounts of pumice clasts, glassy groundmass, coarse phenocrysts and lithoclasts, welding degree), is in agreement with previous experimental observations [23]. eral specimens of Sirok tuff start deteriorating after the 5th cycle, with traces of differential erosion, disintegration, and microcracking, leading to mass decreases up to ~10%. These become even clearer during the following cycles, reaching drops of ~25%-30% for some specimens at the end of the tests. The strength variability observed in Sirok tuffs, dependent on changing textural characteristics (e.g., different amounts of pumice clasts, glassy groundmass, coarse phenocrysts and lithoclasts, welding degree), is in agreement with previous experimental observations [23].

Results of Porosity Estimation
The initial rock physical characterization, carried out through the method of Archimedes buoyancy in water, returned dissimilar porosity characteristics among the tested tuff varieties. In particular, Demjén tuff has the highest total porosity and the lowest bulk density (Table 1), statistically affected by the low standard deviation. This is an index of textural uniformity. On the other hand, Bogács tuff is characterized by the lowest total porosity and the highest bulk density. Among the studied tuffs, Sirok shows the highest standard deviation for total and effective porosity ( Table 1). The degree of pore interconnection is the highest for Demjén, with an effective porosity of 72.6% of the rock volume, and the lowest for Bogács, with an effective porosity representing only 54% of the total porosity. The high pore volume and interconnection in Demjén tuff is associated with a distinctively higher water absorption and salt penetration, which explains the much lower durability to salt crystallization, especially since pore-size distribution is quite similar to Bogács tuff. In addition, the lower tensile strength implies a poorer mechanical performance when in-pore stresses are generated during salt crystallization and dissolution [23]. As the salt crystallization tests go on, considerations can be made only on the strongest specimens of Bogács and Sirok, although no relevant differences between the two tuffs can be highlighted as the cycles progress, with the relevant n' Hg and n He porosity values keeping almost constant trends. This can be explained by the progressive salt accumulation within the rock pores, which balances the ongoing rock weathering, hiding its effects on the porosity. This consideration is strengthened by the positive trend of mass change shown in Figure 3.
Moreover, considerations about the possible changes in pore-size distribution as salt weathering progresses can be presented only about Bogács and Sirok tuffs, supported by the full MIP datasets in Supplementary File S1. Bogács tuff shows no significant changes in pore-size distribution, characterized by capillary pores in the approximate range between 1 and 10 µm, which stays virtually the same during the continuation of the tests. Sirok tuff has a higher concentration of micropores, with broader and multimodal pore-size distribution curves, with peaks mainly in the approximate range between 0.01 and 1 µm; no clear trend of changing pore-size distribution during the experimental tests is noticeable, apart from a certain degree of disturbance (i.e., higher pore-size variability) from the 4th salt crystallization cycle onwards in Table 2. However, since the number of capillary pores (those most involved in the mechanisms of water absorption and movement) is lower in Sirok tuff, the prosecution of the salt tests with a higher number of cycles might eventually highlight the higher durability of this material [23].

Results of IR Cooling Monitoring
The analysis of the thermal images acquired for the untreated tested specimens (Cycle 0) during the cooling monitoring ( Figure 2) shows that the rock cooling proceeds towards the inner rock in a generally uniform way, according to literature accounts [21], thanks to the circular geometry of the specimen' face. Thermal images keep a good definition for the entire monitoring time windows, allowing the good identification of the rock face and sampling its average surface temperature (Figure 2).
According to the estimated temperature values and the calculated cooling rates, the untreated specimens of the Demjén variety were characterized by the highest CRIs, suggesting the fastest cooling among the tested groups, followed by Sirok and Bogács (Figure 4a). For all the cases, the standard deviation is low, thus suggesting the reliability of the considered average values. The highest CRIs are referred to as the 5 min time window, when the temperature difference between the specimen and the laboratory environment is maximum. This is in accordance with Newton's law of cooling and previous literature accounts on this test [19]. Therefore, the cooling speed is maximum at the initial stages of the test and tends to slow down as time passes, as shown by the progressively decreasing angle of the cooling curves (Figure 4b). By analyzing the representative cooling curves reported in Figure 4b, it can be stated that the Demjén variety is characterized by the fastest cooling, considering the common initial stage and the lowest surface temperature reached at the end of the monitoring. This trend is more apparent even focusing on the initial segment, used to calculate CRI 5 , which is steeper for Demjén, thus suggesting the highest cooling rate within the first 5 min of monitoring. Similarly, the less steep curve is related to the Bogács variety, thus suggesting a lower cooling rate, and supporting the average data reported in Figure 4a.

CRI and Rock Porosity on Untreated (Cycle 0) Specimens
To shed light on the statistical correlation between CRI and rock porosity, experimental cooling data referring to Cycle 0 were plotted against both n and n'H2O values (Figure 5). Scatter plots show a positive, linear trend between the variables, proving that the fastest rock cooling is related to the most porous rocks. This outcome, verified herein for By looking at the average CRI trend as the cycles progress (the standard deviation related to each average value is generally <0.1), an initial increase can be reported, with respect to Sirok and Bogács, between Cycle 0 and Cycle 2 (Figure 4c-f). This is about 6%-7% and 5%-12% for Sirok and Bogács, respectively. On the other hand, on average, a CRI decrease of around 1% was found for Demjén (Figure 4c-f). It has to be noted that after the 2nd cycle, no Demjén specimen resisted the crystallization tests, thus the following analysis is referred only to Bogács and Sirok. The related mean CRIs tend to decrease according to an almost constant trend until the 10th cycle. Only the average CRI 5 of Sirok measured between the 4th and the 6th shows a minimal difference (Figure 4c). From the 10th to the 23rd cycle, Sirok's CRIs keep an almost constant value, while Bogács's values mark a slight increase of about 2%-4%.

CRI and Rock Porosity on Untreated (Cycle 0) Specimens
To shed light on the statistical correlation between CRI and rock porosity, experimental cooling data referring to Cycle 0 were plotted against both n and n' H 2 O values ( Figure 5). Scatter plots show a positive, linear trend between the variables, proving that the fastest rock cooling is related to the most porous rocks. This outcome, verified herein for the tested tuffs' total and effective porosity, is in accordance with literature accounts [19][20][21][22]. CRI10, and CRI15 being characterized by a 0.85 R 2 (Figure 5b). Based on these outcomes, the potential of CRI10 for rock porosity prediction is herein strengthened according to Equations (4) and (5). These describe the best trends achieved with the non-standard cylindrical specimens tested for this study, which had a 50 mm diameter and a 30 mm height (indicative volume of 59 cm 3 ). This is significant, since there is no literature experience on such specimen size. The only available literature data of rock specimens with a somewhat similar size, tested by IRTest, are related to 1:1 small rock cylinders with a 54 mm diameter (indicative volume of 124 cm 3 ) [21]. It is, therefore, valuable to compare such data to evaluate their scientific soundness. To this purpose, the experimental CRI10 statistical plots of the Hungarian tuffs were overlapped with those published in the literature [21] related to the 1:1 small cylinders, bearing in mind the volume difference (Figure 5c,d). It can be noted that the rocks tested herein are, indeed, characterized by a higher cooling rate due to their smaller volume. Nevertheless, a balanced trend comparison should be carried out by taking into account only the specimens showing similar porosity values. Therefore, by focusing the analysis only on the greatest porosity ranges (n > 25% and n'H2O ≥ 10%), the statistical trends appear subparallel, and the highest R 2 values belong to the Hungarian tuffs tested for this study in both n and n'H2O cases (Figure 5c,d).
% = 40.13 − 53.45 (5) Figure 5. Statistical relationship between thermal cooling and rock porosity for untreated (Cycle 0) specimens: (a) correlation diagram CRI vs. total porosity; (b) correlation diagram CRI vs. effective porosity; (c) CRI10 vs. total porosity: comparison between data from this study and literature [21] (dotted lines are the best fit of the whole sets, while the dashed black line is the best fit of literature data involving only specimens with n > 25% for a more specific comparison); (d) CRI10 vs. effective porosity: comparison between data from this study and literature [21] (dotted lines are the best fit (c) CRI 10 vs. total porosity: comparison between data from this study and literature [21] (dotted lines are the best fit of the whole sets, while the dashed black line is the best fit of literature data involving only specimens with n > 25% for a more specific comparison); (d) CRI 10 vs. effective porosity: comparison between data from this study and literature [21] (dotted lines are the best fit of the whole sets, while the dashed black line is the best fit of literature data involving only specimens with n' H 2 O ≥ 10% for a more specific comparison).
The best correlation of total porosities is achieved by CRI 10 , with a trend characterized by a 0.74 R 2 against the minimum value of 0.69 belonging to the correlation with CRI 20 (Figure 5a). On the other hand, the best statistical correlations involving n' H 2 O refer to CRI 5 , CRI 10 , and CRI 15 being characterized by a 0.85 R 2 (Figure 5b). Based on these outcomes, the potential of CRI 10 for rock porosity prediction is herein strengthened according to Equations (4) and (5). These describe the best trends achieved with the non-standard cylindrical specimens tested for this study, which had a 50 mm diameter and a 30 mm height (indicative volume of 59 cm 3 ). This is significant, since there is no literature experience on such specimen size. The only available literature data of rock specimens with a somewhat similar size, tested by IRTest, are related to 1:1 small rock cylinders with a 54 mm diameter (indicative volume of 124 cm 3 ) [21].
It is, therefore, valuable to compare such data to evaluate their scientific soundness. To this purpose, the experimental CRI 10 statistical plots of the Hungarian tuffs were overlapped with those published in the literature [21] related to the 1:1 small cylinders, bearing in mind the volume difference (Figure 5c,d). It can be noted that the rocks tested herein are, indeed, characterized by a higher cooling rate due to their smaller volume. Nevertheless, a balanced trend comparison should be carried out by taking into account only the specimens showing similar porosity values. Therefore, by focusing the analysis only on the greatest porosity ranges (n > 25% and n' H 2 O ≥ 10%), the statistical trends appear sub-parallel, and the highest R 2 values belong to the Hungarian tuffs tested for this study in both n and n' H 2 O cases (Figure 5c,d).

CRI and Rock Porosity after Salt Crystallization Cycles
Similarly to what was carried out for the fresh samples, a simple regression analysis between CRI and rock porosity was performed to determine if the statistical relationship persists for weathered rocks (Figure 6a,b). In this case, the CRIs calculated for the rocks belonging to Cycle 2 were taken into account. The results show that the statistical trends described by CRI and n fail to equal the previous ones in terms of statistical goodness (Figure 6a). In fact, more scattered plots were achieved, with maximum R 2 values of 0.53. The worst correlation (R 2 0.23) was found for CRI 5 , while the best one refers to CRI 10 . A similar outcome was outlined for n' H 2 O , although with slightly higher R 2 values, with particular reference to CRI 15 (Figure 6b).
Considering the satisfactory correlations related to the untreated specimens, along with the available literature experience, the poor statistical significance achieved in this case casts doubt on the reliability of the procedure used to estimate rock porosity for the weathered specimens. In fact, the buoyancy method requires the specimen immersion in water for at least 24 h, to achieve saturation. This procedure likely affected the measurements in two ways: (1) weathered rocks tend to gain a looser structure, favoring a certain loss of mass during immersion; (2) salt crystals precipitated within the rock pores could likely be involved in dissolution processes when the rock is saturated. Both circumstances lead to possible dissimilarities between wet and dry specimen masses for the porosity calculation, likely making the so-measured porosity values unrepresentative. Therefore, with the aim of looking for a better statistical relation for such non-standard rock specimens, the n' Hg and n He porosity values, estimated on fragments belonging to 10 selected rock specimens (see Section 3.3), were taken into account. These data were plotted against the CRI values, regardless of the salt crystallization cycle (from Cycle 2 to Cycle 23), and the best fit was achieved for the correlations involving CRI 20 (Figure 6c,d; Equations (6) and (7)).
n Hg = 30 CRI 20 − 12.6 (6) n He = 36.6 CRI 20 − 15.95 (7) Minerals 2023, 13, x FOR PEER REVIEW 13 of 16 Figure 6. Statistical relationship between thermal cooling and rock total porosity, calculated through the buoyancy method, for the specimens subjected to 2 cycles (a); statistical relationship between thermal cooling and rock n'H2O for the specimens subjected to 2 cycles (b); best statistical trend achieved by plotting n'Hg vs. CRI20 (c); best statistical trend achieved by plotting nHe vs. CRI20 (d).

Discussion and Conclusions
The results achieved by this study confirm that the rock cooling process, monitored by IRT through the IRTest, is related to porosity. This is demonstrated by satisfactory linear statistical trends obtained by plotting the rock porosity against the calculated CRIs. It was confirmed that the cooling of rock is directly proportional to the rock volume occupied by voids, concerning both total and effective porosity. In fact, the higher the rock porosity, the faster its cooling. According to the literature background on rock properties and Newton's law of cooling, the cooling velocity of oven-heated rocks is characterized by the highest rates during the first minutes. Then, it tends to slow down as the temperature difference between the rock and the surrounding environment gets smaller. This is confirmed by the trend shown by the cooling curves taken as examples in Figure 4b. The results show also that the cooling rate measured within the first 10 min of cooling (CRI10), which is known in the literature to be a good porosity predictor [19][20][21][22], keeps its prediction potential thanks to the best correlations achieved with n and n'H2O. Its role was strengthened herein by comparing literature data and the experimental results referring to the untreated Hungarian tuffs specimens (Figure 5c,d).
A further consideration is for the specimens' size; in fact, rocks tested for this study had about half of the volume of the specimens from the study taken as Reference [21]. It resulted that the rocks tested herein are characterized by a higher cooling rate due to their smaller volume, underlying that the rock volume plays a key role in the cooling speed, with the highest rates related to lower volumes. Although this study demonstrated that the CRI works fine also for such small rock specimens, there is the need of referring to standardized sample dimensions in the perspective of the potential standardization of such approaches. Figure 6. Statistical relationship between thermal cooling and rock total porosity, calculated through the buoyancy method, for the specimens subjected to 2 cycles (a); statistical relationship between thermal cooling and rock n' H 2 O for the specimens subjected to 2 cycles (b); best statistical trend achieved by plotting n' Hg vs. CRI 20 (c); best statistical trend achieved by plotting n He vs. CRI 20 (d).

Discussion and Conclusions
The results achieved by this study confirm that the rock cooling process, monitored by IRT through the IRTest, is related to porosity. This is demonstrated by satisfactory linear statistical trends obtained by plotting the rock porosity against the calculated CRIs. It was confirmed that the cooling of rock is directly proportional to the rock volume occupied by voids, concerning both total and effective porosity. In fact, the higher the rock porosity, the faster its cooling. According to the literature background on rock properties and Newton's law of cooling, the cooling velocity of oven-heated rocks is characterized by the highest rates during the first minutes. Then, it tends to slow down as the temperature difference between the rock and the surrounding environment gets smaller. This is confirmed by the trend shown by the cooling curves taken as examples in Figure 4b. The results show also that the cooling rate measured within the first 10 min of cooling (CRI 10 ), which is known in the literature to be a good porosity predictor [19][20][21][22], keeps its prediction potential thanks to the best correlations achieved with n and n' H 2 O . Its role was strengthened herein by comparing literature data and the experimental results referring to the untreated Hungarian tuffs specimens (Figure 5c,d).
A further consideration is for the specimens' size; in fact, rocks tested for this study had about half of the volume of the specimens from the study taken as Reference [21]. It resulted that the rocks tested herein are characterized by a higher cooling rate due to their smaller volume, underlying that the rock volume plays a key role in the cooling speed, with the highest rates related to lower volumes. Although this study demonstrated that the CRI works fine also for such small rock specimens, there is the need of referring to standardized sample dimensions in the perspective of the potential standardization of such approaches.
Besides the considerations on the untreated specimens, which are essential to have the IRTest validated even for non-standard-sized rock specimens, interesting discussion points arise from the tests carried out on rocks treated by salt crystallization. From the 2nd to the 23rd cycle, n' Hg and n He were taken as a reference for the pioneering correlation against CRI. Positive linear trends were achieved by considering the 20 min time window (CRI 20 ) for all the selected specimens (Figure 6c,d). This is double the most suitable time window for untreated rocks (10 min) and the reason for that should be searched in the effects that salt crystallization has within the rock texture. This aspect needs to be further investigated, even questioning if an application for the specific IRTest can be posted in terms of crystallization cycles. Figure 4c,f show that CRIs between the 10th and the 23rd cycle remain almost constant for Sirok and are affected by an inverse trend for Bogács. However, the good match between rock cooling and Hg-He porosity suggests that porosity plays a leading role in rock cooling and that IR monitoring is a useful approach for its indirect estimation, even for weathered rocks.
Although needing further experimentation, the application of this method to stone materials used in the Hungarian cultural heritage opens interesting possibilities for the imaging and textural characterization of historical materials. At this stage of development, the technique provides a promising perspective in the frame of stone testing for cultural heritage, although available results refer only to the laboratory setting. Moreover, it is totally non-destructive, since it works with temperatures that do not cause significant changes in structural water content, irreversible strains, or mineralogical transformations. In the frame of technological development, more research is needed to find and set a comparable procedure that could be extended for in situ tests.