Next Article in Journal
The Southern Hemisphere Blocking Index Revisited
Next Article in Special Issue
Effects of Visitor Influx on the Indoor Climate of the Milan Cathedral
Previous Article in Journal
Interdecadal Change in the Relationship between the Winter Siberian High and Summer Tropical Cyclone Genesis Frequency over the Western North Pacific
Previous Article in Special Issue
Non-Invasive-Monitoring Methodology for the Evaluation of Environmental Impacts on Istrian Stone Surfaces in Venice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The European Standard EN 15757 Concerning Specifications for Relative Humidity: Suggested Improvements for Its Revision

by
Dario Camuffo
1,*,
Antonio Della Valle
1 and
Francesca Becherini
2
1
National Research Council of Italy—Institute of Atmospheric Sciences and Climate, Corso Stati Uniti 4, 35127 Padua, Italy
2
National Research Council of Italy—Institute of Polar Sciences, Via Torino 155, 30172 Venice Mestre, Italy
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(9), 1344; https://doi.org/10.3390/atmos13091344
Submission received: 22 July 2022 / Revised: 14 August 2022 / Accepted: 17 August 2022 / Published: 23 August 2022
(This article belongs to the Special Issue Atmospheric Environment and Cultural Heritage Protection)

Abstract

:
The European Standard EN 15757: 2010 ‘Conservation of Cultural Property—Specifications for temperature and relative humidity to limit climate-induced mechanical damage in organic hygroscopic materials’ is a guide specifying the allowed limits of variability of the indoor climate, in particular relative humidity (RH) to preserve cultural heritage objects and collections composed of climate-vulnerable materials. This paper is finalized to provide useful elements to improve the Standard at its next revision, based on focused research. The methodologies and the mathematical tools used are performed on 18 case studies representing different buildings, climates, and use, including heated and unheated buildings, museums, churches, concert halls, archives, and storage rooms. The first aim is to compare the method based on the centred moving average suggested by Annex A of EN15757 with an alternative method based on percentile interpolation to calculate the reference RH values, and in particular the safe band of RH variability, as well as the upper and lower risky bands. It has been found that the two methods provided the same results, but the latter is easier to manage. The second aim is to verify if the duration of the record necessary for the determination of the safe band is really 13 months of measurements as required by the Standard to account for the specific request of the centred moving average with a 30-day time window. This paper demonstrates that the same goal may be reached with a 12-month record, but extracting from the record itself the two periods required by the time window, i.e., the last 15 days of the year will be copied before the start of the record, and the same with the first 15 days after the end. The third aim is to test if the particular choice of the width of the time window is influential on the width of the safe band, and to assess the relationship between the width of the safe band and the width of the time window. The results show that the safe band logarithmically depends on the length of the time window, so it is crucial to respect the 30-day window established by the Standard.

1. Introduction

The European Standard EN 15757: 2010 (This Standard was produced in 2010 and confirmed in 2015 and 2020. In confirmed standards, the use is to cite the year of the first adoption, not the last confirmation).‘Conservation of Cultural Property—Specifications for temperature and relative humidity to limit climate-induced mechanical damage in organic hygroscopic materials’ [1] is a guide specifying the temperature (T) and relative humidity (RH) to preserve cultural heritage by limiting physical damage induced by strain-stress cycles in objects composed of, or containing organic hygroscopic materials, e.g., wood, painting on canvas, books, graphic documents, textiles, bone, ivory, or leather. The aim is to preserve the physical integrity of cultural objects, avoiding either the deterioration of the materials of which they are composed, or critical conditions when different materials are bound together, e.g., in a book the shrinkage of the sheets of paper is different from that of its lambskin cover and twine binding; if the RH drops too much, the tensions between these different elements will tear off the book.
The rationale of the Standard EN 15757 is that over time the considered materials have ‘acclimatised’ to the microclimate conditions where they have been kept for significant periods of time. These conditions are named ‘historical microclimate’. This means that these materials have been adapted to the average values, the seasonal changes, and the everyday fluctuations of the ambient T and RH. The adaptation process is long and causes negative consequences for the materials because seasonal cycles and short-term fluctuations of T and RH cause repeated internal stress, temporary or permanent deformations (e.g., shrinkage or swelling), and create internal micro or macro fractures. The role of these fractures is to act as expansion joints to dissipate internal stress, thus enabling a wider range of acceptable T and RH fluctuations. However, this acclimatisation is achieved at the expense of the integrity of the material which is now weaker and more vulnerable. The repetition of new stress, especially for enhanced RH fluctuations, may induce some internal micro cracks to grow, extend in size, become visible, and disfigure the object.
Section 5 contains the Standard’s explanation for why it is especially focused on RH. The reason is that, in organic hygroscopic materials, the mechanical damage (e.g., permanent deformation, formation of cracks) is more affected by changes in RH than in T. The Standard specifies, however, that in the particular case RH and T have the same priority, i.e., similar impacts, the specifications and the methodologies given for RH hold for T as well. This article adopts the same criterion, i.e., it formally deals with RH, but the same conclusions can be applied to T.
As illustrated in other papers [2,3], this European Standard establishes two basic principles: (i) the need to maintain stable environmental conditions, and (ii) the priority of historical climate. If this Standard is compared with ASHRAE [4,5], one finds that the former priority is considered in both, while only the European Standard considers the historical climate as a priority.
In reality, ASHRAE [4,5] considers selected bands of T and RH, irrespective of the past history of the objects. In addition, ASHRAE introduces the concept of taking advantage of unplanned extreme events, that may be interpreted like a crash test, i.e., suppose an extreme, unplanned RH fluctuation happened by chance, could the object survive this tremendous impact without serious consequences? This unplanned test has been called ‘proofed’ fluctuation [6] and is the largest climate fluctuation that an object has accidentally experienced passing apparently untouched. The underlying idea is that if the object has survived the first time, it can survive all the others as well. However, it has been observed that this is not true and that in the long run the micro cracks grow internally until they reach the surface and appear as macroscopic cracks [2,7,8,9,10,11,12,13].
When the European Standard was drafted (one of the authors is among the drafters), the experts of the working group discussed how to define and calculate the safe and the risky fluctuations. Using the findings of measurements performed by the experts with laboratory tests and several case studies in museums, historical buildings, and churches, as a basis, it was agreed that a safe condition may be reached by excluding the most extreme fluctuations that have the highest potential of risk. It was agreed to exclude the most severe fluctuations up to 14% of the total, i.e., 7% of the lowest and 7% of the highest ones. The RH values falling between two reference percentiles, i.e., the 7-ile and the 93-ile are considered to lie in a safe band (SB) and acceptable; those external to SB lie in a risky area, are dangerous, and should be excluded. This situation determines two risky bands (RB): the lower RB form 0-ile to 7-ile characterised by the most extreme dry conditions, and the upper RB form 93-ile to 100-ile characterised by the most extreme humid conditions.
However, if the RH fluctuations depart less than 10% from the seasonal RH level, all fluctuations are considered acceptable.
This method complies with the World Meteorological Organization (WMO) and the Intergovernmental Panel on Climate Change (IPCC), i.e., “an extreme event is generally defined as the occurrence of a value of a weather or climate variable above (or below) a threshold value near the upper (or lower) ends (‘tails’) of the range of observed values of the variable” [14]. “An extreme weather event is an event that is rare within its statistical reference distribution at a particular place. Definitions of ‘rare’ vary, but an extreme weather event would normally be as rare as, or rarer than, the 10th or 90th percentiles” [15].
In general, for the plot showing SB, the upper and lower RB calculated over a calendar year define the microclimate of a room, with its seasonal variability and daily fluctuations. The method to calculate the target values, i.e., average, upper, and lower limit of SB and RB, as well as their quantitative assessment (i.e., the 7-ile and 93-ile) was not included in the normative part, but shown as an example in an informative Annex A. The example in Annex A was based on a record concerning a church in Poland, occasionally heated with warm air.
A number of researchers have investigated the historical climate of selected case studies and their studies constitute useful examples of application of this standard [13,16,17,18,19,20,21,22,23,24,25,26]. However, the methodology and the calculations required to assess the target values of this standard may be subject to some criticisms and simplifications, as discussed later, and may be improved.
It must be specified that the innovative European Standard EN 15757 was unanimously approved in 2010 and confirmed in 2015 and 2020. The next planned revision will be in 2025. However, this Standard has found some difficulty in practical applications. The Technical Committee 346, ‘Conservation of Cultural Heritage’, of the European Committee for Standardization (CEN) keeps yearly meetings for decision making, planning, and implementation. On the meeting of December 2021, it was decided that the next revision will be useful for ‘an editorial restyling of the Standard in order to improve the editing and comprehensibility, but leaving the content unchanged’. The next revision is an opportunity for a clearer presentation of methods and calculations. The task is assigned to the Working Group 7 (WG7), responsible for indoor climate issues.
This paper has the following aims: (i) to compare the method suggested in Annex A of EN 15757 with an alternative method to calculate the reference values, and in particular SB, and the upper and lower RB; (ii) to verify if the duration of the record necessary for the determination of SB is really 13 months as required by the Standard, or if the same goal can be reached with a 12-month record; (iii) to test if the width of the time window is influential on the width of SB, and to assess the physical relationship between time window and SB width; and (iv) to provide useful elements to WG7 to improve the Standard at its next revision.

