Study and Simulated the Natural Radioactivity (NORM) U-238, Th-232 and K-40 of Igneous and Sedimentary Rocks of Al-Atawilah (Al-Baha) in Saudi Arabia

Bashayer M. Al-Zahrani^1*, Haifa S. Alqannas², Safia H. Hamidalddin²

¹Physics Department, Faculty of Science, Al Baha University, Al Baha, Saudi Arabia.
²Physics Department, Faculty of Science, University of Jeddah, Jeddah, Saudi Arabia.

DOI: 10.4236/wjnst.2020.104015                

Abstract

In this work, gamma-ray spectroscopy based on semiconductor hyper pure germanium (HPGe) detector was used to evaluate the activity concentrations of the natural radionuclides (U-238 (Ra-226), Th-232 and K-40) and the fallout nuclide (Cs-137) for thirty samples of igneous and sedimentary rocks of Al-Atawilah (Al-Baha). The mean values of the activity concentrations of U-238 (Ra-226), Th-232, K-40 and Cs-137 in the igneous samples are found as (11.0, 11.50, 1172.71, 1.47) Bq/Kg respectively. In the sedimentary rocks, the mean values of the activity concentrations of the natural radionuclides (U-238 (Ra-226), Th-232 and K-40) and the fallout nuclide (Cs-137) equal to (12.04, 13.18, 1131.36, 1.60) Bq/Kg respectively. The averages of radiological hazards (Ra_eq, H_ex and I_γ) were calculated and found to be within the UNSCEAR permissible limit values (370 Bq/kg for Ra_eq, and 1 for H_ex and I_γ), except for a slight increase of average value of I_γ in the igneous rock samples (1.36). The results indicate that the dose rate values depend on the kind of rocks (high in some igneous rock samples, and most of sedimentary rock samples have low dose rate). The activities of naturalnuclides were predicted and simulated in T time using a written MATLAB R2020a script based on the average activity concentrations and respective half-lives of U-238 and Th-232 series, and K-40, this is to evaluate the future effects of natural radionuclides on the population and estimate the human inputs in the future.

Keywords

Natural Radioactivity, Gamma-Ray Spectroscopy, Activity Concentration, Igneous and Sedimentary Rocks, Simulated in T Time

1. Introduction

Radiation activity exists everywhere on the surface of Earth and its interior. Uranium-238, thorium-232 (and their progenies) and potassium-40 are most important sources of radiation. These nuclei are found in any type of rocks, especially in igneous and sedimentary rocks. U-238 decays by ejection an alpha particle to generate daughter radionuclide Th-234 which followed by other decays to produce other radionuclides such as Ra-226 and its progeny Rn-222. Similarly, Th-232 disintegrates to produce Rn-228 and followed by other radionuclides [1]. These radionuclides constitute risks by the external exposure to gamma radiation emissions and internally by radon and its progenies. Radon is a human carcinogen, and it is considered the second leading cause of lung cancer [2]. Knowledge of radionuclide distribution is important because it gives helpful information in the observation of natural environmental radioactivity and connected external exposure that resulting from gamma radiation primary based on the geological and geographical conditions and can be seen at various levels in the rocks of each area of the world [3]. The aim of this study is to determine the radionuclides activity concentration of Ra-226, Th-232 and K-40 for igneous and sedimentary rock samples collected from Al-Atawilah (north of Al-Baha region), and to estimate the doses and hazard indices originate from the existence of the natural radionuclides in the surrounding area. Moreover, Ra-226, Th-232 and K-40 activity concentrations were used to simulate and predict the range of decay by using an appropriate program. The results are used to assess the future effect of these radionuclides and evaluate the future of radiation hazards.

2. Methodology

2.1. Sampling and Samples Preparation

Eighteen igneous rock samples and twelve sedimentary rock samples were collected from Al-Atawilah (20.273351˚N, 41.358325˚E), north of Al-Baha region, southeast of Saudi Arabia, Figure 1.

Figure 1. Locations of the igneous and sedimentary rock samples of Al-Atawilah region.

The samples were collected from the study region at depth ≈ 5 cm, by cracking the rock with hammer after removing a thin layer of the mother rock. Then, the collected samples were dried, pulverized and then sieved through less than 1mm-mesh size [4]. All crushed sieved rock samples were filled into Polyethylene Marinelli beakers. The samples were pulverized into fine-grained powder for uniform distribution of radon and its decay products and to avert any accumulation in the top [5]. The weight of each sample was recorded and then hermetically sealed from outside using a thick tape for more than a month so as to establish the secular equilibrium between Ra-226 and Th-232 (and their decay products) [6] [7].

