1. bookVolume 68 (2023): Issue 2 (June 2023)
Journal Details
License
Format
Journal
eISSN
1508-5791
First Published
25 Mar 2014
Publication timeframe
4 times per year
Languages
English
Open Access

Radon exposure in kindergartens in one Bulgarian district

Published Online: 25 Apr 2023
Volume & Issue: Volume 68 (2023) - Issue 2 (June 2023)
Page range: 51 - 56
Received: 15 Nov 2022
Accepted: 16 Jan 2023
Journal Details
License
Format
Journal
eISSN
1508-5791
First Published
25 Mar 2014
Publication timeframe
4 times per year
Languages
English
Introduction

The major part of the radiation dose that humans receive from natural radioactive sources is due to inhalation of radon and its decay products. It has been found that prolonged exposure to high levels of radon (222Rn) increased an individual’s risk of developing lung cancer [1, 2]. Radon is a naturally occurring gas and comes from the decay chain of uranium (238U), which is found in soils and rocks. In the outdoor atmosphere, it presents lower concentrations, but in buildings, the concentration can grow to dangerous levels [3]. The main source of radon in buildings is the geology of the site [2]. On the other hand, the indoor radon concentration (CRn) could vary with the internal surfaces and volumes of the building type and the air change rate [4] as well as with the building factors such as construction, the presence of a ventilation system in the building, etc.

When kindergartens are located in public buildings that carry a relatively high radon exposure risk, this poses a severely problematic situation, because children are particularly sensitive to unhealthy indoor environmental pollutants [5]. Moreover, they spend a large part of their time during the day in kindergarten. Therefore, attention needs to be paid to the monitoring of indoor CRn in such public buildings.

This paper presents a study of the CRn variation in kindergartens in the Bulgarian district of Montana and analyzes the impact of building characteristics, such as the presence of ventilation, implemented energy-saving measures, and sewerage systems, on indoor radon variations. Assessments of the children’s and worker’s doses from exposure to radon and its decay products by applying the measured radon activity concentrations were performed.

Research methodology
Study area

Montana district is located in Northwestern Bulgaria. The district covers an area of 3635.5 km2. Administratively, the district is divided into 11 municipalities: Berkovitsa, Boychinovtsi, Brusartsi, Varshets, Vulchedrum, Georgi Damyanovo, Lom, Medkovets, Montana, Chiprovtsi, and Yakimovo (Fig. 1). The 50 state kindergartens located on the territory of the district have been surveyed. All rooms where children and staff spend their time, located on the ground floor, were measured.

Fig. 1.

Map of Montana district.

Radon measurements

The study was carried out using passive CR-39 track detectors, type RSKS, provided by Radosys Ltd., Budapest, Hungary. The sampling period was for four months from December 2019/January 2020 to April/May 2020, and the exact date of placement and collection of the detector was recorded in the protocol for each kindergarten. The detectors were placed about 1 m from the floor and approximately 0.5 m from the ceiling, away from external windows, walls, doors, and heat appliances, thereby ensuring that they would be exposed to the air but remain out of the reach of children. According to the usual procedure adopted for such requirements, passive detectors were placed on a shelf in the room for the 4-month period of exposure. The total number of detectors provided for the study is 656. The percentage of detector loss is 8% and 602 results for rooms in the kindergartens and nurseries were analyzed. The detectors were calibrated in an accredited laboratory and measurements were made according to the ISO/DIS 11665-4 standard [6, 7].

A detailed questionnaire for each surveyed building was made available to the kindergarten and nursery staff, which they were required to fill in; the same needed the following details to be specified: the exact location and characteristics of the building (presence of a mechanical ventilation system and implemented energy-saving measures; types of sewerage system, heating, type of windows, etc.), and habits of the occupants.

Dose estimation from radon and its decay products

The individual dose was calculated conservatively from the average CRn in kindergartens’ premises using standard and recognized methods employing the ICRP dose conversion factor (DCF) and the estimated exposure time (in hours). When the measured concentration of Rn-222 gas is used in the calculation of the dose, the DCF is 6.7 × 10−6 mSv·Bq−1·h−1·m3, as recommended by the International Commission on Radiological Protection (ICRP), assuming the equilibrium factor F = 0.4 [8].

