RESEARCH PAPER
Andrey Tsapalov1 and Konstantin Kovler1*
1National Building Research Institute, Faculty of Civil and Environmental Engineering, Technion - Israel Institute of Technology, Haifa, Israel
A simple experimental scheme is proposed to study the concentration behavior and mechanism of soil radon transport from beneath a building foundation. In addition to the results of the annual monitoring of soil radon concentration and environmental factors, the calculated results of the annual continuous monitoring of the soil radon entry rate into experimental room with high levels of indoor radon are presented. This room is located in the basement of one of the buildings in Haifa, Israel. The correlation between radon behavior and environmental factors is discussed. It is found that the soil radon concentration beneath the building is maximum in the winter season due to very heavy rains. There is an absence of soil radon entry by only convection mechanism into the experimental room for the summer and autumn seasons; however, the indoor radon concentration remains high (about 500 Bq m-3) due to diffusion of soil radon. A regular soil radon entry into the experimental room by convection is observed in the winter and spring seasons with a maximum, approximately at the beginning of the spring season, and therefore, the indoor radon concentration is twice as high as in the summer and autumn seasons. The obtained results indicate that not only convective process but also diffusion of soil radon plays a significant role in the formation of indoor radon.
Keywords: soil radon; indoor radon; concentration; radon entry rate; diffusion; convection; monitoring
Citation: Journal of the European Radon Association 2022, 3: 7630 http://dx.doi.org/10.35815/radon.v3.7630
Copyright: © 2022 Andrey Tsapalov and Konstantin Kovler. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 International License (http://creativecommons.org/licenses/by-nc/4.0/), permitting all non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Published: 04 March 2022
Competing interests and funding: The authors declare that there is no conflict of interest. This research project has received funding from the European Union’s Horizon 2020 research and innovation programme H2020-MSCA-IF-2017 under Grant agreement No. 792789 (RadonACCURACY).
*Konstantin Kovler, National Building Research Institute, Faculty of Civil and Environmental Engineering, Technion - Israel Institute of Technology, Haifa 3200003, Israel. Email: cvrkost@technion.ac.il
Radon is a dangerous carcinogen (1) constituting 3–14% of lung cancer, as stated by the World Health Organization (WHO) (2). Indoor radon is the main source of human exposure to radiation because about 80% of our life-time is spent inside buildings. According to the recommendations of WHO (3) and Basic Safety Standards (BSS) Directive (4), the annual average indoor radon (AAIR) concentration should be limited to 300 Bq m-3. The national reference levels (RLs) of indoor radon vary in different countries due to differences in regional levels and usually range from 100 to 300 Bq m-3.
Radon is formed as a result of the natural decay of 226Ra, a part of the ‘uranium-radium’ family of radionuclides that are unevenly distributed in the crust of the earth (5). Possessing neither color nor odor, radon is a noble gas, and therefore, can migrate intensively in porous media (as soil) due to molecular diffusion, as well as gas convection. The main source of radon in buildings is soil (5). The flow of soil radon into the atmosphere or into buildings can vary significantly in space and time. This depends on geological, climatic, and weather conditions (6), as well as on the construction of the building foundation and connections of engineering utilities. In contrast to the atmosphere, radon can be accumulated inside building; however, its long-term accumulation is normally prevented by natural or forced ventilation.
The behavior of indoor radon is the result of the two competing and very unstable anthropogenic phenomena: air change and radon entry. Both air change rate (ACR) and radon entry rate (RER) should be considered when designing the buildings, especially in radon priority areas. ACR is an important hygienic parameter that is regulated by building codes and must be ensured at the minimum level through appropriate design decisions (7), considering the particular climate. If the Reference Level (RL) for indoor radon and ACR is known, then using the well-known indoor radon balance equation (5), one can estimate the maximum allowable RER level at which the AAIR concentration will not be exceeded. However, the question arises on how to ensure the calculated (control) RER level is not exceeded using data on the geology (related to soil radon), structure, and materials of the building foundation. For example, in the study by Gulabyants et al. (8), a completely adequate physical algorithm and a solution to this problem are given, which, however, assumes that diffusion is the main process for transport of soil radon into building. At the same time, many scientific reports (9, 10) argue that the decisive transport of soil radon into building is by convection due to the so-called ‘stack effect’ (11), although neither single experimental scheme nor experimental data on RER to prove it have been presented so far. Indeed, there are many publications that provide compelling indirect and model evidence that soil gas convection is responsible for high levels of indoor radon but are we completely sure that there are no buildings in which diffusive transport prevails over the soil radon convection?
In addition, the lack of direct control over ACR due to its complexity (12, 13) does not allow making accurate conclusions about the amplitude (specific values) and temporal variations of RER based only on the results of measurements of indoor radon concentration. In this regard, it is important to develop experimental schemes that allow direct measurements of RER, or at least estimate the magnitudes of parameters, which are closely related to RER.
In contrast to diffusion, the extremely irregular action of convection significantly complicates the behavior of RER, and consequently a creation of an adequate physical model for accurate estimation and reliable control of RER (14, 15). Therefore, when creating mathematical models for the radon transport into building, it would be very important to determine the criteria that allow identifying the mechanism of this transport. Thus, the main goal of this work was to study the patterns of soil radon behavior beneath the building foundation and to define the mechanism of soil radon entry into building using the example of a building with a high concentration of radon in the basement.
Achieving the above goal, first, requires long-term and continuous data collection (at least for 1 year) in relation to the following three main parameters: 1) soil radon concentration beneath the building foundation, 2) soil RER into the building, and 3) radon concentration indoors. Obviously, the organization of monitoring of 1) and 2) requires drilling the building foundation. In addition, the levels of 1) and 3) should be significantly higher than the background levels to ensure satisfactory accuracy of the raw data to be analyzed. Therefore, one of the basement laboratory rooms in the building of the National Building Research Institute (NBRI) located at the campus of the Technion – Israel Institute of Technology (Haifa) was chosen as a research object (see Table 1), where the natural indoor radon concentration, on average, exceeds 500 Bq m-3 (while radon concentrations in the neighbor room are by– two to three times lower). This experimental room (ER1) has a middle size of 6.0 × 5.2 × 3.1 h m (V = 96.7 m3 and S = 31.2 m2) and is rarely visited (most of the time it is closed). Experimental room ER1 is constantly air-conditioned to maintain a stable temperature and humidity without forced outside air inflow. The floors and walls of the room are monolithic structures made of concrete with a thickness of about 20 cm. The concrete foundation slab used as a floor and concrete walls in contact with soil have no entrances of any underground facilities and any defects (which can be visually observed – such as cracks, voids, and openings). The properties of the underlying soil, which obviously serves as a radon source, were not studied because this task was out of the scope of the current project and also because of the difficulty in extracting a sufficient amount of the soil from beneath the foundation slab and from different depths for the analysis. At the same time, it was clear that the main composition of the underlying soil represents mainly carbonate rock containing a relatively minor amount of the quartz sand in its upper part used as a construction fill.
Code | ER1 | ||
Room | AAIR level, Bq∙m-3 | 722 | |
Average temperature (range), °C | 22.5 (20–32) | ||
Average humidity (range), % | 47 (19–65) | ||
Floor | Basement | ||
Type of room (number) | Lab (108) | ||
Occupancy | Rarely | ||
Local system of | Heating | Yes | |
air conditioning | |||
Ventilation | Natural | ||
Room area* | Middle | ||
Building | Location | Technion, Haifa | |
Name | NBRI | ||
Function | Office | ||
Year built | 1970s | ||
Number of stories | 3 | ||
Material of walls | Concrete | ||
Material of floors | Concrete | ||
Central system of | Heating | No | |
air conditioning | No | ||
Ventilation | Natural | ||
Building area** | Large | ||
Asphalt pavement *** | Yes | ||
Foundation type | Slab | ||
Basement | Yes | ||
Geology | Relief | Uneven | |
Soil type | Rocky | ||
*Room area: small (area < 20 m2), middle (20 m2 < area < 60 m2), and large (area > 60 m2). **Building area: small (area < 200 m2), middle (200 m2 < area < 600 m2), and large (area > 600 m2). ***The asphalt or other poorly permeable layers surrounding the building. |
The general measurement scheme in ER1, including the location of the soil radon monitoring points, is shown in Fig. 1.
Fig. 1. Scheme of the soil radon measurement in the ER1 located in the basement.
The space behind the wall is only partially filled with soil, the level of which is approximately 1 m higher than the floor in the basement. Access to this space is not provided; however, it is known that its area in the horizontal section is larger than the floor area of ER1. In the most distant (from the wall of ER1) part of this space, there is a small opening that communicates with the outside air.
Figure 1 shows the location of two sampling points in the floor (FR and FL) and one in the wall (W). The sampling point is a through hole drilled in the structure from which a tube (about 1-m long and with 4 mm inner diameter) comes into the basement room (ER1). The other end of the tube connects to the small plastic box where the radon monitor (Radon Eye Plus2 or Sun Nuclear 1028 XP) is installed. The same tube (outlet) comes out from the opposite side of the box, but is much shorter – only 5 cm long. This scheme was used for continuous measurements of the RER at each of the three measurement points. The soil radon concentration measurements were carried out periodically in the manual sampling mode described in more detail in the next paragraph, through the long tube, which was disconnected from the box for not longer than 5 min. Periodically, for the measurement of pressure differences between the soil gas and indoor air, before manual sampling of the soil radon, the differential pressure meter Testo 512 (with the measurement range of 0–2 hPa and the resolution of 1 Pa) was connected to the long tube.
The indoor radon concentration measurements were carried out in the continuous monitoring mode with a recording period of 1 h using a Radon Eye Plus2 located at a height of about 1.5 m.
To assess the influence of various factors on the behavior of soil radon, the data of synchronous annual monitoring of the following environmental factors were also collected: 1) air temperature, 2) humidity, 3) atmospheric pressure, 4) wind speed, 5) precipitation, 6) insolation (related to cloudiness), 7) tidal forces (Earth’s rotation speed or length of day), and 8) seismic activity. The collection of such data does not require much effort. For example, the nearest professional weather station in Tel-Aviv, located about 100 km from the Haifa, served as a source of meteorological data (http://www.rp5.ru). The data on the tidal forces were taken from The International Earth Rotation and Reference Systems Service (http://www.iers.org). The data on seismic activity in the Haifa area with a radius of about 150–200 km were taken from The Geological Survey of Israel (http://www.earthquake.co.il). In addition, several of the most sensitive earthquakes (with a magnitude of more than 4), which occurred at a greater distance from Haifa, were considered.
To avoid the influence of the measurement method on the behavior of soil radon, it was decided to sample the minimum volumes of soil gas only twice a week in manual mode, because the soil radon concentration under the building is more stable than indoor radon. The soil gas was sampled from the FL (450 mL), FR (450 mL) and W (1,050 mL) points periodically, using a simple original sampler including a syringe and the sorption column (SC-13) containing activated charcoal (13 mL) (Fig. 2). The SC-13 is capable of adsorbing 100% of the radon activity of the soil gas sample of a small volume (<1 L), even at its maximum moisture content (16).
Fig. 2. Equipment for manual sampling of soil radon on the column with activated charcoal.
The concentration of soil radon CS (Bq m-3) is as follows:
where VS is the sample volume (m3) and A is radon activity in charcoal at the time of sampling (Bq).
Radon activity in charcoal was measured using a simple but sensitive detector BDB-13, the description and characteristics of which were discussed by Tsapalov (17). This detector was calibrated using a gamma spectrometer and a fine granular non-emanating sample with a known (certified) specific 226Ra activity.
The radon concentration in the boxes to determine the soil RERs into the boxes were continuously measured at the same points FR, FL, and W (Fig. 1). The continuous radon monitor was located inside the box. Two tubes stem from the box on the opposite sides, one of which was short with an open end, and the other tube (long) was connected to the hole drilled in the concrete floor (or wall). If a convective flow of soil radon occurs through a long tube (while the long tube excludes a diffusion flow into the box), the radon activity concentration in the box increases. However, the increase in radon activity in the box is compensated for by the runoff of radon, which is caused by its natural decay, leakage, and convection of soil gas through the short tube. Thus, the balance of radon activity concentrations in the box C(t) (Bq m-3) is approximately described by the following equation:
where E(t) is the soil RER from the hole in the floor into the box only due to convection (Bq/h), Δt is the registration period of radon concentration in the box equal to 1 h, V is the box free volume, which is equal to 0.0035 m3 with the installed Sun Nuclear 1028 XP continuous radon monitor, or to 0.0030 m3 for the box containing Radon Eye Plus2 device, λR is the radon decay constant (0.00755 1/h), λL is the radon leakage (1/h), and λG(t) is the soil gas change rate in the box (1/h).
In the case of estimating the rate of radon flux in the opposite direction (from the air of the room to the soil), equation (2) is not entirely correct. Nevertheless, extremely low and negative values of E(t) will qualitatively indicate the absence of soil radon entry into room by convection or the existence of a reverse flux. In addition, a quantitative assessment of the reverse radon flux is not of interest in the framework of this study.
The parameter λG(t) is related to the volumetric rate W(t) (m3/h) of the convective flux of soil radon, which enters the box with the know volume V through the long tube:
At the same time, W(t) is related to E(t) and to the current soil radon concentration:
Thus, the equilibrium equation (2) can be rewritten as follows:
The value of the radon leakage from the box is determined using the following simple experiment. The high activity of radon is injected inside the box, and the long tube is closed. Then, a decrease in the radon concentration in the box is recorded until it drops down to 2,000 Bq m-3 in the continuous measurement mode. The combined value of radon decay constant and leakage from the box (λR+λL) = 0.034 1/h is the result of the calculation according to formula (2), where E(t) = 0 and λG(t) = 0. Then the time series of the measured values of radon concentration in the box can be transformed into the time series of the values of E(t) using the following expression derived from Equation 5:
The results of the radon activity concentration measured by different continuous monitors and using the charcoal method in the corresponding ranges were compared. Testing has shown a good convergence of measurement results over a concentration range of 100–10,000 Bq m-3. However, in the higher range Sun Nuclear 1028 XP underestimated the radon concentration relatively to the charcoal method (see above), which is the reference method. Figure 3 shows the experimental points of simultaneous measurements of radon concentration by the Sun Nuclear 1028 XP and the charcoal method, and also provides the correction function to reconstruct the true radon concentration from the measurement results by Sun Nuclear 1028 XP. The results of the radon concentration measurements by Sun Nuclear 1028 XP multiplied by the correction function value represent the ‘true radon concentrations’.
Fig. 3. Correction function for measurement results by Sun Nuclear 1028 XP.
The continuous measurements of indoor radon and radon concentrations in the boxes FR and W were conducted since 12 July, 2019, and the periodical measurements of soil radon concentration in the points FL, FR, and W since 24 October 2019. A bit later, at 01 November 2019 the parallel continuous measurements of radon concentrations in the box FL started. Unfortunately, due to the lockdown of the lab activities caused by the Coronavirus pandemic, the annual monitoring of soil radon concentration is not complete; nevertheless, this practically did not worsen the quality of the experimental observations and the corresponding conclusions.
The results of the annual continuous monitoring of radon concentrations and various environmental factors are shown in Fig. 4.
Fig. 4. Results of continuous monitoring of soil and indoor radon concentration in the ER1 and environmental factors.
As can be seen from Fig. 4, hot weather characterizes the summer season with an average temperature of 25–30°C, without any precipitation, and with an average wind speed of 3–5 m/s. For the winter season, the average temperature was about 15–20°C, with strong gusts of wind being more frequently observed at a slightly higher average speed, in parallel with heavy rainfalls in the form of tropical showers. For this reason, in the winter season, the moisture content of the soil beneath the buildings increased; so on 26 December 2019, it was already impossible to sample a soil gas at the FR point. After the next very heavy rain on 04 January 2019, rainwater penetrated into some of the basements, including ER1. The soil under the foundation of the building became saturated with moisture; therefore, between 05 and 12 January it was impossible to conduct measurements and, accordingly, to determine the concentration of soil radon in the both FR and FL points by manual sampling. The next missing soil radon measurements from 15 March 2015 to 26 April 2020 was due to the strict lockdown of the lab activities caused by the spread of Coronavirus in the country.
The reconstructed results of radon concentration monitoring in the FR and FL boxes using the equation from Fig. 3, including the values of soil RERs into the boxes from the FR and FL points calculated by Equation (6), are shown in Fig. 5. Because the soil radon concentration measurements started after the commencement of the radon concentrations in the box FR, the initial value CS = 50 kBq m-3 was assumed for the calculation of the RER from the FR point in the beginning of the study. The same value of CS was used for the calculation of RER from the FR and FL points during the lockdown period. The modelling shows that the variations of CS in the range 30–70 kBq m-3 do not significantly affect the changes in the calculated values of the RER corresponding to the two periods mentioned above.
Fig. 5. Reconstructed results of radon concentration monitoring in the FR and FL boxes, and calculated soil RER from the FR and FL points.
The average soil radon concentrations in the FR and FL points for the entire measurement period from 24 October, 2019 till 23 July, 2020 (excluding the lockdown period) differed insignificantly and were equal to (55.5 ± 13.6) and (59.1 ± 16.2) kBq m-3, respectively.
The time series of the calculated values of RER from the W point is not given in Fig. 5, because their average value for the period 01 November 2019 to 01 July 2020 was only 0.29 Bq/h. This value is much lower than the RER observed from the FR (4.62 Bq/h) and FL (5.60 Bq/h) points for the same period.
Periodical measurements of the pressure differences between soil gas and indoor air during the year demonstrate that the maximum difference does not exceed 1 Pa in the FR and FL points and 2 Pa in the W point. The maximum values were observed in winter, while in summer the pressure drop fluctuated around zero. The pressure difference of 1–2 Pa can be hardly determined accurately at the device error of only 1 Pa. However, we assume that there is no need for measuring so accurately such a low-pressure difference.
The obtained continuous time series of radon and environmental factors made it possible to carry out a correlation analysis, the results of which are shown in Table 2.
For more convenient comparison and discussion of the obtained data, the results of the annual monitoring and their trends, as well as the calculated and average values of the most important parameters, are shown in Fig. 6.
Fig. 6. The results of the annual monitoring of indoor and soil radon, main environmental factors, and the calculation of time series of RER.
Analysis of the data in Figs. 5 and 6, as well as the correlation coefficients in Table 2, leads to the following observations:
The stable convective transport of soil radon is, however, observed for about 6 months a year during the winter and spring seasons, which corresponds to the period of maximum indoor radon concentration (about 1,000 Bq m-3). During the remaining 6 months, there is no convective transfer of soil radon, with the exception of single, but relatively intense radon entry by convection, the duration of which is not longer than 4 days. However, these sporadic events do not significantly increase indoor radon concentration. Moreover, during the summer and autumn seasons, when there is no convective transport of soil radon, a relatively high level of indoor radon is still observed with an average value of about 500 Bq m-3. Thus, it can be concluded that a sufficiently high level of indoor radon during the summer and autumn seasons is created only due to the diffusion of soil radon into the experimental room.
It should be clarified that until now the term (parameter) ‘soil radon entry rate’ was referred to the rate of radon release only due to convection from specific holes (with a diameter of 4 mm) artificially created by drilling in the floor and wall of the experimental room. Therefore, sometimes a more precise wording was used in the text as ‘soil radon entry rate from the FR point’. However, for the use in the indoor radon balance equation (5), ‘soil radon entry rate’ only refers to the soil radon, which penetrated through the foundation slab (not only due to convection but also due to diffusion) and released from the corresponding area into the room. Since one of the considered parameters does not take into account the diffusion mechanism, the relationship of these parameters cannot have a universal mathematical expression, even when the resistance of the floor slab to the transport of soil radon into the room is zero (i.e. an open soil surface serves as a room floor). Nevertheless, the data presented above quite adequately reflect the nature of temporal variations, as well as the potential of the soil radon entry due to the soil gas convection only. Obviously, the effect of this potential will decrease with increasing resistance of the building structure to the flow of soil radon into the room.
Thus, the proposed experimental scheme and mathematical apparatus make it possible to qualitatively assess the nature of temporal variations in soil RER, as well as the role of soil gas convection in the mechanism of soil radon transport inside the building.
In view of the above, especially mentioned in item 5), the statement that convection is the main mechanism responsible for the transport of soil radon into buildings with an elevated radon concentrations cannot be considered correct. This opinion us shared by the radon standards of the ISO 11665 series (19) mentioning that both these radon transport mechanisms can take a place, while the convective transport is not distinguished as the main one. In addition, it should be borne in mind that unlike convection, the diffusion of soil radon into the building is a regular (constant in time) process, which is due to the always existing concentration gradient between soil and indoor radon (19), the values of which usually differ by more than two orders of magnitude. Therefore, the characterization of the mechanisms of transport of soil radon into the building will look more correct if we say that the main mechanism of transport of radon into the building is diffusion, which can sometimes be enhanced by convection of soil gas.