2. Data and Methods

2.1. Datasets

This paper utilizes 18 case studies representing different buildings, climates, and use, including museums, churches, concert halls, archives, and storage rooms. Some of them are heated and some unheated, as specified in Table 1. The aim is not to compare different systems, but to test mathematical methods on a wide dataset. It must be specified, however, that the used records must be complete, without gaps [19]. As winter is the most critical period for the impact of heating, and the period of severe cold and intensive heating is relatively short, even a short gap in the record may severely affect the results.
The indoor climate of the case studies used in this paper is shown in Figure 1. Even if RH constitutes the key variable, the association with T is very useful to interpret RH and its dynamics for diagnostic purposes. For instance, the plot of T helps to recognize whether a room is heated (i.e., mild temperature (e.g., 12–18 °C) and low RH in winter) or unheated (i.e., low temperature (e.g., 0–6 °C) and unchanged RH in winter), as well its use (e.g., a continuous T line suggests no heating or continuous heating; peaks of T indicate occasional heating). Peaks of T associated with RH drops are typical of sudden central heating (e.g., warm-air heating); RH variability not strictly related to T is mainly due to weather. This helps to distinguish the natural RH variability due to climate, which is accepted by the Standard, from the RH variability due to use, that should be controlled.

2.2. Methods

2.2.1. Datasets Necessary for the Determination of the Safe Band: 12 or 13 Months?

The Standard requires a RH record to document the yearly cycle of RH and the short-term fluctuations. At least a one-year record is requested, or integer multiples of it. The record should include a whole yearly cycle, no matter the starting day.
However, the example reported in Annex A of EN 15757 is based on the centred moving average (CMA) with a 30-day time window, which requires 15 days before and after each reading [27]. In a 365-day record, the truncation at both ends of the record constitutes a problem because in the first and last 14 days of the record the RH readings are widow of those symmetrically distributed for 15 days around them, and it is impossible to calculate CMA. Suppose that we are dealing with a calendar year from 1 January to 31 December 2020 (Figure 2a). This is equivalent to a strip with two ends. To provide the necessary data, the Standard suggests extending the length of the strip, starting the record 15 days earlier (i.e., red rectangle from 16 to 31 December of the previous year 2019) and ending it 15 days later (i.e., red rectangle from 1 to 15 January of the subsequent year 2021).
The conclusion is that a record of 13 months is necessary to calculate CMA over 12 months. Observations continued for this additional month constitute a cost and a serious inconvenience. In addition, the added values recorded in the previous and the subsequent years are always slightly different, and this causes an inconsistency: the SB values calculated for the end of the year are not coincident with those at the beginning. In case the first 15 days, or the last 15 days, of the year depart a bit from the normal climate conditions, by combining them with the moving average, one obtains a mean value, which is better representative than each of the 15-day additional periods taken separately. It must be kept in mind that the year selected for the record is assumed to be, and must be, representative of the normal climate. Consequently, if the winter in the record has been characterized by unusual weather, the anomalous year cannot be taken as a typical reference. If the climate at the beginning of the recorded year does not match with the situation at the end of it, i.e., if the circle does not close, this becomes a serious concern.
We have considered that the additional month is requested by the particular method of calculation but does not add physical information to reach the target values required by the Standard. Therefore, it is advisable to avoid this extra month and perform the same calculations over the calendar year, e.g., from 1 January to 31 December 2020, but with a 12-month record (Figure 2b). The 15-day extension at both ends has been formally performed by repeating the corresponding data at both ends. Before the start of the record on 1 January 2020, we added the ending period of the record from 16 to 31 December 2020 (cyan rectangle). Similarly, after the end of the record on 31 December 2020, we repeated the period 1 to 15 January 2020 (yellow rectangle). This partial repetition at both ends for calculation purposes is acceptable because the record is considered representative of a normal year, so that the days of December of the previous year 2019, as well as those of January of the subsequent year 2021, are consistent with the corresponding days of the 2020 record. This method substantially bends the linear strip of the calendar year and joins the end of December with the beginning of January. In a two-dimensional representation (Figure 2c) this becomes a cyclic year, where all months are in circular sequence. The artifice of joining December with January allows the continuous calculation of CMA over all the months of a calendar year, and needs a 12-month record instead of a 13-month one. In addition, this method combines the two extremes of the year, thus ensuring the consistency between the SB values calculated for the end and the beginning of the year. The same considerations hold for any other selected date to start the climate record.