The effective dose from radon was evaluated for the period of detector exposure (four months – December 2019 to April 2020). Conservatively, the assumption was made that staff work 8 h/day and 5 days/week, giving rise to a total of 720 h over four months; and that children stay in kindergartens for approximately 6 h/day, giving rise to a total of 420 h over four months.

Statistical analysis

IBM SPSS Statistics, v.23 was used for performing the statistical analysis. To assess the impact on radon variation, the measurements were grouped, according to the geographical location of municipalities, into nine groups; and, according to the presence of ventilation, heating, sewage systems, and thermal external wall insulation, into two groups (“yes” and “no”).

Results and discussion
Overall results and effective dose assessment

The statistical parameters of indoor CRn in 602 premises in the measured kindergartens and nurseries are presented in Table 1.

Descriptive statistics of indoor radon concentration

ParameterCRn
Number of rooms  602
Arithmetic mean (Bq·m−3)  125
Standard deviation (Bq·m−3)  135
Median (Bq·m−3)    84
Minimum value (Bq·m−3)    10
Maximum value (Bq·m−3)1439
Coefficient of variation (%)  108
Geometric mean (Bq·m−3)    88
Geometric standard deviation            2.23

The average value of CRn in the premises of the studied kindergartens and nurseries on the territory of the Montana district is 125 Bq·m−3, and the geometric mean (GM) value is 88 Bq·m−3. A review of 63 national and regional indoor radon surveys in kindergartens and schools in Europe, Asia, Africa, and North America has found the average radon arithmetic mean (AM) for all these surveys to be 59 Bq·m−3 [9]. This value is lower than the assessed AM of the CRn for kindergartens and nurseries in the Montana district. The number of premises in which the CRn exceeds the national reference level [10] of the average annual CRn in the air of public buildings and workplaces of 300 Bq·m−3 was 30 or approximately 5% of all measurements made. This percentage is lower than that corresponding to kindergarten premises of the three Visegrad countries, where the rooms with high CRn amounted to 7.5% of the total number measured [11]. The CRn in 12 rooms is >500 Bq·m−3. In these buildings, measures for reducing the concentrations should be taken. The maximum value of CRn was 1439 Bq·m−3, measured in the children’s room, and it is approximately five times higher than the reference level.

The effective dose for children and workers was conservatively estimated for four months (the sampling period from December to April), and the results are presented in Fig. 2.

Fig. 2.

The effective dose for children and workers is estimated for four months by kindergarten-wise presentation.

For four months during the winter period, the staff in five kindergartens are likely to receive effective doses due to radon >0.6 mSv, which is the half of annual average from radon inhalation of 1.26 mSv·a−1, as assessed by UNSCEAR – 2008 [12]. Children in only one kindergarten are likely to receive an effective dose >0.6 mSv.

Summary of results by municipalities

The gathered data and descriptive statistics were analyzed using a classification based on municipalities in the district, and the results are presented in Table 2. The highest AM of the CRn in the premises of kindergartens was assessed in the municipality of Brusartsi (AM = 276 Bq·m−3), and where the highest value of 1439 Bq·m−3 was found. The lowest AM of the CRn was measured in the municipality of Georgi Damyanovo (AM = 70 Bq·m−3). The measured maximum CRn values above the national reference level [10], presented in Table 2, were observed in eight of the municipalities (Berkovitsa, Boychinovtsi, Brusartsi, Vulchedrum, Lom, Medkovets, Montana, and Yakimovo).

Descriptive statistics of indoor radon concentration by municipalities

MunicipalitiesNumber of premisesAM (Bq/m3)SD (Bq/m3)CV (%)Mediana (Bq/m3)Min. (Bq/m3)Max. (Bq/m3)Test Shapiro–Wilk (p)
Berkovitsa123145118.9  8210731  6430.008
Boychinovtsi  27200177.9  8911031  6690.135
Brusartsi  13276426.2154  701314390.152
Vulchedrum122  96  96.4100  7012  8240.165
Varshets  18124  58.2  4712136  2310.120
Georgi Damyanovo  20  70  17.9  25  6445  1130.051
Lom  89  73  62.9  86  5618  4290.011
Medkovets  18246240.4  9814654  8530.162
Montana158130124.6  96  9610  9720.053
Chiprovtsi    6105  69.9  66  7255  2340.153
Yakimovo    8130134.6104  6444  3800.012

