Next Article in Journal
Edge of the Present: A Virtual Reality Tool to Cultivate Future Thinking, Positive Mood and Wellbeing
Next Article in Special Issue
Study of Well Waters from High-Level Natural Radiation Areas in Northern Vietnam
Previous Article in Journal
Study to Investigate the Knowledge of Rare Diseases among Dentists, Orthodontists, Periodontists, Oral Surgeons and Craniomaxillofacial Surgeons
Previous Article in Special Issue
Characteristics of Thoron (220Rn) and Its Progeny in the Indoor Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radiological Assessment of Indoor Radon and Thoron Concentrations and Indoor Radon Map of Dwellings in Mashhad, Iran

1
Institute of Radiochemistry and Radioecology, University of Pannonia, 8200 Veszprém, Hungary
2
Materials and Nuclear Fuel Research School, Nuclear Science and Technology Research Institute, Tehran 11365-8486, Iran
3
Silesian Centre for Environmental Radioactivity, Central Mining Institute, plac Gwarków 1, 40-166 Katowice, Poland
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(1), 141; https://doi.org/10.3390/ijerph18010141
Submission received: 8 December 2020 / Revised: 22 December 2020 / Accepted: 23 December 2020 / Published: 28 December 2020

Abstract

:
A comprehensive study was carried out to measure indoor radon/thoron concentrations in 78 dwellings and soil-gas radon in the city of Mashhad, Iran during two seasons, using two common radon monitoring devices (NRPB and RADUET). In the winter, indoor radon concentrations measured between 75 ± 11 to 376 ± 24 Bq·m−3 (mean: 150 ± 19 Bq m−3), whereas indoor thoron concentrations ranged from below the Lower Limit of Detection (LLD) to 166 ± 10 Bq·m−3 (mean: 66 ± 8 Bq m−3), while radon and thoron concentrations in summer fell between 50 ± 11 and 305 ± 24 Bq·m−3 (mean 115 ± 18 Bq m−3) and from below the LLD to 122 ± 10 Bq m−3 (mean 48 ± 6 Bq·m−3), respectively. The annual average effective dose was estimated to be 3.7 ± 0.5 mSv yr−1. The soil-gas radon concentrations fell within the range from 1.07 ± 0.28 to 8.02 ± 0.65 kBq·m−3 (mean 3.07 ± 1.09 kBq·m−3). Finally, indoor radon maps were generated by ArcGIS software over a grid of 1 × 1 km2 using three different interpolation techniques. In grid cells where no data was observed, the arithmetic mean was used to predict a mean indoor radon concentration. Accordingly, inverse distance weighting (IDW) was proven to be more suitable for predicting mean indoor radon concentrations due to the lower mean absolute error (MAE) and root mean square error (RMSE). Meanwhile, the radiation health risk due to the residential exposure to radon and indoor gamma radiation exposure was also assessed.

1. Introduction

In general, people are exposed to ionizing radiation from various natural and artificial sources. Radon (222Rn), thoron (220Rn), and the progeny of both can be regarded as the largest contributor to the annual effective dose for the public in the world (50% of the total public dose) [1,2]; however, public exposure to ionizing radiation can be higher due to a new dose conversion factor [3]. Exposure to radon and its decay products is the second most common cause of lung cancer after tobacco smoking [4]. The health risks related to radon exposure primarily arise in indoor environments, while outdoor radon levels are generally low. The most important source of indoor radon is soil gas infiltration, and the intensity of this source relies on the composition of the ground, i.e., granite, tile, clay, etc. Soil gas infiltration is produced in mineral grains by the radioactive decay of 226Ra, emanated into the void spaces between the grains, transported by diffusion and advection/convection, and eventually exhaled from the soil into boreholes where it is detected. Moreover, cracks in concrete floors and walls, drainage pipes, connecting parts of buildings, heating, ventilation, and air conditioning ducts are the possible routes through which radon can enter into indoor environments [5].
Measuring indoor radon and thoron concentrations and radon mapping was considered for years and several papers were published on the topic around the world [6,7,8,9,10,11,12,13,14,15,16,17], including in many Iranian cities [18,19,20,21,22,23,24,25,26,27] to increase public awareness of environmental radioactivity and to predict radon-prone areas, which would help authorities with regard to the development of an appropriate strategy to reduce public exposure to radon and thoron. This reduced exposure would increase the quality of life and improve public long-term health. Due to the lack of data concerning indoor radon and thoron concentrations in houses in the city of Mashhad, an attempt was made to measure indoor radon and thoron concentrations in 78 houses to calculate the annual committed effective dose caused by the inhalation of radon and thoron. The measurements were taken during two seasons, summer (July–September 2019) and winter (December 2019–February 2020). The annual average radon and thoron concentrations were estimated by averaging measured concentrations during these periods. Subsequently, a radon map was produced using ArcGIS software and three different interpolation techniques, within a grid with the dimensions of 1 km × 1km. Meanwhile, the soil gas radon concentrations (as the major source of indoor radon) in different districts of Mashhad were also measured using a passive method based on CR-39 detectors during the summer when soil moisture and precipitation is low. Moreover, external exposure rates for terrestrial gamma radiation in Iran from 36 to 130 nGy h−1 with an average of 71 nGy h−1 have been reported [1,2]. In 2015, Sohrabi et al. also measured the indoors and outdoors gamma dose rates for about 1000 houses in 36 cities in Iran and the national mean background outdoor gamma dose rates were reported as being 70.2 nGy h−1 [28]. Therefore, to assess the radiation health risk, the public indoor doses from radon gas and indoor gamma radiation were compared and assessed.

2. Materials and Methods

2.1. Study Area

Mashhad is the second largest metropolis in Iran and is the capital of Razavi Khorasan Province in northeastern Iran. It has an area of 351 km2 and its population is more than 3 million people according to the last census (Statistical Centre of Iran, 2016). It has witnessed rapid growth over the last two decades, mostly as a result of its economic, social, and religious attractions. The city is 985 m above sea level with the geographic coordinates of 36°17′45′′ N, 59°36′43′′ E. Geologically, the Kalaj mountains, which consist of granitic hills covered by silty deposits, are situated to the south of Mashhad, towards the northwest is Kale Ghaemabad that is comprised of more sandy soil, and in all other directions is a plateau with a mix of clay loam and soft sandy soils. Figure 1 shows the location of Mashhad in Iran.