2.2.2. Method of the Informative Annex A: Fluctuations from the Centred Moving Average

The key variables used in the informative Annex A of EN 15757 are summarized in Table 2:
The informative Annex A is conceived to assist users in these calculations step by step. To this aim, Annex makes an example to fix the main concepts, i.e., key variables, and how to perform calculations. The used method, i.e., the centred moving average, is illustrated in this section and in the flowchart (Figure 3, first column).
The first step is to calculate CMA, a tool used to smooth out short-term fluctuations and highlight longer-term trends. The CMA value <RHmo(n, j)> is constituted by the average of all the readings falling within the 30-day window and is attributed to the central value position (n, j). If the record is based on ND readings per day, the time window includes 30 ND +1 readings. If readings are sampled every 5 min, ND = 288; 10 min, ND = 144; 15 min, ND = 96, and so on. After having calculated the average over the selected time window, the software continues repeating the same calculations of the unweighted arithmetic average, by advancing one reading (i.e., lower end from (n, j − 15) to (n + 1, j − 15); central value from (n, j) to (n + 1, j), and upper end from (n, j + 15) to (n + 1, j + 15)), and skipping the last value (i.e., (n, j − 15)).
For every individual reading, RH(n, j), the CMA value <RHmo(n, j)> is calculated, with a 30-day time window. Then the differences, ∆RH(n, j), between every reading and the corresponding mean are calculated, i.e., the most extreme negative value to the most extreme positive one. Fluctuations falling within the interval between the 7-ile and the 93-ile constitute SB, while the external ones RA. The 7-ile represents the lower boundary between SB and RA, while the 93-ile the upper boundary.
∆RH(n, j) = RH(n, j) − <RHmo(n, j)>
These differences constitute the so-called fluctuations. Finally, for each distribution, these fluctuations are ordered from the most extreme negative value to the most extreme positive one.

2.2.3. Alternative Method: Interpolation of the Percentile Distribution

An alternative method of calculating the same values is proposed in this paper and the comparison with the Standard is shown in the flowchart (Figure 3, second column). As the safe band is defined in terms of percentiles, it can be directly obtained from two selected percentiles, i.e., the 7-ile and the 93-ile of the RH(n, j) readings over the calendar year, avoiding the unnecessary calculations of the CMA value <RHmo(n, j)> and all the differences of ∆RH(n; j) from it. Briefly, it is sufficient to plot all the RH(n, j) readings over the calendar year and draw two interpolation lines, i.e., the 7-ile and 93-ile.
In fact, by summing the differences in Equation (1) to the CMA value, one obtains the original RH(n, j) readings:
∆RH(n; j) + <RHmo(n, j)> = [RH(n, j) − <RHmo(n, j)>] + <RHmo(n, j)> = RH(n, j)
i.e., the use of the single values RH(n, j) adopted by the percentile analysis is equivalent to summing the fluctuations to CMA as suggested by the Annex A.
It is evident that Equation (2) holds for any reference used instead of <RHmo(n, j) >, e.g., using the 50-ile instead of CMA.
Using percentiles, the temporal bin to be considered is the monthly time window, e.g., 30 days as in the Annex A. The percentile distribution can be directly calculated and interpolated with specific mathematic tools and the selected time window may be established as an optional input.
Theoretically, it should be noted that mean, mode, and median are coincident in a symmetrical distribution, while they are not in an asymmetrical one, as typically occurs for RH. Therefore, for RH, it is theoretically preferable to deal in terms of median, i.e., 50-ile, instead of mean, i.e., the centred mobile average, and then calculate fluctuations as the difference between the current reading and the 50-ile. The Working Group of the CEN Standard drafters discussed this item but concluded that this is a too subtle theoretical issue for the majority of users, and that it was preferable to use the concept of mean that has the advantage of being largely popular.
In this work, CMA and the method based on the interpolation of percentiles have been applied to the same set of case studies, putting into evidence the 0-ile, 7-ile, 50-ile, 93-ile, and 100-ile that represent the range, the median, and the borders of SB and the two RB. To test if the two methodologies give exactly the same results, the calculation of percentiles has been done using the same 30-day time window around RH(n, j).

2.2.4. Evaluation of the Most Convenient Time Window

The Standard establishes that the most extreme RH fluctuations shall be excluded, starting from the highest and lowest ones, until 7% of them has been excluded from the lower humidity side (i.e., the lower RB) and 7% from the higher humidity side (i.e., the upper RB). If RH readings are represented in a two-dimensional graph, having the calendar year in the abscissa and RH in the ordinate, one obtains a cloud of dots lying on a band that may change average value and width with the season. The central part of the band is more densely populated, and the density decreases approaching the upper and the lower boards. Removing the most extreme 7% of the dots on both sides has the effect of thinning the band. With a percentile computing facility, this is a very simple task.
However, if the issue is solved following the method described in Annex A, it is necessary to calculate all the distances, ∆RH(n; j), from a reference value, that has been suggested (i.e., suggested because the Appendix is only informative) to be the CMA with monthly window, i.e., <RHmo(n, j)>. The monthly window (i.e. 30 days) has been selected to confer stability to the reference value of the calendar day j. In fact, in the normal climate statistics, WMO [28] recommended referring to a period composed of 30 years which implies that the calendar day j is represented by the average of 30 years, e.g., j = 150 is represented by the average of all the days dated 30 May in the 1991–2020 period. As it is extremely rare to dispose of indoor records over a 30-year period, the Standard has considered obtaining a 30-day average by taking advantage of 15 days before and 15 days after the selected reading or the selected calendar day. This choice is finalized to reduce instability, i.e., the random difference between a selected calendar day j and the next one j + 1. The 30-year period is not related to physical issues, e.g., the relaxation time of wooden objects depends on the wooden species, wood thickness, surface treatment, and coatings.
Another crucial issue is related to the time window. In theory, if the RH distribution is homogeneous over the calendar year, always keeping the same distribution (e.g., purely random, Gaussian, Gamma), the SB width is independent of the selected time window, i.e., a horizontal band. However, if the distribution is subject to seasonal cycles, e.g., lower RH values in winter caused by heating, the band will be curved with lower values in the cold season. When the band is curved, by increasing the time window, one meets larger RH ranges and widens SB.
In this work, different options are considered for the time window to investigate if the safe band changes with the time window and how. Finally, we want to verify how critical the choice of one month for the time window is.