A parametric Shapiro–Wilk test was applied to test the hypothesis for normal distribution of the results of CRn by municipalities. The normal distribution of the indoor radon values in the municipalities of Boychinovtsi, Brusartsi, Vulchedrum, Varshets, Georgi Damyanovo, Medkovets, Montana, and Chiprovtsi has been confirmed. The results of the other municipalities in the Montana district do not follow the normal distribution. The Kolmogorov–Smirnov test was applied to test the log-normal distribution of the data, which did not follow the normal distribution. The test confirmed the log normality of distribution (KS, p >0.05).

The AM in the two municipalities of Brusartsi and Medkovets is >200 Bq·m−3, and for Boychinovtsi municipality the average CRn is 200 Bq·m−3. These municipalities are located mainly in mountainous areas. A non-parametric Kruskal–Wallis test (KW, p <0.0001) was applied to compare the data by municipality. A statistically significant difference was found between the surveyed groups, which shows that place position influences the CRn in kindergarten buildings. This circumstance is well known for residential buildings and substantiates the proposition that radon priority areas could be arranged based on municipality [13, 14]. In this way, preventive and corrective measures for reducing indoor CRn could be managed more optimally.

Variation of CRn and building characteristic

The variations of the CRn depend on the various characteristics of the building. The influence of building characteristics on CRn varies from one region to the next [13]. The descriptive statistic by the studied characteristics in the buildings according to the divided groups (presence or not in the building) is presented in Table 3. The normal distribution was assessed using the Shapiro–Wilk test (p >0.05). Not all groups followed a normal distribution.

Descriptive statistics of indoor radon concentration by the common ventilation system, central heating system, thermal external walls insulation, and central sewerage system

Statistic parameterCommon ventilation systemCentral heating systemThermal external walls insulationCentral sewerage system
YesNoYesNoYesNoYesNo
Number  87    515    445  157    282  320  474  128
Median  60      91      87    69    101    68    76  100.5
Arithmetic mean  72.1    134.3    135.9    95.5    152.8  101.1  120.4  143.4
Standard deviation  48.0    143.0    147.3    86.9    161.1  101.8  124.2  169.6
Coefficient of variation    0.67        1.06        1.08      0.91        1.05      1.01      1.03      1.18
p-Value of Shapiro–Wilk  <0.001    <0.001    <0.001  <0.001    <0.001  <0.001  <0.001    <0.001
Minimum  26    10    10  18    10  12  10    13
Maximum3171439143951914398249721439

Figure 3 shows the distribution of average CRn by groups for the presence or absence of ventilation, heating and sewage systems, and thermal external wall insulation. The non-parametric Mann–Whitney test was used to assess the difference between the groups (U, p <0.05), because the value by groups does not follow the normal distribution. A significant difference was found between all the groups (U, p <0.001).

Fig. 3.

The plots of CRn according to the: (a) presence of central heating system, (b) presence of common ventilation system, (c) presence of thermal external walls insulation, and (d) presence of central sewerage system.

The CRn (AM = 72 Bq·m−3) where the common ventilation system is built is approximately twice lower compared to buildings that do not have a ventilation system (AM = 134 Bq·m−3). The results confirm the hypothesis that the presence of a ventilation system in the building reduces the concentration of radon, while the presence of a central heating system increases the CRn in the building. The difference between the temperature in the building and the ground increases, when the whole building is heated. This temperature difference increases the underground pressure, which, in turn, increases the movement of radon gas from the soil toward the building.

Results of CRn from rooms located in buildings with a central sewerage system (AM = 120 Bq·m−3) are lower than those corresponding to usage of a septic tank (AM = 143 Bq·m−3). Radon is likely collected in the septic tank and is transported to the building via a pipeline. The buildings that have built ventilation and sewerage systems have lower CRn, while in a building that has thermal external walls insulation and a common heating system, the CRn is higher.

Conclusion