2.2. Measuring Techniques

Indoor radon and thoron concentrations are measured in the living rooms of houses at ground level. Regarding the recruited of participants, the priority was given to the older houses by selecting 3 to 5 dwellings from each district randomly depending on the size of the residential area. The majority of the houses examined were built 15 to 45 years ago using bricks composed of sand and cement along with cemented floors. The indoor measurements were conducted over a period of 90 days in total during the summer (July–September 2019) and winter (December 2019–February 2020). To determine the indoor radon and thoron concentrations, RADUET, a commercially available passive integrated radon–thoron discriminative detector, was used. These detectors consist of two diffusion chambers with different ventilation rates, and each chamber contains a CR-39 chip with the dimensions of 10 × 10 mm2 (RADUET, Radosys Ltd., Budapest, Hungary) for detecting the alpha particles emitted from radon and thoron as well as their progenies [29,30]. All detectors were hung at a height of 1–2 m above the ground using hard wire and positioned at least 20 cm away from any of the wall surfaces in the living rooms of the houses.
In addition, a solid-state nuclear track detector (SSNTD), CR-39, was used to measure radon concentrations in soil gas. In this regard, a hole was dug in the soil of about 11 cm in diameter and 50–60 cm in depth. Then, a long PVC tube was fixed into the hole with the covered top end of the tube protruding from the ground by about 5 cm. At the bottom of each tube, a NRBP radon dosimeter [31] was placed for a period of 45 days between July and September 2019.
After exposure, all detectors were wrapped in protective aluminum foil and returned for processing at the Institute. In the laboratory, they washed with distilled water and dried and then chemically etched. The etching condition for CR-39 was as follows: Solution: 6.0 M NaOH; Temperature: 90 °C; Time: 3 h. The track densities were counted using an optical transmission microscope and image analysis software. The calibration factors were determined as a result of exposure tests using radon and thoron calibration chambers at the Institute of Radiochemistry and Radioecology of the University of Pannonia, Hungary, and is comprehensively described in [30,31,32].

2.3. Annual Effective Dose, Excess Lifetime Cancer Risk (ELCR) and Lung Cancer Cases (LCC) Associated with Radon/Thoron Exposure

The annual committed effective doses originating from the inhalation of indoor radon or thoron were calculated using the following equation provided by [1]:
E(Rn/Tn) = C(Rn/Tn) × F(Rn/Tn) × t × K(Rn/Tn)
where ERn/Tn denotes the annual committed effective dose from exposure to radon or thoron (mSv yr−1); CRn/Tn stands for the annual average radon or thoron concentrations in houses (Bq m−3); F(Rn/Tn) represents the indoor equilibrium factors for radon or thoron and their respective progenies. The following values were provided by United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) in 2000 [1]: F = 0.40 and F = 0.02 for radon and thoron, respectively; t is the number of hours spent inside annually (7000 h). Also, K(Rn/Tn) denotes the following dose conversion factors recommended by UNSCEAR in 2000 [1]: KRn = 9 nSv and KTn = 40 nSv per unit of integrated radon and thoron concentrations (Bq h m−3), respectively.
In this survey, the average of the radon concentrations in the summer and winter represents the annual radon concentration. Furthermore, between December and February is regarded as the winter season when people tend to close windows because of the cold weather. The total exposure time was 90 days over the two seasons assessed. The Excess Lifetime Cancer Risk (ELCR) per 100,000 people was calculated using the following equation [33]:
ELCR = E(Rn/Tn) × DL × RF
where DL represents the life expectancy, estimated to be 70 years; and RF stands for the risk of fatal cancer per Sievert of 5.5 × 10−2 Sv−1 as recommended by International Commission on Radiological Protection (ICRP) Publication 103. Finally, the Lung Cancer Cases per year per million people (LCC) was estimated by using the risk factor lung cancer induction 18 × 10−6 mSv−1 and calculated by the following equation [34]:
LCC = E(Rn/Tn) × 18 × 10−6.

2.4. Statistical Analysis

The statistical analysis was carried out using IBM SPSS Statistics 21 (Armonk, NY, USA). The Kolmogorov-Smirnov test was applied to test the null hypothesis for the homogeneous distribution of the datasets. The Kruskal–Wallis non-parametric test with the Dunn’s post-hoc analysis was also used to test whether the samples originated from the same distribution based on the comparison of medians.

2.5. Radon Mapping and Cross-Validation

Radon mapping has great economic and social consequences; moreover, a high-resolution, accurate, and statistically powerful radon map is necessary to increase public awareness of environmental radioactivity and influence government policy with the purpose of reducing radon exposure in the general population. Depending on the datasets applied, two types of maps can be used: 1. Indoor Radon Maps which are based on indoor radon measurements (as applied in this study); and 2. Geogenic Radon Maps which are based on geological information [35]. The major merit of indoor radon maps is that radon concentrations are directly measured at the exposure point.
In this study, an indoor radon map was generated by using ArcGIS software version 10.7 (GDi Esri Hungary Ltd., Budapest, Hungary) over a grid with the dimensions of 1 km × 1 km; moreover, three interpolation methods were tested: inverse distance weighting (IDW), ordinary kriging (OK), and Empirical Bayesian kriging (EBK). The arithmetic mean (AM) was used over grid cells with the dimensions of 1 km × 1 km to predict the mean indoor radon concentration on the ground floor of buildings in the grid cells where no data was available. It is important to keep in mind that radon maps are only a probabilistic tool to make policy decisions such as prioritization; they cannot be used to derive radon concentrations for an individual dwelling. By using some indexes, namely the mean absolute error (MAE), root mean square error (RMSE), root mean squared logarithmic error (RMSLE), percentage bias (PB), and coefficient of determination (R2), the accuracy of the different techniques was also examined, as presented in Equations (4)–(8):
MAE = 1 n i = 1 n | Z i X i |
RMSE = 1 n i = 1 n ( Z i X i ) 2
RMSLE = 1 n i = 1 n ( l o g ( Z i + 1 ) l o g ( X i + 1 ) ) 2
PB = 100 i = 1 n ( Z i X i ) i = 1 n X i
R 2 = 1 i = 1 n ( Z i X i ) 2 i = 1 n ( | X i X ¯ | ) 2
where Xi and Zi denote the measured and predicted values in the location, n stands for the number of points in the validation group, and X ¯ represents the mean of Xi. MAE and RMSE are often applied to assess the performance of models. The model fits properly if the aforementioned indicators approach zero when calculated. PB (%) is the mean of the tendency in larger/smaller predicted values than those observed [36]. R2 is the fit line explaining to which degree the model is going to fit to the dataset [37].