3. Results and Discussions

3.1. Comparison between the Two Methods: Centred Moving Average and Percentile Interpolation

The two methods, i.e., the centred moving average (CMA) suggested in Annex A, and the interpolation of the percentile distribution (IPD) have been applied to the selected case studies to investigate if they provide different results and to test which of them is easier to apply. For reasons of space, only four cases are shown in detail.
The CMA method has been calculated as indicated in Annex A that requires only three lines: the mean, 7-ile, and 93-ile. We have added the range, i.e., 0-ile and 100-ile. The lines are reported in Figure 4a, i.e., the red lines represent the contour lines of SB (i.e., 7-ile and 93-ile), while the blue lines the extremes of RB (i.e., 0-ile and 100-ile).
The IPD method has been calculated and reported on the same Figure 4a, but represented in terms of bands, i.e., light grey for SB, and dark grey for RB. Overlapping the two plots (i.e., the plot with bands and the plot with lines), they coincide, giving the appearance of only one plot, with contour lines highlighting the bands. This is obvious because the contours and the bands are given by the same percentiles, i.e., the two methods provide the same results because they calculate in different ways the same items.
In the same figure, we have added the mean (Annex A, green line) and the median (i.e., 50-ile, dotted line). The former is calculated as CMA; the latter is calculated as IPD, i.e., the median of a population included in a bin of 30 days around the central value, and then moving the window (i.e., the bin) to the next value and so on. Although mean and median are different statistical quantities, their values do not reach significant differences: in the plots, they appear very close between them, with some small departures. Instead of using the mean (Annex A), we propose adopting the median (i.e., 50-ile) because it is more consistent with the percentile method, is easier to calculate, and is supplied by the same mathematical tool. The choice of the reference for a central value, either mean or median, is the only small difference.
It is evident that the two methods, i.e., CMA and IPD, do not reach significant differences: i.e., the 0-ile, 7-ile, 93-ile, and 100-ile are exactly the same; the only small difference is the central line, i.e., CMA mean or IPD 50-ile. Therefore, each of the two methods can be equally taken as a target reference.
In the plots, the two external RB, and to a lesser extent SB, have the appearance of rectangles, with 30-day width. This is a consequence of the choice of using a simple mobile average with rectangular window, i.e., the arithmetic average of data all with equal weight 1/NTW, where NTW is the number of data included in the time window [27,29]. When the rectangular window of a mobile average includes the peak of an extreme event, the arithmetic average abruptly increases forming a rectangle, and the rectangle continues until the window with its movement moves forward and leaves out the extreme event. The same happens with the percentile analysis when an extreme event enters a mobile bin. To avoid the issue of extreme events that may generate rectangles, a weighted moving average with a different shape should be preferred, with low weight at the borders and high weight at the centre, either linearly distributed (e.g., triangular), or bell shaped (e.g., Gaussian, cosine, parabolic, Poisson, Hann, Abel, Cauchy, Tukey, Tuckey-Hanning) [29,30]. Another strategy is to apply filters or interpolation tools to smooth edges. However, the aim of this paper is to show and explain the methods given by the Standard. The user is free to consider an edgy or a smooth plot, and adopt any smoothing tool, simply because this topic is not considered in the normative part.
It may be useful to comment on the examples in Figure 4a, by comparing them with their indoor climate in Figure 1. In the case study of the unheated Museum, Udine, with passive climate control, the SB width is basically constant over time and quite horizontally distributed over the calendar year. The case study of the Concert Hall, Padua with continuous basic heating with an occasional extra heating, shows that SB is lower in winter and higher in summer, except in July when the windows are randomly opened for ventilation to resist the hot local climate. In such a case, the abrupt change in thickness of SB is a diagnostic index of incorrect use of doors or windows for summer ventilation. Another typical example is the Church of Rocca Pietore, where the occasional warm-air heating for the Sunday celebrations cause dramatic peaks of temperature and drops of RH, widening very much SB and RB. The peaks of warm-air heating can be recognized by the increase of both bands, but especially RB. The Basilica S. Maria Gloriosa dei Frari, Venice, has a modest local occasional heating for celebrations with a very low impact on the indoor climate, but is a tourist attraction for the important artworks kept inside. The combination of the humid winter climate in Venice and the intense use widens the safe band in winter.
The whole dataset is represented in Figure 4b where a soft, numerical smoothing filter has been applied to smooth the rough edges of the plot that are affected by some random noise for the limited population of readings. We must specify that the filter is not necessary, but has been applied for educational purposes, i.e., to point out that the user should not be concerned to strictly follow the small day-by-day changes, but should have a broad vision of the general trend on a wider timescale. To avoid random noise, it is possible to use various systems. For instance, one may apply this procedure not to the whole dataset, but to a subsample of it, e.g., two or three values per month to obtain smooth plots that show long-term variations, reducing noise. Plots in Figure 4a are without smoothing. Plots in Figure 4b have been obtained by applying the techniques described in Annex A, but with two values per month. Using opportune mathematical tools, in this case a spline interpolation, the resolution of the data has been increased to one value per day. The soft smoothing level of the filter is evident by comparing the four examples in Figure 4a (unfiltered) with the same in Figure 4b (filtered). However, the Standard does not consider these aspects, and the user may either avoid filters or adopt the preferred one.

3.2. Comparison between Different Time Windows