2.2. Measurements Equipment

Measurement of radionuclides activity concentrations in the rock samples were evaluated by using high-resolution gamma-ray spectrometry system includes hyper-purity closed-end coaxial germanium detector. The detector had resolution of 2 keV at 1332.5 keV of Co-60 and 25% relative efficiency and peak to Compton ratio of 50:1. The detector is placed in a cylinder-heavy lead shield to minimize the radiation background values. Activity concentration calculations were based on the establishment of secular equilibrium in the measured sample between ²²⁶Ra and ²³²Th and their progenies of smaller lifetime [8]. ²²⁶Ra activity concentration was evaluated from the peak energies of 295.2 keV and 351.9 keV at ²¹⁴Pb and 609.3 keV, 1120.28 and 1764.5 at ²¹⁴Bi. The activity concentration of ²³²Th was assessed from γ-ray peaks of 338.3 keV, 911.1 keV and 968.96 keV at ²²⁸Ac, 238.6 keV at ²¹²Pb, 727.3 keV at ²¹²Bi and 583.1 keV and 860.5 at ²⁰⁸Tl. The activity concentration of ⁴⁰K was determined directly from the γ-ray line of 1460.8 keV at ⁴⁰Ar, and ¹³⁷Cs activity concentration was evaluated from peak energy of 661.6 keV at ¹³⁷Ba.

2.3. Equations for Calculations

The radionuclide activity concentration (A_c) in the investigated rocks samples were determined in Bq∙kg⁻¹. The activity concentration calculations were carried out using the following formula [9] [10]:

$A_c=N_c/m\beta \epsilon$ (1)

where: N_c is the net peak area per unit time (second), m is sample mass in kg, ε is the detector absolute efficiency at the photo-peak energy and β is the branching ratio of gamma radiation.

Radium equivalent (Ra_eq) is calculated by applying the following equation [11] [12]:

$R{a}_{eq}\left(\text{Bq}\cdot {\text{kg}}^{-1}\right)=C_{Ra}+1.43C_{Th}+0.077C_K$ (2)

where: C_Ra, C_Th and C_K are ²²⁶Ra, ²³²Th and ⁴⁰K activity concentrations, respectively.

External hazard index (H_ex) can be calculated from the following equation [10] [11]:

$H_{ex}=C_{Ra}/370+C_{Th}/259+C_K/4810$ (3)

To keep the radiation hazard insignificant, the calculated value of the external hazard index must be less than unity [13].

Representative level index (I_γ) is determined with following equation [14]:

$I_{\gamma }=\left(C_{Ra}/150\right)+\left(C_{Th}/100\right)+\left(C_K/1500\right)$ (4)

The absorbed dose rate (D_R) in air at 1m overhead the ground level was evaluated from the activity concentrations of the relevant natural radionuclides according to the following equation [15]:

$D_R\left(\text{nGy}\cdot {\text{h}}^{-1}\right)=0.462C_{Ra}+0.604C_{Th}+0.0417C_K$ (5)

2.4. Decay Simulation of Natural Radionuclides

The mean activity concentrations in the rock samples of Al-Atawilah were used to simulate and predict the range of decay. The decay of the radionuclides U-238 (Ra-226), Th-232 and K-40 of the rocks samples were simulated using written MATLAB R2020ascript according to the exponential law of radioactive decay [16]:

$A=A_0{\text{e}}^{-\lambda t}$ (6)

where: A is the change of the radioactive nuclei number with time, A₀ is the initial activity, t is the time and λ is the constant of the decay.

In this study, Forward Different Interpolation Method was applied to reconstruct the activity concentrations of radionuclides. The term e^−λt of the radionuclide decay equation “ $A=A_0{\text{e}}^{-\lambda t}$ ” (n to a 4^th order) was used in Taylor polynomial form. The decay factor e^−λt was approximated to a polynomial form by the following analysis for the fourth order:

$P_n\left(\lambda t\right)=P_n\left(z\right)$ (7)

Since:

$P_n\left(z\right)={\text{e}}^{-\lambda t}={\text{e}}^{-z}$

this yields the polynomial as:

$\begin{array}{l}{\text{e}}^{-z}=P_n\left(z\right)=a_0+a_1\left(z-z_0\right)+a_2\left(z-z_0\right)\left(z-z_1\right)\\ +a_3\left(z-z_0\right)\left(z-z_1\right)\left(z-z_2\right)\\ +a_4\left(z-z_0\right)\left(z-z_1\right)\left(z-z_1\right)\left(z-z_2\right)\left(z-z_3\right)\end{array}$ (8)

where

$a_0=y_0=P_0\left(z_0\right)$ (9)

