Marble decay: towards a measure of marble degradation based on ultrasonic wave velocities and thermal expansion data

Marble as ornamental and dimensional stones as well as in their natural environments show complex weathering phenomena. Physical, chemical, and biological weathering of marble are well documented. The impact of climate change on monuments and historic buildings in terms of modeling and predicting future scenarios requires new approaches to forecast the ongoing decay in the near and far future. Ultrasonic wave velocities are a powerful and sensitive tool for the damage assessment of marble. For a maximum porosity of up to 1%, ultrasonic wave velocities (P-wave velocities) are ranging between 1 km/s and over 6 km/s. Water saturation has an important influence on the magnitude and directional dependence of ultrasonic wave velocities together with the mineralogical composition and the rock fabrics. Ongoing experimental alteration approaches were used to document the state of deterioration using Vp-systematics. In addition, thermal expansion and the residual strain values after applying thermal impacts were used to introduce a new quantitative measure based on experimental length changes and volume changes. To quantify such volume changes, a so-called decay index was proposed. Marbles are sensitive to weathering and have different volume changes under exposure depending on fabric parameters. The volume extension index of marble, based on thermal expansion measurements under dry and water-saturated conditions, is proposed as a decay index for quantifying sample stability and for defining the directions of maximum and minimal dilatation. Such decay index was implemented to different marble types and it was turned out that marbles with the larger decay indexes are more prone to weathering than with smaller ones. The effect of changing climate and, in consequence, different weathering actions can help to calculate or forecast risk numbers based on the Vp data in combination with the proposed decay index especially for marbles.


Introduction
Marbles have been considered to be among the most important building materials since ancient times. Regarded as special stones, they are captivating because of their pureness, formability, and translucence. However, degradation of sculptures, architectural heritage, and façade stones fabricated from marbles can be problematic. There are several approaches to determine weathering rates. One approach could be to carry out experiments in the laboratory or to study natural weathering processes over a very long period. Within the framework of the international program MULTI-ACCESS (Model for multi-pollutant impact and assessment of threshold levels for cultural heritage) coordinated by the Swedish Corrosion Institute in Stockholm, damage functions have been developed for various materials as well as for the Portland limestone (Brimblecombe 2003) The relationship between the released air pollutants, meteorological influencing variables, and the resulting material losses was determined. The worldwide storage tests resulted in the so-called "multi pollutant situations", which yielded damage functions 1 3 395 Page 2 of 34 calculated on the basis of statistical methods. Besides SO 2 and pH values of precipitation, emissions of NO 2 and PM 10 were also considered in the equations. Moreover, papers predicting carbonate stone surface recession for chemical dissolution, induced by changes in precipitation amount and atmospheric CO 2 concentration, have been well known since a very long time. In particular, all these works stressed the increasing importance that climate parameters will have in addition to air pollution, in the decay.
The deterioration of building stones depends in the future mainly on the climate. Many experimental studies have shown that temperature changes, both increases and decreases, induce significant deterioration. Even when the temperature changes are not particularly large, repeated heating and cooling of stone will eventually lead to deterioration over time (e.g., Kessler 1919;Battaglia et al. 1993;Winkler 1994;Siegesmund et al. 2000;Zeisig et al. 2002). This deterioration is most spectacular for marbles, but is also observed for other rocks, like granite, Nagelfluh, and limestone (e.g., Siegesmund and Dürrast 2011). The paper of Bonazza et al. 2009 deals with the prediction of thermal stresses on stone in Europe over the next century. The thermoclastism risk on marbles for the twenty-first century is discussed using a revised damage function. The authors figured out that the susceptibility of marbles derives from the anisotropy in thermal response of calcite, which expands in the vertical c-axis and contracts in the horizontal one. Research carried out so far indicates that marbles are less resistant towards thermoclastism compared to other rock types like, e.g., granites, limestones or sandstones, etc. According to Timoshenko and Goodier (1970), the maximum horizontal thermoelastic stress can be calculated by combining the Young's modulus, the Poisson's ratio, the thermal expansion, and the temperature. Based on this equation, Scheffzük et al. 2004 calculated a thermoelastic stress of around 42MPa for a marble. Considering a strong crystallographic preferred orientation of the calcite crystals, a thermoelastic stress can reach up to between 201 MPa for a temperature difference of 100K. Additionally, Scheffzük et al. 2007 had calculated a residual stress magnitude of 77.5 MPa, i.e., the so-called locked-in stress. Already, Winkler 1994 mentioned that prestressed crystalline marble tends to expand more during the first cycle of heating the locked-in stresses are released. Considering the processes on marble deterioration, it is for sure that the thermal expansion is only important for 1-3 thermal cycles afterwards any residual stresses, i.e., progressive increase in length change occur when water is present (see Siegesmund 2004 or Siegesmund et al. 2008) Thermal degradation in the form of microcracking results in a volumetric expansion of a stone sample. Accordingly, microcracking of a stone sample upon heating or cooling will result in a permanent volumetric expansion of a stone sample when it returns to its ambient temperature. This relative volume change as measured by thermal dilation is thus a measure or index of stone degradation. Thus, the overall objective of this approach is to develop a volumetric extension index (a deterioration index) for the characterization of the degradation potential of widely used types of marbles and to elucidate the influence of fabric parameters on this deterioration index. The index will reflect the relative volume change on the basis of the knowledge of thermal dilatation.
To get an insight into their weathering state, a nondestructive analysis of them is desired. Ultrasonic wave velocity systematics and US tomography are possible tools to receive the desired insight, while a sound interpretation depends on the knowledge regarding the specific material. Systematic analyses were performed on six marble types, which were often used as a heritage stone (Table1, Fig. 1): the Spanish Blanco Macael, the Italian Carrara Gioia and Lasa, the Polish Großkunzendorf and Prieborner, and the less common German Wunsiedel marble which is primarily used regionally. To relate their weathering state to an ultrasonic velocity, a successive artificial weathering has been performed, analyzing all important rock physical changes within the alteration (Menningen 2020). In addition, the thermal expansion systematics were also performed on several Carrara varieties, Sölk, Wachau, Thassos, and Volakas marbles (Table1, Fig.1).