Selected time windows have been considered to investigate how the safe band is affected by this choice. The time windows used are 1 day, 1 week, two weeks, 1 month, 2 months, and 1 year. This choice has been made to test the different impact that the time window may have on sites with higher or lower seasonal cycles. The investigation has shown that the average safe band width <SBW> (i.e., the average of the SB over the calendar year) is highly dependent on the particular choice of the time window; the wider the time window, the wider <SBW>.
In a linear diagram, i.e., with <SBW> in ordinate (Y), and the number of days of which the time window is composed in abscissa (X), and both X and Y have a linear scale, the plots of Y versus X are logarithmic (Figure 5a). Therefore, if a semi-logarithmic diagram is considered, i.e., the axis X with a logarithmic scale, and Y with a linear scale, the plots become straight-lines, as shown in Figure 5b for four case studies in expanded size. The whole set is represented in Figure 5c.
For all the case studies, the relationship between <SBW> and the width X (in days) of the time window is given by the general logarithmic equation
<SBW> = a ln X + b
where ln is the natural logarithm. The coefficient a is representative of the weight of the logarithmic function that represents the slope of the best-fit line in the semi-logarithmic plots in Figure 5a,b. The constant term b [%] is not the intercept with the Y axis (the scale is logarithmic and there is not X = 0) but is the starting value, i.e., the SB obtained with 1-day time window. The coefficients a [%], b [%], are reported in Table 3. In the selected case studies, 1.35 ≤ a ≤ 5.37%; the coefficient b is largely variable, i.e., 0.48 ≤ b ≤ 13.76%.
The last column in Table 3 is the 30-day intercept that corresponds to the value selected by the Standard as a reference. It represents the average value of the RH (%) variability that characterizes SB over the whole calendar year. The lowest values correspond to a thin SB, which indicates that the indoor climate is very stable and provides the best conditions for conservation. On the contrary, a large intercept means a wide SB, high indoor climate variability, and poor conditions for conservation. The observed values of the 30-day intercept lie between 5.1 and 22.9%. Only three cases stay below 10%; more than 50% of them lie between 15 and 23%; and 5 of them exceed 20%.
A comment on the four examples in Figure 5a. The Museum, Udine, with homogeneous SB quite horizontally distributed over the calendar year has very low coefficients a = 1.35% and b =2.38%. The Church of Rocca Pietore, is characterized by an intermediate value of a = 2.37% and a very high value of b = 10.49%. The Concert Hall, Padua, and the Basilica S. Maria Gloriosa dei Frari, Venice, have quite similar values, i.e., high a = 4.72% and 4.26% and intermediate b = 5.52% and 8.34%, respectively.
It is not surprising to find that the best-fit of Equation (3) is a logarithm. The Lognormal is a right-skewed continuous probability distribution introduced by Galton [31], frequently used to represent physical variables that never take negative values, or to represent risk as a function of the time [32,33,34]. This is consistent with the extreme value theory (EVT), where extreme events are those contained in the tail distribution of a given variable [35]. The longer the record, or the wider the time window, the higher the probability of including the most extreme events located in the upper part of the probability distribution.

4. Language and Communication

In this paper, the term ‘safe band’ has been used relying on the standard. However, it may be useful to underline that the term ‘safe’ might be misleadingly interpreted. For instance, if some activity is performed that widens the span between the 7-ile and the 93-ile, this increases the width of SB, but does not correspond to an increase of safety. Conversely, every widening of SB corresponds to a worsening of the conditions for conservation. This becomes clearer when the risky bands are considered: every widening of RB implies a worsening for conservation.
Therefore, instead of ‘safe’ it may be appropriate to speak in terms of ‘accustomed’ band, but only when the natural, unperturbed indoor climate is considered, i.e., non-biased for use, heating, or air conditioning. In the general case, the use of neutral names, e.g., ‘usual variability band’ or ‘band between the 7-ile and the 93-ile’ may be preferable. This is a further aspect that should be considered by WG7 in the forthcoming revision of this Standard.

5. Conclusions

The first important conclusion is that the methodology proposed by the Standard EN 15757 was tested on 18 case studies and was compared with another alternative methodology in view of the next revision of the Standard that has been decided by the European Committee for Standardization. These findings will help to simplify and improve the quality of the standard after having considered the two main difficulties met by most users, i.e., the practical application of centred moving average and the need to monitor over 13 months to get information for a single calendar year.
A relevant conclusion is that the centred moving average suggested in the informative Annex A of the Standard may be substituted with the method of the interpolation of the percentile distribution that gives the same results, but with an easier mathematical approach. It has recognized that the centred moving average requires a lot of unnecessary calculations that may create serious difficulties for most users. On the other hand, the method of interpolation of the percentile distribution is much easier, and specific statistical tools are directly available for the calculation of selected percentiles. In the Standard revision, it will be advisable to present and explain both methods, and leave the user free to choose.
It has been recognized that the adoption of a rectangular window, either for simple (unweighted) arithmetic moving average, or for a moving bin in the case of percentiles, isolated peaks, or drops, do generate rectangles with the same width of the time window. This drawback may be reduced using a weighted moving average or a smoothing filter. Alternatively, the user should be prepared to correct the interpretation of plots.
A relevant simplification is that it is not necessary to monitor for 13 months to characterize a calendar year using centred moving average calculations. A 12-month record is sufficient if one uses the artifice of extending the data set by reporting the first 15 days of the year at the end of the record, and the last 15 days of the year before the beginning of the record. This is equivalent to dispose of a circular, cyclic year, not a single strip broken at both ends. This choice reduces costs and time for preliminary investigations.
A further conclusion is the relevance of the duration of the interval of time selected for the time window. It has been clarified that the safe band increases with the time window, following a logarithmic law. However, the coefficients that determine the equations change case by case, depending on the indoor climate, winter heating, use, and ventilation. Therefore, the impact is different, being smaller in the case of homogeneous RH distributions over the calendar year, and stronger in the case of winter heating, especially when this is made occasionally with warm-air systems that generate sharp peaks of temperature and drops of relative humidity.
A final note concerns language and communication. The term ‘safe band’ might be misleadingly interpreted and it would be preferable to substitute it with a neutral name. Some suggestions have been made, but the most convenient term will be decided by the Working Group 7 at the next Standard revision.

Author Contributions

All authors have equally contributed to the paper. D.C. is member of CEN, TC346 and is coordinator of the Working Group 7 that has drafted the EN 15757: 2010 Standard. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been performed without funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The microclimate data are not publicly available. The ownership of microclimatic data collected and used in this work depends on the particular case study and is subject to specific conditions.

Acknowledgments

We are grateful to the former colleagues Giovanni Sturaro, Antonio Valentino, Emanuela Pagan and Chiara Bertolin, who contributed to the collection of data under contract with the Institute of Atmospheric Sciences and Climate, under various projects, e.g., AER, Sensorgan, Climate for Culture, Unesco for Venice. We are also indebted to the EU funded project Sensorgan (Sensor system to detect harmful environments for pipe organs) and the Project coordinator Carl Johan Bergsten, Gothenburg University, for having kindly shared some records. Finally, we are grateful to Anna Maria Siani, Francesca Frasca and Elena Verticchio, University La Sapienza, Rome, for useful discussions.

