1. Introduction
In recent decades, microclimate analysis has played a key role in the field of preventive conservation, as it allows researchers to investigate which environmental conditions are responsible for the deterioration that was observed on movable and immovable cultural objects. Indeed, understanding the interactions between artwork and surrounding environment means finding solutions for appropriate preventive conservation.
Cultural and natural hypogeal sites constitute a unique and valuable part of world heritage. A list of subterranean UNESCO sites—grouped by cultural, natural, and mixed heritage—is provided by Trofimova and Trofimov [
1]. The term hypogeum refers to different types of subterranean sites: from natural (e.g., caves) to man-made sites (e.g., catacombs, tombs, crypts, foundations, etc.), built above ground and covered by soil or dug into the rocks. Hypogeal sites are characterised by stable temperature and high relative humidity throughout the year. On a yearly basis the indoor mean temperature reflects the outdoor one, whereas the outdoor daily thermal cycle extinguishes a few decimetres below ground level and the seasonal thermal cycle at around 10 m below ground level. These features are affected by the combined effect of the homogenous heat capacity of the envelope and limited indoor-outdoor air exchanges [
2]. Other features of hypogeal sites are low or absent natural light and the abundance of organic and inorganic nutrients as well as soluble salts, which could be transported by percolating water from seepage and soil or by visitors. The stable environmental conditions characterising hypogeal sites are usually not harmful for preservation, but it is fundamental to be aware that this fragile equilibrium could be broken by sudden fresh air infiltrations or exposure to natural/artificial light.
In the last twenty years, several studies have investigated the indoor climate of natural and/or cultural hypogeal sites. A list of recent publications on this topic is given in
Table 1. These works focus on the analysis of measured and/or simulated hygrothermal variables, sometimes along with other environmental data (e.g., irradiance, oxygen, or carbon dioxide concentrations). Some of these studies focus on identifying and preventing the causes of deterioration, e.g., biological and/or crystallisation risks, while others consider the effect of visitors on the sites or, in the case of hypogeal settlements, their energy performance. Nevertheless, none of the listed papers provide a thorough and straightforward approach for the characterisation of microclimate-induced conservation risks.
The crystallisation of salts is one of the frequently occurring factors of the deterioration of porous materials in hypogeal sites [
7,
8,
9]. It is visible at the macro-scale as efflorescences and sub-efflorescences progressing towards the loss of material from superficial to inner layers in terms of powdering, detachments, blistering, exfoliation, and flaking. Crystallisation is influenced by the evaporation of saline solutions, which strongly promotes crystal growth from large to small pores, thereby inducing stresses in the material micro-structure. Several mechanisms are involved, such as crystallisation pressure, thermal expansion, osmotic pressure, and chemical weathering. Common deliquescent salts are sulfates, chlorides, and nitrates, which can often be found in association with cement that is used for restoration interventions [
15], percolation of rainwater from the surrounding soils [
16,
17] and wet or deposition of airborne particles on surfaces [
18]. Price at al. [
19] and López-Arce et al. [
18] proposed strategies to control rising damp and, hence dissolution, mobilisation, and crystallisation of salts. Nevertheless, the natural ventilation of the site should be limited in order to reduce moisture evaporation from walls governing dissolution/precipitation cycles.
The proliferation of biodeteriogens is another detrimental factor for preservation within hypogeal sites [
3,
4,
5]. The presence of biodeteriogens is promoted by the availability of liquid water on surfaces in conjunction with nutrients, natural/artificial light and poor indoor ventilation. Biodeteriogens can form a complex mucous matrix, in which they are embedded, called
biofilm. Biofilms contribute remarkably to rock and building material decay and, depending on the type of biodeteriogen, they usually appear as coloured patinas, which may exert a detrimental aesthetic impact. In Roman hypogeal sites, the main agents that are responsible for biodeterioration are cyanobacteria, i.e., photosynthetic microbes pioneering biofilms by forming a resistant extracellular polymeric substrate. Once the biofilms are formed, heterotrophic biodeteriogen, such as actinobacteria and fungi, can colonise the surface by exploiting the metabolites and biomass that are produced by cyanobacteria [
20,
21]. Bruno et al. [
22] employed an LED lighting system to limit the growth of biofilm by inhibiting photosynthetic activity. Natural ventilation of the site should be favoured in order to reduce spore concentration in the air as well as the moisture evaporation from walls.
This paper presents a comprehensive study that is devoted to the assessment of the microclimate-induced conservation risks inside hypogeal sites. The methodology is applied to the “Mithraeum of the Baths of Caracalla”, a Roman hypogeal site, which suffered from different decay processes and where a long-term microclimate monitoring campaign is still in operation. The case study and the research aims are described in the following paragraphs.
1.1. The Case Study
The “Mithraeum of the Baths of Caracalla” (hereafter called Mithraeum) is one of the most important sites of Mithraism in Rome, as it preserves a valuable fresco of Mithras in the Central Hall. It was built in the early 3rd century AD in an underground corridor of the north-west exedra of the Roman Baths. The Mithraeum, at a depth of 10–15 m, covers a total area of 600 m and has five communicating rooms, where the Central Hall was the main room of the Mithraism cult. The site was built in opus testaceum made of two external parallel terracotta bricks (bipedus spezzati) filled with a layer of lime mortar and inhomogeneous aggregate (earth, tuff, and smashed bricks) for a wall total thickness of 1.5 m.
Several deterioration processes were observed in 2010 in the Central Hall: exfoliation, flaking, and swelling due to salt weathering and green/dark patinas caused by biological colonisation. A microclimate monitoring campaign was started in December 2010 to assess the effect of microclimate on deterioration processes and, thus, to seek solutions for improving site conservation. A maintenance intervention was carried out in 2011 to reduce internal ventilation with the aim of inhibiting efflorescence through the closure of the main openings, i.e., the open lateral windows (“bocche di lupo”) and the roof covering of an adjoining room. The adopted measure visibly reduced the occurrences of efflorescence and sub-efflorescence due to salts crystallisation, but did not inhibit the biological proliferation. Since then, the microclimate monitoring has been ongoing, although not consecutively, over all months of each year, due to a lack of funds and technical issues.
1.2. Research Aims
The aim of this research is to outline a comprehensive approach for the assessment of the microclimate-induced deterioration in cultural hypogeal sites through the definition of a Risk Index. The Risk Index (RI) is targeted at allowing a fast and reliable evaluation of the risk of crystallisation/dissolution of salts and bio-colonisation. The availability of the long-term microclimate data collected in the “Mithraeum of the Baths of Caracalla” provides the opportunity to test the ability of the RI to show the synergism of frequently occurring hazards inside hypogeal sites.
3. Results and Discussion
3.1. Temporal Behaviour of the Microclimate during the Year
The hygrothermal conditions in 2016 are reported in
Figure 3 as being representative of the annual microclimate temporal behaviour inside the Mithraeum. This year was chosen, as it was associated to the most complete (CoI = 1) and continuous (CI = 1) time series that were collected among the years.
On a yearly basis, the indoor mean temperature significantly agreed with the outdoor one: T
= 14.7
C and T
= 15.5
C (
Figure 3a). T
was smoothed out with respect to T
on both a seasonal and daily basis, respectively, showing a time-lag of one month in the T
maximum peak due to conductive heat transfer through soil (high thermal inertia) and a time-lag of few days due to convective heat transfer. The RH values were close to saturation most of the year, ranging between 80% and 100% (
Figure 3b). The median of RH
was 95.5% and the median of RH
= 78.1%, with drops in RH
being related to simultaneous drops in RH
.
3.2. General Microclimate Analysis
The quality of T and RH time series were objectively evaluated before performing the comparison among microclimate data that were collected in three years: 2011 as representative of the microclimate conditions before the intervention, 2013 as representative of the period immediately following the intervention, and 2019 as representative of the most recent microclimate conditions. The preliminary investigation on the microclimate conditions inside the Mithraeum was limited to the time intervals available in all three years from 1st February until 30th April. In 2011, the environmental monitoring was temporarily suspended due to the starting of the maintenance intervention. For the selected periods, the CoI and CI values are close to unity, which means that the time series were both complete (few missing data) and continuous (few intervals with missing data), hence being suitable for studying the hygrothermal behaviour of the site.
The first comparison was carried out investigating the influence of the external weather on the indoor hygrothermal conditions. Looking at the linear regression equations reported on
Figure 4a, the angular coefficient (
m) of T in 2011 was close to unity, while
m was halved in 2013 and 2019. Similarly, the value of
m of the MR linear regressions (
Figure 4b) dropped from 0.9 in 2011 to 0.5 in 2019, which might imply a progressive wetness of the site. This result can be interpreted as a direct consequence of the closure of the openings, as the indoor hygrothermal conditions tended to be less dependent on the outdoor values.
The reduced variability found in T and MR in ambient air since 2013 (
Figure 4) was confirmed by a substantial drop in the intensity of air currents (V) inside the Mithraeum (
Figure 5). In 2011, the V values were well scattered, although a prevailing direction was noticeable, which ranged between 40
and 45
, coming from one of the windows, which was also characterised by the V maximum of 2.3 m/s. In addition to the maintenance intervention, the hygrothermal conditions from October 2012 to January 2013 were affected by temporarily openings to visits. After the reduction of these forcing factors (i.e., natural ventilation and visitors), the dependency on the internal climate on the outdoor weather conditions was visibly reduced in the years 2013 and 2019, with data in the polar plot of air intensity less (2013–2019) scattered in the polar plot of V (
Figure 5) and V always lower than 0.5 m/s.
The temporal behaviour of the hygrothermal conditions measured a few centimeters from the wall (sensor S
in
Figure 1) was analysed while taking the observations of surface temperature (T
), air temperature (T), and mixing ratio (MR) values into account.
The box-and-whiskers plots of T
, T, and MR observations (
Figure 6a) shows that the boxes fully overlapped over the years, meaning that no significant difference was detectable after the closure of the main openings. Indeed, the mean of the medians over the three years was 9.8
C, ranging from 9.4 to 10.2
C for T
, 9.7
C ranging from 9.5 to 9.9
C for T, and 6.9 g/kg ranging from 6.1 to 7.4 g/kg for MR. Conversely, the Interquartile Range (IQR, i.e., the difference between the third and first quartile) shows a decreasing trend year by year: from 10 to 5
C for T
, from 13 to 6
C for T, and from 8 to 4 g/kg for MR.
This effect is also evident in
Figure 6b when comparing
T
and
T
before and after the retrofit. It is important to bear in mind that both of the parameters are based on isobaric cooling until saturation is reached, although T
is related to cooling without changes in MR, whereas T
is linked to cooling due to the evaporation of water. As the values collapsed to the condition where
T
=
T
= 0
C, the risk of condensation associated with both cooling processes greatly increased after the intervention.
Two empirical functions were used to study the temporal behaviour of the daily evaporation rate (E
, Equation (
2)), and the amount of condensed water (CW, Equation (4)) on the walls. As shown in
Figure 7a, in 2011 E
reached up to 0.74 mm/d, whereas in both 2013 and 2019 it was always below 0.10 mm/d, decreasing on average from 0.24 mm/d before the maintenance intervention to 0.02 mm/d after the maintenance intervention.
Figure 7b shows that most of the occurrences of CW with values lower than 1.0 mg/cm
occurred before the intervention (2011), whereas, in 2013 and 2019, CW was higher than 0.5 mg/cm
, reaching values up to 30.0 mg/cm
in 2019. Moreover, the average CW tended to increase, at first, from 0.6 mg/cm
in 2011 to 1.8 mg/cm
in 2013, and then strongly to 11.1 mg/cm
in 2019. Furthermore, in 2019, the CW values tended to raise over time, as P was always close to unity, whereas T
− T
is close to zero.
3.3. Specific Risk Assessment
The conservation risk assessment was performed on the time series of the years 2011, 2013, and 2019, while taking the hygrothermal conditions in proximity to the walls into account (sensor S
in
Figure 1). Mechanical and biological risk were investigated, as described in
Section 2. It should be noticed that, in 2011, data were only collected from January to the end of April, as the maintenance intervention started in May.
The risk related to crystallisation was estimated by plotting T
and RH
on phase diagrams showing the critical RH for the deliquescence salts ordered by increasing harmfulness, i.e., chlorides, nitrates, and sulfates (
Figure 8). In 2011 (first row of the 3 × 3 matrix), RH
values spanned from 40% to saturation, which means that chlorides, nitrates, and sulfates could be present in both solid and dissolved forms. These conditions could be associated with high risk of crystallisation for all of the deliquescent salts considered, except for antarcticite. The risk associated with sodium chloride (halite) and potassium chloride (sylvite) could have produced negligible damage to the stone, limited to aesthetically-disturbing efflorescences on the surface [
32]. Concerning mirabilite, the main conservation risk could be related to the formation of sub-efflorescence within the stone pores [
35], threatening the conservation of porous materials because of its high crystallisation pressure [
15].
After the intervention, the RH values always exceeded 60% due to the reduced ventilation in the Mithraeum. In 2013 (second row of the matrix), nitromagnesite and nitrocalcite could have only existed as dissolved ions, thus excluding the risk of mechanical damage caused by their precipitation. In 2019 (third row of the matrix), the RH values were always higher than 80% and only sylvite, nitrokalite, and the sulfates might have formed crystals. Nevertheless, the risk associated with the crystallisation of mirabilite might still remain high, as mechanical damage is governed by the frequency of the alternation between dissolution and precipitation.
The risk of bio-colonisation by mould fungi was assessed by plotting the collected hygrothermal data on Sedlbauer isopleths describing spore germination (
Figure 9, left panels) and mycelial growth (
Figure 9, right panels) on substrates characterised by biologically adverse recyclable building materials with considerable soiling (Substrates category I). The biological risk greatly increased year by year, especially during winter and spring, as most data lied considerably beyond the LIM (Lowest Isopleth for Mould). Moreover, in summer and autumn of 2019 the risk of spore germination rose to eight days, whereas the mycelial growth is always higher than 1 mm/d if compared to data that were collected in 2013.
3.4. Application of the Risk Index
The Risk Index (RI) for hypogeal sites was calculated on a yearly basis (
Figure 10, upper panels) and in each season (
Figure 10, lower panels) to synthetically characterise the biological proliferation and crystallisation/dissolution risks inside the Mithraeum.
In 2011 (
Figure 10, left panels), the hygrothermal conditions promoted both the conservation risks, with RI ranging between 20% (nitrates) and 70% (germination and growth of spores) on a yearly basis. Deterioration seemed to be mainly governed by the air moisture content close to the wall surface. In spring, the increase in moisture content determined an enhanced biological risk and phase changes of sulfates, as well as an inhibition of phase changes of chlorides; on the contrary, nitrates did not seem to be affected by the moisture seasonal cycle. In 2013 and 2019, the biological risk significantly increased with respect to the condition in 2011 (RI = 100%), as the closure of the main openings limited air movements (
Figure 5) and RH higher than 80% most of the time. On the contrary, the crystallisation risk diminished following a seasonal cycle where the most hazardous period was winter. In 2013 (
Figure 10, mid panels), the risk of phase change of chlorides was reduced as compared to 2011 and limited to winter, when the site was still open to visits. Conversely, the risk of phase change of nitrates and sulfates increased, especially in winter and spring, showing a seasonal cycle with minimum in summer. This finding might be caused by the increased evaporation that occurred in winter 2013, as shown in
Figure 7a. In this respect, crystallisation/dissolution risk seemed to be governed by temperature. In 2019 (
Figure 10, right panels), the biological risk still remained high throughout the year (RI = 100%). As in 2013, the risk that was related to chlorides was low and limited to winter and spring. The risk of phase change of nitrates strongly decreased from 25% in 2013 to 5% in 2019, showing a low seasonal trend. The risk of crystallisation/dissolution alternation of sulfates was comparable with that in 2013, except for spring.
The visualisation of T-RH on scatter plots with critical RH curves provided a meaningful risk assessment, which cannot be fully interpreted without detailed knowledge of the damage potential of individual degradation agents. In contrast, RI made a fast and reliable identification of the risks as well as an assessment of the severity levels possible. Indeed, the advantage of the application of RI was the ability to synthesise the study of long-term time series, providing a global characterisation of deterioration on both a yearly and a seasonal basis. In this perspective, RI might be used by heritage site managers as a practical parameter to effectively compare different individual risks throughout the years and after maintenance interventions. There is some uncertainty regarding the assumptions made in this study, as microclimate-induced deterioration in the Mithraeum could have been different from those most frequently occurring in other Roman hypogeal sites. However, the results can be generally applicable in the assessment of the conservation risks. The experimental verification of crystallisation/dissolution and bio-colonisation risks will be planned in future analytical campaigns to encompass the above limitation.
4. Conclusions
A comprehensive assessment of the microclimate-induced deterioration was performed inside the ”Mithraeum of the Baths of Caracalla” by evaluating the most common conservation risks in hypogeal sites, i.e., crystallisation/dissolution of deliquescent salts and mould fungi proliferation. The analysis of long-term microclimate series provided an accurate characterisation of the site over the years, clearly showing the effect of a maintenance intervention on the hygrothermal conditions. The investigation highlighted a reduced variability of temperature and relative humidity due to the limited internal air flows. The closure of the main openings has reduced the moisture exchanges with the outdoors, leading to an increase in the dampness conditions. The conservation risk close to wall surfaces was assessed on hygrothermal scatter plots with the critical RH curves of three deliquescent salts systems and of mould fungi proliferation. These diagrams provide an evaluation of the crystallisation and biological risks in response to the hygrothermal conditions. Such data visualisation cannot allow easy interpretation, especially in the case of crystallisation risk, as the mechanical deterioration is related not only to crystalline growth, but also to crystallisation and dissolution cycles of each saline solution around its critical RH curve.
The Risk Index (RI) provided a quantitative measure of the individual risks and made it possible to analyse their synergism, supporting the interpretation of diagrams of both deliquescent salts phase change and mould fungi proliferation. A synthetic visualisation of the overall risks was proposed by means of radar and histogram plots of RI, respectively, on a yearly and seasonal basis. The Risk Index proved to be effective in diagnosing and following up on the microclimate-induced risks in the hypogeal sites subjected to specific interventions. Although the Risk Index here was specifically designed for cultural hypogeal sites with long-term time series, RI can be adjusted to fit different aims and requirements. The evaluation of crystallisation/dissolution risk was limited to typical occurrences of deliquescent salts in Roman hypogeal sites, but it could be extended to other salt systems that are associated with different environmental conditions. The same applies to the biological risk, which focuses on spore germination and mould growth, but might also be conveniently extended to cyanobacteria and colonisation by insects.
Some final suggestions for improving the risk assessment procedure for future research perspectives are proposed. First, a database of microclimate-induced deterioration functions for cultural materials could be implemented in view of the application of RI to different environments and various conservation needs. Subsequently, a standalone executable of the RI algorithm can be developed to better support historical site managers and other stakeholders in tailoring preventive conservation strategies.