3. Results and Discussion

3.1. Activity Measurements

Within this study, the radon and thoron concentrations in 78 houses were surveyed over a total exposure time of 90 days (45 days in both the summer and the winter) by using a RADUET detector. The frequency distribution of the indoor radon and thoron for the 78 houses assessed in Mashhad over the two seasons are shown in Figure 2. In addition, a comparison among normal distribution and log-normal distribution of data is also illustrated in Figure 3. The indoor 222Rn and 220Rn concentrations in the winter ranged from 75 ± 11 to 376 ± 24 Bq·m−3 with a mean value of 150 ± 19 Bq·m−3 and from below the Lower Limit of Detection (LLD) to 166 ± 10 Bq·m−3 with a mean value of 66 ± 8 Bq m−3, respectively. In the case of the summer, the indoor 222Rn and 220Rn concentrations ranged from 50 ± 11 to 305 ± 24 Bq·m−3 with a mean value of 114 ± 18 Bq·m−3 and from below the LLD to 122 ± 10 Bq·m−3 with a mean value of 48 ± 6 Bq m−3, respectively. In addition, the annual average indoor 222Rn and 220Rn concentrations in the studied areas were 132 ± 19 Bq·m−3 and 58 ± 7 Bq m−3, respectively. The main source of indoor 222Rn originates from soil gas infiltration, building materials, and ventilation [4]. Meanwhile, during cold winters, residents use natural gas and close all vents, causing the radon accumulation found in houses.
World organizations such as ICRP, WHO, and U.S. Environmental Protection Agency (EPA) have recommended various guidelines for radon exposure [1,4,33,38,39]. The annual average indoor radon concentration is below the recommendation values (300 Bq m−3) provided by the ICRP in 2010. The results concerning the annual average radon concentration exceed the action level (100 Bq m−3) recommended by the WHO in 2009. When compared to the worldwide geometric mean (GM) of 37 Bq·m−3 (geometric standard deviation (GSD) = 2.2) reported by UNSCEAR in 2000, the indoor radon in the city of Mashhad is almost 4 times (139.68 Bq m−3) and approximately 3 times (105.8 Bq m−3) higher than the world average in the winter and summer, respectively. It was also found that during the winter and summer, the indoor radon concentrations in 31% and 20% of the dwellings were higher than the reference level of 148 Bq·m−3 recommended by the US EPA in 2003.
The graph in Figure 4 shows the correlation between indoor 222Rn and 220Rn concentrations for the dwellings examined in Mashhad. Regarding the relationship between radon and thoron concentrations, no clear and strong correlation between was observed and thoron concentrations could not be predicted from widely available information concerning radon. However, indoor radon and thoron concentration might directly depend on the activity of 226Ra and 232Th (228Th) in building materials, ground is the main entry path of radon at dwellings; therefore, it could say that the content of both 222Rn and 220Rn depends on the building materials and soil composition.
As previously discussed, the main source of indoor 222Rn originates from soil gas infiltration, 222Rn concentrations in the soil gas of different districts in Mashhad were measured by using a passive method based on CR-39 detectors in the summer when soil moisture and precipitation are low. In order to determine soil gas radon concentrations, only 36 NRPB dosimeters were retrieved from where they were set up, while the remaining 6 dosimeters were considered lost. Figure 5 shows a histogram of soil gas radon concentrations in Mashhad during the summer. The soil gas radon concentrations recorded in the studied area fell within the range of 1.07 ± 0.28 to 8.02 ± 0.65 kBq·m−3 with a mean value of 3.07 ± 1.09 kBq m−3. As is shown in Figure 4, the activity concentrations of 222Rn vary from location to location, possibly because of the physic geological properties of the types of soil studied, topographic differences, as well as geomorphology and meteorological conditions of the region. The average radon concentrations in both soil gas and indoor environments are approximately the minimum and maximum values in the same region, respectively. Moreover, the correlation between indoor radon and soil gas radon concentrations for the districts studied is shown in Figure 6. The correlation analysis yielded a positive correlation (R2 = 0.361) between average indoor radon and soil-gas radon concentrations.
The normality distribution of data was checked using the Kolmogorov–Smirnov test. Considering the normality assumption in the null hypothesis of the Kolmogorov–Smirnov test, the probability value (p-value) in all tests was less than 1%; therefore, the normality distribution of radon and thoron concentrations in any of the following subfactors was rejected. In this study, by applying the Kruskal–Wallis nonparametric test with Dunn’s post hoc analysis, the null hypothesis, due to the absence of a statistically significant difference in the average gas concentration, was rejected; therefore, the season and type of gas affect the gas concentration (p-value < 0.05). The difference in radon concentrations between well-ventilated and poorly ventilated dwellings was statistically significant (p < 0.05). It was assumed that houses with natural ventilation are poorly ventilated and houses mechanical ventilation systems are well-ventilated houses. The finding indicates that the radon concentration is lower for well-ventilated dwellings compared to poorly ventilated ones. The results of this study are consistent with others that have been conducted concerning this topic [40,41,42]. Because the level of indoor radon concentration depends on the degree of indoor ventilation, moreover, in well-ventilated dwellings, radon can easily escape and does not accumulate inside, meaning indoor radon concentrations are less high in well-ventilated dwellings compared to in poorly ventilated ones [40].
Table 1 shows descriptive statistics that resulted from the measurement of the indoor radon and thoron concentrations in the 78 houses studied in Mashhad during the two seasons considered. Furthermore, the results reveal a seasonal variation in indoor radon and thoron concentrations, which were higher in the winter than in the summer. This is because the doors and windows of dwellings remain closed most of the time in the winter compared to in the summer, hence ventilation is poorer in the winter. The ratio of winter to summer concerning indoor radon and thoron concentrations was also established for all 78 dwellings studied. This ratio of indoor radon concentrations ranged from 1.23 to 1.48 with an average value of 1.31. With regard to the indoor thoron concentration, the average of this ratio was similar, at 1.36. The reason of heterogeneous behavior of seasonal variations in 222Rn and 220Rn concentrations might be that the source of 220Rn is mainly limited to the concentration of 232Th (228Th) in building materials, while in case of radon, the ground’s concentration is additionally considered. Therefore, in summer due to a high air exchange rate, e.g., using a ventilator or opening windows, the concentration of both 220Rn and 222Rn goes down, while in winter since the air exchange rate is lower than summer, the concentration of 220Rn and 222Rn build up but as the source of indoor 222Rn is both ground and building material rather than the only source of 220Rn as building materials, the seasonal change of indoor thoron concentration is less than that of indoor radon.
A comparison of radon concentration in the soil gas under investigation with those reported in other countries is also given in Table 2. It can clearly be seen that the radon concentration in soil samples from the Sri Ganganagar district and the northern state of Rajasthan in India, the city of Najaf in Iraq, and Yemen are in close agreement with the present work. It can be concluded that the soil in Mashhad is suitable for construction without posing any health hazards.