Conflicts of Interest

The authors declare no conflict of interest. The periodic revision and improvement of Standards is part of the normal task of CEN members and WG coordinators.

References

  1. EN 15757 Standard; Conservation of Cultural Heritage—Specifications for Temperature and Relative Humidity to Limit Climate-Induced Mechanical Damage. European Committee for Standardization (CEN): Brussels, Belgium, 2010.
  2. Camuffo, D. Microclimate for Cultural Heritage—Measurement, Risk Assessment, Conservation, Restoration and Maintenance of Indoor and Outdoor Monuments, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  3. Camuffo, D. The European Standards on Environment and Cultural Heritage, and the Activity of WG4 Working Group at UNI. In La Normativa Tecnica Italiana Ed Europea Per Il Patrimonio Culturale. Atlanti Della Conservazione Del Patrimonio Culturale; Fassina, V., Ed.; Nardini: Florence, Italy, 2022; pp. 61–75. [Google Scholar]
  4. ASHRAE. Museums, Galleries, Archives and Libraries. In ASHRAE Handbook—HVAC Applications; American Society of Heating, Refrigerating, and Air-Conditioning Engineers: Atlanta, GA, USA, 2007; Chapter 21; pp. 21.1–21.23. [Google Scholar]
  5. ASHRAE. Museums, Galleries, Archives and Libraries. In ASHRAE Handbook—HVAC Applications; American Society of Heating, Refrigerating, and Air-Conditioning Engineers: Atlanta, GA, USA, 2011; Chapter 23; pp. 23.1–23.23. [Google Scholar]
  6. Michalski, S. The Ideal Climate, Risk Management, the ASHRAE Chapter, Proofed Fluctuations, and Toward a Full Risk Analysis Model. In Proceedings of the Contribution to the Experts’ Roundtable on Sustainable Climate Management Strategies, Tenerife, Spain, 23–27 April 2007; The Paul Getty Institute: Los Angeles, CA, USA, 2009. [Google Scholar]
  7. Bratasz, Ł.; Kozłowski, R.; Camuffo, D.; Pagan, E. Impact of Indoor Heating on Painted Wood. Monitoring the Altarpiece in the Church of Santa Maria Maddalena in Rocca Pietore, Italy. Stud. Conserv. 2007, 52, 199–210. [Google Scholar] [CrossRef]
  8. Camuffo, D.; Pagan, E.; Schellen, H.; Limpens-Neilen, D.; Kozlowski, R.; Bratasz, L.; Rissanen, S.; Van Grieken, R.; Spolnik, Z.; Bencs, L.; et al. Church Heating and Preservation of the Cultural Heritage: A Practical Guide to the Pros and Cons of Various Heating Systems; Electa Mondadori: Milano, Italy, 2007. [Google Scholar]
  9. Camuffo, D.; Pagan, E.; Rissanen, S.; Bratasz, Ł.; Kozłowski, R.; Camuffo, M.; della Valle, A. An Advanced Church Heating System Favourable to Artworks: A Contribution to European Standardisation. J. Cult. Herit. 2010, 11, 205–219. [Google Scholar] [CrossRef]
  10. Bratasz, Ł. Acceptable and Non-Acceptable Microclimate Variability: The Case of Wood. In Basic Environmental Mechanisms Affecting Cultural Heritage—Understanding Deterioration Mechanisms for Conservation Purposes; Camuffo, D., Fassina, V., Havermans, J., Eds.; COST Action D42 “Enviart”; Nardini: Florence, Italy, 2010; Chapter 4; pp. 49–58. [Google Scholar]
  11. Jakieła, S.; Bratasz, Ł.; Kozłowski, R. Acoustic Emission for Tracing the Evolution of Damage in Wooden Objects. Stud. Conserv. 2007, 52, 101–109. [Google Scholar] [CrossRef]
  12. Lukomski, M.; Strojecki, M.; Pretzel, B.; Blades, N.; Beltran, V.L.; Freeman, A. Acoustic Emission Monitoring of Micro-Damage in Wooden Art Objects to Assess Climate Management Strategies. Insight-Non-Destr. Test. Cond. Monit. 2017, 59, 256–264. [Google Scholar] [CrossRef]
  13. Califano, A.; Baiesi, M.; Bertolin, C. Analysing the Main Standards for Climate-Induced Mechanical Risk in Heritage Wooden Structures: The Case of the Ringebu and Heddal Stave Churches (Norway). Atmosphere 2022, 13, 791. [Google Scholar] [CrossRef]
  14. WMO. Guidelines on the Definition and Monitoring of Extreme Weather and Climate Events—TT-DEWCE WMO 14 April 2016; World Meteorological Organization: Geneva, Switzerland, 2016. [Google Scholar]
  15. IPCC. Climate Change 2014: Synthesis Report; Mach, K.J., Meyer, L.A., Pachauri, R.K., Planton, S., von Stechow, C., Eds.; Annex II: Glossary. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  16. Camuffo, D.; Bertolin, C.; Bonazzi, A.; Campana, F.; Merlo, C. Past, Present and Future Effects of Climate Change on a Wooden Inlay Bookcase Cabinet: A New Methodology Inspired by the Novel European Standard EN 15757:2010. J. Cult. Herit. 2014, 15, 26–35. [Google Scholar] [CrossRef]
  17. Bertolin, C.; Camuffo, D.; Bighignoli, I. Past Reconstruction and Future Forecast of Domains of Indoor Relative Humidity 