Rock materials, macroscopic description, and decay phenomena
The investigated marbles originated from 15 different locations. They are well known and widespread used as building stones (Blanco Macael, Rosa Estremoz, Carrara (several varieties), Großkunzendorf, Prieborn, Wunsiedel, Lasa, Thassos, Volakas, Sölk, Wachau) in the past and also in our days. Mineralogically all samples belong to calcitic and dolomitic marbles with quartz, biotite, muscovite, phlogopite, and ore minerals as accessories. The most important characteristics of the fabric elements are summarized in Table 1.
The average grain size which was calculated from microscopical observations or from the digitized images of the grain boundaries varies from fine-grained to medium-grain size 50 μm up to 3 mm. All samples are also characterized by a large variation considering the shape fabric examples are given in Figs. 1,2,3,4,5,6,and 7. The grain boundary orientation given for the XY-, XZ-, and YZ-planes may differ from a more or less random distribution to a distinct planar or linear shape fabric. An example for a random distribution is the Arabella marble, while most of the other marbles (e.g., Sölk, Thassos and Carrara) exhibit a shape preferred orientation (Fig. 1). A more detailed discussion on the fabric is given for Blanco Macael, Carrara, Lasa, Großkunzendorf, Prieborn, and Wunsiedel marbles, because they are more or less representative for the whole collection and exhibit fabric characteristics which are important for the understanding of the observed petrophysical properties and degradation phenomena. The shape of the grain aggregates can generally be classified between equigranular-polygonal, inequigranular-polygonal, or seriate-interlobate in the sense of Moore (1970). Since a dolomitic composition is discussed as a main factor for the durability of marbles, dolomitic and  Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture calcitic marbles were introduced. Carrara (calcite) and Thassos (dolomite) are fine-grained marbles with more or less straight grain boundaries, a so-called foam structure. In contrast, Thassos and Sölk are coarser grained, while the grain boundaries can be recognized as equigranular-interlobate and inequigranular-interlobate. The much higher irregularity of the lobate grain boundaries is clearly documented in Thassos and Sölk exhibiting two sets of conjugate planar fabrics with an opening angle of around 50° which are slightly asymmetrically disposed to the foliation. These fabrics can easily be observed in the thin sections and in the corresponding grain boundary surfaces (Fig.1). In summary, the grain boundary geometry or the grain interlocking for all marbles ranges between equigranular-polygonal (Rosa Estremoz) up to seriate-interlobate (Großkunzendorf) with a slight tendency to be inequigranular-amoeboid (Wachau). This implies that the metamorphic and deformation history over geologically long time spans controlling the development of rock fabrics was quite different for all selected samples. Recrystallization processes, for instance, may control the grain size and the configuration of grain boundaries. Evidence for grain boundary migration recrystallization is the presence of highly irregular grain boundaries. The driving force is the difference in the dislocation density resulting Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture to a bulging of the grain boundaries into crystals with the higher dislocation density (Gottstein and Mecking 1985). In other cases, subgrain recrystallization may lead to an equilibrium fabric of polygonal crystals with interfacial angles of approximately 120° (see Carrara), a so-called foam structure. From the mechanical point of view, it seems to be clear that a decrease in grain boundary energy is correlated with the decrease in strength of a polycrystalline material. Beside a lot of other processes (for a review, see Skrotzki 1995), the deformation may be responsible for the development of a lattice preferred orientation (here referred as texture) which has also a significant control on the physical properties and their directional dependence (e.g., Siegesmund and Dahms 1995). All samples were analyzed with respect to their textures. Leiss and Ullemeyer (1999) discussed the fundamental texture types of calcite and dolomite which are found in nature. In summary, they can be described by a rotating single crystal with the c-axis or a-axis as the rotation axis, respectively. These c-axis and a-axis fiber types can overlap forming intermediate texture types. The recalculated c-axis and a-axes pole figures of Blanco Macael,Rosa Estremoz,Carrara,Großkunzendorf,Prieborn,Wunsiedel and Lasa,are illustrated in Figs. 2,3,4,5,6,and 7. Usually, a single maximum which is in most cases elongated to an Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture Page 7 of 34 395 oval-shaped pattern can be observed. According to the crystallography of calcite, the (110) poles are arranged on a great circle around the (001) pole density maximum. The maximum intensity is highly variable and at maximum with 5.9 m.r.d. (multiples of random distribution) of (001) indicates a strong texture and, consequently, a pronounced anisotropy of the physical properties must be expected. To characterize the textures of all 18 samples, the tensor shape T calculated from the pole figure tensor is applied to describe the different c-axis concentrations (Jelinek 1981). T calculated from the pole figure tensor characterizes the texture differences very well. It covers the range from T = -0.859 (Wachau) to T = 0.881 (Rosa Estremoz) indicating a well-pronounced cluster-like and a moderate girdle-like shape of the intensity distribution as the extreme cases. Extremely planar c-axes distributions (i.e., complete girdles) are not observed. The maximum intensity obtained from a simplified texture reproduction with L = 2 ranges from 1.4 to 5.9. It more or less concurs with the shape parameter T.

Petrophysical properties
The alteration of marble leads to a significant porosity increase and change of the pore space properties. Generally, the porosity has a direct and indirect effect on most physical rock properties and is the most important rock Main fabric characteristics of the Großkunzendorf marble. Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture parameter (Siegesmund and Dürrast 2014). An increasing porosity has often an unfavorable influence on the weathering behavior, especially if the pore radii shift to the pore classes of 0.1-1.0 µm, the so-called capillary pores (Klopfer and Homann 1985). Those pore sizes change the water transport mechanisms of crystalline marble significantly and increase the weathering process. The quarry fresh effective porosity of crystalline marble analyzed by hydrostatic weighing is very low ranging about 0.1-0.46 % (Koch 2006), 0.2-0.3 % (Rüdrich 2003;Rüdrich et al. 2013), and up to 0.57 % (Poschlod 1990). The alteration of marble causes an increase of this porosity which is still very low and ranging between 0.4 and 0.7% up to 2.33 % reported for a Carrara marble by Rüdrich et al. 2013, but effective regarding the resulting modifications of the rock properties. Therefore, the change of porosity within the stepwise artificial alteration was determined for 7 samples ( Table 2). The measured porosities (φ) range between 0.12 % and 0.41 % for the quarry fresh marble samples, while the varieties Prieborn with φ=0.44 % as well as Carrara Gioia with φ =0.37 % were already weathered at their 'fresh' state. Over the weathering course, the Rosa Estremoz and the Wunsiedel marble Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture Fig. 7 Main fabric characteristics of the Wunsiedel marble. Characteristic line drawing of a microscope image, grain size [relative frequency of specific area class has been calculated with respect to the total area (grey bars; A total area) and total number of grains (black bars; N total number of grains)], shape preferred orientation and texture shows a different behavior within the whole artificial alteration. The porosity increase is low and even the last cycle did not show any significant increasing effect, but this marble shows a comparatively high standard deviation (0.14-0.19), which increases after each weathering step. This reflects the macro-and microscopic visible heterogeneity and shows that the crack growth is not uniform for this marble and depends on specific conditions.

