Analysis of Ventilation Efficiency as Simultaneous Control of Radon and Carbon Dioxide Levels in Indoor Air Applying Transient Modelling

Abstract

The impact of ventilation efficiency on radon (²²²Rn) and carbon dioxide (CO₂) concentrations in the indoor air of a residential building was studied by applying transient data analysis within the CONTAM 3.4 program. Continuous measurements of ²²²Rn and CO₂ concentrations, together with basic meteorological parameters, were carried out in an apartment (floor area about 27 m²) located in Ljubljana, Slovenia. Throughout the experiment (October 3–15), frequent ventilation (several times per day), poor ventilation (once to twice per day) and no ventilation scenarios were applied, and the exact ventilation and occupancy schedule were recorded. Based on the measurements, a transient simulation of ²²²Rn and CO₂ concentrations was performed for six sets of scenarios, where the design ventilation rate (DVR) varied based on the ventilation requirements and recommendations. On the days of frequent ventilation, a moderate correlation between the measured and simulated concentrations (r = 0.62 for ²²²Rn, r = 0.55 for CO₂) was found. The results of the simulation indicated the following optimal DVRs: (i) 36.6 m³ h⁻¹ (0.5 air changes per hour, ACH) to ensure a CO₂ concentration below 1000 ppm and a ²²²Rn concentration below 100 Bq m⁻³; and (ii) 46.9 m³ h⁻¹ (0.7 ACH) to ensure a CO₂ concentration below 800 ppm. These levels are the most compatible with the 5C_Cat I (category I of indoor environmental quality, defined by EN 16798-1:2019) scenario, which resulted in concentrations of 656 ± 121 ppm for CO₂ and 57 ± 13 Bq m⁻³ for ²²²Rn. The approach presented is applicable to various types of residential buildings with high overcrowding rates, where a sufficient amount of air volume to achieve category I indoor environmental quality has to be provided. Lower CO₂ and ²²²Rn concentrations indoors minimise health risk, which is especially important for protecting sensitive and fragile occupants.

1. Introduction

A built environment is defined as a four-dimensional human-made space that ranges from indoor to outdoor and provides the setting for human activity [1,2]. As a primary health determinant, it accounts for almost 20% of all deaths in the WHO European Region that are attributable to a degraded urban environment and housing-related inequalities, among which poor air quality presents a major contribution [3]. However, indoor air is often more seriously polluted than outdoor air, even in the largest and most industrialised cities [4]. It may contain over 900 chemicals, including particles and biological materials with potential health effects [5].

The building is an individual component of the built environment that contributes positively or negatively to both built and natural environments [6]. The design of buildings, either residential or non-residential, should follow the morphology of their engineering design [7], which will be defined and shaped by the context of the human-environment relationships [6]. The main interface between the indoor and outdoor environment is the building envelope (i.e., with transparent and nontransparent parts of the walls, floor and roof), which enables a continuous transfer of heat, mass and information through a medium (solid, fluid or gas) [7]. Human interventions in the built environment must be sustainable and not cause environmental degradation, to prevent negative impacts upon the occupants’ health and well-being. Among them, ventilation is essential to ensure the breathable air is healthy, by diluting pollutants originating in the building and removing them [8].

Two significant pollutants regulated by international and national legislation in the built environment are radon (²²²Rn) and carbon dioxide (CO₂). Both can accumulate in the indoor air but are most often reduced rapidly with proper ventilation. The indoor/outdoor ratio (I/O) may vary for ²²²Rn from approximately 2 to over 100 [9], and for CO₂ from 2 to 10 [10,11].

Radon (²²²Rn, Rn) is a radioactive noble gas that accumulates in the indoor air of insufficiently ventilated buildings and may increase lung cancer risk. It is ranked as the second most common cause of lung cancer immediately after smoking [9,12]. ²²²Rn is primarily formed by the α-transformation of radium (²²⁶Ra) in the earth’s crust, from where it migrates towards the surface via diffusion and advection and exhales in the atmosphere. In general, outdoor radon concentrations are low (about 10 Bq m⁻³) [13], depending on the geological characteristics of the terrain and the atmospheric mixing state [14,15]. On the other hand, indoor radon concentrations are usually higher by one order of magnitude (up to several 100 Bq m⁻³) or even more (up to several 1000 Bq m⁻³). There are four possible sources of radon entering a building: (1) the soil beneath the building, if the building envelope is leaky in contact with the ground; (2) construction products containing radium (e.g., fly-ash bricks); (3) tap water, insofar as it is obtained from groundwater sources, such as springs, wells and boreholes, which generally have higher radon levels than surface waters (rivers, lakes and reservoirs); and (4) natural gas released into the air via combustion. By far, the most important is the first source. In buildings with elevated indoor radon concentrations, only mitigation measures can adequately reduce radon entry into the building. Otherwise, if the radon concentration is close to or slightly exceeds the reference value, adequate and regular ventilation (natural, mechanical or hybrid) can significantly reduce the radon levels.

CO₂ in indoor air is a metabolic product, a bio-effluent. It is the crucial indicator of room ventilation and a well-established measure of good indoor air quality (IAQ) [16,17]. Typical background CO₂ concentration in outdoor ambient air is 350 to 450 ppm [18], where the dominant factor for the emissions is fuel combustion. The indoor concentrations of CO₂ depend on the occupancy load, the room size, and the qualitative and quantitative ventilation characteristics. A range of 600 to 800 ppm of CO₂ provides reliable indoor air quality, with an upper limit of 1000 ppm. Concentrations above 1000 ppm can lead to an increase in absenteeism, lower attendance, and reduced productivity. The maximum workplace concentration over 8 h is 5000 ppm, and the critical, only short-term, exposure concentration range is 6000 to 30,000 ppm. The effects of different CO₂ concentrations [18] are an increased breathing frequency, headache (3–8%), nausea, vomiting, loss of consciousness (>10%), rapid loss of consciousness, and death (>20%).

Among the quantitative aspects of ventilation, the crucial parameter is the design ventilation rate (DVR), defined by legal requirements and/or recommendations [16,19,20,21,22,23,24,25]. The DVR can be determined as the amount of fresh air: (i) per floor area; (ii) per room volume; (iii) per occupant; and (iv) for a specific contaminant. The final selection of the DVR is, therefore, often left to the designer, who, intending to achieve the lowest possible ventilation heat losses, favours lower DVR values [26]. Persley [27] highlighted the problem that many practitioners and researchers claim a building has good IAQ because it complies with the 1000 ppm CO₂ limit set in the standard. Therefore, identifying relevant CO₂ concentrations that correspond to ventilation rate requirements must consider the building type and its occupancy, as well as other contaminants [27]. This aspect is essential, especially for the control and prevention of radon entry. For example, a Decree on the national radon programme [28] defines that ventilation is the primary measure required in all buildings where indoor radon concentrations are right below the 300 Bq m⁻³ (i.e., the reference level of the annual average indoor radon concentration in living and working spaces). In other buildings with radon concentrations above the reference level, it is necessary to set up an active radon ventilation system and seal the structural assemblies in contact with the ground to prevent radon entry.

The tightening of the building energy efficiency requirements, especially after 2010 [29,30], has been reflected in increased building airtightness as well as decreased DVRs in several ongoing construction projects [26]. As a result, such engineering measures might be related to a deterioration in indoor environmental quality. IAQ in energy-efficient residential and non-residential buildings has already been analysed by many authors [26,31,32,33,34,35,36,37,38,39]. The problem of increased indoor radon concentrations in renovated residential buildings has also been highlighted in several studies [32,35,36,37]. As reported, the concentration of ²²²Rn was increased from 17.5 to 49.6 Bq m⁻³ (11–33%) [36] and for 32 Bq m⁻³ (20%) [37] right after the building energy renovation. Similarly, the surveys of CO₂ in new and renovated residential [39] and non-residential buildings [33,38] showed an increase in CO₂ concentrations during occupancy to 2500 ppm [33,38] and 3000 ppm [39] (approximately 5–6 times, if the initial CO₂ concentration is about 500 ppm). The increased CO₂ concentrations were associated with lower ventilation rates, particularly in younger dwellings [40] that are naturally [40,41,42,43] or mechanically ventilated [43,44].

To provide an in-depth analysis of ventilation efficiency, some authors have included simulations of the selected indoor air pollutants in their studies and compared them to measurements. Concerning ²²²Rn, García-Tobar [45] proposed a methodology for estimating radon levels in a naturally and mechanically ventilated dwelling in a radon-prone area by using the CONTAM program. Further, García-Tobar [46] analysed the weather factors on indoor radon concentration in a new multistorey building in a radon-prone area. In the next study, García-Tobar [47] used CONTAM and computational fluid dynamics (CFD) transient simulations, including weather effects. Several authors have performed a transient simulation of CO₂ concentrations in residential buildings and compared them to measurements. Szczepanik-Ścisło and Flaga-Maryańczyk [44] focused on a bedroom in a passive house. Using the CONTAM tool, the influence of occupancy schedules and the ventilation efficiency on the CO₂ concentration was analysed over 10 days. According to the literature review, there has been an increased focus on the relationship between ²²²Rn or CO₂ and ventilation efficiency. However, to characterise IAQ and the effectiveness of ventilation, it is crucial to identify the relevant ²²²Rn and CO₂ levels simultaneously with those corresponding to ventilation rate requirements.