This survey represents a study on indoor radon in kindergartens and nurseries in Montana district, Bulgaria. The average value of CRn in the premises is 125 Bq·m−3, and the geometric mean value is 88 Bq·m−3. The number of premises in which the CRn exceeds the national reference level of the average annual CRn of 300 Bq·m−3 is 30 or approximately 5% of all measured cases. To reduce the high CRn, appropriate corrective measures should be implemented. Conservatively evaluated doses of workers in five kindergartens for four months are above those assessed by the UNSCEAR – 2008 annual average from a radon inhalation of 1.26 mSv·a−1 [12].

The difference in the measurements by municipalities shows that the geographical location affects the indoor CRn in buildings.

The study found that the central heating system and the thermal insulation of the external walls increased the CRn, while the general ventilation system and the central sewage system reduced the CRn in the studied kindergarten buildings. These findings could be used in planning activities for radon prevention and corrective measures for public buildings in the Montana region.

Fig. 1.

Map of Montana district.
Map of Montana district.

Fig. 2.

The effective dose for children and workers is estimated for four months by kindergarten-wise presentation.
The effective dose for children and workers is estimated for four months by kindergarten-wise presentation.

Fig. 3.

The plots of CRn according to the: (a) presence of central heating system, (b) presence of common ventilation system, (c) presence of thermal external walls insulation, and (d) presence of central sewerage system.
The plots of CRn according to the: (a) presence of central heating system, (b) presence of common ventilation system, (c) presence of thermal external walls insulation, and (d) presence of central sewerage system.

Descriptive statistics of indoor radon concentration by the common ventilation system, central heating system, thermal external walls insulation, and central sewerage system

Statistic parameter Common ventilation system Central heating system Thermal external walls insulation Central sewerage system
Yes No Yes No Yes No Yes No
Number   87     515     445   157     282   320   474   128
Median   60       91       87     69     101     68     76   100.5
Arithmetic mean   72.1     134.3     135.9     95.5     152.8   101.1   120.4   143.4
Standard deviation   48.0     143.0     147.3     86.9     161.1   101.8   124.2   169.6
Coefficient of variation     0.67         1.06         1.08       0.91         1.05       1.01       1.03       1.18
p-Value of Shapiro–Wilk   <0.001     <0.001     <0.001   <0.001     <0.001   <0.001   <0.001     <0.001
Minimum   26     10     10   18     10   12   10     13
Maximum 317 1439 1439 519 1439 824 972 1439

Descriptive statistics of indoor radon concentration

Parameter CRn
Number of rooms   602
Arithmetic mean (Bq·m−3)   125
Standard deviation (Bq·m−3)   135
Median (Bq·m−3)     84
Minimum value (Bq·m−3)     10
Maximum value (Bq·m−3) 1439
Coefficient of variation (%)   108
Geometric mean (Bq·m−3)     88
Geometric standard deviation             2.23

Descriptive statistics of indoor radon concentration by municipalities

Municipalities Number of premises AM (Bq/m3) SD (Bq/m3) CV (%) Mediana (Bq/m3) Min. (Bq/m3) Max. (Bq/m3) Test Shapiro–Wilk (p)
Berkovitsa 123 145 118.9   82 107 31   643 0.008
Boychinovtsi   27 200 177.9   89 110 31   669 0.135
Brusartsi   13 276 426.2 154   70 13 1439 0.152
Vulchedrum 122   96   96.4 100   70 12   824 0.165
Varshets   18 124   58.2   47 121 36   231 0.120
Georgi Damyanovo   20   70   17.9   25   64 45   113 0.051
Lom   89   73   62.9   86   56 18   429 0.011
Medkovets   18 246 240.4   98 146 54   853 0.162
Montana 158 130 124.6   96   96 10   972 0.053
Chiprovtsi     6 105   69.9   66   72 55   234 0.153
Yakimovo     8 130 134.6 104   64 44   380 0.012

United Nations Scientific Committee on the Effects of Atomic Radiation. (2009). Effects of ionizing radiation: United Nations Scientific Committee on the Effects of Atomic Radiation 2006 Report to the General Assembly, with scientific annexes. Vol. II. Annex E: Sources-to-effects assessment for radon in homes and workplace. New York: UN. https://www.unscear.org/docs/publications/2006/UNSCEAR_2006_Annex-E-CORR.pdf. United Nations Scientific Committee on the Effects of Atomic Radiation . ( 2009 ). Effects of ionizing radiation: United Nations Scientific Committee on the Effects of Atomic Radiation 2006 Report to the General Assembly, with scientific annexes. Vol. II. Annex E: Sources-to-effects assessment for radon in homes and workplace . New York : UN . https://www.unscear.org/docs/publications/2006/UNSCEAR_2006_Annex-E-CORR.pdf . Search in Google Scholar