Capillary water uptake
Almost all weathering processes as freeze-thaw cycles, solution processes, salt-related weathering, and thermohygric processes are determined by the presence of water. The pore space and pore system relates to the amount of water absorption, the nature of water transport, and the size of the surface area available for chemical interactions (Siegesmund and Dürrast 2014). Capillary water uptake is an important parameter for materials with pore diameters within the range of the so-called capillary pores. To determine the transformation within the artificial alteration, the capillary water uptake was analyzed for all marble varieties. The capillary water absorption coefficient [kg/m 2 *√h] given in Table 3 was analyzed for one cubic sample in the unweathered state and after heating cycles up to 200 °C. The w value is given for 3 directions, as well as the water absorption (Wcap) in weight and volume percent. The porosity can also be found in Table 2. The given water absorption value in weight percent is given as mean value of all analyzed cubes under the specific conditions (10 samples at the fresh state, 5 samples for the 200 °C state). The w value can be used to classify rocks (Snethlage 2008) into three types (1) slightly absorbing [w < 0.5 (kg/m 2 *√h), (2) medium absorbing (w 0.5-3.0 kg/m 2 *√h), and (3) highly absorbing [w > 3.0 (kg/m 2 *√h)]. According to this classification, all analyzed marble varieties are slightly absorbing regardless if they are quarry fresh [w values between 0.02 and 0.17 (kg/m 2 *√h) or artificially weathered (0.09-0.45 [kg/m 2 *√h], but all of them increase according to the weathering intensity. Only the strongly weathered Großkunzendorf shows a slightly higher w value of 0.55 [kg/m 2 *√h] in X-direction. The thermal-induced deterioration of all marble types shows an increasing capillary water uptake capability, which is connected to the increasing porosity (see Siegesmund and Dürrast, 2014). Regardless of the weathering state, it should be noted that the capillary water uptake occurs relatively fast and is finished within the first 60-180 min. The saturation coefficient S describes how much of the total pore space is accessible to water absorption and helps to determine the frost resistance. Siegesmund and Dürrast 2014 describe that for metamorphic rocks, the saturation coefficient has limited significance due to the very small water absorption values under atmospheric and under vacuum conditions. They are influenced by the measurement procedure, the standard deviation of different samples, and the resolution limit. According to the classification of Hirschwald 1912, S values below 0.75 stand for a weather and frost resistant rock, while values below 0.9 indicate a rock which is not frost and weather resistant. For any rock type with a saturation coefficient between 0.75 and 0.9, further investigations are recommended. The analyzed marble varieties represent the entire spectrum of the S value with 0.57 for the Blanco Macael up to 0.93 for the Carrara Gioia, whereby most of the samples are assigned to the group of rocks, where further investigations are recommended.

Ultrasonic measurements
The correlation of the P-wave velocity determined by ultrasonic measurements and the weathering state of marble was confirmed by many scientific analyses (Ahmad et al. 2009;Babacan and Gelisli 2015;Boudani et al. 2015;Dürrast et al. 1999;Fleischer 2002;Köhler 1991Köhler , 2014Lindner et al. 1999;Menningen et al. 2018a, 2018b, 2019, Rüdrich 2003Rüdrich et al. 2013;Siegesmund et al. 2000a;Siegesmund et al. 2004b;Weiss et al. 1999;Weiss et al. 2002a). Whereas the P-wave velocity of a quarry fresh marble is relatively high [6-7 (km/s)], it decreases with increasing weathering state. This velocity decrease is connected to the formation of microcracks throughout the weathering process, which has a huge influence on the velocity reduction (e.g., Weiss et al. 2002a). Köhler 1991, established a velocity-based classification for the weathering state of Carrara marble which helps to define the weathering state of objects made of marble. As those measurements allow a non-destructive evaluation of the material, they are used for quality assurance as well as for preservation purposes in the context of heritage conservation.

Transmission measurements
The distribution of elastic waves depends on the mineralogical composition, the fabric as well as the type and arrangement of cracks. Furthermore, the crack filling medium, water or air, has a significant impact (e.g., Siegesmund 1996;Weiss et al. 2002a). To assign a specific weathering state to an ultrasonic velocity, the systematic alteration of the marble samples was connected to repeated ultrasonic measurements. All cubic samples were analyzed after each temperature cycle under dry and wet conditions. In Table 4, the mean dry ultrasonic velocities are given in three directions as well as the amount of anisotropy (%). The variation of the maximum and minimum velocity is used to calculate the anisotropy A(%) = (Vpmax−Vpmin)/Vpmax*100 (Siegesmund and Dürrast 2014). The Blanco Macael shows the strongest and most distinct velocity decrease after almost every weathering cycle with a total reduction down to 1.60 km/s which is a reduction of 65.58% in comparison to the unweathered velocity in Z-direction. The strongest impact is caused by the heating cycle up to 200 °C, with an average reduction of 27.4% in comparison to the 150 °C cycle. The average velocity reduction per weathering step is 0.51 km/s (±0.23) and combined with a steadily increasing anisotropy up to 30.7%. In comparison, the Lasa marble shows an alternating intensity of the velocity reduction. Here, every second weathering step shows a higher reduction and reaches its strongest influence at the combination of 200 °C heating and freezing. That this combination has the biggest impact on the velocities can also be observed for the Großkunzendorf, Wunsiedel and the Lasa marble (see Fig. 8). The amount of anisotropy increases with each weathering step. While the Carrara Gioia shows almost no changes and the anisotropy ranges between 5.24 and 8.84 %, the Wunsiedel marble shows the highest initial anisotropy of 17.36% and increases successive up to 34.05%. It can be observed that the main velocity reduction for each sample takes place in Z-direction, whereas the Großkunzendorf shows an equalization after the second heating-freeze cycle. This may indicate a foliation parallel preferred crack formation.

Influence of humidity and water saturation
The weathering state of marble can be determined by ultrasonic measures, as long as the material is dry. The effect of an increasing porosity caused by cracks is reduced as soon as water penetrates the pore space and increases the velocities. This increase can reach almost intrinsic values for a fully saturated sample. The ultrasonic velocities for the water-saturated samples after each heating cycle are given in Table 5. There, the mean values, out of 4-9 cubes and 5 measurements per direction, are displayed. Independent of the weathering state, all cubes and varieties show a significant increase of the compressional wave velocity, covering almost the span of velocities for a calcite single crystal [5.7 km/s parallel to the c-axis and 7.7 km/s parallel to the a-axes according to Dandekar (1968)]. Until a stress of 150 °C, the velocities range between 5.73 and 6.75 km/s, with a slightly fluctuation trend. After the next heating step up to 200 °C, these velocities are slowly but clearly reduced for all marble varieties, while the Großkunzendorf still keeps higher velocities. After the combination of heat and freezing temperatures, none of the marbles show higher velocities than 6.03 km/s. Rüdrich 2003 observed the same for strongly weathered Carrara and Großkunzendorf marble and named the lower compressibility of water and/or the incomplete water saturation as the reason. As all samples show the same behavior, the lower water compressibility seems to be more plausible, as the high number of cracks and, therefore, higher water proportion reduce the overall velocity. Regardless of this, the directional dependence of the wave velocity can still be observed for all samples and is consistent with the dry samples.

Intrinsic velocity distribution
As already described above, the water saturation of deteriorated marble can cause almost intrinsic velocity distributions. The intrinsic velocity distribution represents the crackfree compressive wave velocity, exclusively caused by the rock-forming minerals and their spatial arrangement. The P-wave velocities determined for a calcite single crystal are strongly anisotropic and differ with regard to the measured direction. A pronounced texture caused by the arrangement of those single crystals leads to an overprint of the anisotropic properties of calcite on the whole rock properties. Based on the textural data, it was possible to calculate the polycrystal compressional wave velocity for each sample and to quantify their directional dependence, known as intrinsic anisotropy (Table 6). The calculated velocities range between 6.96 and 6.32 km/s, while the anisotropy varies between 2.94 and 9.2%. As shown, the intrinsic anisotropy depends on the intensity of the texture (Fig. 11). The Großkunzendorf, with a low texture of 1.57 mrd, shows only a low intrinsic anisotropy. On the other hand, the Prieborner with a strong texture of 3.3 mrd has a high intrinsic anisotropy of 9.2%.

Ultrasonic wave velocities measured on spherical samples
To determine the spatial distribution of the velocity decrease, spherical samples were used to perform measurements in more than 100 directions. Those measurements were repeatedly executed for the same sample under dry and fully water-saturated conditions after each weathering step.

Crack-related velocity reduction
While the intrinsic velocity distribution is only influenced by rock-forming minerals, the measured velocity will always be lower. This velocity reduction is caused by cracks forming throughout the exhumation, exploitation process, and the further processing of the marble, and finally, it is affected by weathering processes. As shown, the alteration of marble significantly reduces the ultrasonic velocities due to the formation of cracks. To determine the specific crack evolution within the weathering process, spherical samples were analyzed after each weathering step. The knowledge about an existing preferential orientation of microcracks and their influence on the intrinsic velocity distribution is of particular importance with regard to ultrasonic tomography (Rüdrich 2003).
Therefore, the spherical samples were analyzed under dry conditions and compared to the almost intrinsic, watersaturated sample. The calculated ∆-Vp-Plot is the difference between the water saturated and dry sample and allows an approximately determination of the preferred orientation. In Figs. 12 and 13, the results for all varieties after their alteration up to 200° C are given.
The comparison of the water saturated, weathered, and calculated ∆-Vp-Plot shows that the main crack pattern caused by artificial weathering increases the intrinsic anisotropy (Fig. 14). The comparable distribution pattern of the water-saturated and the calculated pole figure points on this conclusion. It can be observed that former high ultrasonic velocities show the lowest velocity reduction, and on the other hand, the former low velocities show the highest reduction. While this general observation can be made, there are also several modifications between the different varieties. The samples Carrara Gioia and Wunsiedel show a crack pattern that is more comparable to the weathered distribution than to the water saturated. The pole figures of the varieties Blanco Macael, Lasa, and Großkunzendorf show more or less the same distribution pattern, independent if they are water saturated or weathered. Therefore, the calculated distribution pattern is comparable to them. In contrast, the results for the Rosa Estremoz differ. Here, we have a strong increase of the reduction parallel to the Z-direction, while the lowest increase can be observed in Y-direction. All anisotropies increase due to the weathering process, while the anisotropy of the ∆-Vp-Plot distinguishes between comparable anisotropies, as were found for the water-saturated sample (Blanco Macael), medium intensities, ranging between the intrinsic, water saturated, low anisotropy, and the high anisotropy of the weathered sample. The calculated Wunsiedel pole figure shows an anisotropy which is higher (54%) than the weathered anisotropy of 50%. This marble shows a strong foliation parallel preferred grain shape orientation, influencing the crack distribution pattern.
The fresh, unweathered sample shows velocities between 4.9 and 4.1 km/s, while the velocity distribution is primarily dominated by the texture, with a velocity minimum in Z-direction, slightly elongated in Y-direction and a maximum velocity in X-direction, slightly dipping into the Y-direction. Compared to the unweathered cubic samples, Fig. 11 Correlation of texture (mrd) and anisotropy (%) determined for seven analyzed marble varieties the sphere has lower velocities, which can be explained by the sample preparation process. The ∆-Vp-Plot indicates a crack distribution that increases the former intrinsic anisotropy. The first alteration step, heating up to 60 °C and cooling down to 4 °C seem to activate different mechanisms. The dry sample shows two velocity maxima with velocities between 4.4 and 4.5 km/s, following the axes belt, while the velocity minimum is still oriented perpendicular to the Z-direction. The ∆-Vp-Plot shows a more diffuse crack distribution, with two minima caused by the velocity maxima and two maxima with 2.8 and 2.7 km/s. This pattern could be caused by the crystallographic r-planes, which show sub maxima in this area. Even if the temperature was only increased up to 60 °C, the velocities decrease down to 4.5-3.6 km/s. The temperature step up to 90 °C decreases the maximum velocity down to 4.0 km/s and the minimum velocity down to 3.3 km/s, while the calculated difference between the maximum and minimum reduction decreases and lies between 2.8 and 2.3 km/s. The velocity maximum of 4.0 km/s has now almost the same position oriented in X-direction, while the lower velocities are rearranged in Z-direction. The crack pattern shows that a small channel in Y-Z-direction appears, decorated with spots of reductions between 2.8 and 2.6 km/s. Within the next heating cycle up to 120 °C, the main velocity reduction can be found in Z-direction. While the maximum velocity decreases from 4.0 to 3.9 km/s, the lower velocities decrease from 3.3 down to 2.7 km/s. This evolution can be observed by an increasing anisotropy, reaching values of 30.7%. The crack distribution is concentrated on the areas of the highest and lowest velocities, following the texture induced velocity distribution, while the maxima are still arranged as two points within the Z-Y-direction. While the next heating cycle causes only moderate velocity reductions, the ∆-Vp-pattern becomes more equal to the intrinsic distribution. The crack pattern develops more uniform, as the anisotropic influence has less preferred orientations than before. The last heating step was performed with an additional freeze-thaw cycle, which can be recognized due to the significant velocity reduction down to 1.9-1.4 km/s. As described before, the main crack pattern increases the intrinsic anisotropy. This exemplaric crack evolution could be observed for almost all analyzed samples. With increasing weathering state, the velocity distribution adapts to the intrinsic distribution and amplifies it.

Relationship of velocity and porosity
That a decreasing compressive wave velocities correlates with a porosity increase is known and was described by different authors (Mahmutoğlu 2017;Rüdrich 2003;Rüdrich et al. 2013;Siegesmund et al. 2010;Weiss et al. 2002a). The analyzed marble varieties were examined with regard to this relationship. In Fig. 15, the mean value for the different weathering states are given for the Blanco Macael (BM), Lasa (LA), Carrara Gioia (GI), and the Großkunzendorf (GK). The results for the Wunsiedel are not inserted, as they show a huge standard deviation and complicate the readability of the graph. Each weathering state is labeled with It can be observed that each weathering cycle causes a low porosity increase, starting around 0.1-0.2% and 4-6 km/s and increases to almost 0.9% and velocities between 1 and 3 km/s. The Blanco Macael shows the biggest porosity increase by heating up to 200 °C (V), while the following cycle (VI) causes a low further porosity increase. The combination of heating and freezing (VI) has the biggest influence on the porosity of the GI and GK, while this effect is less distinct for the Lasa. The strongest impact can be stated for the Großkunzendorf marble; here, every weathering step causes an increase of porosity, much higher than the one before, with a low range within each weathering group. Whereas the Carrara Gioia shows a clear influence of the combination heating and freeze (VI), the antecedent cycles cause relatively small changes between 4.1 and 2.1 km/s and a porosity between 0.32 and 0.59%. The relationship of the compressive wave velocity and the porosity was analyzed by different authors. Köhler 1991 developed a correlation function between velocity and porosity which was based on experimental data, determined by numerous measurements on Carrara marble in different weathering states. Weiss et al., 2001 presented a model hypothesis based on the mathematical approach of O'Connell and Budiansky 1974 and include the principle that an existing porosity is formed by different ellipsoidal cracks. His approach considers possible crack geometries as pore aspect ratio defined by a length relationship. Rüdrich et al. 2013 used both approaches to test a data set of 123 marble samples. Their data ranged between 7.0 and 2.0 km/s with very low porosities based on very flat cracks with an aspect ratio of 0.005, which could be correlated to the approach of O'Connell and Budiansky 1974. Two strongly weathered samples showed a modification of the crack geometry based on corrosion, what leads to a porosity increase, but has no distinct effect on the ultrasonic velocities. In those cases, Rüdrich et al. 2013 received a correlation to the empirical correlation function of Köhler 1991. They concluded that purely mechanical crack formation follows the model approach of O'Connell and Budiansky 1974 and thereby requires a very small porosity increase with a strong velocity decrease. As soon as the fabric is affected by disintegration, solution processes can become active and modify the microcracks. They assume that in those cases, a strong porosity increase correlates to a relatively small velocity reduction. To see how the determined data should be interpreted within this context, Fig. 16 shows the average compressional wave velocities for all analyzed marble varieties within the artificial alteration under dry conditions. The correlation of velocity and porosity after Köhler 1991 is given in red, the theoretical predictions according to the models of O'Connell and Budiansky 1974 are given as black hatched lines. According to the results of Rüdrich 2003, Rüdrich et al. 2013, and Weiss et al. 2001, all samples follow the model approach of O'Connell and Budiansky 1974 with an aspect ratio of 0.005. However, the Wunsiedel marble, with a partially high weathering resistance, follows the porosity-velocity function of Köhler 1991. The biggest difference between those marbles is that the Wunsiedel shows a relatively low weathering sensitivity which is associated to locally occurring macro-cracks. These foliation parallel macro-cracks, causing the high standard deviation (see Table 4), result in a comparable high porosity, combined with a high average velocity. Those results confirm the analyses of Rüdrich et al. 2013 in particular with regard to the mechanical crack formation. The high deterioration degree does not exclusively cause a high porosity, even if the, for example, strongly deteriorated Blanco Macael shows clear signs of total granular disintegration.

Experimental tomography
For an in-depth knowledge with regard to the development of a tomogram throughout an alteration, phantoms made of Blanco Macael, Carrara Gioia, Lasa, Großkunzendorf, and Wunsiedel marble were manufactured. Those five samples were prepared according to the known spatial orientation of the fabric elements and transferred to the reference coordinate system. The phantoms were made out of the same block as all the other specimens to ensure comparable conditions, to determine the weathering behavior by ultrasonic tomography with well-known material properties during the course of artificial weathering (Fig. 17). A methodical recording of the weathering behavior needs constant test conditions. To ensure these conditions, the ultrasonic transducer was fixed on the surface. This required the use of numerous ultrasonic transducers, which have been developed with a frequency of 1 MHz. The transducers were connected to a multiplex which allowed the switching through the different probe positions (Fig. 18). With this equipment, a repeatability under stable conditions was ensured. Figure 19 shows the modeled tomograms of the marble phantoms in their previous state without consideration of the anisotropy and on the right considering the final percentage of anisotropy. It can be observed that for all samples, a certain degree of anisotropy already has to be considered under unweathered conditions. This correlates very well with the data obtained from the systematic measurements of the cubes and spherical samples (see Tables 4 and 5). Only the unweathered sample of Großkunzendorf (Fig. 19), needs no consideration of anisotropy. This decision was made on basis of the distribution patterns for the root-mean-square (RMS) residuals, which are provided by the software GeoTomCG.
The velocities range between 3.8 km/s for the Carrara Gioia and 6.5 km/s for the Wunsiedel marble. The orientation of all tomograms is the same, while the direction of the lowest velocities correlates with the Z-direction. The anisotropy ranges between 0% (Großkunzendorf) and 20% (Wunsiedel), causing relatively homogeneous tomograms. The irregularities that occur at the edges of the tomograms are triggered by the receiver positions and an increased number of rays and are not caused by the marble itself.
While all tomograms show elongated zones of higher velocities tending to the Z-direction, there is one area with lower velocities tending to the X-direction, which can be observed for the Wunsiedel marble. In this area, macroscopically visible healed cracks can be observed, influencing the velocity.
After the marble phantoms were recorded tomographically in their previous state, they were artificially alterated by thermal stress and then measured again. For all marble varieties, a higher degree of anisotropy results from an increasing weathering state, which corresponds with the systematic observations of the other samples. In Figs. 20 and 21, the tomogram evolution is given for the Blanco Macael and the Carrara Gioia. Exemplaric for the other samples, the Blanco Macael is depicted as a case where a high anisotropy (50%) needs to be considered for the tomographic inversion. The Carrara Gioia still has a high anisotropy of 20%, but was depicted to show the effect of the application of higher anisotropies than necessary. Furthermore, it is shown that by taking the anisotropy into account, previously occurring zoning is put into perspective. This is an important aspect that leads to misinterpretations if not considered. In particular, Fig. 20 shows the large transformation for the velocity distribution of the Blanco Macael after recalculation of the weathered tomogram with 50% anisotropy (Fig. 20). The previously occurring low-velocity channels are relativized and bond to a considerably more homogeneous distribution, which also shows higher velocities than in the isotropic representation. A comparable evolution can be observed for the Carrara Gioia with an increasing homogeneity until the necessary amount of anisotropy was chosen with 20%. As soon as this percentage becomes exceeded, a new channel of higher velocities appears (Fig. 21). The shown tomograms demonstrate that the anisotropy has a large influence and definitely needs to be determined and well chosen as otherwise misinterpretations are possible.

Dynamic Young's modulus
The dynamic Young's modulus (Edyn) was determined based on the ultrasonic velocities for the cubes (Table 4) with a frequency of 350 kHz. The Poisson's ratio was calculated and ranged between 0.296 and 0.303, while the shear wave data were determined on basis of the modeled textural velocities. The Edyn varies between 25.2 ± 2.9 kN/mm 2 for the Carrara Gioia and 71 ± 1.6 kN/mm 2 for the Großkunzendorf in the unweathered state (Fig. 22). Each weathering step causes a clear reduction of the Young's modulus and decreases to values between 6 ± 0.3 kN/mm 2 for the Carrara Gioia and 25 ± 4.6 kN/mm 2 for the Wunsiedel. The highest reduction can be observed for the Blanco Macael with 85.91% and the lowest for the Wunsiedel marble with 62.3% and is for all samples almost linear. As the dynamic Young's modulus describes the elasticity of a material and its resistance to deformation, all marble types show a strong decrease within this property.  . 19 Tomograms investigated for the dummies in the unweathered state. The isotropic tomograms are given in the left column, while the anisotropic tomograms are shown in the right column by using the GeoTom code (GeoTomLLC 2020). The amount of anisotropy is given below the tomograms. Velocities are given in km/s

Thermal dilatation
The specific weathering behavior of marble is based on the thermal properties of calcite and dolomite. Main influencing parameter is the lattice preferred orientation where, according to a pronounced texture, the single crystal properties influence the properties of the whole rock. The calcite single crystal properties are not only anisotropic regarding their ultrasonic wave velocities, but also for the thermal dilatation coefficient α. While heating leads to an expansion parallel to the c-axis, it causes contraction parallel to the a-axes. A strong texture leads to an expansion of the marble sample according to the preferred orientation of the c-axis and a less intense expansion or contraction according to the preferred orientation of the a-axes. As the expansion and contraction affect the grain-to-grain contacts even within a small temperature range (Battaglia et al. 1993) and may cause an irreversible length change (Luque et al. 2010;Rüdrich et al. 2002;Rüdrich 2003;Siegesmund et al. 2000;Tschegg et al. 1999;Zeisig et al. 2002), the application of thermal cycles is a good instrument to determine the weathering vulnerability of marble. The thermal dilatation, ε (mm/m), as a function of temperature describes the expansion behavior during thermal exposure and was analyzed simultaneously for six specimens using a pushrod dilatometer. For these investigations, six specimens in X-, Y-and Z-direction as well as in XY-, XZ-, and YZdirection, according to the reference system, were measured for each marble type. The samples were heated from 20 to 90 °C and subsequently cooled down to 20 °C. To ensure the equilibration of the samples, heating and cooling was performed with a velocity of 1 °C/min, while the maximum and minimum temperature was hold for 6 h. This cycle was performed for four times and used to calculate the thermal dilatation coefficient α. The results are plotted in Fig. 23 and shown for all investigated marble samples in context with their texture. It can be stated that the thermal expansion for all marble varieties shows a clear directional dependency according to the crystallographic preferred orientation, represented by the pole figures of the c-and a-axis distribution. Each sample shows the highest expansion in Z-direction, perpendicular to the foliation plane, which correlates with the direction of the highest intensity of the c-axis maxima, while the lowest expansions can be found along the a-axes in X-direction. All marbles show a more or less pronounced anisotropy of expansion in relation to the spatial position of the test specimens taken, which is primarily reflected in the main directions (X, Y, and Z). The minor directions show a significantly less pronounced anisotropy, depending on how close the specimen location is to the respective aand c-axis maxima. The Wunsiedel marble shows the most obvious variation, with different expansion rates in all six sample directions. Basically, in the first cycle, a successive increase of the slope until the maximum heating temperature can be observed for all samples and directions. Within this cycle, the slope of the hysteresis curve is flatter until a Fig. 23 Thermal dilatation ε (mm/m) shown as a loop of hysteresis with directional dependence (left graphs and middle) and pole figures of the lattice preferred orientation for c-and a-axis (right plots), for the specimens temperature of 40-60°. After this part, the curve starts to become steeper. As soon as the cooling phase is started, the curve declines almost linear until lower temperatures of about 40 °C have been achieved. In this temperature range, the temperature continues to decrease, while the contraction does not progress. This phenomenon can be clearly observed for the Blanco Macael and the Carrara Gioia and is less pronounced for the Großkunzendorf and the Lasa marble. The Wunsiedel shows the same phenomenon, while it is only weakly pronounced. The latter explain the flatter increase in a lower temperature range with a buffering effect due to already existing microcracks. With increasing temperature, this effect declines, and the coefficient increases. The result of this cycle is an irreversible length change; therefore, this phase seems to be the crack initializing. For the second-tofourth cycle, the slope of the hysteresis curve is more linear, and the expansion and contraction shows almost the same temperature-dilatation rate. According to this, the first cycle shows the highest α-coefficients (Table 7), defined by the gradient of the slope, whereas they are lower for the second cycle and remain approximately constant for the following cycle. The smallest value in the first cycle was determined for the Wunsiedel (3.82 × 10 −6 K −1 in X-direction) and the highest for Blanco Macael (16.23 × 10 −6 K −1 in Z-direction). Cycles 2-4 show a range between 2.87 × 10 −6 K −1 (X-direction Wunsiedel) and 14.2 × 10 −6 K −1 (Z-direction) determined for the Blanco Macael and Carrara Gioia. The α-coefficients of the other directions are somewhere between the values of the Z-and X-directions, which can also be observed by the expansion intensity. The variation intensity in every direction for one sample represents an anisotropic rock property and shows with Aα 77.81 % [Aα= αmax−αmin)/αmax * 100 (%)] the highest anisotropy for the Blanco Macael, whereas the lowest anisotropy was measured for the Carrara Gioia with Aα 40.93 % ( Table 7).
The directional dependence of the thermal dilatation is strongly influenced by the texture and the resulting amount of anisotropy can be correlated with the texture (see Fig. 24). As stronger the intensity of the c-axis maximum, as higher the anisotropy. These findings agree with the results by Koch 2006, who analyzed an almost linear trend of the function anisotropy/texture. While there is a correlation between the amount of anisotropy and the texture, the amount of maximum expansion is independent from the texture intensity.
A particularly weathering relevant aspect of the graphics in Fig. 23 is the measure of the residual strain, the irreversible length change after the heating cycle. Here, all samples show a certain measure of the residual strain (εRS), whereby the Wunsiedel marble shows the lowest residual strain and the Großkunzendorf the highest. As described before, the first cycle leads to an irreversible length change, which hardly increases under dry conditions. Comparable data were obtained by Koch and Siegesmund 2004, who . 24 Correlation of the anisotropy for the thermal expansion and the texture also described that the same cycles under hygrothermal conditions can, however, reactivate and increase this irreversible expansion. Therefore, the samples were subsequently thermohygrically weathered for five further cycles. For this purpose, the samples were again cyclically heated from 20 to 90 °C and cooled down under implementation of water into the system at 20 °C (Fig. 25) shows the results of the irreversible residual strain (εRS) under thermal (red highlighted) and thermohygric (blue highlighted) load for all marbles in the six directions. The graphs show that under dry conditions (highlighted in red), the greatest increase of εRS occurs in the first heating cycle and is then significantly lower or does not further increase. An exception is the Großkunzendorf, whose curve progression potentially indicates a further increase of εRS. In contrast, the Wunsiedel shows a very low degree of residual strain and also exclusively after the first heating cycle. Afterwards, no further increase can be observed. All varieties show an anisotropy of residual strain under dry conditions, which is most clearly developed for the Großkunzendorf and Blanco Macael. The thermohygric load (highlighted in blue) leads for all varieties to a further increase of the residual strain. It reaches a maximum value of 0.56 mm/m for the Blanco Macael in the Z-direction and at 0.1 mm/m in the X-direction of the Wunsiedel marble. The increase in residual strain also seems to have not yet been completed for the varieties Großkunzendorf, Blanco Macael, and Lasa, and therefore, higher values of εRS can probably be reached. Apart from the varieties Gioia and Wunsiedel, all marbles show the highest residual strain in Z-direction and a more (Großkunzendorf, Blanco Macael) or less pronounced directional dependence (Gioia, Lasa). Since the residual strain is also influenced by other parameters, such as grain shape and grain boundaries, a simple correlation of the values with the results of the texture analysis is not possible. According to the four types of thermal dilatation (Siegesmund and Dürrast 2014), all tested samples belong to Type IV with an anisotropic α value and residual strain εRS.
Thermal degradation in the form of microcracking results in a volumetric expansion of a stone sample. Indeed, bowing of a stone façade panel most likely results from asymmetric microcracking on one side of the panel. Accordingly, microcracking of a stone sample upon heating or cooling will result in a permanent volumetric expansion of a stone sample when it returns to its ambient temperature. This relative volume change as measured by thermal dilation is thus a measure or index of stone degradation. Thus, the overall objective of this approach is to develop a volumetric extension index (a deterioration index) for the characterization of the degradation potential. The index will reflect the relative volume change on the basis of the knowledge of thermal dilatation.

Determination of volume expansion index
Quantitatively, a continuous surface can be described by the following series of Shushakova (2014): where summation over repeated indexes is implied (the Einstein summation convention). The central symmetry implies that only even term will be present. Accordingly, up to fourth-order terms, the continuous surface is described by The matrix a (2) ij has to be symmetric and therefore, given by six independent elements, which can be found if at least six values of F(x 1 , x 2 , x 3 ) are known from the dilatation experiments. In general, the second-order surface describing the shape change is complex surface, but as seen below, is well approximated as an ellipsoid.
In most practical instances, the deformation surface for marble is an ellipsoid to very high accuracy due to the axial symmetry of the deformation with respect to the crystallographic c-axis. Accordingly, the dilatation experiment data are analyzed by the ellipsoidal approximation. The symmetrical second rank tensor a (2) ij is given by the following matrix: Its values can be found from the following set of equations: Thus, the matrix a (2) ij is defined as and a (2) 12 = a (2) 21 , a (2) 13 = a (2) 31 , where x , y , z , xy , xz , yz are values of residual strain measured in these directions (see Fig. 26).
When the matrix elements are found, it is possible to determine matrix eigenvalues and eigenvectors, by solving the equation: Each 3×3 matrix has 3 eigenvectors and 3 corresponding eigenvalues. The relative volume change of the finite deformation is given by the three eigenvalues according to For infinitesimal deformations, the volume change can be approximated by This value of the relative volume change represents a decay index. A larger value of this index corresponds to greater instability of the material and to its consequent failure upon heating or cooling.
The eigenvalues are also useful in another very important sense, and notably, the smallest eigenvalue is equal to the smallest dilatational strain and the largest eigenvalue is equal to the largest dilatational strain. Additionally, the directions (6) a (2) ij y i = j y j . (7)

Fig. 26
The choice of cut directions for thermal expansion measurements. X-direction corresponds to lineation and Z-direction is perpendicular to a foliation of sample of corresponding eigenvectors correspond to the directions of minimal and maximal dilatation, respectively. As will be seen below, the samples directions do not necessarily coincide with the directions of minimal and maximal dilatation. The sample choices for the thermal expansion experiments are illustrated in Fig. 26.
To verify the approach of decay, three marble cubes were cut from the same piece of Carrara marble (See Fig  27, Table 8). From each cube, six samples (X-, Y-, Z-, XZ-, XY-, and YZ-direction) were cut for the thermal dilatation measurements.

Results and discussion: verification of the approach
Five repetitive heating-cooling cycles under dry conditions and afterwards ten cycles under water-saturated conditions were carried out (see Koch and Siegesmund 2004). The thermal dilatation is plotted for each of three cubes in Fig. 28. The directional dependence under dry conditions is not pronounced, i.e., they show more or less an isotropic thermal behavior. The largest residual strain of 0.17 mm/m is observed in YZ-direction for cube 1, XZ-direction for cube 2 exhibits 0.14 mm/m expansion, and the thermal dilatation of 0.11 mm/m is observed in XZ-direction for cube 3. Under wet conditions, the thermal behavior has totally changed (see Fig. 28).
From Fig 28, the effect of the heterogeneity is evident. Accordingly, the total relative volume change is proposed. The results under dry and wet conditions are presented in Table 8. The decay index is comparable for the three cubes after the five dry cycles, while a difference in values is observed after wet cycles. Such a difference could be caused by the heterogeneity, i.e. influence preferred orientation and pre-existing microcracks on thermal response of marble. The presence of water has significant impact, as well (Sage et al. 1988, Siegesmund et al. 2000. To compare the direction of the maximum and minimum thermal dilatation and its coincidence after the thermal treatment under wet conditions for the three cubes, the pole figure projections of the second rank tensors were plotted and are shown in Fig. 29. All three samples show large values for the relative volume change (Table 8) and the quite same direction of maximum and minimum dilatation (Fig. 29). In sum, this marble type (Carrara) has a high risk associated with failure while exposed under higher temperature and moisture. The studies of Weiss et al. 2002Weiss et al. , 2003Weiss et al. , 2004Koch and Siegesmund 2004;Siegesmund et al. 2008;Shushakova et al. 2011;2012;Weiss et a. 2004;Siegesmund & Dahms 2014 clearly demonstrated that this Carrara marble exhibits rapidly deterioration phenomena. Thus, the decay index is a good quantification tool to predict the thermal behavior of marble, and such an approach was applied for the characterization of the investigated marbles.

Decay indexes for the investigated marbles
Six independent directions for each sample (Fig. 26) were cut for the thermal expansion measurements. Four cycles under dry conditions and seven cycles under wet conditions were carried out. The resulted residual strains in X-, Y-, and Z-directions after the dry and wet cycles for the 12 marbles are listed in Tables 9 and 10, respectively. Decay indexes were calculated for each of the 12 marble samples on the basis of thermal expansion measurements after the dry and wet cycles. The results are given in Table 11.
Decay indexes are quite different under dry and watersaturated conditions. After the dry cycles, marbles such as Arabescato Altissimo, Sölk, Thassos, and Wachau show small deterioration indexes (less than 0.01 % volume expansion). A larger value of the decay index is observed for a majority of the samples after the wet cycles. A significant increase is evidenced for the Carrara (Arabescato Altissimo) marble (from 0.05 to 3.46). Only the Wachau marble does not exhibit any significant change between dry and wet conditions (from 0.06 to 0.09).  To elucidate maximum and minimum thermal dilatation, plots of residual strain tensor after dry and wet conditions were created for each marble sample. Two examples, the Arabescato Altissimo (AA) marble and the Sölk (SK) marble, are illustrated in Fig. 30. The direction of maximum and minimum thermal expansion changes with the presence of moisture. The AA marble with the weak texture exhibits a significant directional change, while the SK marble with strong texture shows only a slight difference. Some of the marbles with a strong texture exhibit a directional change, as well. Eigenvectors corresponding to the largest and the smallest dilatation could be found and defined in the reference frame of the marble.

Summary and conclusions
Although marble has a very simple mineralogical composition, i.e., calcite or dolomite as rock-forming minerals, the weathering features appear in a large variety of forms. In the case of the physical weathering of marbles, the phenomena may range from superficial disintegration to complete loss of cohesion along grain boundaries due to dilatancy, i.e., the total decay of the material. The decay of marbles is sometimes the most spectacular deterioration feature, although the loss in cohesion and the loss in strength are already evident without the well-known bowing of marble panels. Some of the mechanisms that have been suggested to explain the decay of marble include marble type and its fabrics, stress relief, thermal expansion, and expansion due to moisture and temperature (thermohygric properties) or the building physics and the chemical and or biological action. Some following conclusions can be drawn: 1. The investigations have shown that marble decay causes a significant decrease in the strength properties (Young's modulus, etc.), while the porosity as well as the pore radii distribution change. The decrease of the mechanical properties with the increase in deterioration is significant but non-linear. Strength loss versus the increase in porosity relationship is representative for the investigated marbles. 2. A conclusive result was obtained by combining the effect of strength loss and its relationship to the porosity and the compressional wave velocities. Several samples representing a broad variation in decay phe- Fig. 28 The residual strain as a function of the numbers of heating-cooling cycles (20-90-20 °C) under dry (from 1 to 5 cycles) and water-saturated (from 6 to 15 cycle) conditions for the three cubes in X-,Y-, Z-, XY-YZ-, and XZ-directions nomena were used to constrain the loss in cohesion of this sample because the Vp velocities can be applied to characterize the Young's modulus of the marble, and thus, the derivation of the materials quality or their heterogeneities, respectively.

Fig. 29
Plots of residual strain tensors after ten wet cycles: a cube 1, b cube 2, and c cube 3 in the reference frame of sample 2 3. Ultrasonic wave velocities (Vp) is also a measure for the decrease in strength. The reduction of the velocities with increasing state of weathering is correlated with an increasing microcrack density and microcrack width. In addition, the crack type and their distribution in the marble have a significant control on the velocity reduction. The microcrack formation is also crystallographically controlled, i.e., parallel to the cleavage planes or twins. The loading direction for example may control the initiation of intragranular cracks, while in other cases, the shape fabric controls the opening of grain boundaries. Therefore, coarse-grained marbles with interlocking grain boundaries may show completely different microcrack patterns compared to finegrained marbles with polygonal grain boundaries. Consequently, the Vp velocities may differ significantly. If any lattice preferred orientation of the calcite and dolomite crystals exists, the directional dependence of the Vp velocities must be considered when discussing the state of weathering. 4. Ultrasonic wave velocity measurements are a powerful and sensitive tool for the damage assessment of marble. Experimental data reveal that the state of preservation of a marble is clearly documented by compressional wave velocities. For the maximum porosity of around 2.5%, velocities determined on dry samples range on the order of 1.5 km/s. Model calculations reveal that the velocity reduction is caused by cracks with an extreme aspect ratio of about 0.005 or even less. 5. Water saturation has an important influence on the magnitude and directional dependence of ultrasonic velocities. Hence, it is essential to gather sufficient information on the state of water saturation of an object made from marble to quantify the state of deterioration, since water as a pore fluid significantly increases the velocities. This observation is a key constrain when on-site inspections are performed. 6. For the marble, a close relationship between quantified rock fabrics and the mechanical properties was also demonstrated as a function of an increasing state of weathering. In summary, the mechanical and physical properties decrease with weathering, while the porosity and the pore radii distribution increase. 7. The modal composition is an important constrain for the thermal properties of a marble. Freeze-Thaw cycles as well as heating cycles clearly control the weathering of marbles documented by a drastic change in the Youngs' modulus. 8. Based on a systematic study, a quantitative measure as a deterioration index was constrained and introduced on the basis of thermal expansion measurements. This index will reflect the relative volume change. 9. On the basis of six independent measured thermal expansion directions, it was possible to find the values of the smallest and largest possible dilatations. This is also of practical importance if a special cutting direction can help to overcome or hamper the bowing behavior of marbles. 10. Deterioration index was investigated respective to grain size, texture, and other fabric parameters. The tendency how deterioration index varies with changing properties of marble was found. 11. Assess the approach of describing marble's stability using deterioration index on the basis of studies of durability of the investigated marbles, while the loss in cohesion along microcracks and grain boundaries will reduce the mechanical and elastic properties of the rocks. 12. The quantified deterioration index seems to be a very important tool to define risk numbers when for example the overall situation of outdoor sculptures should be forecasted. In this case, the question of the application of winter shelters can be helpful to overcome the missing basic data by modeling approaches in terms of the decay index. 13. Observations of climate change over the past hundred years have induced the research sector to invest greater resources in forecasting its future evolution and impacts. Efforts to focus on improving knowledge of future climate effects on the risk assessment of highranking cultural heritage by quantifying damage processes are of great importance.