$a_1=\left(y_1-y_0\right)/h=\Delta y_0/h$ (10)

$a_2=\left(y_2-2y_1+y_0\right)/2h^2={\Delta }^2{y}_0/2h^2$ (11)

$a_3=\left(y_3-3y_1+3y_1-y_0\right)/2!h^3$ (12)

$a_4=\left(y_4-4y_3+6y_2-4y_1+y_0\right)/4!h^4$ (13)

The coefficients a, b, c and e were calculated for the equation (8), and it was used with MATLAB R2020a to simulate the decay of the radionuclides ²³⁸U, ²³²Th and ⁴⁰K using their half-lives [17].

$P_n\left(z\right)=a{z}^4+b{z}^3+c{z}^2+dz+e$ (14)

3. Results and Discussion

3.1. Activity Concentrations

Analytical results for the samples have been applied to evaluate the activity concentration of ²²⁶Ra, ²³²Th and ⁴⁰K, and the artificial radionuclide (¹³⁷Cs) in Bq∙kg⁻¹ together with their total uncertainties. The results are presented in Table 1. In the igneous rock samples, the activity concentrations of ²²⁶Ra ranged from 5.86 to 17.50 Bq∙kg⁻¹ with average value of 11.00 Bq∙kg⁻¹, ²³²Th from 21.42 to 5.12 Bq∙kg⁻¹ with average value of 11.50 Bq∙kg⁻¹. For ⁴⁰K activity concentrations varied from 328.08 to 4854.00 Bq∙kg⁻¹ with average value of 1172.71 Bq∙kg⁻¹. The concentrations of the fallout nuclide ¹³⁷Cs varied from 5.63 to 0.90 Bq∙kg⁻¹ and the average value is 1.47 Bq∙kg⁻¹. ²²⁶Ra activity concentrations in the sedimentary rock samples varied from 5.23 to 54.85 Bq∙kg⁻¹ with average value of 12.04 Bq∙kg⁻¹. The highest value of ²³²Th is 61.26 Bq∙kg⁻¹ and lowest value is 5.15 Bq∙kg⁻¹, the average value is 13.18 Bq∙kg⁻¹. ⁴⁰K activity concentrations ranged from 110 Bq∙kg⁻¹ to 6070.75 Bq∙kg⁻¹ with mean value of 1131.36 Bq∙kg⁻¹. The maximum value of ¹³⁷Cs is 4 and the lowest value is 0.84 Bq∙kg⁻¹, and the average value is 1.60 Bq∙kg⁻¹. ¹³⁷Cs radionuclide is found in the most of the studied igneous and sedimentary rocks samples. Radium and Thorium activity concentrations of the studied igneous and sedimentary rocks varied depending on the types of the rocks. However, the average values of ²²⁶Ra and ²³²Th activity concentrations are less than the world average (50 Bq∙kg⁻¹) reported in UNSCEAR 2008. Potassium activity concentration (C_K) values are higher than recommended value in all igneous rock samples, except for two samples. The mean value of ⁴⁰K the igneous rock samples is 2.35 times greater than the acceptable average value of 500 Bq∙kg⁻¹ that recommended by UNSCEAR 2008. ⁴⁰K is the most important radionuclides identified, it is the extremely abundant natural radionuclide in all igneous rock samples under investigation. The average contribution of ⁴⁰K in the igneous rock samples is 98.1% while the average contributions of ²²⁶Ra and ²³²Th are 0.92% and 0.96%, respectively.

⁴⁰K activity concentration are measured higher than the acceptable value in all the sedimentary rock samples in this study except for four samples. The calculated mean value of ⁴⁰K is 2.2 times higher than the recommended value reported in UNSCEAR 2008 (500 Bq∙kg⁻¹). ⁴⁰K radionuclide is the most significant radionuclide had measured, it is most abundant natural radionuclide in all sedimentary rock samples under studied. The mean contribution of ⁴⁰K in the sedimentary rock samples is 97.82% while the mean contributions of ²²⁶Ra and ²³²Th are 1.04% and 1.14%, respectively.

Table 1. Values of activity concentrations, radium equivalent, external hazard index, γ-ray representative level index and absorbed dose rate, and their minimum, maximum and mean values of the rock samples.

^aSD: standard deviation; ^bND: not determined.

3.2. Radiation Hazards

Table 1 summarized the estimated values of radiation hazard indices (radium equivalent (Ra_eq), external hazard index (H_ex), gamma ray representative level index (I_γ) and absorbed dose rate (D_R)) for all igneous and sedimentary rock samples, with the average recommended values reported by UNSCEAR 2008 (370 Bq/kg for Ra_eq, and 1 for H_ex and I_γ). In the igneous rock samples, (Ra_eq) values are ranged from 42.53 to 406.73 Bq∙kg⁻¹, and the mean value is 163.83 Bq∙kg⁻¹. (H_ex) are ranged between 0.11 and 1.10 with average value of 0.44. The lowest value of (I_γ) is 0.33, while the highest value is 3.46. And the mean value is 1.36. The estimated total absorbed dose rate for the igneous rock samples are varied from 28.57 to 216.86 nGy∙h⁻¹, and the mean value is 86 nGy∙h⁻¹. The proportions of radionuclides contributions to the total absorbed dose rate in the igneous rock samples are 5.84% of ²²⁶Ra, 8.13% of ²³²Th and 86.04% of ⁴⁰K.

In the sedimentary rock samples, the lowermost value of (Ra_eq) equal to 23.45 Bq∙kg⁻¹ and the uppermost value is 610 Bq∙kg⁻¹, and the average value is 118.01 Bq∙kg⁻¹. (H_ex) values are varied between 0.06 and 1.64 respectively, with mean value of 0.31. The lowest value of (I_γ) is 0.18, where the highest value equal to 5.02, and the mean value of (I_γ) in the sedimentary rock samples is 0.96. The total absorbed dose rate values are ranged between 11.37 and 315.5 nGy∙h⁻¹, with a calculated mean value equal to 60.70 nGy∙h⁻¹, and the ratios of radionuclides contributions to the total absorbed dose rate in the sedimentary rock samples are 9.16%, 13.18% and 77.71% of ²²⁶Ra ²³²Th and ⁴⁰K, respectively.

3.3. Natural Radionuclides Decay Simulation of the Rock Samples

The decay of ²²⁶Ra (²³⁸U equivalent), ²³²Th and ⁴⁰K of the rock samples were predicted over 10³ years, 10⁶ years, 10⁹ years and finally over 10¹⁰ years using MATLAB R2020a script based on their present mean activity concentrations. Figures 2-5 display the exponential decay graphs of the natural radionuclides. The mean activity concentrations for all rocks samples under investigation are 11.35 Bq∙kg⁻¹, 12.17 Bq∙kg⁻¹ and 1516.17 Bq∙kg⁻¹ for ²²⁶Ra, ²³²Th and ⁴⁰K, respectively. The graphics have been zoomed in to clarify ²³⁸U (²²⁶Ra) and ²³²Th decay curves. The following observations from the simulation outcomes can be indicated as:

1) According to the Equation (6), the decay curves were expected to give an exponential graph, the lines approach zero if the background radiation is ignored.

2) The decay will be almost constant in the several next years in the study area, as the decay was not observed in Figure 2 and Figure 3, this due to the extremely long half-lives of the terrestrial radionuclides, thus the decay that the radionuclides will undergo during thousands or millions years will be insignificant.

3) From Figure 4 and Figure 5, a significant collapse was observed in the decay curve of ⁴⁰K compared to ²³²Th and ²³⁸U (²²⁶Ra), this due to its shorter half-life compared to ²³⁸U and ²³²Th half-lives.

Figure 2. The decay simulation of the natural radionuclides of rock samples over ×10³ years.

Figure 3. The decay simulation of the natural radionuclides of rock samples over ×10⁶ years.

Figure 4. The decay simulation of the natural radionuclides of rock samples over ×10⁹ years.

Figure 5. The decay simulation of the natural radionuclides of rock samples over ×10¹⁰ years.

4) This study evaluates the future effect of natural radionuclides on the population of this region and any increase in the concentration of radioactivity should be due to human inputs.

4. Conclusion

Gamma-ray spectroscopy of hyper-purity Germanium (HPGe) detector is a good experimental tool for studying levels of the radioactivity in various environmental samples such as rocks. The mean activity concentrations of ²²⁶Ra, ²³²Th and ⁴⁰K in the igneous rock samples are 11, 11.5, 1172.71 Bq∙kg⁻¹, respectively. For the sedimentary rock samples, the activity concentrations are found to be 12.04, 13.18, 1131.36 Bq∙kg⁻¹ for ²²⁶Ra, ²³²Th and ⁴⁰K, respectively. Fallout nuclide (¹³⁷Cs) was found in the most of the rock samples under investigation, the low mean values of ¹³⁷Cs are not of radiologically significant. The calculated mean values of radium equivalent (Ra_eq), external hazard index (H_ex), representative level index (I_γ) and absorbed dose rate (D_R) are within the suggested limit values, except for a slight increase of I_γ in igneous rock samples. This study also evaluated the radioactivity levels in the future and their effects on the population of the study region and any increase in the concentration of radioactivity should be due to human inputs. The results of the present study can help us to understand the distribution of natural radionuclides in the environment of Al-Baha region and provide a main map of radioactivity levels in Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Omar, M., Hamzah M.S. and Wood, A.K. (2008) Radioactive Disequilibrium and Total Activity Concentration of NORM Waste. Journal of Nuclear Science and Technology, 5, 47-56.
[2]	Alharbi, W.R. and Abbady, A.G. (2013) Measurements of Radon Concentrations in Soil and the Extent of Their Impact on the Environment from Al-Qassim, Saudi Arabia. Natural Science, 5, 93-98. https://doi.org/10.4236/ns.2013.51015
[3]	Akkurt, I. and Günoglu, K. (2014) Natural Radioactivity Measurements and Radiation Dose Estimation in Some Sedimentary Rock Samples in Turkey. Science and Technology of Nuclear Installations, 2014, Article ID: 950978. https://doi.org/10.1155/2014/950978
[4]	Iaea, I. (1989) Measurement of Radionuclides in Food and the Environment. International Atomic Energy Agency, Technical Report Series, No. 295.
[5]	Kerur, B.R., Rajeshwari, T., Sharanabasappa, S.A., Narayani, K., Rekha, A. and Hanumaiah, B. (2010) Radioactivity Levels of Rocks in North Karnataka, India. Indian Journal of Pure & Applied Physics, 48, 809-812.
[6]	Ahmad, N., Jaafar, M. and Alsaffar, M. (2015) Natural Radioactivity in Virgin and Agricultural Soil and its Environmental Implications in Sungai Petani, Kedah, Malaysia. Pollution, 1, 305-313.
[7]	Zubair, M. (2015) Measurement of Natural Radioactivity in Rock Samples Using Gamma Ray Spectrometry. Radiation Protection and Environment, 38, 11-13. https://doi.org/10.4103/0972-0464.162820
[8]	Achola, S.O. (2009) Radioactivity and Elemental Analysis of Carbonatite Rocks from Parts of Gwasi Area, South Western Kenya. Doctoral Dissertation, University of Nairobi, Nairobi.
[9]	Al-Zahrani, J. (2017) Estimation of Natural Radioactivity in Local and Imported Polished Granite Used as Building Materials in Saudi Arabia. Journal of Radiation Research and Applied Sciences, 10, 241-245. https://doi.org/10.1016/j.jrras.2017.05.001
[10]	Younis, H., Qureshi, A.A., Manzoor, S. and Anees, M. (2018) Measurement of Radioactivity in the Granites of Pakistan: A Review. Health Physics, 115, 760-768. https://doi.org/10.1097/HP.0000000000000917
[11]	Beretka, J. and Matthew, P. (1985) Natural Radioactivity of Australian Building Materials, Industrial Wastes and By-Products. Health Physics, 48, 87-95. https://doi.org/10.1097/00004032-198501000-00007
[12]	Rangaswamy, D., Srilatha, M., Ningappa, C., Srinivasa, E. and Sannappa, J. (2016) Measurement of Natural Radioactivity and Radiation Hazards Assessment in Rock Samples of Ramanagara and Tumkur District, Karnataka, India. Environmental Earth Science, 75, Article No. 373. https://doi.org/10.1007/s12665-015-5195-8
[13]	Bavarnegin, E., Moghaddam, M.V. and Fathabadi, N. (2013) Natural Radionuclide and Radiological Assessment of Building Materials in High Background Radiation Areas of Ramsar, Iran. Journal of Medical Physics/Association of Medical Physicists of India, 38, 93-97. https://doi.org/10.4103/0971-6203.111325
[14]	Nuclear Energy Agency (1979) Exposure to Radiation from the Natural Radioactive in Building Materials: Report. OECD, Paris.
[15]	United Nations Scientific Committee on the Effects of Atomic Radiation (2000) UNSCEAR 2000 Report to the General Assembly, with Scientific Annexes. UNSCEAR, New York.
[16]	Krane, K.S. (2019) Modern Physics. John Wiley & Sons, Hoboken.
[17]	Doyi, I., Essumang, D., Dampare, S., Duah, D. and Ahwireng, A. (2017) Evaluation of Radionuclides and Decay Simulation in a Terrestrial Environment for Health Risk Assessment. Scientific reports, 7, Article No. 16537. https://doi.org/10.1038/s41598-017-16659-w
[18]	United Nations Scientific Committee on the Effects of Atomic Radiation (2008) UNSCEAR 2008 Report to the General Assembly, with Scientific Annexes. UNSCEAR, New York.