Zeeb, H., Shannoun, F., & World Health Organization. (2009). Handbook on indoor radon: A public health perspective. Geneva, Switzerland: WHO. https://experts.umn.edu/en/publications/who-hand-book-on-indoor-radon-a-public-health-perspective. Zeeb H. Shannoun F. World Health Organization . ( 2009 ). Handbook on indoor radon: A public health perspective . Geneva, Switzerland : WHO . https://experts.umn.edu/en/publications/who-hand-book-on-indoor-radon-a-public-health-perspective . Search in Google Scholar

Pervin, S., Yeasmin, S., Khandaker, M. U., & Begum, A. (2022). Radon concentrations in indoor and outdoor environments of atomic energy center Dhaka, Bangladesh and concomitant health hazards. Front. Nucl. Eng., 1, 901818. https://doi.org/10.3389/fnuen.2022.901818. Pervin S. Yeasmin S. Khandaker M. U. Begum A. ( 2022 ). Radon concentrations in indoor and outdoor environments of atomic energy center Dhaka, Bangladesh and concomitant health hazards . Front. Nucl. Eng. , 1 , 901818 . https://doi.org/10.3389/fnuen.2022.901818 . Search in Google Scholar

Laughlin, J. M. (2012). Radon: Past, present, and future. Rom. J. Phys., 58, S5–S13. https://rjp.nipne.ro/2013_58_Suppl/0005_0013.pdf. Laughlin J. M. ( 2012 ). Radon: Past, present, and future . Rom. J. Phys. , 58 , S5 S13 . https://rjp.nipne.ro/2013_58_Suppl/0005_0013.pdf . Search in Google Scholar

International Commission on Radiological Protection. (1992). Protection against radon-222 at home and work. Ann. ICRP, 23, 1–38. https://doi.org/10.1080/09553009414551371. International Commission on Radiological Protection . ( 1992 ). Protection against radon-222 at home and work . Ann. ICRP , 23 , 1 38 . https://doi.org/10.1080/09553009414551371 . Search in Google Scholar

International Organization for Standardization. (2012). Measurement of radioactivity in the environment – Air: radon-222 – Part 4: Integrated measurement method for determining average activity concentration using passive sampling and delayed analysis. ISO 11665-4:2020(en). https://www.iso.org/obp/ui/#iso:std:iso:11665:-4:ed-2:v1:en. International Organization for Standardization . ( 2012 ). Measurement of radioactivity in the environment – Air: radon-222 – Part 4: Integrated measurement method for determining average activity concentration using passive sampling and delayed analysis . ISO 11665-4:2020(en) . https://www.iso.org/obp/ui/#iso:std:iso:11665:-4:ed-2:v1:en . Search in Google Scholar

International Atomic Energy Agency. (2019). Design and conduct of indoor radon surveys. Vienna: IAEA. (Safety Reports Series no. 98). https://www-pub.iaea.org/MTCD/Publications/PDF/PUB1848_web.pdf. International Atomic Energy Agency . ( 2019 ). Design and conduct of indoor radon surveys . Vienna : IAEA . ( Safety Reports Series no. 98 ). https://www-pub.iaea.org/MTCD/Publications/PDF/PUB1848_web.pdf . Search in Google Scholar

International Commission on Radiological Protection. (2017). Occupational intakes of radionuclides: Part 3. ICRP Publication 137. Ann. ICRP, 46(3/4). https://journals.sagepub.com/doi/pdf/10.1177/ANIB_46_3-4. International Commission on Radiological Protection . ( 2017 ). Occupational intakes of radionuclides: Part 3. ICRP Publication 137 . Ann. ICRP , 46 ( 3/4 ). https://journals.sagepub.com/doi/pdf/10.1177/ANIB_46_3-4 . Search in Google Scholar

Zhukovsky, M., Vasilyev, A., Onishchenko, A., & Yarmoshenko, I. (2018). Review of indoor radon concentrations in schools and kindergartens. Radiat. Prot. Dosim., 181, 6–10. DOI: 10.1093/rpd/ncy092. Zhukovsky M. Vasilyev A. Onishchenko A. Yarmoshenko I. ( 2018 ). Review of indoor radon concentrations in schools and kindergartens . Radiat. Prot. Dosim. , 181 , 6 10 . DOI: 10.1093/rpd/ncy092 . Open DOISearch in Google Scholar

Republic of Bulgaria Nuclear Regulatory Agency. (2020). Regulation on radiation protection, adopted by CM Decree no. 20/14.02.2018, amended by CM no. 455/22.12.2020. https://bnra.bg/media/2022/02/regulation-on-radiation-protection.pdf. Republic of Bulgaria Nuclear Regulatory Agency . ( 2020 ). Regulation on radiation protection, adopted by CM Decree no. 20/14.02.2018, amended by CM no. 455/22.12.2020 . https://bnra.bg/media/2022/02/regulation-on-radiation-protection.pdf . Search in Google Scholar

Mullerova, M., Kozak, K., Kovacs, T., Smetanova, I., Csordas, A., Grzadziel, D., Holy, K., Mazur, J., Moravcsik, A., Neznal, Martin, & Neznal, Matej. (2016). Indoor radon survey in Visegrad countries. Appl. Radiat. Isot., 110, 124–128. DOI: 10.1016/j.apradiso.2016.01.010. Mullerova M. Kozak K. Kovacs T. Smetanova I. Csordas A. Grzadziel D. Holy K. Mazur J. Moravcsik A. Neznal Martin Neznal Matej . ( 2016 ). Indoor radon survey in Visegrad countries . Appl. Radiat. Isot. , 110 , 124 128 . DOI: 10.1016/j.apradiso.2016.01.010 . Open DOISearch in Google Scholar

United Nations Scientific Committee on the Effects of Atomic Radiation. (2010). Sources and effects of ionizing radiation. United Nations Scientific Committee on the Effects of Atomic Radiation UNSCEAR 2008 Report to the General Assembly with scientific annexes. Vol. I. New York: UN. https://www.unscear.org/unscear/uploads/documents/unscear-reports/UNSCEAR_2008_Report_Vol.I-CORR.pdf. United Nations Scientific Committee on the Effects of Atomic Radiation . ( 2010 ). Sources and effects of ionizing radiation. United Nations Scientific Committee on the Effects of Atomic Radiation UNSCEAR 2008 Report to the General Assembly with scientific annexes. Vol. I . New York : UN . https://www.unscear.org/unscear/uploads/documents/unscear-reports/UNSCEAR_2008_Report_Vol.I-CORR.pdf . Search in Google Scholar

Ivanova, K., Stojanovska, Z., Kunovska, B., Chobanova, N., Badulin, V., & Benderev, A. (2019). Analysis of the spatial variation of indoor radon concentrations (national survey in Bulgaria). Environ. Sci. Pollut. Res. Int., 26(7), 6971–6979. DOI: 10.1007/s11356-019-04163-9. Ivanova K. Stojanovska Z. Kunovska B. Chobanova N. Badulin V. Benderev A. ( 2019 ). Analysis of the spatial variation of indoor radon concentrations (national survey in Bulgaria) . Environ. Sci. Pollut. Res. Int. , 26 ( 7 ), 6971 6979 . DOI: 10.1007/s11356-019-04163-9 . Open DOISearch in Google Scholar

Ivanova, K., Stojanovska, Z., Cenova, M., & Kunovska, B. (2017). Building-specific factors affecting indoor radon concentration variations in different regions in Bulgaria. Air Qual. Atmos. Health, 10(9), 1151–1161. DOI: 10.1007/s11869-017-0501-0. Ivanova K. Stojanovska Z. Cenova M. Kunovska B. ( 2017 ). Building-specific factors affecting indoor radon concentration variations in different regions in Bulgaria . Air Qual. Atmos. Health , 10 ( 9 ), 1151 1161 . DOI: 10.1007/s11869-017-0501-0 . Open DOISearch in Google Scholar

Recommended articles from Trend MD