3.2. Radiation Dose and Risk Assessment

In this study, the average of the radon concentrations in the summer and winter was assumed to be the annual average radon concentration. The corresponding annual effective dose from the inhalation of radon and thoron was calculated as 3.7 ± 0.5 mSv yr−1. The committed effective dose from indoor radon and thoron was found to vary from 2.11 to 9.73 mSv yr−1 with a mean value of 4.22 mSv yr−1 for the winter, and 1.51 to 7.92 mSv yr−1 with a mean value of 3.14 mSv yr−1 for the summer. It can be seen that the committed effective doses were higher in the winter compared to the summer. The excess lifetime cancer risk and lung cancer cases per year per million people were also calculated to be 14.1 and 65.4, respectively.
According to the WHO in 2009, the risk of lung cancer increases by 16% per a 100 Bq·m−3 increase in radon concentration (long-time average) [4]. Therefore, the dose-response relation is linear, i.e., the risk of lung cancer is proportional to radon exposure [4]. Nevertheless, based on a report produced by the Ministry of Health’s Center for Disease Management in Iran, cancer is the third most significant cause of death after road traffic accidents and cardiovascular mortality. In 2019, Roshandel et al. reported that the age-standardized rates (ASR) of lung cancer were 127 and 52.1 per 100,000 Iranian males and females, respectively [52]. In the case of Razavi Khorasan Province, the ASR was 121.2 and 54.0 per 100,000 Iranian males and females, respectively. Therefore, the annual average excess risk due to radon inhalation in Mashhad is 14/100,000, i.e., less than the age-standardized death rate from cancer in Mashhad. Hence, indoor radon exposure is responsible for approximately 12% of lung cancer deaths in this city, which is close to the estimates by WHO in 2009 of the worldwide proportion of lung cancer due to radon (3–14%).
Given the annual average excess risk values by comparing local radiological risks with national cancer incidence data, it can be concluded that the local risks are raised but are not necessarily representative of the city as a whole. This radiation risk assessment should be considered with caution as the radon measurements are not sufficiently representative of the investigated area; moreover, calculations made using ICRP data only provide a broad overview of the risk and comparison with the national cancer incidence rate. Therefore, extensive measurements are needed for a reliable comparison.
It is also essential to measure the amount of natural radiation in each area as this can determine the suitability of the environment for a healthy lifestyle. The indoor and outdoor gamma exposure rates in the air 1 m above the ground from terrestrial radionuclides and cosmic rays in Mashhad are 155.73 ± 13.92 nGy h−1 and 126.15 ± 15.66 nGy h−1, respectively [28]. Using a determined conversion factor as 0.7 Sv Gy−1, converting the absorbed doses to effective doses [1], the annual indoor and outdoor effective dose rates of the public from gamma exposure were found to be 0.95 ± 0.08 and 0.77 ± 0.09 mSv yr−1. Therefore, by comparing these values with the corresponding annual effective doses from the inhalation of radon and thoron (3.7 ± 0.5 mSv yr−1), it could be concluded that most of the doses received indoors in the dwellings studied in Mashhad city are from the inhalation of radon and thoron (about 79 % of the total dose). A comparison of the indoor radon concentration and radiation risk assessment under investigation, with those reported in other Iranian cities also provided in Table 3.

3.3. Spatial Distribution Map of Indoor Radon Concentrations

The Distance Weighting (IDW) and Ordinary Kriging (OK) techniques, known as Kriging techniques, depend on the distance between two points, namely those of observation and estimation in the interpolation. IDW weighted the contribution of the observed points on the estimated interpolation with regard to this distance alone. On the other hand, OK also considers the correlation between the points and forms an initial function, i.e., covariance or variogram, which can iteratively be updated.
The spatial distribution map of indoor radon concentrations in Mashhad dwellings were plotted in Figure 7 by various interpolation techniques, e.g., IDW, OK, and Empirical Bayesian Kriging (EBK) over a grid with the dimensions of 1 km × 1 km using ArcGIS software version 10.7. Accordingly, radon concentrations were lower than standard values in eastern residential areas and were higher in central as well as southern districts. Nevertheless, when the spatial autocorrelation between cells was considered, predictions about radon concentrations using different methods range from 65 to 260 Bq m−3. These values may be more realistic and similar to average values found in some dwellings in the region.
Moreover, the accuracies of the various techniques applied according to five indicators are given in Table 4. IDW, which predicts unknown values using known values concerning their distance, was proven to be more suitable for predicting mean indoor radon concentrations over grids with the dimensions of 1 km × 1 km (i.e., arithmetic mean, ground floor), due to the lower MAE and RMSLE values of 28.159 and 0.01210, respectively, in addition to a lower bias, 20.069 to be exact. However, all mentioned models have a tendency to overestimate bias (PB > 0). In addition, the model with the higher R2 is IDW, which indicates that this model fits the data better.

4. Conclusions