Fluctuations Calculated According to EN 15757:2010. Energy Build. 2015, 102, 197–206. [Google Scholar] [CrossRef]
  18. Taylor, T. Preservation of Cultural Heritage: The Design of Low-Energy Archival Storage. In Advanced Technologies for Sustainable Systems. Lecture Notes in Networks and Systems; Bahei-El-Din, Y., Hassan, M., Eds.; Springer: Berlin, Germany, 2017; Volume 4. [Google Scholar]
  19. Frasca, F.; Siani, A.M.; Casale, G.R.; Pedone, M.; Bratasz, L.; Strojecki, M.; Mleczkowska, A. Assessment of Indoor Climate of Mogiła Abbey in Kraków (Poland) and the Application of the Analogues Method to Predict Microclimate Indoor Conditions. Environ. Sci. Pollut. Res. 2017, 24, 13895–13907. [Google Scholar] [CrossRef]
  20. Leijonhufvud, G.; Broström, T. Standardizing the Indoor Climate in Historic Buildings: Opportunities, Challenges and Ways Forward. J. Archit. Conserv. 2018, 24, 3–18. [Google Scholar] [CrossRef] [Green Version]
  21. Żaba, A.; Marchacz, M. Historical Climate of the Historic Church of St. George at Ostropa. E3S Web Conf. 2018, 49, 00138. [Google Scholar] [CrossRef] [Green Version]
  22. Díaz-Arellano, I.; Zarzo, M.; García-Diego, F.J.; Perles, A. A Methodology for the Multi-Point Characterization of Short-Term Temperature Fluctuations in Complex Microclimates Based on the European Standard EN 15757:2010: Application to the Archaeological Museum of L’Almoina (Valencia, Spain). Sensors 2021, 22, 7754. [Google Scholar] [CrossRef] [PubMed]
  23. Manfriani, C.; Gualdani, G.; Goli, G.; Carlson, B.; Certo, A.R.; Mazzanti, P.; Fioravanti, M. The Contribution of IoT to the Implementation of Preventive Conservation According to European Standards: The Case Study of the “Cannone” Violin and Its Historical Copy. Sustainability 2021, 13, 1900. [Google Scholar] [CrossRef]
  24. Verticchio, E.; Frasca, F.; Bertolin, C.; Siani, A.M. Climate-Induced Risk for the Preservation of Paper Collections: Comparative Study Among three Historic Libraries in Italy. Build. Environ. 2021, 206, 108394. [Google Scholar] [CrossRef]
  25. Fabbri, K. Historic Climate in Heritage Building and Standard 15757: Proposal for a Common Nomenclature. Climate 2022, 10, 4. [Google Scholar] [CrossRef]
  26. Verticchio, E.; Frasca, F.; Cavalieri, P.; Tedonio, L.; Fugaro, D.; Siani, A.M. Conservation Risks for Paper Collections Induced by the Microclimate in the Repository of the Alessandrina Library in Rome (Italy). Herit. Sci. 2022, 10, 80. [Google Scholar] [CrossRef] [PubMed]
  27. Burt, J.E.; Barber, G.M.; Rigby, D.L. Elementary Statistics for Geographers, 3rd ed.; Guilford Press: New York, NY, USA, 2009. [Google Scholar]
  28. WMO. Guidelines on the Calculation of Climate Normals—WMO No. 1203; World Meteorological Organization: Geneva, Switzerland, 2017. [Google Scholar]
  29. Wei, W.W.S. Time Series Analysis; Addison-Wesley: Redwood City, CA, USA, 1990. [Google Scholar]
  30. Poularikas, A.D. Windows. In The Handbook of Formulas and Tables for Signal Processing; Poularikas, A.D., Ed.; CRC Press LLC: Boca Raton, FL, USA, 1999; Chapter 7; pp. 131–148. [Google Scholar]
  31. Galton, F. The Geometric Mean in Vital and Social Statistics. Proc. R. Soc. Lond. 1879, 29, 365–367. [Google Scholar]
  32. Hart, P.E. Lognormal Distribution. In Econometrics; Eatwell, J., Milgate, M., Newman, P., Eds.; Palgrave Macmillan: London, UK, 1990. [Google Scholar]
  33. Papalexiou, S.M.; Koutsoyiannis, D.; Makropoulos, C. How Extreme is Extreme? An Assessment of Daily Rainfall Distribution Tails. Hydrol. Earth Syst. Sci. 2013, 17, 851–862. [Google Scholar] [CrossRef] [Green Version]
  34. Colombo, D.; Abreu, D.; Ramos Martins, M. Application of Markovian Models in Reliability and Availability Analysis: Advanced Topics. In Safety and Reliability Modeling and Its Applications; Phamand, H., Ram, M., Eds.; Elsevier: Dordrecht, The Netherlands, 2021; pp. 91–160. [Google Scholar]
  35. Coles, S.G. An Introduction to Statistical Modeling of Extreme Values; Springer: New York, NY, USA, 2001. [Google Scholar]
