The contamination from radioactive decay of radon and thoron gases among different sources of artificial and anthropogenic background radiation accounts for about 55% of the exposure. These radioactive components attach to airborne particles and may be retained in the lungs(UNSCEAR, 2000; ICRP, 2014; Mehra et al., 2016). Exhalation and emanation cycles cause radon and thoron gases to transit from their respective isotopes (found in soil or building materials) to the interior environment, resulting in radiation exposure to the public (Butterweck et al. 2002). Multiple seasonal variance data of radon, thoron, and other related contributing factors are limited in India, and these variables have been given lesser attention in investigations of radiological elements. Thoron’s impact on actual radiation dose is smaller; therefore, most studies have neglected its estimations. Experiments in several countries have revealed that thoron and its daughter 212Pb, which has a half-life of 10.6 hours, may not aggregate to deadly amounts in the air (Guo and Cheng, 2005; Ramola, 2016).
Dose rate due to natural gamma levels in air ranged from 0.526 to1.139 mSv y-1 having a mean value of 0.756 mSv y-1. The outcomes show that the mass exhalation rate varies significantly at different locations in the study area. The surface exhalation rate in the studied samples varied from 0.14 to 94.65 mBq kg-1 h-1with a mean value of 20.65 mBq kg-1 h-1 (Table 1). 4.2% samples have a higher value of mass exhalation rate than the global average value i.e. 57 mBq kg-1 h-1. The variation in exhalation depends on several factors such as soil type, parent nuclei content, emanation factor and radon diffusion coefficient. Mass exhalation rate was higher at Sotha village ( 94.65 mBq kg-1 h-1) and Kirmara village (86.27 mBq kg-1 h-1), while minimum at Nangthala village (0.14 mBq kg-1 h-1) (Table 1). The mass exhalation rate may be higher due to the presence of radon gas in the soil, which may be due to the deposition of fly ash from a nearby thermal Power Plant. Flyash from the power plant continuously gets mixed with the soil of the nearby region by different physical processes, which may increase the radium level in the soil (Chauhan et al. 2001). Due to radium contamination in soil, radon concentration gets higher in the soil. Radon is the decay product of radium, so there is a significant positive correlation between them which has been reported elsewhere (Pisapak et al. 2017; Chauhan et al. 2001). Uranium may also migrate from deeper soil horizons to the O-horizon with the help of phosphate fertilizers and some other micronutrients, which may increase the radium concentration in the soil. It may be another plausible reason for a higher concentration of mass exhalation rate than average in the study area.
The mass exhalation rate in the study area has been compared with available literature. It was found to be higher in the study area than in Northern Rajasthan; Karnal, Haryana; Mohali, Punjab; Kolhapur, Maharashtra [Duggal et al. (2015), Devi et al. (2020), Mehta et al. (2015), Raste et al. (2020) ]. In comparison, it was found to be almost equal to Amritsar and Tarn-taran, Punjab, and Himachal Pradesh [Kaur et al. (2018), Parminder Singh et al. (2015)]. The mass exhalation rate in the study area is lesser than in Kumaun, Himalayas; Uttarakhand and Faridabad, Haryana; central Haryana [Kumar et al. (2021), Anamika et al., (2020), Singh et al., (2020), Panghal et al. (2018)] (Table 3).
The study area’s Thoron surface exhalation rate ranged from 46.42 to 619.88 mBq m-2 h-1, with an arithmetic mean 161.73 mBq m-2 h-1 (Table 1), which is much lower than the global average value i. e. 3600 Bq m-2 h-1. All studied values were below the global average value. An extensive deviation in thoron surface exhalation rate may be due to the geology and ecology of the study area. The sampling sites may have unrelated geometries and soil structures that influence the thoron exhalation rate of the soil sample. The maximum Thoron surface exhalation rate was observed in a sample collected at Bhojraj village (619.88 mBq m-2 h-1) while the lowest in the sample collected from village Adampur (46.42 mBq m-2 h-1) (Table 1). The elevated concentration of 220Rn could be because of significant thorium activity at the sampling locations. The adjoining region of the sampling locations is the Aravali hills mining area (Tosham and Khanak). Surface exhalation rate contributes to indoor thoron concentration. In this region, the soil is widely used as a construction material so that a higher exhalation rate can contribute to a higher annual dose rate for the residents (Kumar et al. 2014). The selection of soil for the construction of houses is an essential parameter for reducing the radiation level in dwellings.
The surface exhalation rate obtained in this study has been compared with available literature. It was found to be higher than in Hamirpur, Himachal Pradesh; Ambala, Haryana; Mohali, Punjab; Khasi Hills, Meghalaya [Singh et al. (2015), Mehta and Shikha et al. (2016), Mehta et al. (2015), Pyngrope et al. (2022)] while it is almost equal to Kurukshetra, Haryana [Chauhan and Chakarvarti et al. (2002)]. The values are lesser than Kamaun, Himalaya; Uttarakhand, Faridabad, Haryana; central Haryana [Kumar et al. (2021), Anamika et al. (2020), Singh et al. (2020), Panghal et al. (2018] respectively (Table 3).
Figure 2(A) indicates that data distribution of overall gamma dose rate in the air is not consistent, while it is almost consistent for radon mass exhalation rate (Fig 3(A)). The histograms for the gamma dose rate and radon mass exhalation rate are slightly right-skewed and unimodal, which is clear from the thick tail for lower values as compared to the thin tail for higher values. Although skewed to the right, the histogram for thoron surface exhalation rate was distorted unimodal. The adverse impact of outlier values results in distorted unimodal distribution. The skewness and kurtosis values of gamma dose rate, mass exhalation rate and surface exhalation rate (Table 2) indicate that histogram tails are shorter and thinner, and central peaks are frequently lower and not sharp. This distribution is called platykurtic, where several figures are larger than the mode. It may be because some outlier values are on the higher side. The tail section may act as an outlier for the statistical model, and outliers have an unfavourable impact on the model’s presentation, particularly in regression-based models. This shape signifies that several figures are larger than the mode, possibly outliers. This is also supported by Box Plot (Figure 2(B), 3(B), 4(B)). The normal Q-Q plot (Figure 2(C)) showed that data of gamma dose rate in air is scattered almost symmetrical around the reference line with discontinuity. Asymmetrical radon mass exhalation rate and thoron surface exhalation rate distributions in the study area are also depicted by the normal Q-Q plot [Figure 3(C), 4(C)] where the scatter plot crosses the reference line two times.
The distribution of radon mass exhalation rate and thoron surface exhalation rate of soil samples vs outdoor gamma dose rate is given in Figure. 5. It is evident from the 3D scatter plot that there is not much effect on the dose rate with increasing concentrations of Rn222 and Rn220. The dose level changes most of the time arbitrarily. As a result, a weak positive correlation was seen between the outdoor gamma dose rate and the Rn222 mass exhalation rate, with R2 = 0.01 and R2 = 0.01 for the outdoor gamma dose rate and Rn220 surface exhalation rate. Thus an insignificant relationship was formed between the gamma dose rate from decay products and the radon/thoron exhalation rates.