Our study focuses on the ventilation efficiency of a residential building. The primary purpose was to use a transient analysis of ²²²Rn and CO₂ concentrations simultaneously for the first time using the CONTAM 3.4 program [48]. Methodologically, our research was divided into four steps. In the first step, a ventilation zone based on an actual apartment was modelled. In the second step, measurements of ²²²Rn and CO₂ concentrations, together with basic meteorological parameters (air temperature, relative air humidity, barometric pressure), were conducted, and an accurate schedule of window opening was recorded. Based on the measurements, a model validation was carried out in the third step. In the fourth step, six sets of scenarios were critically analysed, defined by legal requirements and recommendations for the ventilation of residential buildings. Based on the findings, recommendations with practical benefits for constructions and renovations were developed, especially those where more efficient ventilation is sufficient as a radon protection measure.

2. Materials and Methods

2.1. Study Design

The study was conducted according to the steps below, which are explained in detail in the following sub-chapters.

• Selecting the measurement location for indoor (an apartment) and outdoor measurements (meteorological and air quality station);

• Defining the ventilation zone in the apartment;

• Determining the schedule for the ventilation of the apartment;

• Conducting the measurements of ²²²Rn and CO₂ concentrations and selected meteorological parameters (T–air temperature, RH–relative air humidity, P–barometric pressure);

• Simulating measured ²²²Rn and CO₂ concentrations in the air of the apartment by using the CONTAM 3.4 [48] program;

• Validating the model;

• Verifying six ventilation scenarios for ²²²Rn and CO₂ concentrations in the apartment.

2.2. Selection of Locations for Indoor and Outdoor Measurements

The study was conducted in Ljubljana (299 m above sea level, a.s.l.), the capital of Slovenia, located in the Ljubljana Basin in the central part of the county. It is characterised by a continental climate (Koppen–Geiger classification Cfb [49]) with an average minimum daily temperature of 5 °C and a maximum of 17 °C in October (time of measurements).

Two locations (one indoors and one outdoors) were selected for the measurements (at a distance of approximately 3 km from the city centre and approximately 2 km from each other):

• Indoor air measurements: A small apartment in an apartment building, part of a larger settlement in the city;

• Outdoor air measurements: The central meteorological and air quality station at the Environment Agency of Slovenia (ARSO).

2.3. Ventilation Zone

The ventilation zone was modelled according to the dimensions of an actual apartment in the apartment building. The building is a part of a larger settlement of apartment buildings and terraced houses built in 2002. In the basement, below the entire surface of the settlement, there is a garage with parking lots, which has a mechanical ventilation system installed. The apartment has a net size of 4.51 m × 6.33 m (26.6 m² of net floor area, A_u), with a height of 2.60 m (69.3 m³ of conditioned volume, V_e). It faces east and is located on the 3rd floor (Figure 1) of a three-storey apartment building. The exterior wall assembly consists of reinforced concrete (16 cm) and facade plaster. The apartment is naturally ventilated by two French doors, with dimensions of 2.25 m × 2.70 m. Additional ventilation is possible through the kitchen hood and bathroom fan. Heating is based on a gas central heating boiler. The geometry of the ventilation zone with the position and dimensions of the openings is consistent with the actual apartment. The occupational load is 1.

2.4. Ventilation Schedule

The ventilation schedule of the apartment (with the day of the week, date, absence of occupant, and ventilation duration) is presented in Table 1. Throughout the experiment, only one door, the same French door, was open in full-screen mode (on the left side from the entrance). During the periods without ventilation, all of the French doors were closed, and the door to the bathroom was open (Figure 1). In addition, the kitchen hood and bathroom fan were not used during the measurement period.

In the first part of the measurement period, October 3–8, the schedule for the window opening (i.e., frequency and ventilation duration) was adjusted to maintain the CO₂ concentrations below 1000 ppm. In the second part, October 9–10 (weekend), on Saturday, the dwelling was not ventilated and on Sunday, the previous ventilation regime was applied. In the last part, October 11–15, the ventilation was minimised to twice per day (Monday) and once per day (Tuesday and Thursday), with no ventilation on Wednesday.

2.5. Measurements

The measurements were conducted in the period 3–15 October 2021, and all presented data are reported in local time (LST = UTC + 2 h). A standardised protocol for characterising IAQ in residential buildings was followed [16,17,19,25,50]. In the apartment, the instrument for continuous measurement of the selected parameters was placed in the respiratory zone (living zone) at the height of 1.1 m above the floor; 3 m from the external window and wall, door and radiator; and 0.8 m from the internal wall (Figure 1) [16,17,19,25,50].

The selection of instruments was based on the expected radon (²²²Rn) concentrations in indoor and outdoor air and the requirements of our radon laboratory [51], accredited according to ISO/IEC 17025 [52]. Both devices were operated continuously in a diffusion mode with a frequency of once per hour.

Indoor air: radon C_Rn-in [Bq m⁻³] and carbon dioxide $C_{{\mathrm{CO}}_2}$ [ppm] concentrations, room air temperature T_in [°C] and relative air humidity RH_in [%] were measured with the Radon Scout Professional device (Sarad). The Radon Scout Professional monitor operates in the range from 0 Bq m⁻³ to 2 MBq m⁻³ with the sensitivity to Rn > 2.5 cpm/(kBq m⁻³). The integration interval of the data should be adjusted to the concentration range. If the expected radon concentrations are of the order of the reference level of 300 Bq m⁻³ or below, an interval of 60 min should be used. The sensor for CO₂ operates in the range of 400 to 5000 ppm [53]. The integrated CO₂ sensor uses the non-dispersive infrared (NDIR) operational principle.

Outdoor air: radon C_Rn-out [Bq m⁻³] concentration, temperature T_out [°C], relative humidity RH_out [%], and pressure P_out [hPa] were measured with the AlphaGUARD (Bertin Instruments) monitor, placed into a Stevenson screen at a height of 1.5 m above the ground. The instrument operates in the range from 2 Bq m⁻³ to 2 MBq m⁻³, and the efficiency of the detector is 1 cpm at 20 Bq m⁻³ [54].

2.6. Simulation

The simulation was based on the CONTAM 3.4 program [48]. This is a multizone analysis program, designed to analyse the IAQ in relation to the selected contaminants, ventilation rates, and the effectiveness of ventilation. According to the net dimensions of our test apartment, one ventilation zone was modelled. Openings (French doors, interior door) for natural ventilation were considered as airflow paths in our model.

Conservation of mass was applied to the zone, leading to a set of nonlinear algebraic equations that must be solved interactively. The detailed calculation protocol is presented in the CONTAM user guide [48]. The selected type for our analysis was transient and followed all of the required steps presented in the work by García-Tobar [45,46].

The input data in our model are as follows:

• i.. Airflow paths: one-way flow using power law for French door and two-way flow for the indoor door (type of model); orifice area data for French door and one opening for the interior door (selected formula); 13,500 cm² for French door and 20,000 cm² for the interior door (cross-sectional data); 1.3111 cm for French door (hydraulic diameter); 30 for French door (Transition Reynolds number); 0.78 for French door and 0.78 for the interior door (discharge coefficient); 0.5 for French door and 0.5 for the interior door (flow exponent). The program enables a simultaneous mass balance of air in the ventilation zone to determine zonal pressures and airflow rates through each airflow path.

• ii.. Measured data in outdoor air (hourly weather data, [55]): radon concentration C_Rn-out [Bq m⁻³], temperature T_out [°C], relative humidity RH_out [%], pressure P_out [hPa], and wind speed v_w [m s⁻¹].

• iii.. Measured data in indoor air: radon concentration C_Rn-in [Bq m⁻³], carbon dioxide concentration $C_{{\mathrm{CO}}_2}$ [ppm], temperature T_in [°C], and relative humidity RH_in [%].

• iv.. Default data: the radon generation rate [Bq h⁻¹] was determined for every hour according to the methodology defined in Dovjak et al. [26]. The CO₂ metabolic emission rate is 0.0027 dm³ s⁻¹ during sleeping and 0.0038 dm³ s⁻¹ during light activity [56]. The outdoor CO₂ concentration is 400 ppm. Uncontrolled ventilation is 0.1 air changes per hour, ACH (6.9 m³ h⁻¹).

• v.. Defined schedules: the ventilation schedule of the apartment and the presence of the occupant were determined according to the records (Table 1).

• vi.. Defined type of calculation: transient calculation of airflows and concentrations of ²²²Rn and CO₂. The ²²²Rn and CO₂ concentrations were determined from predefined indoor and outdoor sources. The main characteristics of ²²²Rn are an atomic weight of 222 kg kmol⁻¹, a diffusion coefficient in the air of 5.91 mm² s⁻¹, and a half-life of its α-transformation of 3.8 days [45]. The main characteristics of CO₂ are an atomic weight of 44 kg kmol⁻¹ and a diffusion coefficient in the air of 20 mm² s⁻¹. Airflow and contaminants information are then used to determine the ²²²Rn and CO₂ concentrations within the zone.

2.7. Ventilation Scenarios

The simulation was performed for 6 different sets of ventilation scenarios, where the DVR was changed according to the legal requirements and recommendations (Table 2). Scenarios 1, 2, 3 and 4 are based on the requirements of the rules relating to the ventilation and air conditioning of buildings [19]. Scenarios 5-I, 5-II, 5-III, 5-IV are based on the recommendations of the standard SIST EN 16798-1: 2019 [25], where all four categories of indoor environment quality (I–IV) were considered and applied to residential buildings. Scenario 6 is based on the Proposal of Rules for efficient use [22] and the Proposal of TSG-1-004: 2021 [23].