To estimate the impact of indoor radon and thoron on residentials as well as develop and implement the most economical method to reduce radon exposure using a radon map, this paper presents the measured indoor radon and thoron concentrations in 78 dwellings as well as soil gas radon concentrations in different districts of Mashhad, Iran during summer and winter. As the average of the radon concentrations in the summer and winter were assumed to be the annual average radon concentration in this study, the annual average indoor radon and thoron concentrations were calculated as being 132 ± 19 and 58 ± 7 Bq m−3, respectively. Soil gas radon concentrations also ranged from 1.07 ± 0.28 to 8.02 ± 0.65 kBq·m−3 with a mean value of 3.07 ± 1.09 kBq·m−3 during the summer.
The corresponding annual effective dose from the inhalation of radon and thoron was calculated as being 3.7 ± 0.5 mSv yr−1. Subsequently, the excess lifetime cancer risk was calculated as 14.13. Hence, exposure to indoor radon is responsible for approximately 12% of lung cancer deaths in Mashhad, which is close to the WHO estimates of the worldwide proportion of lung cancer due to radon (3–14%). By comparing the annual indoor effective dose rate from gamma exposure with the annual effective dose from the inhalation of radon and thoron, it was concluded that most of the dose received inside the dwellings studied in Mashhad, approximately 79% of the total dose, originates from the inhalation of radon and thoron.
Since high-risk areas can be recognized on radon maps, which are useful for targeting landlords and the building industry, an indoor radon map was generated by using ArcGIS software over a grid with the dimensions of 1 km × 1 km using three interpolation techniques. The arithmetic mean was used over the grid cells to predict a mean indoor radon concentration on the ground-floor level of buildings in the grid cells where no data was available. The IDW technique was proven to be most suitable one for predicting mean indoor radon concentrations over grids with the dimensions of 1 km × 1 km.
In addition to the results and given the significant health impacts of radon and thoron, it is hoped that both radon and thorn gases will be studied more seriously in Iran and that these techniques as well as complementary procedures will be used to minimize its concentration. It is recommended that radon gas concentrations should be measured in all regions of the country by numerous devices supplied by the Atomic Energy Organization of Iran (AEOI). As a result, it would be possible to compile a radon map of Iran to estimate the concentration and number of radon-induced incidences of cancer, as well as decide how to distribute the population. This would aid to reduce the number of cases of lung cancer and other radon-induced human health problems.

Author Contributions