Figure 1. Indoor climate, i.e., temperature and relative humidity of the selected case studies.
Figure 1. Indoor climate, i.e., temperature and relative humidity of the selected case studies.
Atmosphere 13 01344 g001
Figure 2. (a) The method proposed by EN 15757 requires a 13-month record: i.e., 12 months of the calendar year (white strip) plus two periods of 15 days at both extremes (red rectangles). (b) The method proposed in this paper, i.e., 12-month record and repetition of the first 15 days (yellow rectangle) and last 15 days (cyan rectangle) at both ends of the year for calculation purposes. R for ‘repeated’. (c) A cyclic year, where all months are in circular sequence with December connected to January.
Figure 2. (a) The method proposed by EN 15757 requires a 13-month record: i.e., 12 months of the calendar year (white strip) plus two periods of 15 days at both extremes (red rectangles). (b) The method proposed in this paper, i.e., 12-month record and repetition of the first 15 days (yellow rectangle) and last 15 days (cyan rectangle) at both ends of the year for calculation purposes. R for ‘repeated’. (c) A cyclic year, where all months are in circular sequence with December connected to January.
Atmosphere 13 01344 g002
Figure 3. Flowcharts showing the method of the centred moving average used in the Annex A of EN15757, and the method based on the interpolation of percentiles to obtain the safe band (SB).
Figure 3. Flowcharts showing the method of the centred moving average used in the Annex A of EN15757, and the method based on the interpolation of percentiles to obtain the safe band (SB).
Atmosphere 13 01344 g003
Figure 4. (a) Comparison between the CMA method suggested in the Annex A (represented with contour lines, i.e.,: blue: 0-ile and 100-ile, red: 7-ile, 93-ile; green: mean), and the IPD method (represented with bands, i.e., light grey: safe band; dark grey: risky band; dotted line: median, i.e., 50-ile). Lines and bands are calculated as moving average, or moving bins, with a 30-day time window. (b) Safe band (light grey) and risky band (dark grey) of all the case studies to which smoothing has been applied.
Figure 4. (a) Comparison between the CMA method suggested in the Annex A (represented with contour lines, i.e.,: blue: 0-ile and 100-ile, red: 7-ile, 93-ile; green: mean), and the IPD method (represented with bands, i.e., light grey: safe band; dark grey: risky band; dotted line: median, i.e., 50-ile). Lines and bands are calculated as moving average, or moving bins, with a 30-day time window. (b) Safe band (light grey) and risky band (dark grey) of all the case studies to which smoothing has been applied.
Atmosphere 13 01344 g004
Figure 5. The safe band width is logarithmically related to the time window. (a) The example of a case study in a linear coordinate paper. The plot is logarithmic. (b) In semi-logarithmic paper (i.e., days in logarithmic scale) the plots become straight-lines. The four case studies used as detailed examples. (c) The whole set of case studies.
Figure 5. The safe band width is logarithmically related to the time window. (a) The example of a case study in a linear coordinate paper. The plot is logarithmic. (b) In semi-logarithmic paper (i.e., days in logarithmic scale) the plots become straight-lines. The four case studies used as detailed examples. (c) The whole set of case studies.
Atmosphere 13 01344 g005
Table 1. Characterization of the case studies.
Table 1. Characterization of the case studies.
NoNameLocationUseClimate Control Period
1Uffizi Gallery °Florence (IT)MuseumHVAC during opening hours; continuous RH control1998
2Ala Ponzone °Cremona (IT)MuseumHVAC, humidity control2011
3Museo Vescovile °Udine (IT)Museumpassive climate control; filtered light and shutter control2015–2016
4Ca’ Granda °Milano (IT)Archive & book storage with restricted accesspassive climate control2011–2012
5Liviano °Padua (IT)Monumental concert hallcontinuous basic heating + occasional extra heating 2002–2003
6S. Maria Gloriosa dei Frari °Venice (IT)Church; tourist attractionmodest local heating2009–2010
7S. Maria Maddalena °Rocca Pietore (IT)Churchoccasional winter heating 2002–2003
8S. Maria Maggiore °Rome (IT)Church; tourist attractionno HVAC 1996–1997
9St Andrew the Apostle °Olkusz (PO)Churchoccasional winter heating 2007–2008
10Madonna di Campagna *Valtellina (IT)Church with concertsno HVAC2003–2004
11S. Maria Collemaggio *L’Aquila (IT)Church with concertsno HVAC2003–2004
12Walloon Church
(Waalse Kerk) *
Amsterdam (NL)Church with concertscontinuous winter heating2003–2004
13St Willibrord
(Groene Kerk) *
Oegstgeest (NL)Church with concertscontinuous winter heating2003–2004
14St Jacob *Hamburg (GE)Church with concertsoccasional winter heating2003–2004
15St Michael *Leuven (BE)Church with concertscontinuous winter heating 2004–2005
16Notre Dame du Sablon *Brussels (BE)Church with concertsuneven winter heating2004–2005
17Old Choir
Monastic Church °
Padua (IT)Church with restricted accessno HVAC; almost unperturbed indoor climate2011–2012
18Church basement °Venice (IT)Basement with restricted accessno HVAC; unperturbed indoor climate2020
The case studies with ° have been monitored by the Institute of Atmospheric Sciences and Climate. Those with * have been shared by Carl Johan Bergsten, Gothenburg University, during the common project Sensorgan. All records are unpublished.
Table 2. Key variables used in the informative Annex A of EN 15757.
Table 2. Key variables used in the informative Annex A of EN 15757.
VariableDefinitionExample
n sequential number of the reading sampled in the calendar day jn = 1 is the first reading of the day j, starting from midnight; n = 2 is the second reading and so on
jsequential number of the calendar day (1≤ j ≤ 365)j = 45 represents the 45th day of the year, i.e., 14 February
RH(n, j)individual readings of the record RH(n, j)the nth value of RH recorded in the day j of the calendar year
<RHmo(n, j)>average of RH calculated with the centred moving average, using a monthly window <RHmo(n, 45)> for j= 45 (i.e., 14 February). It represents the average RH around the reading n, from j − 15 = 30 (i.e., 30 January) to j + 15 = 60 (i.e., 1 March)
∆RH(n, j)fluctuations defined as the difference RH(n, j) − <RHmo(n, j)>difference between every individual reading RH(n, j) and the moving average <RHmo(n, j)> centred on it
<RHyr>yearly average of RHaverage from 1 January to 31 December
Table 3. Coefficients a, b of the logarithmic Equation (3) representing the dependence of the average safe band on the selected time window duration. The last column is the 30-day intercept.
Table 3. Coefficients a, b of the logarithmic Equation (3) representing the dependence of the average safe band on the selected time window duration. The last column is the 30-day intercept.
NoCase Studyab30-d Intercept
1Museum, Florence (IT)2.7013.7622.9
2Museum, Cremona (IT)1.360.485.1
3Museum, Udine (IT)1.352.387.0
4Archive and book storage, Milan (IT)1.862.538.9
5Concert hall, Padua (IT)4.725.5221.6
6Church, Venice (IT)4.268.3422.8
7Church, Rocca Pietore (IT)2.3710.4918.6
8Church, Rome, (IT)2.448.7217.0
9Church, Olkusz (PL)3.185.1215.9
10Church + concerts, Ponte in Valtellina (IT)3.527.3119.3
11Church, L’Aquila (IT)3.473.1615.0
12Church + concerts, Amsterdam (NL)2.871.3311.1
13Church + concerts, Oegstgeest (NL)2.416.3614.6
14Church + concerts, Hamburg (DE)3.193.1214.0
15Church + concerts, Leuven (BE)3.142.4913.2
16Church + concerts, Brussels (BE)5.373.4721.7
17Monastic Church, Padua (IT)2.654.2713.3
18Church basement, Venice (IT)4.028.7622.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Camuffo, D.; Della Valle, A.; Becherini, F. The European Standard EN 15757 Concerning Specifications for Relative Humidity: Suggested Improvements for Its Revision. Atmosphere 2022, 13, 1344. https://doi.org/10.3390/atmos13091344

AMA Style

Camuffo D, Della Valle A, Becherini F. The European Standard EN 15757 Concerning Specifications for Relative Humidity: Suggested Improvements for Its Revision. Atmosphere. 2022; 13(9):1344. https://doi.org/10.3390/atmos13091344

Chicago/Turabian Style

Camuffo, Dario, Antonio Della Valle, and Francesca Becherini. 2022. "The European Standard EN 15757 Concerning Specifications for Relative Humidity: Suggested Improvements for Its Revision" Atmosphere 13, no. 9: 1344. https://doi.org/10.3390/atmos13091344

APA Style

Camuffo, D., Della Valle, A., & Becherini, F. (2022). The European Standard EN 15757 Concerning Specifications for Relative Humidity: Suggested Improvements for Its Revision. Atmosphere, 13(9), 1344. https://doi.org/10.3390/atmos13091344

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop