- Open Access
Microclimate monitoring in the Carcer Tullianum: temporal and spatial correlation and gradients evidenced by multivariate analysis; first campaign
- Giovanni Visco1Email author,
- Susanne Heidi Plattner†1,
- Patrizia Fortini†2,
- Serena Di Giovanni†1 and
- Maria Pia Sammartino†1
© Visco et al; licensee BioMed Central Ltd. 2012
- Published: 2 May 2012
Too often microclimate studies in the field of cultural heritage are published without any or scarce information on sampling design, sensors (type, number, position) and instrument validation. Lacking of this fundamental information does not allow an open discussion in the scientific community. This work aims to be an invitation for a different approach.
Three main parameters (temperature, humidity, luminance) were monitored in a selected part of a complex construction by an inexpensive self-assembled system along some horizontal and vertical vectors. All data was then processed and analyse by chemometric methods. Some measurements of oxygen, carbon monoxide and dioxide and pressure were also performed.
Correlation of some indoor and outdoor data was shown for temperature and humidity. In case of outdoor changes the indoor environment reacted with a certain delay which is position-dependent and more evident for humidity data. The two observed rooms (Carcer and Tullianum) behave differently and the hypogean one is less influenced by the outdoor environment. Instrument validation before and after the campaign, allows to consider detected variations as significant.
The fundamental importance of Sampling Design and of instrument validation before and after the monitoring campaign was enhanced. The choice of two main and two minor vectors allowed detection of different behaviour for the two rooms, also permitting to detect for both rooms a trend towards a spontaneous microclimate necessary for a conservation project. In the next campaign we will focus on the choice of the best sampling frequency to use more sophisticated statistical methods.
- Control Chart
- Barometric Pressure
- Measurement Campaign
- Stone Column
- Matrix Layer
Microclimate never can be univocally defined in space; should one consider the indoor part of a building or restrict it to a single room?
Obviously, correct spatial definition of a microclimate should also consider all surfaces (walls, statues, coins, etc.); their climatic conditions are indeed influenced by that of the surrounding environment, but depend also on the material.
Very often the definition is referred to a building or at least to a partially isolated area.
For a building a ‘mesoclimate’ can be defined, which is influenced by the building itself, its structure, geometry, shape and materials and by the surroundings.
For our purpose we define a macroclimate, characterised by locally measured meteor-climatic values, which are not influenced by the building itself .
In the case of historic buildings, one or more centuries old, both the meso- and the microclimate should have reached an equilibrium condition, that can include small natural variation; however higher change and or fluctuation of climatic parameters might occur due to restorations, excavations, reinforcements (sometimes with improper materials) or to the presence of humans (visitors, restorers) making measurement campaigns more complex. For conservation such fluctuations mean that, for example, a wall structure of stone, ceramics, glass, frescos, metals and also organic material which were conserved for centuries or millennia due to equilibrium conditions are put in danger in a very short time range compared to the sites life-time.
The present study concerns with a hypogean, which is a very particular part of a building and its investigation can be difficult; one reason is that - while identification of the spontaneous condition of the research environment is particularly important and fundamental for a conservation project - it is impossible to define a priori for a hypogean the time necessary to reach constancy or at least spontaneous periodicity of fluctuation, referring to Banham’s definitions .
Once macro- and microclimate are sufficiently defined, one can compare both in order to identify their correlation type, e.g. “Conservative Mode”, where the walls or the burying soil exclude or strongly slow down any exchange; “Selective Mode”, where some parameters are filtered by the barrier while others are not; “Regenerative Mode”, where an independent internal microclimate is created, which instead leads to retroactive mechanisms (typical for modern living environments with high waste of energy). For a hypogean the external construction often functions as filter for the macroclimate and might guarantee good microclimatic and air quality conditions. Two significant parameters, for conservation, are not studied in this research; the humidity distribution inside the walls  and the presence of microorganisms like fungi, bacteria, yeast, moss and also anaerobic organisms , both related to microclimatic conditions and necessary in a conservation project.
But what are good microclimatic conservation conditions? Ideal or limit values of microclimate parameters are defined in some official Italian publications  as well as in scientific ones , but knowledge on these values is still poor. Given that none of these values derive from an accurate (long-lasting?) experimental study with double blind technique as in the case of new pharmaceuticals and instead they are often the result of intensive expert-arguing [personal communication by one component of Normal commission], there is an urgent need for many, accurate, reproducible, comparable works, in order to avoid that one day microclimate parameters are imposed by law . Unfortunately many papers on microclimate studies of indoor environments, especially in the field of cultural heritage conservation, could be discussed critically, because they do not state the type, number or position of the sensors used. Often the Design of Experiment (DoE), the result and instrument validation is not described. Generally few sensors, sometimes only one, are applied and time series obtained, which are then analysed by means of control cards and/or univariate graphs, but surely there is a correlation between temperature and humidity so a multivariate approach is necessary .
In this preliminary study we propose a, we think, sufficiently rigorous procedure, but for a too short observation period, where15 inexpensive but accurate data loggers were used for monitoring three environmental parameters (temperature, humidity, lighting) every 10 minutes along some vectors in a hypogean constituted by two rooms one above the other; chemometry was used to extract the most significant information.
The position of the building in the ancient Forum Romanum is shown in fig.A in additional file 1 (all figures in additional file 1 have a capital letter as index); it was probably built in the V. or VI. century b.C., above a water source.
Basing on the archaeologists maps the dimensions of the two rooms were obtained; for the Tullianum wall thickness of about 1.7 m, an internal surface of ca. 25 m2 and a volume of ca. 56 m3 were found. For the Carcer walls thickness varies from 1.7 m (external) to 1 m (versus other indoor environments), surface is approximately 60 m2 and volume ca. 300 m3.
Luckily, the entire complex is not equipped with a heating, ventilation and air-conditioning system (HVAC), so, to ensure against anthropogenic variables, it was sufficient to close both rooms to public for the entire measurement campaign. The pump used to drain water from the Tullianum was set on automatic mode in order to avoid that the water level raised up some cm so ensuring against sensor diving.
For microclimate monitoring we used 13 Hobo U12-012 by Onset USA for measurements of humidity (RH%), temperature (°C) and Luminance (foot-candles or lumens/ft2).
A EB20-THP, by Ebro GmbH, Germany was also used as equipped with certified sensors for humidity (RH%), temperature (°C) and pressure (mbar).
In the preliminary coherence study (see paragraph Sensors calibration) two VWR cod. 61027-056, certified Hg thermometers AA grade, were used for controlling and calibrating of data loggers.
A LabQuest data station with sensors for CO2, O2 and temperature measurements was used as well (CO2-bta, O2-bta and STS-bta by Vernier Tech., USA).
Two Lascar data loggers (UK) were used for testing (EL-usb-CO for CO in air and EL-usb2 for RH%, °C).
In order to compare micro- and macroclimate, meteorological (macroclimate) data was obtained from "Osservatorio Meteorologico del Collegio Romano", active since 1782 and located in the centre of ancient Rome and on the top of a 63m (above sea level) tower, situated in 41.900°N, 12.480°E. Tullianum instead is situated in 41.893°N, 12.484°E.
Data elaboration software were mainly free ones, e.g. Past, Datalab, Gnumeric, WinIDAMS, Libre Office and some commercial ones like MVSP and Lotus 123.
To read data from data-logger producer software were used (available on brand web site). Only one software was not free and we apologize for this.
Our work focused on the monitoring of two principle parameters, temperature (°C) and relative humidity (RH%).
As above said, an other parameter strongly influencing the conservation state of an object is bound to lighting. In a running research we are considering the effect of both the spectra of the light source and its wavelength-weighted power while in the present study only the last was monitored; illuminance was measured due to its greatest significance in studies regarding frescoes, paintings and any object in general.
15 data loggers were employed as passive samplers to monitor the three above said parameters; they were easily fixed also along vertical vectors (weight 46g) without disturbing the spontaneous environmental air flux.
Comparison of the values of the 3 data loggers, according to the producer
drift, one years
rise time, min
Even if active samplers may be used at a higher frequency with respect to passive ones, we choose the last as, in our opinion, for a long monitoring period the adopted frequency allows to obtain significant information and, overall, the artificial air flux of known flow rate created by a pump and oriented on the active sensor , involves a self-heating especially in restricted indoor environments. Really, the self heating of a thermistor is a problem well known to electrical engineers, even a small 2W pump motor can modify the microclimate and thus measured values, when employed for longer periods.
In this first campaign an experimental data station with some sensors was also used in order to measure gas concentrations (e.g. CO2 and O2) in air. Further an attempt was made to evaluate the CO level as an indicator for contamination of the hypogean with external, urban, air. In order to identify the "spontaneous" values, the two environments were not entered during the 15 days measurement campaign.
In this work the choice of vectors direction and position for sensor placement was based on Sampling Design (SD) using an expert-sampling among all authors for vector position; on each vector a systematic-sampling with even spaced distances including the extremes, was applied.
Position of the two main vectors is vertical (trespassing the two rooms, 6 sensor points) and horizontal (in the basement, 6 sensor points); the minor two vectors are one along a horizontal axis in the upper room (3 sensor points) from the door to the altar and one in the basement (3 sensor points) from the stairs to the opposite wall (fig.1).
The long horizontal vector is designed at about 50 mm from the floor and at the cross of the two vectors, the concentrations of CO, CO2 and O2 were measured too.
In the time domain we used a systematic-sampling with a reading every 10 minutes, but the choice of the proper sampling frequency was based on a preliminary 24h measurement campaign, using the faster available sensor, at 2 readings each minute.
Obviously immediately after the campaign one notices that some additional data logger somewhere would have been useful, leaving space for further campaigns and making work even more complex.
Sensors control before and after measurement campaigns is a minimum requirement which is generally not stated in the microclimate studies on Cultural Heritage Sites.
Especially when inexpensive sensors are used problems could arise due to lacking of coherence of the single measurements; the same data-loggers are often available in a certified version, but as our aim was also the assembling of a inexpensive measurement system, the price of these latter was obviously too high. So, at the same cost we preferred a double number of data points buying more data-loggers and to verify by our own the sensors quality.
To do this, in a preliminary measurement campaign all 15 sensors where placed for a week in a hypogean ambient similar to the one to be studied, near to each other and around two certified instruments (the EB20-THP and VWR Hg thermometer); in this way we checked the sensors rising time, accuracy, hysteresis, and coherence before starting our measurement campaign. This allows to estimate the sensors response to environmental variations, e.g. a coherence of 0.5 °C allows us to regard a variation of 0.8 °C between two rooms or along a vertical axis as significant.
Once the sensors coherence was confirmed they were reprogrammed and positioned in the Tullianum for our study purpose.
At the end of our work all sensors were re-checked in the same way and in the same place and results are shown in additional file 1 as fig.B (temperature) and fig.C (humidity); the graphs show that no damages or calibration lacks occurred to the sensors during the measurement campaign because data are coherent considering the accuracy limit of the data-loggers (table 1). A Vernier data station was also used equipped with temperature, CO2 and O2 sensors (fig.D in additional file 1, CO trend is also reported); temperature data showed an offset of ca.0.5 °C and was therefore excluded even if used by the data station to apply the right correction factor to data coming from the two other sensors.
The obtained data was put in a 3D-matrix composed by 16 columns, with sensors as variables, 2150 rows for each sampling time and 7 layers for temperature, RH, luminance, barometric pressure, CO2, O2, CO. Some layer columns are empty as not all data loggers can detect all parameters. For example the barometric pressure layer has only one column (Ebro sensor).
The above described matrix was then studied by Exploratory Data Analysis in order to identify any correlation of indoor and outdoor measured values as well as to identify the rooms spontaneous values.
Data of each measured parameter (layer) was plotted in monovariate graphs. Normally Run-Sequence plots of parameters (temperature, humidity, light intensity..) versus time or a function of time (here minutes from midnight of the first day).
In these graphs prefixed limits can be added and Shewhart Charts are obtained (also known as Control Charts), which enhance the environments characteristics.
The successive step in EDA is data analysis by Box & Whisker Notched Box Plots of a single parameter (a matrix layer), where variables (sensors) are plotted on the X-axis and raw measured distribution on the Y-axis.
This powerful analytical tool shows, here in the McGill  version, in the median domain, the distribution, spread, skew and percentile values 5% and 95% by whiskers; no outliers are shown. In our case only one matrix layer was plotted at a time so avoiding problems which can arise when variables with different scaling are plotted together.
Successively the three Scatter Plot Matrixes were designed showing up all possible Var-Var graphs between all variables (sensors), starting from raw data, in order to enhance correlations as well as distribution type of each variable on the diagonal.
On the diagonal of the Scatter Plot Matrix a histogram of ten bars is shown for each variable, which roughly allows identification of the distribution form; the values on histogram show the height of the major bar as the total is normalised to 1.
This important graph is often substituted by a correlation matrix, using the r2 Pearson coefficient; this latter works only on linear data, which is not always the case, especially in the cultural heritage sector (please see the complex column of Hobo10, on the altar, on the next figure).
In temperature and humidity layers there is also a column with macroclimate (outdoor) data; correlation of data in this column and those of the others is possible and plausible. However macroclimate data has a 30 minutes sampling interval; therefore a new matrix was constructed in which the indoor data was downsampled in order to allow better comparison of indoor and outdoor data.
Evaluation of the delay between internal and external barometric pressure studying the lag correlation.
parameters \ lag (min)
Pearson product-moment correlation coefficient, sample
square of the Pearson correlation coefficient
assess of Std Err of Pearson's r
variance of post-fit residue
residue sum of squares (RSS), sample
root mean square deviation (RMSD), sample
standard error of the Y estimate
If one plots in graphs parameters which best describe correlation, in terms of quality, here rows enhanced in blue, it is possible to compare two or more correlation, which differ for instance for the data amount (many validation techniques like LeaveOneOut or Bootstrap are based on these descriptors).
A building isolation can be evaluated also by Loss/Gain for Transmission graphs, where internal and external data are shown and compared to a perfect correlation which is shown as graph diagonal. Any sign above the diagonal is a gain of the building compared to the outdoor environment, e.g. can be due to a heating system for a temperature graph. By using labels for single data points these graphs allow also evaluation of hysteresis, spontaneous behaviour, changes with time. The sensors accuracy is fundamental in these graphs.
A temperature-humidity correlation study was the next analysis step. Thus two columns of two overlaying layers were compared and for each point variation of both parameters as well as the 98 % confidence ellipses were obtained.
All these data exploration methods together allow evaluation of the buildings isolation from the outdoor environment, as well as enhancing internal gradients; further spontaneous values, not normal behaviour and sudden changes, which are the most dangerous events for a cultural heritage site, can be evaluated.
These simple graphs already enhance:
– a more or less evident correlation between indoor and outdoor (macroclimate) data.
– a certain delay, more evident for humidity data, for indoor changes with respect to outdoor changes.
– as expected, the two rooms behave differently; the hypogean is less influenced by macroclimate and show a lower variability for both the measured parameters.
– a tendency versus a spontaneous microclimate, more evident for humidity so putting in evidence the need to extend in time the next measurement campaign.
– finally, the coherence among sensors was checked before and after the campaign, (examples of the latter in fig.B and fig.C in additional file 1), allowing to consider all variations shown in graphs of fig.2 and fig.3 as significant.
– careful graph observation allows stating that for these two different rooms the attempt to draw univariate Shewhart or Control Charts and define limit values to be respected by a HVAC would be confusing and stupid.
Study of the values distribution
For remarking of each distribution and mutual behaviour Box&Whisker plots can be used starting from raw data.
In fig.5b, the humidity distribution evidences the same two main data groups well separated; one characterised by wider ranging values and a lower median (points 1, 2, 4 and 10) and another one including all remaining data points. Sensor 6, in this case, is clearly included in the Tullianum group even if, together with sensor 5, located on the stairs connecting the rooms at about 2 m from the lower floor, shows a similar median value enough low with respect to all the others. The presence of the subgroup including sensors 11, 12 and 13 is also better evidenced; such sensors are positioned on the main horizontal vector, inside a tunnel probably coming versus the Cloaca Maxima and connected to the room through a metallic door that was left open during the campaign.
Correlation among points
Analysis of the structure’s isolation from the external environment cannot be performed by a simple scatterplot, instead one has to calculate all correlation parameters  among two columns of the same layer.
Before starting lag-correlation analysis it seems useful to identify time-windows hidden in time-series (repetition of the same values with a certain frequency due to cyclic phenomena such as night-day or seasonal variations). A number of software allow detection, expensive ones and free ones like Past .
Time series analysis are performed in the same way as Spectral Analysis, starting with Fourier Transformation which allows shifting from the time domain to the frequency domain. For temperature and RH data series of sensors 2 (top of Carcer) and 8 (top of Tullianum), a time window of 1430 minutes is detected for the first peak of the FFT spectra.
Usually in combination with spectral analysis Sinusoidal Regression (often called Sinusoidal Curve Fitting, SCF) is used, where a sum of n sinusoids with specified periods is used to fit a known list of time dependent data. Past algorithm of SCF is based on a least-squares criterion and singular value decomposition.
This can be useful to check and model periodic phenomena in time series, as in our case. Further it is interesting to compare the results of FFT and SCF in order to obtain result validation and to avoid overfitting.
Expressing the identified time window in hours, one finds a values close to 24 (1440 minutes), which confirms influence of the external environment on our indoor ones.
To check fast change in values (temp. RH, lux) we must use a sampling frequency more than two times faster than the change under study (Nyquist–Shannon sampling theorem), but this is difficult with a sensor with rise time of 6 or 10 minutes, one other big problem in microclimate monitoring often silenced.
The meteorological data used had a 30 minutes sampling interval, thus our measured microclimate data had to be reduced by downsampling as described in methods section.
For correlation analysis of meteorological data (outdoor, X axis) and measured microclimate data (indoor, Y axis) Ordinary Least Square was used; data of the Y axis was then shifted by 0, 30, 60, 90 and 120 minutes and obtained parameters registered in new matrixes, like the one in table 2.
Correlation coefficients were calculated for all available data without using a time window in order to avoid any anomaly. If number of data is much higher (10-30 times) than the width of time window this techniques is sufficiently robust to avoid dependence from some not normal data (e.g. presence of a small heater on day 12 or a window open all night on day 18).
In contrast to other works our first approach was data validation and results coherence; correlation analysis of barometric pressure (mbar) of Tullianum data and meteorological data was performed first. As no HVAC is present in the Tullianum, barometric pressures must be equal, see fig.F in additional file 1.
In table 2 obtained coefficients are shown. A value of r2=0.9665 is observed, which is very good considering the distance of the two measurement points (850m) and the different instrument types and cost.
Further, lag-correlation is maximum for t=0 while all blue values worsen the situation, introducing a delay on the values of the Y axis, which are those measured in the Tullianum.
Computing all correlation parameters listed in table 2 for all space points, for temperature and humidity data, graphs containing the data referring to all vectors can be draft.
Humidity always shows the same delay with a r2=0.061, thus X and Y must be considered independent.
A simple visual comparison is not sufficient and certainly the common analysis of delays between minima and/or maxima (if detectable) is not useful in order to obtain a correlogramm.
The lag correlation is already evident in figs.4a and 4b, i.e. the line plots of the raw data obtained with a 1440 minutes time window for temperature data on day 5 and for RH data on day 4 respectively (days 4 and 5 were chosen due to strongly varying values of parameters).
The next step in our analysis on the building was to show up r2 values and slopes of the linear correlation of each main axe obtained by comparing outdoor and our indoor sampling points using a proper delay prior identified as above described.
Fig.10b shows lower values for both r2 and slope and a less evident trend for both parameters relative to the series of the shorter horizontal Carcer vector.
A flat trend and a further decrease for both slope (close to zero) and r2 values can be seen in fig.10c and d relative to the two horizontal Tullianum vectors, remarking that the ambient has to be considered an isolated one.
For all the vectors, lower slope and r2 values are evidenced for humidity with respect to temperature; such behaviour is consistent with a greater influence of the Tullianum humidity with respect to the macroclima.
Rooms isolation can be estimated as well by a Loss/Gain representation; firstly a daily mean is calculated for both the outdoor (meteorological) and the indoor data (the latter for both the Carcer and the Tullianum separately, considering all the relative sensors) then correlation graphs were drawn for temperature and humidity data.
To understand at the best the above behaviour, it is useful to consider that the position of the building entrance, that is the only point in direct communication with outdoor through an open-work door, is east-south-east oriented, so, sunlight can enter from about 8 until about 14. Further, it must be pointed out that Tullianum is a classical hypogean but also the Carcer is partially underground; as a consequence, illumination level is low everywhere as can be seen in fig.E in additional file 1 where only the upper sensor on the vertical vector (Hobo2) seems to be stimulated. On the other hand the underground position of the Carcer allows a lower influence of eventual external wind through the open-work door. Finally, all the walls are very thick so ensuring a good insulation.
An expected opposite trend, i.e. value are above the outdoor ones, is shown in fig.11b for humidity with the exception of data relative to Carcer on days 17, 21, 23, 24 and 29. The same flatting trend on variations, above seen for temperature, is evidenced with maximum value of about 25, 18 and 10 % for outdoor, Carcer and Tullianum respectively. Really the presence of groundwater ascending in the Tullianum and the huge walls funded directly on the water layer cause high internal humidity levels, like the one shown in the final part of graph in fig.3; any reasonable conservation concept should not interfere on this.
Column centering is an important data pretreatment because it leads to scatterplots which show directly the swing of the values on the axes; for fig.12 we obtain about 14 °C for temperature and about 70% for humidity with a negative correlation. Further outliers, trends and leverage are easily enhanced while centering does not modify distribution structure; instead it can lead to different results in the successive EDA phase .
In order to enhance the study of the distribution form for temperature and humidity data for each sensor, one can compute bar charts with class distribution. Choice of class width is not a trivial task, should consider values well above the instruments resolution level and producing a number of classes ranging from 10 to 30. Here we chose a width of 0.5 °C and of 2 % for temperature and humidity respectively, but different choices can change "the shape".
The complete series of figures are available in additional file 1, from fig.G-a to fig.V-c.
Among this huge amount of data, one should focus on data from sensor 2 (fig.13b), located at the top of the vertical vector (see fig.1); a wider ranging value for both parameters is evident with respect to the previous sensor (spread equal to 6.4 °C and 37%RH) and a no clear correlation could be seen.
Interesting is also point 5 (located in the lower part of the inner stair, see fig.14a; this figure shows a narrower and high sloping ellipsis with respect to the previous two; it must be also noted an enough crowded queue, out of the ellipsis, where a humidity variation of about 5% corresponds to an almost constant temperature. Humidity spread results equal to 20% so being no so different from that of sensor 1; an other similarity with such sensor is the positive correlation while temperature spread results significatively lower and equal to 2.3 °C. Looking for a similarity due to the position, it can be noted that, with respect to Carcer, sensor 1 is the nearest to outdoor and sensor 5 is the nearest to Tullianum; it is evident that macroclima have the maximum influence on the first while the effect of the very humid Tullianum is higher on the last. In fig.13c, relative to sensor 6, one can observe a near spherical shape for the ellipsis and a queue, similar to the one discussed above; in fig.14b, relative to sensor 8, a narrower shape along y axis results from the grouping of points inside the ellipsis and the external queue is more similar to the one in fig.13c. The two sensors are located above and below the hole connecting the two rooms and the above trends evidence the influence of the high humidity inside Tullianum on that of the Carcer. Basing on fig.14c it can be stated that, as expected, no influence from the macroclima is evidenced by sensor 11, located in the farest position inside the tunnel toward Cloaca Maxima; in particular, humidity is practically 100% for most of the points (see fig.R-c in additional file 1).
The class distributions for all the sensors (fig.G to fig.V, letters b and c in additional file 1) show a decreasing number of classes in the order meteo>Carcer>Tullianum with decreasing temperature and increasing humidity in the same order.
Referring to the Tullianum the main exchange with outdoor seems to take place along the stairs and not through the hole in the ceiling. Really figs.5(a-b) evidence significative differences between sensors 6 and 8 that are located, at only about half meter each, above and behind the grid connecting the two rooms; this can be explained with a scarce air flow through the grid and/or with a bigger influence of the Tullianum microclima on sensor 6. An air flow from the tunnel coming toward Cloaca Maxima, together with a water layer often present on the floor, seem to be the main responsible for the Tullianum microclima; anyway, also an air flow from the stairs plays a role as humidity never reaches 100% that, on the contrary, is the almost constant value evidenced in the tunnel. More, the low macroclima influence is also evident in figs.2 and 3.
In fig.18 the same trend is evidenced for the Carcer except for sensor 2 that, as already seen, is located in the position most influenced by the macroclima. Such figure also reveals that spontaneous values flat toward values more closer to those of Tullianum probably as a results of the good insulation.
We cannot give the right reply! Really this paper can give an input to perform a monitoring at the best but time was too short and, as a consequence, no seasonal but only daily variations were obtained; further, it was performed not in real conditions that foresee the entrance of visitors with all the consequent problems. Probably a longer monitoring would allow to obtain more accurate values for the spontaneous conditions an surely in such conditions the bigger damage due to variations can be eliminated, but who can ensure that spontaneous conditions would be the best for all the materials present in the building? Really in the building are present different tuffs, frescoes, stone column …; each material has different chemical-physical characteristics so being differently influenced by the environment. Lastly but not in importance: what conditions can ensure against damage from biological source? The presence of visitors allows microorganism entering the building, increase the carbon dioxide concentration and so on, but the absence of visitors cannot ensure that no biological damage can occur. More, in our case we have a sure capillary rising that can allow presence of salts that must be taken into account.
Probably is exact the statement of F.G. Ojea is exact "la conservación de las pinturas exige un mantenimiento de la temperatura y humedad a niveles lo más constantes posibles, y de un orden aproximado al que tenía el monumento antes de ser abierto, condiciones óptimas para su conservación, y que, obviamente, si fueran desfavorables no permitirían que el monumento llegara hasta nosotros”  and it is applicable to our hypogeum.
What must be done? A correct conservation project foresees a perfect knowledge of all the materials that must be saved and of all the surrounding to find a compromise between the best environmental conditions for all. In our case, surely frescoes are the most precious artworks and also those more sensitive to damage; surely the microenvironment must take into account this. A compromise must be adopted to allow visitors to enter the building but surely microclima cannot be projected for their well-being, as stated in some book.
In practice the right procedure is running as the Sovrintendence asked for cooperation to “La Sapienza” University to carry out analysis of constituent materials, map of salts and microclimate monitoring.
For the last a new campaign is running in the same seasonal period to evaluate the influence of visitors. Other rooms will be monitored to evaluate the reciprocal influence and seasonal campaign will follow.
Also citing by reference  "…not only is the relationship among all the variables at any one time important, but so is the entire past history of the trajectories of all these variables." confirm the use of multivariate analysis as necessary .
We may conclude this work stating that careful conservation of an entire ambient, here hypogean, declared cultural heritage, should foresee a triple action:
a) identification of real conditions microclimate over a long time-scale, with season changes and lasting possibly years;
b) comparison of the identified conditions with those indicated in official publications (see  for a list) in order to plan possible site tailored interventions, preferring passive methods;
c) publication of the done work in international journals in order to allow wide knowledge diffusion; work description should include details on measurement procedure AND RAW DATA to allow eventual further data elaboration by of other researchers.
This article has been published as part of Chemistry Central Journal Volume 6 Supplement 2, 2012: Proceedings of CMA4CH 2010: Application of Multivariate Analysis and Chemometry to Cultural Heritage and Environment. The full contents of the supplement are available online at http://journal.chemistrycentral.com/supplements/6/S2.
- Nilsson LO: Interaction between microclimate and concrete a prerequisite for deterioration. Constr Build Mater. 1996, 10 (5): 301-308. 10.1016/0950-0618(95)00046-1.View ArticleGoogle Scholar
- Banham R: The Architecture of the Well-tempered Environment. 1969, Architectural Press, London, 13: 978-0226036984Google Scholar
- Hedenblad G: Moisture permeability of mature concrete, cement mortar and cement paste. PhD dissertation. 1993, Lund University, Sweden, http://www.dissertations.se/dissertation/aa1ef97c3b/accessed at feb 2011,Google Scholar
- Cataldo R, Leucci G, Siviero S, Pagiotti R, Angelini P: Diagnostic of the conservation state in the crypt of the Abbey of Montecorona: biological, microclimatic and geophysical evaluations. J Geophys Eng. 2009, 6 (3): 205-220. 10.1088/1742-2132/6/3/001.View ArticleGoogle Scholar
- . List of the UNI-Normal available at: http://webstore.uni.com/unistore/public/searchproducts?action=search&usecache=false&committee=BENI+CULTURALI+-+NORMAL, accessed Feb 2011.
- Delogu S: Qualità microclimatica e consumi energetici nei musei. Thesis dissertation. 2010, Politecnico di Torino, ItalyGoogle Scholar
- Strada M, Carbonari A, Peron F, Porciani L, Romagnoni P: The microclimate analysis of tezone ‘105’ of the Venetian Arsenale. J Cult Herit. 2002, 3 (1): 89-92. 10.1016/S1296-2074(02)01164-0.View ArticleGoogle Scholar
- García-Diego F-J, Zarzo M: Microclimate monitoring by multivariate statistical control: The renaissance frescoes of the Cathedral of Valencia (Spain). J Cult Herit. 2010, 11 (3): 339-344. 10.1016/j.culher.2009.06.002.View ArticleGoogle Scholar
- Karner DB, Lombardi L, Marra F, Fortini P, Renne PR: Age of Ancient Monuments by Means of Building Stone Provenance: a Case Study of the Tullianum, Rome, Italy. J Archaeol Sci. 2001, 28 (4): 387-393. 10.1006/jasc.2000.0567.View ArticleGoogle Scholar
- Fortini P: Carcer tullianum il carcere mamertino al Foro romano. 1998, Mondadori Electa, 8843566237Google Scholar
- Ryhl-Svendsen M: Luftschadestoffe in Museen, Eine Einführung in Wirkungsweise, Monitoring und Kontrolle. Restauro. 2001, 8: 613-619.Google Scholar
- McGill R, Tukey JW, Larsen WA: Variations of Box Plots. Am Stat. 1978, 32 (1): 12-16. 10.2307/2683468.Google Scholar
- Pearson K: On lines and planes of closest fit to systems of points in space. Philos Mag. 1901, 2 (6): 559-572.View ArticleGoogle Scholar
- Hammer Ø, Harper DAT, Ryan PD: PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol Electron. 2001, 4 (1): 9-19.Google Scholar
- Cadima J, Jolliffe I: On Relationships Between Uncentred And Column-Centred Principal Component Analysis. Pak J Statist. 2009, 25 (4): 473-503.Google Scholar
- Montenegro Rúa EJ, Blanco-Rotea R, Benavides R, Portela C: Santa Eulalia De Bóveda. 2008, Xunta de Galicia, 20-978-84-453-4606-8Google Scholar
- Nomikos P, MacGregor JF: Monitoring batch processes using multiway principal component analysis. AIChE J. 1994, 40 (8): 1361-1375. 10.1002/aic.690400809.View ArticleGoogle Scholar
- La Gennusa M, Lascari G, Rizzo G, Scaccianoce G: Conflicting needs of the thermal indoor environment of museums: In search of a practical compromise. J Cult Herit. 2008, 9 (2): 125-134. 10.1016/j.culher.2007.08.003.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.