Conceptualization, M.A., A.S., and T.K.; methodology, M.A., A.S., and T.K.; software, M.A. and M.I.; validation, M.A., A.S., T.K., and S.C.; formal analysis, M.A., A.S.; investigation, M.A., A.S., and T.K.; resources, T.K.; data curation, T.K. and S.C.; writing—original draft preparation, M.A.; writing—review and editing, M.A., A.S., S.C., and T.K.; visualization, M.A. and A.S.; supervision, T.K.; project administration, T.K.; funding acquisition, T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the TKP2020-IKA-07 project financed under the 2020-4.1.1-TKP2020 Thematic Excellence Programme by the National Research, Development and Innovation Fund of Hungary and by the Hungarian National Research OKTA grant No. K128805 and K128818.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Not applicable for studies not involving humans.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United Nations Scientific Committee on the Effects of Atomic Radiation. UNSCEAR 2000 Report, Sources 473 and Effects of Ionizing Radiation; Annex B Exposures from Natural Radiation Sources; UNSCEAR: New York, NY, USA, 2000; Volume 1. [Google Scholar]
  2. United Nations Scientific Committee on the Effects of Atomic Radiation. UNSCEAR 2008 Report, Sources and Effects of Ionizing Radiation; Annex B Exposures of the Public and Workers from Various Sources of Radiation; UNSCEAR: New York, NY, USA, 2010; Volume 1. [Google Scholar]
  3. Shahrokhi, A. Application of the European Basic Safety Standards Directive in Underground Mines: A Comprehensive Radioecology Study in a Hungarian Manganese Mine. Ph.D. Thesis, Pannon Egyetem, Veszprem, Hungary, 2018. [Google Scholar]
  4. World Health Organization. WHO Handbook on Indoor Radon: A Public Health Perspective; WHO: Geneva, Switzerland, 2009. [Google Scholar]
  5. United States Environmental Protection Agency. EPA Protocol for Conducting Radon and Radon Decay Product Measurements in Multifamily Buildings; Environmental Protection Agency: Washington, DC, USA, 2017.
  6. Dubois, G. An Overview of Radon Surveys in Europe; European Commission, Ed.; EUR 21892 EN; Publications Office of the European Union: Luxembourg, 2005. [Google Scholar]
  7. Rahman, S.; Mati, N.; Ghauri, B.M. Seasonal indoor radon concentration in the North West Frontier Province and federally administered tribal areas—Pakistan. Radiat. Meas. 2007, 42, 1715–1722. [Google Scholar] [CrossRef]
  8. Rafique, M.; Rahman, S.U.; Mahmood, T.; Rahman, S.; Matiullah. Assessment of seasonal variation of indoor radon level in dwellings of some districts of Azad Kashmir, Pakistan. Indoor Built Environ. 2011, 20, 354–361. [Google Scholar] [CrossRef]
  9. Kim, Y.; Chang, B.U.; Park, H.M.; Kim, C.K.; Tokonami, S. National radon survey in Korea. Radiat. Prot. Dosim. 2011, 146, 6–10. [Google Scholar] [CrossRef]
  10. Szeiler, G.; Somlai, J.; Ishikawa, T.; Omori, Y.; Mishra, R.; Sapra, B.K.; Mayya, Y.S.; Tokonami, S.; Csordás, A.; Kovács, T. Preliminary results from an indoor radon thoron survey in Hungary. Radiat. Prot. Dosim. 2012, 152, 243–246. [Google Scholar] [CrossRef]
  11. Sainz-Fernandez, C.; Fernandez-Villar, A.; Fuente-Merino, I.; Gutierrez-Villanueva, J.L.; Martin-Matarranz, J.L.; Garcia-Talavera, M.; Casal-Ordas, S.; Quindós-Poncela, L.S. The Spanish indoor radon mapping strategy. Radiat. Prot. Dosim. 2014, 162, 58–62. [Google Scholar] [CrossRef]
  12. Shahrokhi, A.; Shokraee, F.; Reza, A.; Rahimi, H. Health risk assessment of household exposure to indoor radon in association with the dwelling’s age. J. Radiat. Prot. Res. 2015, 40, 155–161. [Google Scholar] [CrossRef] [Green Version]
  13. Singh, P.; Singh, P.; Singh, S.; Sahoo, B.K.; Sapra, B.K.; Bajwa, B.S. A study of indoor radon, thoron and their progeny measurement in Tosham region Haryana, India. J. Radiat. Res. Appl. Sci. 2015, 8, 226–233. [Google Scholar] [CrossRef] [Green Version]
  14. Cinelli, G.; Tositti, L.; Capaccioni, B.; Brattich, E.; Mostacci, D. Soil gas radon assessment and development of a radon risk map in Bolsena, Central Italy. Environ. Geochem. Health 2015, 37, 305–319. [Google Scholar] [CrossRef] [Green Version]
  15. Elío, J.; Crowley, Q.; Scanlon, R.; Hodgson, J.; Long, S. Logistic regression model for detecting radon prone areas in Ireland. Sci. Total Environ. 2017, 599, 1317–1329. [Google Scholar] [CrossRef]
  16. Hoffmann, M.; Aliyev, C.S.; Feyzullayev, A.A.; Baghirli, R.J.; Veliyeva, F.F.; Pampuri, L.; Valsangiacomo, C.; Tollefsen, T.; Cinelli, G. First map of residential indoor radon measurements in Azerbaijan. Radiat. Prot. Dosim. 2017, 175, 186–193. [Google Scholar] [CrossRef]
  17. Park, J.H.; Lee, C.M.; Lee, H.Y.; Kang, D.R. Estimation of seasonal correction factors for indoor radon concentrations in Korea. Int. J. Environ. Res. Public Health 2018, 15, 2251. [Google Scholar] [CrossRef] [Green Version]
  18. Sohrabi, M.; Babapouran, M. New public dose assessment from internal and external exposures in low and elevated-level natural radiation areas of Ramsar, Iran. Int. Congr. Ser. 2005, 1276, 169–174. [Google Scholar] [CrossRef]
  19. Hadad, K.; Doulatdar, R.; Mehdizadeh, S. Indoor radon monitoring in northern Iran using passive and active measurements. J. Environ. Radioact. 2007, 95, 39–52. [Google Scholar] [CrossRef]
  20. Hadad, K.; Mokhtari, J. Indoor radon variations in central Iran and its geostatistical map. Atmos. Environ. 2015, 102, 220–227. [Google Scholar] [CrossRef]
  21. Bouzarjomehri, F.; Ehrampoosh, M.H. Radon level in dwellings basement of Yazd-Iran. Iran. J. Radiat. Res. 2008, 6, 141–144. [Google Scholar]
  22. Mowlavi, A.A.; Fornasier, M.R.; Binesh, A.; De Denaro, M. Indoor radon measurement and effective dose assessment of 150 apartments in Mashhad, Iran. Environ. Monit Assess. 2012, 184, 1085–1088. [Google Scholar] [CrossRef]
  23. Fahiminia, M.; Fard, R.F.; Ardani, R.; Naddafi, K.; Hassanvand, M.S.; Mohammadbeigi, A. Indoor radon measurements in residential dwellings in Qom, Iran. Int. J. Radiat Res. 2016, 14, 331–339. [Google Scholar] [CrossRef] [Green Version]
  24. Yarahmadi, M.; Shahsavani, A.; Mahmoudian, M.H.; Shamsedini, N.; Rastkari, N.; Kermani, M. Estimation of the residential radon levels and the annual effective dose in dwellings of Shiraz, Iran, in 2015. Electron. Phys. 2016, 8, 2497–2505. [Google Scholar] [CrossRef]
  25. Shahbazi Sehrani, M.; Boudaqpoor, S.; Mirmohammadi, M. Measurement of indoor radon gas concentration and assessment of health risk in Tehran, Iran. Int. J. Environ. Sci. Technol. 2019, 16, 2619–2626. [Google Scholar] [CrossRef]
  26. Mirbag, A.; Shokati Poursani, A. Indoor radon measurement in residential/commercial buildings in Isfahan city. J. Air Pollut. Health 2019, 3, 209–218. [Google Scholar] [CrossRef]
  27. Shahrokhi, A.; Adelikhah, M.; Chalupnik, S.; Kocsis, E.; Toth-Bodrogi, E.; Kovács, T. Radioactivity of building materials in Mahallat, Iran–an area exposed to a high level of natural background radiation–attenuation of external radiation doses. Materiales de Construcción 2020, 70, e233. [Google Scholar] [CrossRef]
  28. Sohrabi, M.; Roositalab, J.; Mohammadi, J. Public effective doses from environmental natural gamma exposures indoors and outdoors in Iran. Radiat. Prot. Dosim. 2015, 167, 633–641. [Google Scholar] [CrossRef] [PubMed]
  29. Tokonami, S.; Takahashi, H.; Kobayashi, Y.; Zhuo, W.; Hulber, E. Up-to-date radon-thoron discriminative detector for a large-scale survey. Rev. Sci. Instrum. 2005, 76, 113505. [Google Scholar] [CrossRef]
  30. Adelikhah, M.E.; Shahrokhi, A.; Chalupnik, S.; Toth-Bodrogi, E.; Kovács, T. High level of natural ionizing radiation at a thermal bath in Dehloran, Iran. Heliyon 2020, 6, e04297. [Google Scholar] [CrossRef]
  31. Shahrokhi, A.; Nagy, E.; Csordás, A.; Somlai, J.; Kovács, T. Distribution of indoor radon concentrations between selected Hungarian thermal baths. Nukleonika 2016, 61, 333–336. [Google Scholar] [CrossRef] [Green Version]
  32. Csordás, A.; Fábián, F.; Shahrokhi, A.; Somlai, J.; Kovács, T. Calibration of CR-39-based thoron progeny device. Radiat. Prot. Dosim. 2014, 160, 169–172. [Google Scholar]
  33. International Commission on Radiological Protection. Lung Cancer Risk from radon and Progeny and Statement on Radon; Annals of the ICRP: New York, NY, USA, 2010. [Google Scholar]
  34. Özen, S.; Celik, N.; Dursun, E.; Taskın, H. Indoor and outdoor radon measurements at lung cancer patients’ homes in the dwellings of Rize Province in Turkey. Environ. Geochem. Health 2018, 40, 1111–1125. [Google Scholar] [CrossRef]
  35. Bossew, P.; Tollefsen, T.; Gruber, V.; De Cort, M. The European radon mapping project. In Proceedings of the IX Latin American IRPA Regional Congress on Radiation Protection and Safety, IRP, Rio de Janeiro, Brazil, 14–19 April 2013. [Google Scholar]
  36. Janik, M.; Bossew, P.; Kurihara, O. Machine learning methods as a tool to analyse incomplete or irregularly sampled radon time series data. Sci. Total Environ. 2018, 630, 1155–1167. [Google Scholar] [CrossRef]
  37. Alexander, D.L.; Tropsha, A.; Winkler, D.A. Beware of R2: Simple, unambiguous assessment of the prediction accuracy of QSAR and QSPR models. J. Chem. Inf. Model. 2015, 55, 1316–1322. [Google Scholar] [CrossRef] [Green Version]
  38. International Commission on Radiological Protection. Protection Against Radon-222 at Home and at Work; Annals of the ICRP: New York, NY, USA, 1993. [Google Scholar]
  39. United States Environmental Protection Agency. Assessment of Risk from Radon in Homes; Environmental Protection Agency: Washington, DC, USA, 2003.
  40. Mehra, R.; Badhan, K.; Kansal, S.; Sonkawade, R.G. Assessment of seasonal indoor radon concentration in dwellings of Western Haryana. Radiat. Meas. 2011, 46, 1803–1806. [Google Scholar] [CrossRef]
  41. Dey, G.K.; Das, P.K. Estimation of radon concentration in dwellings in and around Guwahati. J. Earth Syst. Sci. 2012, 121, 237–240. [Google Scholar] [CrossRef]
  42. Kapdan, E.; Altinsoy, N. A comparative study of indoor radon concentrations between dwellings and schools. Radiat Phys. Chem. 2012, 81, 383–386. [Google Scholar] [CrossRef]
  43. Cosma, C.; Cucoş-Dinu, A.; Papp, B.; Begy, R.; Sainz, C. Soil and building material as main sources of indoor radon in Băiţa-Ştei radon prone area (Romania). J. Environ. Radioact. 2013, 116, 174–179. [Google Scholar] [CrossRef]
  44. Kunovska, B.; Ivanova, K.; Stojanovska, Z.; Vuchkov, D.; Zaneva, N. Measurements of radon concentration in soil gas of urban areas, Bulgaria. Rom. J. Phys. 2013, 58, S172–S179. [Google Scholar]
  45. Szabó, K.Z.; Jordan, G.; Horváth, Á.; Szabó, C. Mapping the geogenic radon potential: Methodology and spatial analysis for central Hungary. J. Environ. Radioact. 2014, 129, 107–120. [Google Scholar] [CrossRef]
  46. Hasan, A.K.; Subber, A.R.; Shaltakh, A.R. Measurement of radon concentration in soil gas using RAD7 in the environs of Al-Najaf Al-Ashraf City-Iraq. Adv. Appl. Sci. Res. 2011, 2, 273–278. [Google Scholar]
  47. Mittal, S.; Rani, A.; Mehra, R. Estimation of radon concentration in soil and groundwater samples of Northern Rajasthan, India. J. Radiat. Res. Appl. Sci. 2016, 9, 125–130. [Google Scholar] [CrossRef] [Green Version]
  48. Kikaj, D.; Jeran, Z.; Bahtijari, M.; Stegnar, P. Radon in soil gas in Kosovo. J. Environ. Radioact. 2016, 164, 245–252. [Google Scholar] [CrossRef]
  49. Vaupotič, J.; Gregorič, A.; Kobal, I.; Žvab, P.; Kozak, K.; Mazur, J.; Kochowska, E.; Grządziel, D. Radon concentration in soil gas and radon exhalation rate at the Ravne Fault in NW Slovenia. Nat. Hazard. Earth Sys. 2010, 10, 895–899. [Google Scholar] [CrossRef]
  50. Duggal, V.; Rani, A.; Mehra, R. Measurement of soil-gas radon in some areas of northern Rajasthan, India. J. Earth Syst Sci. 2014, 123, 1241–1247. [Google Scholar] [CrossRef]
  51. Al Mugahed, M.; Bentayeb, F. Studying of Radon Gas Concentrations in Soil Qaa Al-Hakel Agricultural Area, Ibb, Yemen. Mater. Today Proc. 2019, 13, 525–529. [Google Scholar] [CrossRef]
  52. Roshandel, G.; Ghanbari-Motlagh, A.; Partovipour, E.; Salavati, F.; Hasanpour-Heidari, S.; Mohammadi, G.; Khoshaabi, M.; Sadjadi, A.; Davanlou, M.; Tavangar, S.M.; et al. Cancer incidence in Iran in 2014: Results of the Iranian National Population-based Cancer Registry. Cancer Epidemiol. 2019, 61, 50–58. [Google Scholar] [CrossRef]
Figure 1. Location of the city of Mashhad in Razavi Khorasan Province, Iran.
Figure 1. Location of the city of Mashhad in Razavi Khorasan Province, Iran.
Ijerph 18 00141 g001
Figure 2. Normal Distribution of indoor (A) radon and (B) thoron concentrations (Bq m−3) over the two seasons at the ground level of dwellings examined in Mashhad.
Figure 2. Normal Distribution of indoor (A) radon and (B) thoron concentrations (Bq m−3) over the two seasons at the ground level of dwellings examined in Mashhad.
Ijerph 18 00141 g002
Figure 3. Log-Normal Distribution of indoor (A) radon and (B) thoron concentrations (Bq m−3) over the two seasons at the ground level of dwellings examined in Mashhad.
Figure 3. Log-Normal Distribution of indoor (A) radon and (B) thoron concentrations (Bq m−3) over the two seasons at the ground level of dwellings examined in Mashhad.
Ijerph 18 00141 g003
Figure 4. Correlation between the indoor radon and thoron concentrations of the 78 houses examined in Mashhad over two seasons.
Figure 4. Correlation between the indoor radon and thoron concentrations of the 78 houses examined in Mashhad over two seasons.
Ijerph 18 00141 g004
Figure 5. The soil-gas radon concentrations (kBq m−3) in the studied area.
Figure 5. The soil-gas radon concentrations (kBq m−3) in the studied area.
Ijerph 18 00141 g005
Figure 6. The correlation between indoor radon and soil-gas radon concentrations.
Figure 6. The correlation between indoor radon and soil-gas radon concentrations.
Ijerph 18 00141 g006
Figure 7. Predicted Indoor Radon Map of Mashhad dwellings over a grid with the dimensions of 1 km × 1 km using the (A) Inverse distance weighting, (B) Empirical Bayesian Kriging, and (C) Ordinary Kriging interpolation techniques.
Figure 7. Predicted Indoor Radon Map of Mashhad dwellings over a grid with the dimensions of 1 km × 1 km using the (A) Inverse distance weighting, (B) Empirical Bayesian Kriging, and (C) Ordinary Kriging interpolation techniques.
Ijerph 18 00141 g007
Table 1. Basic Statistics of indoor and soil gas 222Rn and 220Rn concentrations in samples from Mashhad.
Table 1. Basic Statistics of indoor and soil gas 222Rn and 220Rn concentrations in samples from Mashhad.
SeasonParameterA.M 1G.M 2S.D 3MinMaxWinter/Summer Ratio (Mean)
222Rn220Rn
SummerIndoor air radon (Bq m−3)115.02105.850.6450.8305.21.311.36
Indoor air thoron (Bq m−3)48.7337.427.95<LLD122.5
Soil-gas radon (kBq m−3)3.072.711.6211.0788.021
WinterIndoor air radon (Bq m−3)150.3139.6862.7475.3376.6
Indoor air thoron (Bq m−3)66.1749.4134.24<LLD166.3
1 A.M = Arithmetical mean, 2 G.M = Geometrical mean, 3 S.D = Standard deviation.
Table 2. Comparison of soil-gas radon concentrations under investigation with those in other countries using different methods and sampling depths.
Table 2. Comparison of soil-gas radon concentrations under investigation with those in other countries using different methods and sampling depths.
RegionRadon in Soil-Gas (KBq m−3)Measurement MethodSampling Depth (cm)Reference
Bǎita-Stei, Romania5.5–512Lucas Cell40–80[43]
Bolsena, Italy7–176RAD 760–70[14]
Bulgaria3–97AlphaGuard100[44]
Hungary1–47.1RAD 780[45]
Najaf, Iraq0.009–9.29RAD 75–60[46]
Rajasthan, India0.94–10.05RAD 7100[47]
Sharr-Korabi, Kosovo0.295–32SSNTDs (CR-39)80[48]
Slovenia0.9–32.9AlphaGuard100[49]
Sri Ganganagar, India0.9–10.10RAD 710–100[50]
Yemen0.15–13.56SSNTDs (CR-39)0–150[51]
Mashhad, Iran1.07–8.02SSNTDs (CR-39)50–60Present study
Table 3. Studies on indoor radon concentration (Bq m−3) and radiation health risk in various Iranian cities.
Table 3. Studies on indoor radon concentration (Bq m−3) and radiation health risk in various Iranian cities.
RegionNumber of DwellingMean Radon Concentration (SD 1)Mean Effective Dose (mSv yr−1)ELCR2LCC3 × 10−6Excessive Rate (%)Reference
Isfahan5128.57 (39.38)0.722.7 × 10−112.964% > 100 Bq m−3[26]
Lahijan400163 (57)3.431.3 × 10−261.74In most dwellings > 100 Bq m−3[19]
Mashhad14831.9(0.25–3.78)--5.3% of apartments > 100 Bq m−3[22]
Qom12395.832.419.2 × 10−343.3824.3% > 100 Bq m−3[23]
Ramsar500Autumn: 355
Winter: 476
Autumn: 8.95
Winter: 12
3.44 × 10−2
4.6 × 10−2
161.11
216
-[18]
Shiraz18557.6 (33.06)1.455.6 × 10−326.15.4% > 100 Bq m−3[24]
Tehran301042.621 × 10−247.1638% > 100 Bq m−3[25]
Yazd84137.4 (149.5)3.461.3 × 10−262.2830% of basements> 148 Bq m−3[21]
Mashhad78Summer: 115(51)
Winter: 150 (62)
Summer: 3.1
Winter: 4.2
(12.3 × 10−3)
(15.9 × 10−3)
(56.7)
(74.1)
Summer: 20% > 148 Bq m−3
Winter: 31% > 148 Bq m−3
Current study
1 SD= Standard deviation, 2 ELCR= Excess Lifetime Cancer Risk, 3 LCC= Lung Cancer Cases.
Table 4. Summary of the cross-validation results.
Table 4. Summary of the cross-validation results.
MethodMAE 1RMSE 2RMSLE 3PB 4R2
Inverse distance weighting28.15934.9310.0121020.0690.234
Empirical Bayesian Kriging28.23535.1480.0121820.1230.224
Ordinary Kriging28.42436.3640.0134620.2680.169
1 MAE = mean absolute error, 2 RMSE = root mean square error, 3 RMSLE = root mean squared logarithmic error, 4 PB= percentage bias.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adelikhah, M.; Shahrokhi, A.; Imani, M.; Chalupnik, S.; Kovács, T. Radiological Assessment of Indoor Radon and Thoron Concentrations and Indoor Radon Map of Dwellings in Mashhad, Iran. Int. J. Environ. Res. Public Health 2021, 18, 141. https://doi.org/10.3390/ijerph18010141

AMA Style

Adelikhah M, Shahrokhi A, Imani M, Chalupnik S, Kovács T. Radiological Assessment of Indoor Radon and Thoron Concentrations and Indoor Radon Map of Dwellings in Mashhad, Iran. International Journal of Environmental Research and Public Health. 2021; 18(1):141. https://doi.org/10.3390/ijerph18010141

Chicago/Turabian Style

Adelikhah, Mohammademad, Amin Shahrokhi, Morteza Imani, Stanislaw Chalupnik, and Tibor Kovács. 2021. "Radiological Assessment of Indoor Radon and Thoron Concentrations and Indoor Radon Map of Dwellings in Mashhad, Iran" International Journal of Environmental Research and Public Health 18, no. 1: 141. https://doi.org/10.3390/ijerph18010141

APA Style

Adelikhah, M., Shahrokhi, A., Imani, M., Chalupnik, S., & Kovács, T. (2021). Radiological Assessment of Indoor Radon and Thoron Concentrations and Indoor Radon Map of Dwellings in Mashhad, Iran. International Journal of Environmental Research and Public Health, 18(1), 141. https://doi.org/10.3390/ijerph18010141

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

Article Metrics

Back to TopTop