The calculated concentrations of ²²²Rn and CO₂ for all of the variants were compared with the legal requirements and recommendations presented in Table 3.

So far, the Federation of European Heating, Ventilation and Air Conditioning Associations (REHVA) has also prepared the ventilation guidelines to prevent the spread of SARS-CoV-2 in workplaces [57]; the guidelines for residential buildings have not yet been prepared.

3. Results

3.1. Results of Measured ²²²Rn and CO₂ Concentrations and Meteorological Parameters

The results of the measurements are presented in Figure 2 for the entire period, 3–15 October 2021. Figure 2a shows the outdoor radon concentration (C_Rn-out) and air temperature (T_out); Figure 2b shows the indoor radon concentration (C_Rn-in) and the temperature difference between the indoor and outdoor air (∆T, T_in − T_out); and Figure 2c shows the indoor carbon dioxide concentration ($C_{{\mathrm{CO}}_2}$).

Outdoor ²²²Rn concentrations (Figure 2a) range from 3.3 to 39 Bq m⁻³ with the average and standard deviation of 13.7 ± 7.0 Bq m⁻³. A typical daily run, with the highest concentrations in the early morning and the lowest in the afternoon, is not always pronounced. The outdoor temperature (range 4.1–24.4 °C and average 12.5 ± 4.1 °C) decreases rapidly from the beginning to the end of the measurement period. It rarely drops below 14 °C in the first days, hovers around 14 °C in the next two days (7–8 October), and is mostly below 14 °C in the last days (9–15 October). The correlation of C_Rn-out with T_out is weakly negative (r = 0.34), and of C_Rn-out with the pressure time gradient (∆P/∆t) in the hourly scale is very weakly negative (r = 0.09). A high correlation was not expected because outdoor ²²²Rn concentration is a sum of local (exhalation from the ground) and synoptic (remote) sources [15]. The contribution of each source was not sought because, in this study, only the outdoor ²²²Rn concentration during the ventilation of the apartment was needed for the simulations. Figure 2a does not show the relative air humidity (range 46–94% and average 78 ± 12%).

Due to a relatively low indoor ²²²Rn concentration (range 5–151 Bq m⁻³ and average 57 ± 30 Bq m⁻³) and the lower sensitivity of the instrument (the average error of a single measurement is ±32%), the hourly values fluctuate significantly (Figure 2b). A longer integration time of the measurements (e.g., 3 h) would give a smoother curve, but less information about the decrease in ²²²Rn concentration due to the ventilation of the apartment. The indoor ²²²Rn concentration and the temperature difference T_in − T_out show a weak positive correlation (r = 0.32) for the entire measurements.

In the first part of the measurements (3–8 October), when the door was opened to maintain the CO₂ concentration below 1000 ppm, the ²²²Rn concentration also remained below 100 Bq m⁻³ for most of the time (range 5–149 Bq m⁻³ and average 46 ± 23 Bq m⁻³). In the second part (9–10 October, weekend), no ventilation on Saturday and the previous ventilation regime on Sunday were applied, and the following indoor ²²²Rn concentrations were obtained: Saturday (9 October) 32–141 Bq m⁻³ (80 ± 23 Bq m⁻³) and Sunday (10 October) 5–91 Bq m⁻³ (42 ± 23 Bq m⁻³). In the last part (11, 12, 14 October), when the ventilation was minimised to once or twice per day, ²²²Rn concentrations in the range of 14–123 Bq m⁻³ and an average of 66 ± Bq m⁻³ were obtained. On 13 October, the apartment was not ventilated and ²²²Rn concentration in the range of 36–151 Bq m⁻³ (93 ± 32 Bq m⁻³) was obtained, which is similar to October 9 when the apartment was also not ventilated (Figure 2b). The average ²²²Rn concentration in the non-ventilated apartment (87 Bq m⁻³) dropped by about 25% (66 Bq m⁻³) when ventilated once to twice per day, and by about 50% (45 Bq m⁻³) when ventilated more frequently.

The CO₂ concentrations (Figure 2c) range from 400 to 2340 ppm with the average and standard deviation of 1010 ± 490 ppm for the entire period of measurements. Similar to indoor ²²²Rn concentration, indoor CO₂ concentration also fluctuates according to the frequency of the ventilation (Figure 2c). In the first part (3–8 October), when the apartment was ventilated three to five times per day (except Friday, when it was not occupied for 9.5 h) to keep the CO₂ concentration below 1000 ppm, the $C_{{\mathrm{CO}}_2}$ was 420–1490 ppm (average 759 ± 222 ppm). In the second part (9–10 October, weekend), under closed conditions on Saturday and frequent ventilation on Sunday, the following CO₂ concentrations were observed: 1435–2220 ppm (average 1860 ± 260 ppm) on Saturday (9 October) and 410–2070 ppm (average 812 ± 338 ppm) on Sunday (October 10). In the last part (11, 12, 14 October), during minimal ventilation (once to twice per day), CO₂ concentration in the range of 400–2045 ppm (average 1033 ± 424 ppm) was obtained. On 13 October, when no ventilation was performed, the range of 1120–2340 (average 1800 ± 311 ppm) was recorded, similar to 9 October.

The indoor air temperature (range 17.7–25.9 °C, average 22.4 ± 1.6 °C) and indoor air humidity (range 42–59%, average 50 ± 5%) are not presented in Figure 2. Although they are influenced by ventilation, they have less impact on the results as the temperature difference between indoor and outdoor air.

3.2. Comparison of Measured and Simulated ²²²Rn and CO₂ Concentrations

Figure 3a shows a comparison of the measured and simulated datasets of ²²²Rn concentrations (C_Rn-m, C_Rn-s) and a similar trend of both curves is observed. In the first part (3–8 October), when the apartment was ventilated several times per day, a moderate correlation between the measured and simulated datasets is obtained (r = 0.62). In the second part (9–10 October, weekend), there was no ventilation on Saturday (exhibiting moderate correlation, r = 0.59) and frequent ventilation on Sunday (with slight correlation, r = 0.32). In the last part (11, 12, 14 October), the ventilation was done once or twice per day (moderate correlation, r = 0.68) and 13 October was without ventilation (moderate correlation, r = 0.47). The differences between C_Rn-m and C_Rn-s are as follows: 12 ± 20 Bq m⁻³ (3–8 October); 1 Bq m⁻³ ± 20 Bq m⁻³ (9 October); 21 ± 22 Bq m⁻³ (10 October); 9 ± 19 Bq m⁻³ (11 October); 14 ± 33 Bq m⁻³ (12 October); 14 ± 33 Bq m⁻³ (14 October); and 79 ± 29 Bq m⁻³ (13 October).

The measured and simulated concentrations of CO₂ ($C_{{\mathrm{CO}}_2-\mathrm{m}}{, C}_{{\mathrm{CO}}_2-\mathrm{s}}$) are shown in Figure 3b. During the days of frequent ventilation (3–8 October) a moderate correlation is obtained (r = 0.55). In a closed condition with no ventilation (9 October), a very high correlation (r = 0.94) is observed, and in well ventilated (seven times) conditions (10 October) a moderate correlation (r = 0.69) is observed. In a poor ventilated condition (11, 12, 14 October, a moderate correlation (r = 0.55) is observed, and in a no ventilated condition (13 October), very weak negative correlation (r = −0.06) is observed. The difference between $C_{{\mathrm{CO}}_2-\mathrm{m}}$ and $C_{{\mathrm{CO}}_2-\mathrm{s}}$ is as follows: 151 ± 110 ppm (3–8 October, without night time, when the highest discrepancy was observed); 252 ± 164 ppm (9 October); 107 ± 134 ppm (10 October); 107 ± 93 ppm (11 October); 51 ± 181 ppm (12 October); 837 ± 278 ppm (14 October); and 281 ± 375 ppm (13 October). The exact ventilation data referred to in the above data are summarised in Table 2.

3.3. The Influence of Required and Recommended DVRs on Simulated ²²²Rn and CO₂ Concentrations

The simulated concentrations are shown in Figure 4a for ²²²Rn (C_Rn-s), and in Figure 4b for CO₂ ($C_{{\mathrm{CO}}_2-\mathrm{s}}$) by varying the DVRs in the apartment (Table 3) within six sets of scenarios. Deviations (expressed in h and % of simulated time, 288 h in total) from the limit values of 100 Bq m⁻³ for ²²²Rn concentration [25,59], and 1000 ppm [60] and 800 ppm for CO₂ concentration [17], are presented in Table 4.

The best scenario is represented by the DVR in case 5C_Cat I (46.9 m³ h⁻¹ (0.7 ACH)), recommended by EN 16798-1 [25]. For this case, the simulated ²²²Rn and CO₂ concentrations were below the limit values (100 Bq m⁻³, 1000 and 800 ppm) for the entire simulation (the total 288 h). The worst scenarios are represented by the DVRs in the 1st and 2nd case (13.9 m³ h⁻¹ and 15.0 m³ h⁻¹ (0.2 ACH)), required by the rules relating to the ventilation and air conditioning of buildings [19]. In the case of 13.9 m³ h⁻¹ (1st case), the simulated CO₂ concentrations exceeded 1000 ppm 64% of the time (185 h), and 800 ppm 82% of the time (237 h); the simulated ²²²Rn concentration exceeded 100 Bq m⁻³ 4% of the time (10 h). In the case of 15.0 m³ h⁻¹ (0.2 ACH), simulated CO₂ concentration exceeded 1000 ppm 61% of the time (176 h), and 800 ppm 93% of the time (267 h); the simulated Rn concentrations exceeded 100 Bq m⁻³ 2% of the time (6 h). A similar deviation was also found for cases 5A_Cat III and 5B_Cat III, both representing categories III of indoor environmental quality (IEQ) defined by EN 16798-1 [25]. The DVRs in the scenarios (i.e., 3, 6, 5A_ Cat I, 5B_Cat I, 5C_Cat II) resulted in a simulated ²²²Rn concentration below the limit values of 100 Bq m⁻³ and a CO₂ concentration below 1000 ppm, but not below 800 ppm.

In the last step of the study, the optimal DVRs were simulated. As can be seen in Figure 5a, they permanently assure a concentration of ²²²Rn below 100 Bq m⁻³, and concentrations of CO₂ below 1000 ppm (Figure 5b), and below 800 ppm (Figure 5c).

4. Discussion

The deterioration of IAQ in residential buildings is a subject of numerous studies worldwide. One of the main features of energy-efficient buildings is increased airtightness, which leads to lower air leakage through the building envelope [61]. A controlled infiltration rate adjacent to overlooked building ventilation might result in elevated indoor air pollutant concentrations [62]. Therefore, it is not surprising that recent research has highlighted the dependence of indoor air pollutant concentrations on ventilation, either based on measured or modelled data or as a synthesis of both [63].

Dealing with ²²²Rn and CO₂, the authors study the effects of ventilation on the measured and simulated concentrations separately, either on ²²²Rn or CO₂. In our research, we evaluated the results of measurements and simulations of both pollutants simultaneously (Figure 3).

In the indoor air of the test apartment in our study, the measured average ²²²Rn concentration of 57 ± 30 Bq m⁻³ (range 5–151 Bq m⁻³) is about 4-fold higher than in outdoor air (13.7 ± 7.0 Bq m⁻³, 3.3–39 Bq m⁻³) (Figure 2a,b). The primary indoor ²²²Rn source is assumed ²²²Rn diffusion from building materials [64,65]. The apartment has one outdoor wall, two walls border the staircase, and one adjoins the next apartment. A minor source of ²²²Rn could be attributed to infiltration through the walls from the apartment next to and below, and less to infiltration from the staircase, where the window is open most of the time. Similar to our results, the average measured ²²²Rn concentrations obtained by García-Tobar [45] were: 62 ± 17 Bq m⁻³ (dwelling A) and 77 ± 20 Bq m⁻³ (dwelling B), with the highest values in the hallway and bathroom (115 Bq m⁻³ in dwelling A, and 150 Bq m⁻³ in dwelling B), mainly due to lower ventilation rates. Similar ²²²Rn concentrations were also obtained in the following study (62 ± 5 Bq m⁻³ in dwelling A, 78 ± 6 Bq m⁻³ in dwelling B) in which García-Tobar [46] analysed the weather factors on indoor ²²²Rn concentration. The most significant impacts on indoor Rn concentrations are reported to be wind speed and wind direction, followed by air temperature and barometric pressure. On the contrary, in our study, the dominant factor that affected ²²²Rn and CO₂ concentrations in the non-ventilation period was the difference between the indoor and outdoor air temperature. The other parameters, such as the wind speed and direction, were not analysed.

In our study, the indoor CO₂ concentration reflects the resident presence. At the time of the experiment, only one person was present, according to the records in Table 2. Considering that during the 12-day measurements (Figure 2c), the apartment was intensively ventilated for 6 days, poorly ventilated for 4 days and not ventilated for 2 days, the average concentration of 1010 ± 490 is relatively high (range 400–2340 ppm). Similar findings have been reported for naturally [40,41,42,43] and mechanically ventilated buildings [43,44]. A cross-sectional study in 79 Greenlandic dwellings in the town of Sisimiut found that 73% of bedrooms were insufficiently ventilated ($C_{{\mathrm{CO}}_2}$ > 1000 ppm), that younger dwellings (built after 1990) had poorer IAQ than older dwellings, and that children’s bedrooms (2000–4000 ppm) had higher CO₂ concentrations than adults’ bedrooms [40].

Our results showed that the actual ventilation rate (controlled by ventilation; uncontrolled by infiltration, breakthroughs, shafts) significantly influences the accuracy of the simulation. When the accurate occupancy schedule and ventilation rates were included in the model, the measured and simulated data of ²²²Rn and CO₂ concentrations were well-matched (Figure 3).

A comparison of the measured and simulated time series of ²²²Rn concentrations (Figure 3a) shows a moderate correlation of frequent ventilation days (r = 0.62), sparse ventilation days (r = 0.68), and no ventilation days (r = 0.59 and r = 0.47). The slight correlation (r = 0.32) on 10 October, when the maximum number of ventilations was performed, revealed a too low ²²²Rn monitor sensitivity for this concentration range, which resulted in the unreliable ²²²Rn generation rates used in the simulation. In general, the time series of simulated ²²²Rn concentration is underestimated, except during the last 3 days (12–14 October), when it is significantly overestimated. In the last 3 days, the outside air temperature dropped sharply overnight (Figure 2a), leading to a temperature difference ∆T between the outdoor and indoor air (Figure 2b) sufficiently high to trigger the so-called ‘chimney effect’, a natural draft of air through the chimney, and thus decreased the indoor ²²²Rn concentration. This type of (uncontrolled) ventilation is not considered in our model and, therefore, leads to the overestimation in the simulation in this period (Figure 3a). In the study by García-Tobar [45], ²²²Rn data are compared for: (i) average measured and simulated values (62/66 Bq m⁻³ in dwelling A, and 77/74 Bq m⁻³ in dwelling B); and (ii) the highest measured and simulated values (115/110 Bq m⁻³ in dwelling A, and 150/110 Bq m⁻³ in dwelling B). When the flow rate of the fan in the bathroom was doubled (from 2 to 4 dm³/(s m⁻²)), the ²²²Rn concentration was reduced by 50% (from 110 to 55 Bq m⁻³) and the following measured and simulated average ²²²Rn concentrations obtained: 62 ± 5 Bq m⁻³/75 ± 3 Bq m⁻³ (dwelling A) and 78 ± 6 Bq m⁻³/83 ± 4 Bq m⁻³ (dwelling B), with the Pearson correlation coefficients of 0.341 (dwelling A) and 0.198 (dwelling B).

The measured and simulated concentrations of CO₂ in our study (Figure 3b) show moderate correlation with the same correlation coefficient (r = 0.55) during the days of frequent ventilation and poor ventilation. A very high correlation (r = 0.94) was obtained on 9 October under no ventilation, and the resident was present all day. On the other hand, under a no ventilation condition, on 13 October, a very weak negative correlation (r = −0.06) was noticed, most probably due to the chimney effect not considered in our simulation (uncontrolled ventilation). The curve of $C_{{\mathrm{CO}}_2-\mathrm{s}}$ is underestimated during the whole dataset, and the discrepancy is more pronounced in the last part. The reason for the sudden overestimation on the last day has not been understood yet. A similar study carried out by Szczepanik-Ścisło and Flaga-Maryańczyk [44] also performed a CONTAM 3.2 simulation, where the influence of occupancy schedules and the ventilation efficiency on the CO₂ concentration was analysed over 10 days. In the first case, when the simulations were conducted at the minimum ventilation rate and the door to the wardrobe was open, the simulation data were mainly higher than the measurement (300–800 ppm). In the second case, with the minimum ventilation rate and a closed door, the difference between the measurement and the simulation data was 600–1000 ppm. In the following case, with medium ventilation and both doors open, the measurement data were higher than those of the simulation (100–600 ppm) in the first and the last part of the period, and lower (approximately 400 ppm) in the middle part of the period. In the last case, which included the occupancy schedule and the real ventilation rate, the simulated values were 100–300 ppm lower than the measured ones. Regarding the results, the authors proved that the CONTAM 3.2 program was able to recreate the conditions of CO₂ inside the analysed room, especially if the occupancy schedules and real ventilation rate (with all openings) were precisely considered in the simulation.

Studies that simultaneously address ²²²Rn and CO₂ are rare, especially on the critical analysis of the DVR. Our study analysed the influence of required and recommended DVRs on CO₂ and ²²²Rn concentrations with six sets of scenarios. The best scenario is the DVR in case 5C_Cat I (46.9 m³ h⁻¹ (0.7 ACH)), recommended by EN 16798-1 [25], which resulted in the lowest CO₂ (656 ± 121 ppm) and ²²²Rn concentrations (57 ± 13 Bq m⁻³). The worst scenarios are the DVRs in the 1st and 2nd case (13.9 m³ h⁻¹ and 15.0 m³ h⁻¹ (0.2 ACH)), required by the rules relating to the ventilation and air conditioning of buildings [19], as well as the DVRs in cases 5A_Cat III and 5B_Cat III (15.0 m³ h⁻¹; 14.4 m³ h⁻¹ (0.2 ACH)), defined by EN 16798-1 [25]. For example, scenario two with 15 m³ h⁻¹ resulted in CO₂ (1159 ± 291 ppm) and ²²²Rn concentration (59 ± 21 Bq m⁻³), and scenario 5B_Cat III with 14.4 m³ h⁻¹ (0.2 ACH) resulted in CO₂ (1188 ± 300 ppm) and ²²²Rn concentration (61 ± 21 Bq m⁻³). This was the case for the scenarios defined for the III and IV categories of IEQ. 5A_Cat III and 5B_Cat III resulted in concentrations exceeding the limits for C_Rn-s (6 h above and 8 h above 100 Bq m⁻³) and $C_{{\mathrm{CO}}_2-\mathrm{s}}$ (1159 ± 291 ppm and 1188 ± 300 ppm). Moreover, category III of IEQ is also defined by the national Proposal of Rules on efficient use [22] and Proposal of TSG-1-004:2021 [23].

Particularly in residential buildings with a lower net floor area per occupant, and therefore, lower DVRs (if defined per floor area), higher CO₂ emissions result. Several authors have pointed out the same problem. Bekö et al. [42] inspected Danish homes and found that in 57% of new dwellings, the ventilation rate is lower than the minimum required 0.5 h⁻¹, in 32% of bedrooms, an average $C_{{\mathrm{CO}}_2}$ < 1000 ppm during the measured nights; in 23% of rooms, at least 20 min during the night a $C_{{\mathrm{CO}}_2}$ > 2000 ppm; and in 6% of rooms a $C_{{\mathrm{CO}}_2}$ > 3000 ppm, which is similar to the findings of Kotol et al. [40]. The lower ventilation rate problem was also addressed in a detached passive house in the Silesian region of Poland [44]. The stated reason for the increased CO₂ concentration in the bedroom (peak 1800 ppm) was the minimised ventilation rate to reduce the noise of the ventilation system. Additionally, Sekkhar and Goh [43] conducted a study in naturally and mechanically ventilated and air-conditioned bedrooms in Singapore’s hot and humid climate. Higher CO₂ concentrations in their study were related to the use of split-system air-conditioning units that only recirculate air and do not provide fresh outdoor air. Increased DVRs, however, resulted not only in the reduction in CO₂ but also in the reduction in ²²²Rn concentrations. In the García-Tobar [45] study, doubling the fan’s flow rate in the bathroom reduced the ²²²Rn concentration by 50%. The significant reduction in ²²²Rn concentration due to increasing the DVR was also the case in our study.

This problem needs to be understood in the broader context of housing policy. For example, in Slovenia, according to the Statistical Office of the Republic of Slovenia [66] for 2018, 38% of housing consisted of small apartments with a net floor area of less than 50 m² (i.e., studio, one- and two-bedroom apartments). Two-thirds of the population of Slovenia reside in one- or two-bedroom apartment buildings, and, therefore, the percentage of people living in underoccupied apartments is 30.4% (in 2018). Compared to the EU-27, the percentage is slightly higher, at 33.0%, and the most exposed populations are young people and children [67]. Smaller apartments are, thus, the most critical in terms of CO₂, particularly in the presence of overcrowding. However, following the Energy Performance of Buildings Directive [29,30], mechanical ventilation with heat recovery has become an unavoidable element of a nearly zero energy building (NZEB); and according to the current construction practice [26], they often operate on too low DVRs. Consequentially, poor IAQ might be in contradiction to the objectives of the Resolution on the National Housing Programme 2015–2025 [68] for quality and functional housing. ²²²Rn can further impair IAQ, especially in low-floor apartments directly connected to the ground. On upper floors, ²²²Rn usually does not pose any significant health threat. However, more frequent ventilation, which lowers the CO₂ concentration, also lowers the ²²²Rn concentration and thus, further reduces the health risk.

Based on the results of our research, optimal DVRs are proposed that can be regulated by demand-controlled ventilation using sensory information [69] or as an automatic control system [26], which may also include other health risks of chemical or biological origin [70]. As simulated in our study (Figure 5), to ensure CO₂ below 1000 ppm and ²²²Rn below 100 Bq m⁻³, permanent ventilation of at least 36.6 m³ h⁻¹ (0.5 ACH) is required. To ensure CO₂ below 800 ppm, the DVR must always be at least 46.9 m³ h⁻¹ (0.7 ACH).

The simultaneous analysis of measured and simulated ²²²Rn and CO₂ data highlights the benefits in data evaluations compared to studies of either ²²²Rn or CO₂, as in our previous work [10,11,26]. However, the most beneficial evaluation was the transient analysis on measured and simulated ²²²Rn and CO₂ datasets, which we used for the first time. There are some further steps foreseen in our future research, which might improve the transient analysis of the datasets: (i) to add meteorological parameters, e.g., wind speed and direction, precipitation, etc.; (ii) to provide detailed information on the infiltration rate (which was considered in a simplified manner and defined as a constant value, as the primary purpose was to examine the quantitative requirements for DVRs with the controlled ventilation); and (iii) to prolong the measurement duration to all yearly seasons (in the present study, simulations were based on a short period, in which natural ventilation was executed to examine only the characteristic of high, moderate and no ventilation periods). In particular, our future work will upgrade the methodology for characterizing the IAQ with measurements and simulations in other buildings types, such as family houses and non-residential buildings. The effectiveness of natural and mechanical ventilation will also be evaluated, and the optimal DVRs examined.

5. Conclusions

Based on the measurements, ²²²Rn and CO₂ concentrations in a small apartment were simulated by using the CONTAM 3.4 program and six sets of the DVR scenarios modelled. The optimal DVRs, which permanently assure ²²²Rn concentrations below 100 Bq m⁻³, and CO₂ concentrations below 1000 ppm and 800 ppm, were sought. The main findings of the research are as follows:

• A comparison of measured and simulated time series of ²²²Rn and CO₂ concentrations shows a moderate correlation (r = 0.62 for ²²²Rn and 0.55 for CO₂) during the days of frequent ventilation, which was our main focus in the study.

• A critical analysis of six sets of ventilation scenarios showed that the optimal DVR values were those defined as the maximum amounts of fresh air, determined per floor area and per person, and applied for category I of indoor environmental quality (for test apartment: 5C_Cat I with 46.9 m³ h⁻¹ (0.7 ACH) that resulted in 656 ± 121 ppm, 57 ± 13 Bq m⁻³). Lower DVR values, especially those defined for categories III or IV of IEQ (5A_Cat III with 15.0 m³ h⁻¹ and 5B_Cat III with 14.4 m³ h⁻¹ (0.2 ACH)), resulted in CO₂ and ²²²Rn concentrations above the limit values (CO₂: 1159 ± 291 ppm and 1188 ± 300 ppm; ²²²Rn: 6 h above and 8 h above 100 Bq m⁻³), which can present a problem for buildings located in Zone 3 areas.

• To increase the accuracy of our simulation, a more extended time series of measured data is needed, which should include all seasons of the year.

• Although the measured and simulated data matched relatively well, uncontrolled air infiltration through the building envelope needs to be further studied and defined to improve the model.

• The approach presented in our study can be applied to various building types to determine the optimal DVR values for ventilation. However, special attention should be paid to small apartments which, in the EU, have a high overcrowding rate of 33.0% (in Slovenia 30.4%). Accordingly, to protect sensitive and fragile occupants, a sufficient amount of fresh air volume for category I of indoor environmental quality has to be provided in terms of CO₂. In addition, by lowering the CO₂ concentration, the ²²²Rn concentration is also reduced, thus minimising the health risk.

• Our findings might be implemented in national legislation and the existing construction practice, which will result in safer and healthier indoor environments.

Author Contributions

Conceptualization, M.D. and J.V.; methodology, M.D. and J.V.; measurements, J.V.; software, O.V. and M.D.; validation, M.D., J.V., O.V.; writing—original draft, M.D. and J.V.; writing—review and editing, M.D., J.V., O.V. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the financial support from the Slovenian Research Agency (research core funding No. P2-0158, Structural engineering and building physics; No. P1-0143, Cycling of substances in the environment, mass balances, modelling of environmental processes and risk assessment and research project No. J1-1716; Sources, transport and fate of persistent air pollutants in the environment of Slovenia–STRAP). They also acknowledge the financial support of project No. C3330-19-952056, Development of research infrastructure for the international competitiveness of the Slovenian RRI space–RI-SI-EPOS. The operation is co-financed by the Republic of Slovenia, the Ministry of Education, Science and Sport, and the European Union from the European Regional Development Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACH	Air changes per hour
A u	Net floor area
CO2	Carbon dioxide
$C_{{\mathrm{CO}}_2}$	Carbon dioxide concentration
$C_{{\mathrm{CO}}_2-\mathrm{m}}$	Measured carbon dioxide concentration
$C_{{\mathrm{CO}}_2-\mathrm{s}}$	Simulated carbon dioxide concentration
C Rn	Radon activity concentration (in text generally without activity)
C Rn-in	Indoor radon concentration
C Rn-out	Outdoor radon concentration
C Rn-m	Measured radon concentration
C Rn-s	Simulated radon concentration
DVR	Design ventilation rate
IAQ	Indoor air quality
IEQ	Indoor environmental quality
I/O	Indoor/outdoor ratio
NOX	Nitrogen oxides
O3	Ozone
PM2.5	Particulate matter with an aerodynamic diameter smaller than 2.5 μm
P out	Outdoor barometric pressure
RH in	Indoor relative air humidity
RH out	Outdoor relative air humidity
222Rn (Rn)	Radon isotope
T out	Indoor air temperature
T in	Outdoor air temperature
t	Time spent in dwelling
V e	Conditioned volume
vw	Wind speed
∆T	Temperature difference between indoor and outdoor air

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Floor plan of the tested apartment ventilation zone (red dot indicates the location of the instruments).

Figure 2. Results of the measurements in the period 3–15 October 2021: (a) radon concentration (CRn-out) and air temperature (Tout) outdoors; (b) the radon concentration (CRn-in) indoors and the temperature difference between the indoor and outdoor air (∆T = Tin − Tout); (c) the carbon dioxide concentration ([Forumla omitted. See PDF.]) indoors. The blue regions in (b,c) indicate the ventilation periods, and the blue lines in (b) the Rn limit according to WHO recommendations [59] and the CO2 limit according to [60], respectively. The solid lines indicate midnight, and the broken lines indicate noon in the gridlines.

Figure 3. Measured and simulated concentrations of: (a) radon CRn-m and CRn-s [Bq m−3]; and (b) carbon dioxide [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] [ppm] in the period 3–15 October 2021. The solid lines indicate midnight and the broken lines indicate noon in the gridlines.

Figure 4. Simulated concentrations of: (a) radon, CRn-s [Bq m−3]; and (b) carbon dioxide ([Forumla omitted. See PDF.]) [ppm] in the apartment for 6 sets of scenarios in the period October 3–15, 2021. The blue lines indicate the Rn limit according to WHO recommendations [59] and the CO2 limits according to [60] and [17], respectively. The solid lines indicate midnight and the broken lines indicate noon in the gridlines.

Figure 5. Optimal DVR in the apartment (3–15 October 2021), based on the criteria: (a) CRn-s < 100 Bq m−3; (b) [Forumla omitted. See PDF.] < 1000 ppm; and (c) [Forumla omitted. See PDF.] < 800 ppm. In the gridlines, the solid line indicates midnight and the broken line is noon.

Table 1

Ventilation schedule of the apartment (except in the periods of absence, one person was present).

Day of the Week, Date	Absence Time Start–End	Absence Duration [h]	Ventilation Time Start–End	Ventilation Duration [h]
Sunday, 03.10.2021			09:40–10:0012:40–13:2015:30–16:5017:55–18:5020:00–21:15	0.330.671.330.921.25
Monday, 04.10.2021	12:00–19:00	7.00	00:30–01:5008:00–08:3511:05–11:5519:30–20:25	1.330.580.830.92
Tuesday, 05.10.2021	8:30–17:00	8.50	00:58–01:3306:20–06:4518:40–19:3021:00–21:40	0.580.420.830.67
Wednesday, 06.10.2021	10:00–18:15	8.25	06:50–07:2009:35–09:5520:05–20:50	0.500.330.75
Thursday, 07.10.2021	8:50–12:3012:40–15:15	3.672.58	00:35–01:3507:15–07:4108:35–08:4618:55–19:5022:50–23:20	1.000.430.180.920.50
Friday, 08.10.2021	10:30–20:00	9.50	06:35–07:1009:51–10:00	0.580.15
Saturday, 09.10.2021	–	–	–	–
Sunday, 10.10.2021			01:03–02:4209:30–10:0011:50–12:3013:00–13:4016:32–17:0521:06–21:2623:28–23:55	1.650.500.670.670.550.330.45
Monday, 11.10.2021			06:40–06:55	0.25
			09:30–09:48	0.30
Tuesday, 12.10.2021	10:00–18:20	8.33	07:30–08:00	0.50
Wednesday, 13.10.2021	11:15–15:55	4.67	–	–
Thursday, 14.10.2021	17:30–19:25	1.92	07:15–07:35	0.33

Table 2

List of scenarios with the required and/or recommended design ventilation rate (DVR).

Scenario	Level of Obligation	Required, Recommended DVR			Reference
		Descriptive Criterion	Quantitative CriterionGeneral	Quantitative CriterionTest Apartment
1	Requirement	Minimal air changes per hour (ACH) in the absence of occupants to remove building emissions and prevent harm (can be considered in the 24 h cycle)	0.20 h−1	13.9 m3 h−1 (0.2 ACH)	[19]
2	Requirement	Minimal outdoor air intake	15.0 m3 h−1 person−1	15.0 m3 h−1 (0.2 ACH)	[19]
3=6	Requirement	Minimal ACH	0.50 h−1	34.6 m3 h−1 (0.5 ACH)	[19,22,23]
4	Requirement	Minimal volume of air per floor surface area (without consideration of other sources)	1.50 m3 h−1 m−2	40.0 m3 h−1 (0.6 ACH)	[19]
5A:Cat I-III	Recommendation	Ventilation rate per person and per m2 floor area	Cat I: 12.6 m3 h−1 person−1 + 0.9 m3 h−1 m−2Cat II: 9.0 m3 h−1 person−1 + 0.54 m3 h−1 m−2Cat III: 5.4 m3 h−1 person−1 + 0.36 m3 h−1 m−2	36.6 m3 h−1 (0.5 ACH)23.4 m3 h−1 (0.3 ACH)15.0 m3 h−1 (0.2 ACH)	[25]
5B:Cat I-III	Recommendation	Ventilation rate per person	Cat I: 36.0 m3 h−1 person−1Cat II: 25.2 m3 h−1 person−1Cat III: 14.4 m3 h−1 person−1	36.0 m3 h−1 (0.5 ACH)25.2 m3 h−1 (0.4 ACH)14.4 m3 h−1 (0.2 ACH)	[25]
5C:Cat I-IV	Recommendation	Ventilation rate per m2 floor area with infiltration	Cat I: 1.76 m3 h−1 m−2Cat II: 1.51 m3 h−1 m−2Cat III: 1.26 m3 h−1 m−2Cat IV: 0.83 m3 h−1 m−2	46.9 m3 h−1 (0.7 ACH)40.2 m3 h−1 (0.6 ACH)33.6 m3 h−1 (0.5 ACH)22.1 m3 h−1 (0.3 ACH)	[25]

Table 3

Requirements and recommendations for the concentrations in indoor air of: (a) radon (²²²Rn) [19,25,28,58,59]; and (b) carbon dioxide (CO₂) [17,19,60].

Obligatory Level	Required, Recommended Concentration	Reference
(a) 222Rn
Requirement: the permissible value of Rn in indoor air	400 Bq m−3	[19]
Requirement: the reference level of the average annual concentration of radon in closed living and working spaces	300 Bq m−3	[28]
Recommendation: WHO guideline value	100 Bq m−3	[25,59]
Recommendation: WELL Building Standard. The following conditions are met in projects with regularly occupied spaces at or below grade: radon less than 4 pCi/L in the lowest occupied level	4 pCi L−1 (148 Bq m−3)	[58]
(b) CO2
Requirement: the permissible value of CO2 in indoor air	1667 ppm	[19]
Recommendation: for the design and assessment of energy performance in buildings	Cat I: 350 ppm aCat II: 500 ppm aCat III: 800 ppm aCat IV: <800 ppm a	[17]
Recommendation:	Max: 2500 ppmRecommended: 1000 ppm	[60]

Note: ^a value above outdoor background concentration. Cat I: presents a high level of expectation and is recommended for spaces occupied by very sensitive and fragile persons with special requirements such as disabled, sick, very young children, and elderly persons; Cat II: normal level of expectation, should be used for new buildings and renovations; Cat III: acceptable, moderate level of expectation, may be used for existing buildings; Cat IV: values outside the criteria for the above categories. The last category should only be accepted for a limited part of the year [17].

Table 4

Deviations of simulated C_Rn-s [Bq m⁻³] and $C_{{\mathrm{CO}}_2-\mathrm{s}}$ [ppm] from the limit values [17,25,59,60] for 6 sets of scenarios.

Scenario	DVR	Duration of $C_{{\mathrm{CO}}_2-\mathrm{s}}$above 1000 ppm		Duration of $C_{{\mathrm{CO}}_2-\mathrm{s}}$ above 800 ppm		Duration of CRn-s above 100 Bq m−3
		[h]	[%]	[h]	[%]	[h]	[%]
1	13.9 m3 h−1 (0.2 ACH)	185	64	237	82	10	4
2	15.0 m3 h−1 (0.2 ACH)	176	61	267	93	6	2
3=6	34.6 m3 h−1 (0.5 ACH)	0	0	93	32	0	0
4	40.0 m3 h−1 (0.6 ACH)	0	0	60	21	0	0
5A_Cat I-III	36.6 m3 h−1 (0.5 ACH)	0	0	83	29	0	0
	23.4 m3 h−1 (0.3 ACH)	87	30	169	59	2	1
	15.0 m3 h−1 (0.2 ACH)	176	61	218	76	6	2
5B_Cat I-III	36.0 m3 h−1 (0.5 ACH)	0	0	81	28	0	0
	25.2 m3 h−1 (0.4 ACH)	61	21	159	55	1	0.4
	14.4 m3 h−1 (0.2 ACH)	188	65	226	79	8	3
5C_Cat I-IV	46.9 m3 h−1 (0.7 ACH)	0	0	0	0	0	0
	40.2 m3 h−1 (0.6 ACH)	0	0	61	21	0	0
	33.6 m3 h−1 (0.5 ACH)	0	0	133	46	0.5	0.2
	22.1 m3 h−1 (0.3 ACH)	117	41	163	57	2	0.7

References

1. Roof, K.; Oleru, N. Public health: Seattle and King County’s push for the built environment. J. Environ. Health; 2008; 71, pp. 24-27. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18724501]

2. Dovjak, M.; Kukec, A. Creating Healthy and Sustainable Buildings: An Assessment of Health Risk Factors; Springer Open: Cham, Switzerland, 2019; pp. 1-42.

3. WHO (World Health Organisation). An Estimated 12.6 Million Deaths Each Year Are Attributable to Unhealthy Environments. 2016; Available online: https://www.who.int/news/item/15-03-2016-an-estimated-12-6-million-deaths-each-year-are-attributable-to-unhealthy-environments (accessed on 10 October 2021).

4. EPA (Environmental Protection Agency). Indoor Air Quality (IAQ). 2021; Available online: https://www.epa.gov/indoor-air-quality-iaq (accessed on 10 October 2021).

5. EC (European Commission). Health & Consumer Protection Directorate-General, Opinion on Risk Assessment on Indoor Air Quality. 2007; Available online: http://ec.europa.eu/health/ph_risk/committees/04_scher/docs/scher_o_055.pdf (accessed on 10 October 2021).

6. Bartuska, T.J. Introduction: Definition, design, and development of the built environment. The Built Environment: A Collaborative Inquiry into Design and Planning; McClure, W.R.; Bartuska, T.J. Wiley: New Jersey, NJ, USA, 2007; pp. 3-14.

7. Krainer, A. Toward Smart Buildings; European Commission: London, UK, 1993; 84p.

8. Awbi, H.B. Ventilation of Buildings; 2nd ed. Taylor & Francis: New York, NY, USA, 2003; 522p.

9. WHO (World Health Organisation). Radon and Health. 2021; Available online: https://www.who.int/news-room/fact-sheets/detail/radon-and-health (accessed on 10 October 2021).

10. Dovjak, M.; Slobodnik, J.; Krainer, A. Det as a collateral damage of present day extensive renovations. Stroj. Vestn.; 2019; 65, pp. 31-40.

11. Dovjak, M.; Slobodnik, J.; Krainer, A. Consequences of energy renovation on indoor air quality in kindergartens. Build. Simul.; 2020; 13, pp. 691-708. [DOI: https://dx.doi.org/10.1007/s12273-020-0613-6]

12. EPA. United States Environmental Protection Agency Radon. 2020; Available online: https://www.epa.gov/radon (accessed on 10 October 2021).

13. UNSCEAR (United Nations Scientific Committee on the Effects of Atomic Radiation). Sources and Effects of Ionizing Radiation. UNSCEAR 2006 Report to the General Assembly, with Scientific Annexes; United Nations Publications: New York, NY, USA, 2008.

14. Kikaj, D.; Vaupotič, J.; Chambers, S.D. Identifying persistent temperature inversion events in a subalpine basin using radon-222. Atmos. Meas. Tech.; 2019; 12, pp. 4455-4477. [DOI: https://dx.doi.org/10.5194/amt-12-4455-2019]

15. Kikaj, D.; Chambers, S.D.; Kobal, M.; Crawford, J.; Vaupotič, J. Characterizing atmospheric controls on winter urban pollution in a topographic basin setting using Radon-222. Atmos. Res.; 2020; 237, 1104838. [DOI: https://dx.doi.org/10.1016/j.atmosres.2019.104838]

16. ANSI/ASHRAE (American Society of Heating, Refrigerating and Air-Conditioning Engineers Inc.). ANSI/ASHRAE Standard 62.1. Ventilation for Acceptable Indoor Air Quality; ASHARE: Atlanta, GA, USA, 2019.

17. ISO (International Organisation for Standardization). EN 15251. Indoor Environmental Input Parameters for Design and Assessment of Energy Performance of Buildings—Addressing Indoor Air Quality, Thermal Environment, Lighting and Acoustics; ISO: Geneva, Switzerland, 2007.

18. Küçükhüseyin, Ö. CO₂ monitoring and indoor air quality. REHVA J.; 2021; 1, pp. 45-59. Available online: https://www.rehva.eu/rehva-journal/chapter/co2-monitoring-and-indoor-air-quality (accessed on 10 October 2021).

19. OJ RS (Official Journal of the Republic of Slovenia). Rules on the Ventilation and Air-Conditioning of Building. OJ RS, No 42/02 with Amending; Office of the Government of the RS: Ljubljana, Slovenia, 2002.

20. OJ RS (Official Journal of the Republic of Slovenia). Rules on Efficient Use of Energy in Buildings with a Technical Guideline. OJ RS, No 52/10 with Amending; Office of the Government of the RS: Ljubljana, Slovenia, 2010.

21. MOP RS (Ministry of the Environment and Spatial Planning, Republic of Slovenia). TSG-1-004. Technical Guideline, Efficient Use of Energy; Ministry of the Environment and Spatial Planning, Republic of Slovenia: Ljubljana, Slovenia, 2010.

22. MOP RS (Ministry of the Environment and Spatial Planning, Republic of Slovenia). Proposal on Rules on Efficient Use of Energy in Buildings with a Technical Guideline; Ministry of the Environment and Spatial Planning, Republic of Slovenia: Ljubljana, Slovenia, 2021.

23. MOP RS (Ministry of the Environment and Spatial Planning, Republic of Slovenia). Proposal of Technical Guideline on the Field of Building Construction TSG-1-004: 2021 Efficient Use of Energy in Buildings; Ministry of the Environment and Spatial Planning, Republic of Slovenia: Ljubljana, Slovenia, 2021.

24. ISO (International Organisation for Standardization). CR 1752. Ventilation for Buildings—Design Criteria for the Indoor Environment; ISO: Geneva, Switzerland, 1999.

25. ISO (International Organisation for Standardization). EN 16798-1. Energy Performance of Buildings. Ventilation for Buildings. Indoor Environmental Input Parameters for Design and Assessment of Energy Performance of Buildings Addressing Indoor Air Quality, Thermal Environment, Lighting and Acoustics. Module M1-6; ISO: Geneva, Switzerland, 2019.

26. Dovjak, M.; Virant, B.; Krainer, A.; Šijanec-Zavrl, M.; Vaupotič, J. Determination of optimal ventilation rates in educational environment in terms of radon dosimetry. Int. J. Hyg. Environ. Health; 2021; 234, 113742. [DOI: https://dx.doi.org/10.1016/j.ijheh.2021.113742]

27. Persily, A.K. Quit Blaming ASHRAE Standard 62.1 for 1000 ppm CO₂. Proceedings of the 16th International Conference on Indoor Air Quality and Climate (INDOOR AIR 2020); Online, 1 November 2020; International Society of Indoor Air Quality and Climate (ISIAQ): Herndon, VA, USA, 2020.

28. OJ RS (Official Journal of the Republic of Slovenia). Decree on the National Radon Programme. OJ RS, No 18/18 with Amending; Office of the Government of the RS: Ljubljana, Slovenia, 2018.

29. OJ EU (Official Journal of the European Union). Directive 2010/31/EU of the European Parliament and of the Council of 19 May 2010 on the Energy Performance of Buildings (Recast); The European Parliament and the Council of the European Union: Brussels, Belgium, 2010; Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2010:153:0013:0035:EN:PDF (accessed on 10 October 2021).

30. OJ EU (Official Journal of the European Union). Directive (EU) 2018/844 of the European Parliament and of the Council of 30 May 2018 Amending Directive 2010/31/EU on the Energy Performance of Buildings and Directive 2012/27/EU on Energy Efficiency; The European Parliament and the Council of the European Union: Brussels, Belgium, 2018; Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32018L0844&from=EN (accessed on 10 October 2021).

31. Jiránek, M.; Kačmaříková, V. Dealing with the increased radon concentration in thermally retrofitted buildings. Radiat. Prot. Dosim.; 2014; 160, pp. 43-47. [DOI: https://dx.doi.org/10.1093/rpd/ncu104]

32. Cucoş, A.; Dicu, T.; Cosma, C. Indoor radon exposure in energy-efficient houses from Romania. Rom. J. Phys.; 2015; 60, pp. 1574-1580.

33. Hou, Y.; Liu, J.; Li, J. Investigation of indoor air quality in primary school classrooms. Procedia Eng.; 2015; 121, pp. 830-837. [DOI: https://dx.doi.org/10.1016/j.proeng.2015.09.037]

34. Földváry, V.; Bekö, G.; Langer, S.; Arrhenius, K.; Petráš, D. Effect of energy renovation on indoor air quality in multifamily residential buildings in Slovakia. Build. Environ.; 2017; 122, pp. 363-372. [DOI: https://dx.doi.org/10.1016/j.buildenv.2017.06.009]

35. Vasilyev, A.; Yarmoshenko, I.; Zhukovsky, M. Radon safety in terms of energy efficiency classification of buildings. IOP Conf. Ser. Earth Environ. Sci.; 2017; 72, 012020. [DOI: https://dx.doi.org/10.1088/1755-1315/72/1/012020]

36. Pampuri, L.; Caputo, P.; Valsangiacomo, C. Effects of buildings’ refurbishment on indoor air quality. Results of a wide survey on radon concentrations before and after energy retrofit interventions. Sustain. Cities Soc.; 2018; 42, pp. 100-106. [DOI: https://dx.doi.org/10.1016/j.scs.2018.07.007]

37. Yang, S.; Pernot, J.G.; Jörin, C.H.; Niculita-Hirzel, H.; Perret, V.; Licina, D. Radon investigation in 650 energy efficient dwellings in western Switzerland: Impact of energy renovation and building characteristics. Atmosphere; 2019; 10, 777. [DOI: https://dx.doi.org/10.3390/atmos10120777]

38. Bogdanovica, S.; Zemitis, J.; Bogdanovics, R. The effect of CO₂ concentration on children’s well-being during the process of learning. Energies; 2020; 13, 6099. [DOI: https://dx.doi.org/10.3390/en13226099]

39. Jung, C.; Awad, J. The improvement of indoor air quality in residential buildings in Dubai, UAE. Buildings; 2021; 11, 250. [DOI: https://dx.doi.org/10.3390/buildings11060250]

40. Kotol, M.; Rode, C.; Clausen, G.; Nielsen, T.R. Indoor environment in bedrooms in 79 Greenlandic households. Build. Environ.; 2014; 81, pp. 29-36. [DOI: https://dx.doi.org/10.1016/j.buildenv.2014.05.016]

41. Gładyszewska-Fiedoruk, K. Pomiary stężenia dwutlenku węgla w sypialniach domku jednorodzinnego. Ciepłownictwo Ogrzew. Went.; 2008; 6, pp. 32-34.

42. Bekö, G.; Lund, T.; Nors, F.; Toftum, J.; Clausen, G. Ventilation rates in the bedrooms of 500 Danish children. Build. Environ.; 2010; 45, pp. 2289-2295. [DOI: https://dx.doi.org/10.1016/j.buildenv.2010.04.014]

43. Sekhar, S.C.; Goh, S.E. Thermal comfort and IAQ characteristics of naturally/mechanically ventilated and air-conditioned bedrooms in a hot and humid climate. Build. Environ.; 2011; 46, pp. 1905-1916. [DOI: https://dx.doi.org/10.1016/j.buildenv.2011.03.012]

44. Szczepanik-Ścisło, N.; Flaga-Maryańczyk, A. Measurements and simulation of CO₂ concentration in a bedroom of a passive house. Tech. Trans.; 2018; 9, pp. 163-180.

45. García-Tobar, J. A comparative study of indoor radon levels between two similar dwellings using CONTAM Software. Environments; 2018; 5, 59. [DOI: https://dx.doi.org/10.3390/environments5050059]

46. García-Tobar, J. Weather-dependent modelling of the indoor radon concentration in two dwellings using CONTAM. Indoor Built Environ.; 2019; 28, pp. 1341-1349. [DOI: https://dx.doi.org/10.1177/1420326X19841119]

47. García-Tobar, J. Study of radon propagation in a dwelling using the CFD modelling capabilities of CONTAM. Phys. J.; 2020; 5, pp. 72-79.

48. NIST (National Institute of Standards and Technology). CONTAM User Guide and Program Documentation Version 3.4; Dols, W.S.; Polidoro, B.J. NIST: Gaithersburg, MD, USA, 2021.

49. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World map of the Köppen-Geiger climate classification updated. Meteorol. Z.; 2019; 15, pp. 259-263. [DOI: https://dx.doi.org/10.1127/0941-2948/2006/0130]

50. EPA (Environmntal Protection Agency). A Standardized EPA Protocol for Characterizing Indoor Air Quality in Large Office Buildings; EPA, Environmntal Protection Agency: Washington, DC, USA, 2003; Available online: https://www.epa.gov/sites/default/files/2016-04/documents/standardized_iaq_base_protocol_for_characterizing_iaq_in_large_buildings-pdf.pdf (accessed on 20 November 2021).

51. IJS (Jožef Stefan Institute). Manual, Kontinuirne Meritve Koncentracije Aktivnosti Radona-222 in Njegovih Potomcev z Aktivnimi Merilniki (Continuous Measurements of Radon-222 and its Daughter Activity Concentrations by Active Monitors); ELME-DN-23, Version April 2019 IJS: Ljubljana, Slovenia, 2019; (In Slovene)

52. ISO (International Organisation for Standardization). ISO/IEC 17025. General Requirements for the Competence of Testing and Calibration Laboratories; ISO: Geneva, Switzerland, 2017.

53. SARAD GmbH. Manual, Radon Scout Professional, Professional Radon Monitor/Personal Dosimeter, Version March 2019; SARAD GmbH: Dresden, Germany, 2019.

54. Saphymo GmbH. AlphaGUARD, Portable Radon Monitor, Use Manual, 08/2012; Saphymo GmbH: Frankfurt, Germany, 2012.

55. ARSO MOP (Slovenian Environmental Agency, Ministry of the Environment and Spatial Planning). Meteorological Data, Ljubljana, Slovenia; Slovenian Environmental Agency, Ministry of the Environment and Spatial Planning: Ljubljana, Slovenia, 2021; Available online: http://meteo.arso.gov.si (accessed on 10 October 2021).

56. Persily, A.; de Jonge, L. Carbon dioxide generation rates for building occupants. Indoor Air; 2017; 27, pp. 868-879. [DOI: https://dx.doi.org/10.1111/ina.12383]

57. REHVA (Federation of European Heating, Ventilation and Air Conditioning Associations). REHVA COVID-19 Guidance Document How to Operate HVAC and Other Building Service Systems to Prevent the Spread of the Coronavirus (SARS-CoV-2) Disease (COVID-19) in Workplaces; REHVA: Brussels, Belgium, 2021; Available online: https://www.rehva.eu/fileadmin/user_upload/REHVA_COVID-19_guidance_document_V4_23112020_V2.pdf (accessed on 10 October 2021).

58. Delos Living LLC. WELL Building Standard; Delos Living LLC: New York, NY, USA, 2021; pp. 22-60.

59. WHO (World Health Organization). WHO Handbook on Indoor Radon: A Public Health Perspective; Zeeb, H.; Shannoun, F. World Health Organization: Geneva, Switzerland, 2009; Available online: https://www.who.int/ionizing_radiation/env/9789241547673/en (accessed on 10 October 2021).

60. ANSI/ASHRAE (American Society of Heating, Refrigerating and Air-Conditioning Engineers Inc.). ANSI/ASHRAE Standard 62.1. Ventilation for Acceptable Indoor Air Quality; ASHARE: Atlanta, GA, USA, 2010.

61. Kraus, M.; Kubeková, D. Airtightness of energy efficient buildings. GSTF J. Eng. Technol.; 2013; 2, pp. 74-80. [DOI: https://dx.doi.org/10.5176/2251-3701_2.2.74]

62. Howieson, S.G.; Sharpe, T.R.; Farren, P. Building tight—Ventilating right? How are new air tightness standards affecting indoor air quality in dwellings?. Build. Serv. Eng. Res. Technol.; 2013; 35, pp. 475-487. [DOI: https://dx.doi.org/10.1177/0143624413510307]

63. Ye, W.; Zhang, X.; Gao, J.; Cao, G.; Zhou, X.; Su, X. Indoor air pollutants, ventilation rate determinants and potential control strategies in Chinese dwellings: A literature review. Sci. Total Environ.; 2017; 586, pp. 696-729. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2017.02.047]

64. Zhukovsky, M.V.; Vasilyev, A.V. Mechanisms and sources of radon entry in buildings constructed with modern technologies. Radiat. Protect. Dosim.; 2014; 160, pp. 48-52. [DOI: https://dx.doi.org/10.1093/rpd/ncu111]

65. Yarmoshenko, I.; Malinovsky, G.; Vasilyev, A.; Onishchenko, A. Model of radon entry and accumulation in multi-flat energy-efficient buildings. J. Environ. Chem. Eng.; 2021; 9, 105444. [DOI: https://dx.doi.org/10.1016/j.jece.2021.105444]

66. SURS (Statistical Office of the Republic of Slovenia). Dwellings, Slovenia, 1 January 2018; Statistical Office of the Republic of Slovenia: Ljubljana, Slovenia, 2021; Available online: https://www.stat.si/StatWeb/en/news/Index/8160 (accessed on 10 October 2021).

67. EC (European Commission)Eurostat. Is Your Home too Crowded?; European Commission: Brussels, Belgium, 2020; Available online: https://ec.europa.eu/eurostat/web/products-eurostat-news/-/ddn-20200422-1 (accessed on 10 October 2021).

68. MOP RS (Ministry of the Environment and Spatial Planning, Republic of Slovenia). Resolution on the National Housing Programme 2015–2025; Ministry of the Environment and Spatial Planning, Republic of Slovenia: Ljubljana, Slovenia, 2016; Available online: https://www.gov.si/assets/ministrstva/MOP/Publikacije/d42acebd4d/resolucija_nsp_2015_2025.pdf (accessed on 10 October 2021).

69. Batog, P.; Badura, M. Dynamic of changes in carbon dioxide concentration in bedrooms. Procedia Eng.; 2013; 57, pp. 175-182. [DOI: https://dx.doi.org/10.1016/j.proeng.2013.04.025]

70. Fujiyoshi, S.; Tanaka, D.; Maruyama, F. Transmission of airborne bacteria across built environments and its measurement standards: A Review. Front. Microbiol.; 2017; 8, 2336. [DOI: https://dx.doi.org/10.3389/fmicb.2017.02336] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29238327]