Radionuclide concentrations in sand samples from riverbanks of Muzaffarabad, Azad Kashmir

NUCLEAR ENERGY SCIENCE AND ENGINEERING

Radionuclide concentrations in sand samples from riverbanks of Muzaffarabad, Azad Kashmir

Abdul Razzaq Khan，

Muhammad Rafique，

Abdul Jabbar，

Saeed Ur Rahman，

Muhammad Ikram Shahzad，

Muhammad Ejaiz Khan，

Mulaika Yasin

Nuclear Science and Techniques

Vol.29, No.7

Article number 93

Published in print 01 Jul 2018

Available online 28 May 2018

DOI：10.1007/s41365-018-0442-9

41800

This paper presents the results of a radiological risk assessment arising from the presence of naturally occurring radionuclides in sand samples from three riverbanks in Muzaffarabad. The mean values obtained for ²³²Th, ²²⁶Ra, and ⁴⁰K were found to be 44.58 ± 3.34, 48.25 ± 1.77, and 239.92 ± 22.73 Bq.kg^-1, respectively. To assess the uniformity of exposure, the radium equivalent activity (Ra_eq) was calculated and was found to be 130.47 ± 8.29 Bq.kg^-1. The current reported value for Ra_eq is lower than the maximum permissible value, that is, 370 Bq.kg^-1, and equivalent to a gamma dose of 1.5 mSv y^-1. To investigate the possible contribution to health risks of alpha particle exposure, the radon exhalation rate (RER) from the sand samples was determined. The mean RER for all the samples was found to be 335 mBq.m^-2 h^-1. About 43% of the samples were found to have an indoor excess lifetime cancer risk value slightly higher than recommended safety limit of 1, as proposed by the ICRP. A normalized parameter, the equivalent multiplicative factor (EMF), was defined and used to compare the current results with the findings of studies performed in other countries. Our findings are relevant to both human health and future environmental radiation monitoring.

radiological riskradon exhalation rateabsorbed dose ratehazard indexannual effective dose

1. INTRODUCTION

If building materials, either of natural origin or formed from industrial by-products, have elevated levels of radioactivity, they can contribute to radiation exposure of the population. Almost all building materials with a rock or soil origin, contain various amounts of radionuclides. Long-lived radionuclides with half-lives comparable to the age of the earth, e.g., those derived from either ²³⁸U or ²³²Th series and ⁴⁰K, can still be found today in the Earth’s crust. ²²⁶Ra has the greatest radiological importance in the ²³⁸U decay series, causing it to be referenced very often instead of ²³⁸U [1]. The consequences of radionuclides being present in building materials is two-fold, that is, external irradiation of the body from direct gamma radiation and internal irradiation of the lung tissues from ²²²Rn, ²²⁰Rn (thoron) and its alpha-emitting progenies [2]. Human exposure arising from direct gamma radiation, called external exposure, is considered to be the largest external dose contributor to the world’s population [3,4]. Building materials with elevated levels of ²²⁶Ra, ²³²Th, and ⁴⁰K may be responsible for external doses.

On the other hand, internal exposure is mainly caused by alpha-emitting radionuclides. Radon (²²²Rn), thoron (²²⁰Rn) and their short-lived progenies are responsible for internal exposure. Since uranium is present in all types of building materials in trace amounts, with radon being its sixth decay product in series, the presence of radon within indoor environments is regarded as being natural. Although the main source of indoor radon (IR) is the underlying soil, in some cases the building materials may contribute significantly to the IR levels. As mentioned in the European Commission Radiation Protection 112 report, typical excess IR levels due to building materials are about 10–20 Bq.m^-3, although this value may escalate and become greater than 1000 Bq.m^-3 in some rare circumstances. In addition to adding radon to the indoor environment, building materials are also believed to be the most significant source of indoor thoron (IT). However, the importance of IT is limited due to its short half-life, but its concentration solely depends upon the level of thorium in building materials [5]. While keeping the importance of radiological hazards in mind, as they relate to the presence of radionuclides in building materials, many researchers across the globe have conducted studies to identify and mitigate this problem [6-25]. Regarding the building materials themselves, sand is an integral part of many materials and it is usually formed by the weathering of rocks. Broadly, it is classified as pit sand, river sand, and sea sand. Sand from rivers in Muzaffarabad is commonly used for building. It is obtained from the banks or bed of the rivers. It consists of fine rounded grains.

The radon exhalation rate plays an important role in characterizing building materials and their contribution to indoor radon levels [26]. Within mineral grains ²²²Rn is produced from the α–particle decay of ²²⁶Ra. Since there are very few pore spaces, a large amount of radon remains within the mineral grains and only a small fraction of it can escape. During the radon emanation process, the radon escapes from the grains of the material in to the pores [27,28]. The fraction of radon escaping from the pore spaces is called the "emanation coefficient" (EC). Radon emanation depends upon different parameters including the ²²⁶Ra distribution and amount within the grains, the water content within the pore spaces, the porosity, and the grain size [29,30]. For different materials, the value of EC ranges from 0.05 to 0.7. In the case of sand EC ranges from 0.15 to 0.3 [31]. The fraction of radon reaching the surface of a building material is subject to diffusion and convective flow processes. This liberated quantity of radon per unit surface area of building material per unit time is called the "radon exhalation rate" (RER) [32].

The rationale of the current study was to recognize and compute the significant gamma emitting radionuclide’s concentration in sand samples collected from three rivers in Muzaffarabad. The study of the radon exhalation rate and gamma spectroscopy of the sand samples can allow us to assess the radiological risk due to the inhalation of airborne radon emanating from the sand, building material, and external human exposure to gamma radiation originating from naturally occurring radionuclides. The radionuclide concentration is found using gamma ray spectroscopy while the radon exhalation rate is determined using the closed "CAN" technique [13,15,32].

2. Experimental Procedures

2.1 Sample Collection

The samples of sand were collected from the banks of the Neelum, Jhelum, and Khunhar rivers. A total of thirty sand samples were collected from along a 40-km stretch of each river. Ten sand samples were collected from each river. The study area thus spanned four districts, namely, Muzaffarabad, Hattian Bala, Neelum, and Manshera (KPK). To collect a sand sample from a site, the upper sand layer was removed to access the sand that is local to that site. (The upper layers are subject to movement as a result of the action of water.) Then, the polythene bag was filled with sand and sealed after being indexed. Every sample bag was marked with a sample code. In this way, a total of 30 samples were collected, each from a different location. To find the coordinates of the sampling locations, a global positioning system (GPS) was used.

2.2 Sample Preparation

Sand samples, weighing 200 g, were pretreated in the Solid State Nuclear Tracks Laboratory of the University of Azad Jammu and Kashmir. Initially, the samples were dried by exposing them to sunlight. The sand samples were sieved through a sieve with a mesh size of less than 0.16 mm. Samples were then placed in an oven at 110 °C for 4–6 hours. After the removal of the moisture, the samples were then packed in a cylindrical polyethylene container. This container was airtight, so that secular equilibrium between the parent and daughter radionuclides were achieved within one month. At this point, each sample was placed in the detector to measure gamma ray activity by using an HPGe detector.

2.3 Gamma spectrometric analysis

A PC-based high-resolution gamma spectrometry system was used for the radiometric analysis of the heat-treated sand samples [33]. This system used a high-purity germanium (HPGe) coaxial detector with an efficiency of 30%, relative to a NaI(Tl) detector. To reduce the effect of background noise, the detector was shielded by 15-cm lead with an inner lining of 3-mm copper and 4-mm tin. [12,18,21,34,35]. The calibration of the system was performed experimentally using IAEA soil-326. The minimum detectable activity (MDA) for each of the radio nuclides ²²⁶Ra, ²³²Th, ⁴⁰K, and ¹³⁷Cs were determined using the following equation:

Minimum Detectable Activity (M D A) = \frac{4.66 {(Continum Counts+ Back Ground Peak Counts)}^{\frac{1}{2}}}{Sample Mass (kg) \times Efficiency \times Live time(s) \times Yield}

(1)

where, MDA is expressed in units of Bq.Kg^-1and 4.66 is the statistical coverage factor. The MDA values for radionuclides ²²⁶Ra, ²³²Th, ¹³⁷Cs, and ⁴⁰K were found to be 3.60, 2.25, 1.35, and 6.70 Bq.Kg^-1, respectively.

The activity concentrations due to ²²⁶Ra, ²³²Th, ¹³⁷Cs, and ⁴⁰K were calculated using the relationship proposed by Fredj et al. [36],

A_{S} (Bq/sample) ​ ​ = \frac{(Net peak counts for sample -Background counts)}{counting time×detection efficiency×emission probability of-ray}

(2)

3. Determination of Gamma Hazard Indices and Dose Rate

Hazard indices, i.e., the radium equivalent activity, Ra_eq, internal, external hazard indices (H_in, H_ex), gamma index (I_γ), alpha index (I_α), and the representative level index I_r, were calculated to assess the radiation hazards associated with the individual sand samples. In addition, the absorbed dose rate, $\dot{D}$ , in the air as well as the annual effective dose were estimated. The formulae for each of said indices and the dose rate (from Eq. 3–10) are listed in Tables 1 and 2 [37–40,2,4,5].

Formulae for determination of hazard indices

Hazard Index		Characteristics		Reference
Radium equivalent activity (Ra_eq)
	$R a eq .= A Ra+(1 .43 \times A Th)+(0 .077 \times A K)$ (3)	The criterion Ra_eq⪯ 370 Bq.kg^-1 should be met to guarantee that the external dose rate (D) remains less than or equal to 1.5 mGy y⁻¹.		[37]
External Hazard Index (H_ex)
	$H_{ex} = \frac{A_{Ra}}{370} + \frac{A_{Th}}{259} + \frac{A_{K}}{4810}$ (4)	To minimize radiation hazards, the External Hazard Index must be less than unity i.e., H_ex⪯ 1		[37]
Internal Hazard Index (H_in)
	$H_{in} = \frac{A_{Ra}}{185} + \frac{A_{Th}}{259} + \frac{A_{K}}{4810}$ (5)	For minimizing radiation hazards, Internal Hazard Index must be less than unity for i.e., H_in⪯ 1		[37]
Gamma Index (I_γ)
	$I_{γ} = \frac{A_{Ra}}{300} + \frac{A_{Th}}{200} + \frac{A_{K}}{3000}$ (6)	For bulk material		[5]
			Iγ ⪯ 0.5 ⇒ D⪯ 0.3 mSv y^-1	[38]
			Iγ ⪯ 1.0 ⇒ D ⪯ 1.0 mSv y^-1
		For restricted material (Tiles etc.)
			Iγ ⪯ 2 ⇒ D ⪯ 0.3 mSv y^-1
			Iγ ⪯ 6 ⇒ D ⪯ 1.0 mSv y^-1
Alpha Index (I_α)
	$I_{α} = \frac{A_{Ra}}{200}$ (7)	Exemption level = 100 Bq.kg^-1		[38,39]
		Recommended upper level = 200 Bq.kg^-1.
Representative level index (I_r)
	$I_{r} = \frac{A_{Ra}}{150} + \frac{A_{Th}}{100} + \frac{A_{K}}{1500}$ (8)	To keep radiation, hazard insignificant, I_r must be less than unity.		[2]

A_Ra, A_Th, and A_K represent the activity concentrations of ²²⁶Ra, ²³²Th and ⁴⁰K respectively, in Bq.kg^-1.

Dose rate, annual effective dose due to radionuclides ²³²Th, ²²⁶ Ra and ⁴⁰K

Hazard Index		Characteristics	Reference
Absorbed gamma dose rate in air (nGy h^-1)
	$\dot{D} = (0.462 \times A_{Ra}) + (0.604 \times A_{Th}) + (0.0417 \times A_{K})$ (9)	Assumption: All decay products of ²²⁶Ra and ²³²Th are in radioactive equilibrium with their precursors.	[4]
Annual effective dose eq. E (mSv y^-1)
	$E = Q \times T \times O \times \dot{D} \times 10^{- 6}$ (10)	Q=0.7 Sv. Gy^–1	[40]
"Q" is a conversion coefficient "T" is the time in hours,		T= 8760 hours in one year $\dot{D}$ is in nGy h^-1
" $\dot{D}$ " is the dose rate, "O" is the occupancy factor.

Since radionuclides are not uniformly distributed in building materials, a quantity known as the radium equivalent activity, Ra_eq, is defined to represent an equivalent uniform exposure to radiation due to radium, thorium, and potassium, in units of Bq/Kg. External hazards are defined for the external γ- radiation dose from building materials, whilst an internal hazard index is estimated to determine the radiation hazard to respiratory organs due to exposure to radon and its short-lived progenies. The gamma index (I_γ) and alpha index (I_α) are defined for the excess external gamma and excess α-radiation exposures, respectively. From a radiation safety perspective, a criterion defining an upper level for each of these parameters has been suggested. These criteria are shown in Table 1 in the 2^nd column, titled "Characteristics." The absorbed gamma dose rate in air and the annual effective dose due to radium, thorium, and potassium radionuclides are listed in Table 2.

4. Radon exhalation rate (RER) determination

The RER value was determined for thirty sand samples using the 'CAN’ technique [13,15,32]. To this end, all the samples were dried for 4 h in an oven at 110 °C. Samples, weighting 200 g, were put into plastic cans with a volume of $8.55 \times 10^{3}$

cm³. Sheets of CR-39 detectors were cut into small strips measuring 2 × 2 cm² and placed inside the NRPB detector holder. Closing the holder prevents all radon decay products, as well as aerosol particles, from entering the active volume of the NRPB detector holder. That is, only radon can enter the sensitive volume and, when it decays, the resulting α-particles interact with the CR-39 detectors. The CR-39 detectors were tacked to the at upper inner surface of cans 25 cm from the samples. The thirty cans were then hermetically sealed and left for 30 days. After exposure to radon for 30 days, the CR-39 detectors were retrieved and etched in 6N NaOH at 80 °C. The etched detectors were cleaned with fresh tap water and then the tracks were counted with the help of an optical microscope. The measured track densities were corrected for background contributions.

The calibration of the system was carried out at the PINSTECH Laboratory using the same method as that used by Matiullah [41]. Uranium ore samples of a known grade (to act as a source) were placed in plastic containers with a volume of 8.55 × 10³ cm³ and the CR-39 detectors, in NRPB-type holders, were exposed to the uranium samples for four weeks. After retrieval, the detectors were etched, as described above, to determine the track densities. A calibration factor of 2.713 tracks cm⁻² h⁻¹/kBq.m⁻³ was obtained and subsequently used to find the radon concentrations.

4.1 RER Measurement Methodology

The CR-39 detectors within the NRPB dosimeters were exposed for a period of 30 days, resulting in a variable radon exposure. The experimental setup for the measurement of the RER is shown in Figure 1.

Figure. 1

(Color online) Experimental setup for measuring radon exhalation rate from sand samples [13,15, 32]

The effective detector exposure time to radon was calculated using the following relationship [32].

T_{effective} = Total exposure time in Days (t) - ​ Radon ​ Mean Life Time (τ) [1 - e^{- λ t}]

(11)

where, $τ$ (radon mean life time) = 5.5 days, $λ$ is ²²²Rn decay constant.

Durrani and Ilic stated that this type of correction is needed only for a closed system [31].

The RER was determined using the following relationship [13,15,32].

{(R E R)}_{0} = \frac{C (t) [ω + λ V / A]}{[1 - e^{- (\frac{ω A}{V} + λ) t}]}

(12)

Equation (12) was corrected to account for the back-diffusion factor. The radon exhalation rate under back diffusion conditions is given by Eq. 13 [32],

R E R = {(R E R)}_{0} - ω C ​ = {(R E R)}_{0} - ε λ Z_{0}

(13)

where the surface area of the sample was denoted by $A$ and, within the closed chamber, the volume of the void spaces was denoted by $V$ . The radon accumulation time within chamber is denoted by $t$ . Here, w=ελZ₀ was used for the back-diffusion constant depending upon the given material type. $Z_{o}$ was the thickness of the sample. Parameter $C (t)$ represents the concentration of ²²²Rn, that has to be exhaled.

5. Results and Discussion

Ten samples were collected from banks of each of the Kunhar, Neelum, and Jhelum rivers and were then analyzed to determine the gamma ray spectroscopy and radon exhalation rate. The obtained activity concentration (Bq.kg^-1) of ²³²Th, ²²⁶Ra, and ⁴⁰K, on a dry weight basis, and the radon exhalation rate, for 30 samples of sand are presented in Figures 2-8.

Figure. 2

(Color online) Activity Concentration of ²²⁶Ra, ²³²Th, and ⁴⁰K in sand samples collected from the banks of the Jhelum, Kunhar, and Neelum rivers.

Figure. 8

Annual effective dose (mSv/y) of sand samples

5.1 Determination of ²²⁶Ra, ²³²Th, and ⁴⁰K activity levels

For the sand samples collected from the Jhelum river, the activity concentration due to ²²⁶Ra varied from 13.75 ± 1.55 to 59.87 ± 1.77 should be avoided ^-1with a mean value of 40.82 ± 1.71 Bq.kg^-1. Similarly, the activity concentrations of ²³²Th and ⁴⁰K varied from 26.48 ± 3.26 to 64.39 ± 3.35 with a mean value of 38.08 ± 3.28 Bq.kg^-1 and 118.99 ± 22.12 to 265.86 ± 23.28Bq.kg^-1 with a mean value of 192.87 ± 22.41 Bq.kg^-1, respectively. The radon exhalation rate varied from 139.89 to 411.04 mBq.m^-2 h^-1with a mean value of 289.92 mBq.m^-2 h^-1.

Similarly, for those sand samples collected from the Kunhar river, the activity concentration due to ²²⁶Ra varied from 37.95 ± 1.8 to 69.19 ± 1.83 Bq.kg^-1with a mean value of 54.15 ± 1.81 Bq.kg^-1. Similarly, the activity concentrations of ²³²Th and ⁴⁰K varied from 25.09 ± 3.41 to 57.89 ± 3.48 Bq.kg^-1with a mean value of 44.55 ± 3.36 Bq.kg^-1and 192.31 ± 22.85 to 385.8 ± 23.47 Bq.kg^-1with a mean value of 273.01 ± 23.03Bq.kg^-1 respectively. The radon exhalation rate varied from 261.94 to 504.31 mBq.m^-2 h^-1 with a mean value of 375.24 mBq.m^-2 h^-1.

For the sand samples collected from the Neelum river, the activity concentration due to ²²⁶Ra varied from 35.04 ± 1.61 to 61.78 ± 1.83 Bq.kg^-1with a mean value of 49.78 ± 1.77 Bq.kg^-1. Similarly, the activity concentrations of ²³²Th and ⁴⁰K varied from 32.95 ± 3.27 to 66.55 ± 3.75Bq.kg^-1 with a mean value of 51.14 ± 3.38 Bq.kg^-1and 132.17 ± 22.48 to 402.59 ± 25.14Bq.kg^-1 with a mean value of 253.9 ± 22.74 Bq.kg^-1. The radon exhalation rate varied from 227.97 to 417.38 mBq.m^-2 h^-1 with a mean value of 340.18 mBq.m^-2 h^-1. Figure 2 shows the average values for ²²⁶Ra, ²³²Th, and ⁴⁰K in sand samples collected from the banks of the Jhelum, Kunhar, and Neelum rivers. Figure 3 shows the average values for ²²⁶Ra, ²³²Th, and ⁴⁰K in all three rivers.

Figure. 3

(Color online) Average values of ²²⁶Ra, ²³²Th, and ⁴⁰K in sand samples.

5.2 Radium Equivalent Activities (Ra_eq) and Radiological Indices

The Ra_eq value of the sand samples collected from the banks of the Jhelum river varied from 60.81 ± 7.96 to 168.49 ± 8.23 Bq.kg^-1with a mean value of 110.07 ± 8.13 Bq.kg^-1 (see Figure 4).

Figure. 4

Radium equivalent activity (Ra_eq) of sand samples

The external (H_ex) and internal hazard (H_in) indices varied from 0.164 ± 0.0215 to 0.455 ± 0.0222 Bq.kg^-1with a mean value of 0.297 ± 0.0219 Bq.kg^-1and 0.201 ± 0.026 to 0.617 ± 0.027 Bq.kg^-1with a mean value of 0.408 ± 0.0266 Bq.kg^-1, respectively (see Figure 5). The gamma (Iγ) and alpha (I_α) indices varied from 0.218 ± 0.029 to 0.593 ± 0.029Bq.kg^-1with a mean value of 0.391 ± 0.029 Bq.kg^-1and 0.069 ± 0.008 to 0.299 ± 0.009 Bq.kg^-1with a mean value of 0.204 ± 0.009 Bq/kg, respectively. The representative level index (I_r) varied from 0.436 ± 0.052 to 1.186 ± 0.054 Bq.kg^-1with a mean value of 0.781 ± 0.053 Bq.kg^-1(see Figure 6).

Figure. 5

External and internal hazard index

Figure. 6

Alpha, gamma, and representative level index for sand samples

For the sand samples collected from the banks of the Kunhar river, the radium equivalent activity was found to have a minimum value of 97.14 Bq/kg and a maximum value of 168.16 Bq/kg with an average value of 138.88 Bq/kg.

The values of the external (H_ex) and internal (H_in) hazard indices ranged from 0.26 ± 0.023 to 0.45 ± 0.025 with a mean value of 0.375 ± 0.022 and from 0.365 ± 0.28 to 0.640 ± 0.030 with a mean value of 0.521 ± 0.275, respectively. The gamma (Iγ) and alpha (I_α) indices varied from 0.353 ± 0.0312 to 0.591 ± 0.033 with a mean value of 0.559 ± 0.031 and from 0.190 ± 0.009 to 0.346 ± 0.009 with a mean vale of 0.270 ± 0.009. The value of the representative index ranged between 0.706 ± 0.057 and 1.182 ± 0.061 with an average value of 0.988 ± 0.056.

The value of Ra_eq for the sand samples collected from the Neelum river was found to be in a range of 107.876 ± 8.19 to 179.153 ± 9.22 Bq/kg, with an average value of 142.468 ± 8.358 Bq/kg (see Figure 4).

The values of the external (H_ex) and internal (H_in) hazard indices in the sand samples taken from the river Neelum varied from 0.29 ± 0.022 to 0.48 ± 0.024 with a mean value of 0.38 ± 0.022 and from 0.39 ± 0.025 to 0.63 ± 0.030 with a mean value of 0.52 ± 0.027, respectively. The estimated values of the gamma (Iγ) and alpha (I_α) indices for the Neelum river varied from 0.38 ± 0.029 to 0.64 ± 0.034 with a mean value of 0.51 ± 0.030 and from 0.18 ± 0.008 to 0.31 ± 0.008 with a mean value of 0.25 ± 0.009. For this river, the value of the representative index varied from 0.75 ± 0.055 to 1.29 ± 0.061 with a mean value of 1.013 ± 0.055.

The overall value of Ra_eq for all of the 30 studied samples ranged from 60.81 ± 7.67 to 179.15 ± 9.22 Bq/kg with a mean value of 130.47 ± 8.29 Bq/kg (see Figure 4).

The values of the external (H_ex) and internal (H_in) hazard indices for all the sand samples from all three rivers varied from 0.164 ± 0.021 to 0.484 ± 0.025 with a mean value of 0.352 ± 0.022 and from 0.201 ± 0.025 to 0.64 ± 0.03 with a mean value of 0.483 ± 0.027, respectively. The overall estimated ranges of the gamma (Iγ) and alpha (I_α) indices were found to be 0.218 ± 0.028, 0.644 ± 0.034 with a mean value of 0.464 ± 0.03 and 0.069 ± 0.008, and 0.346 ± 0.01 with a mean value of 0.241 ± 0.009. For all the samples, the value of the representative index ranged from 0.436 ± 0.051 to 1.287 ± 0.061 with a mean value of 0.927 ± 0.055.

For all the sand samples, the ⁴⁰K content was very much dependent on the location, with the lowest value being 118.9 Bq.kg^–1, that being obtained for the sand sample collected from the bank of the Jhelum river at site SJ3. The maximum value of 402.59 Bq.kg^–1was obtained for the sample collected from site SN3 on the bank of the Neelum river. The higher value for ⁴⁰K could perhaps be due to the occurrence of K-feldspar within the mineral matrix of the sand deposits in the study area.

5.3 Estimation of Absorbed Dose Rate (ADR) and Total Annual Effective Dose (AED)

For those sand samples collected from the banks of the Jhelum river, the absorbed gamma dose rate in air (nGyh^-1) at a point 1 m from the ground, and the annual effective dose equivalent E (mSv y^-1) varied from 27.32 ± 3.63 to 75.51 ± 3.74nGy h^-1with a mean value of 49.88 ± 3.71 nGyh^-1and 0.168 ± 0.022 to 0.463 ± 0.023 mSv y^-1 with a mean value of 0.306 ± 0.023 mSv y^-1, respectively (see Figures 7 and 8).

Figure. 7

Absorbed gamma dose rate in air (nGyh^-1) at point 1 m from the ground

The gamma ray absorbed dose rate in air at a point 1 m above the ground was also calculated for the samples collected from the Kunhar river, and was found to range from 45.314 ± 3.917 to 75.831 ± 4.716 nGyh^-1 with an average value of 63.311 ± 3.829 nGyh^-1. The value of the annual effective equivalent dose E varied from 0.278 ± 0.0240 to 0.465 ± 0.0256 mSv y^-1 whereas the mean value of this dose for the Kunhar river was found to be 0.388 ± 0.023 mSv y^-1.

The estimated value of the gamma ray absorbed dose rate in air at a point 1 m above the ground for the samples collected from the banks of the Neelum river was found to range from 48.78 ± 3.74 to 81.46 ± 4.20 nGyh^-1 with an average value of 64.48 ± 3.81 nGyh^-1. The value of the annual effective equivalent dose E was found to be in a range of 0.29 ± 0.023 to 0.49 ± 0.026 mSv y^-1 with a mean value of 0.39 ± 0.023 mSv y^-1.

The total gamma ray absorbed dose rate (GDR) in air at a point 1 m above the ground for all the studied samples was found to range from 27.324 ± 3.501 to 81.466 ± 4.2 nGyh^-1 with an average value of 59.222 ± 3.782 nGyh^-1. The total annual effective equivalent dose E for the three rivers varied from 0.168 ± 0.021 to 0.5 ± 0.026 mSv y^-1 with a mean value 0.363 ± 0.023 mSv y^-1.

5.4 Excess Life Time Cancer Risk (ELCR)

The gamma-ray-induced cancer risk was calculated, based on the values obtained for the annual effective dose equivalent (AEDE) and other factors, using Eq. 14.

\begin{array}{l} Excess lifetime cancer risk (ELCR) \\ = AEDE \times Average duration of life (DL) \\ \times Risk Factor (RF) \end{array}

(14)

where, AEDE, DL, and RF are the annual effective dose equivalent, duration of life (66 years) [42], and risk factor (Sv^-1), together giving the fatal cancer risk per Sievert. For low-dose background radiation which is considered to produce a stochastic effect, ICRP 60 uses values of 0.05 for public exposure [43,44].

The ELCR defines the probability of developing cancer once over an entire lifetime for a given exposure level. ELCR due to gamma radiation exposure from sand samples varied from 0.55 × 10^-3 to 1.65 × 10^-3 with an average value of 1.20 × 10^-3. ELCR for indoor exposure was found to range from 0.44 × 10^-3 to 1.32 × 10^-3 with an average value of 0.96 × 10^-3. For outdoor exposure, ELCR varied from 0.11 × 10^-3 to 0.33 × 10^-3 with an average value of 0.24 × 10^-3.

The outdoor values of the excess lifetime cancer risk are within the recommended safety limit of 1 [43] while about 43% of the indoor excess lifetime cancer risk is slightly higher than 1. Hence, dwellers of buildings constructed using the sand addressed in the study area are at risk, to a certain degree, of developing cancer over their lifetimes.

5.5 Annual gonadal equivalent dose (AGED)

Different organs of human body, e.g., the gonads, bone marrow, and the bone surface cells are radio-sensitive. Hence, for any level of radiation exposure, these organs are vulnerable [4]. The annual gonadal dose equivalent (AGED) was calculated using the specific activities of ²²⁶Ra, ²³²Th, and ⁴⁰K according to Eq. 15 [45].

AGED ({μSv y}^{- 1}) = 3.09 C_{Ra} + 4.18 C_{Th} + 0.314 C_{K}

(15)

where, C_Ra, C_Th, and C_K are the specific activity concentrations of radium, thorium, and potassium, respectively. The average values of the annual gonadal equivalent dose varied from 0.191 ± 0.025 to 0.568 ± 0.029 mSv y^-1with an average value of 0.411 ± 0.026 mSv y^-1. The average value of the AGDE, resulting from gamma exposure from the sand samples, is higher than the global average value of 0.298 mSv·y⁻¹ reported for soil [39].

5.6 Correlation between radon exhalation rate (RER) and radium content

The correlation between the radon exhalation rate (RER) and radium content in every sand sample was investigated (see Figure 9). The linear regression mathematical model (Y = A + B × X) for the RER and ²²⁶Ra content correlation is given as;

Figure. 9

Linear regression analysis of RER with ²²⁶Ra.

^{226} Ra activity concentration (Bq/kg) = 0 .14867 \times Value of RER - 1 .569 83

(16)

The linear correlation analysis between the RER and ²²⁶Ra activity gives a slope of 0.14867 and an R² value of 0.93. The coefficient of determination (R²) is a statistical measure expressing how close the data are to a fitted regression line. Wherever the radium content of a sand sample is relatively high, the RER is also high for that sample, as would be expected since ²²²Rn is a decay product of radium, and there is a good correlation between them. A positive correlation of R² = 0.93 has been observed between the radium content and the radon exhalation rate.

5.7 Geo-Spatial Analysis

As the measurements were representative of the samples collected from a specific location, the acquisition of information about the unsampled locations would have to rely on geostatistical methods [46]. To represent the spatial distributions of the natural radionuclides concentration present in the sand samples of the study area, the Kriging interpolation method was used. The Kriging method (named after D. Krige (Krige 1966)), was further developed by Matheron (Matheron 1970) and is now regarded as being the best linear unbiased estimator (BLUE) [47]. The Kriging formula is given by Eq. 17 [46];

Z (S_{0}) = \sum_{i = 1}^{N} λ_{i} Z (S_{i})

(17)

where, $Z (S_{i})$ is the estimated value at the ith place, $λ_{i}$ is the unknown weight for the measured value at the ith location (ith), and $S_{0}$ is the measurement location.

For the current study, the ²²⁶Ra, ²³²Th, and ⁴⁰K activity concentrations are shown in Figs. 10, 11 and 12. The GDR for the current study area is shown in Figure 13.

Figure. 10

(Color online) 2-D contour map of the ²²⁶Ra activity concentrations obtained in the study

Figure. 11

(Color online) 2-D contour map for ²³²Th activity concentrations obtained in the study

Figure. 12

(Color online) 2-D contour map for ⁴⁰K activity concentrations obtained in the study

Figure. 13

(Color online) 2-D Contour map of gamma dose rate in study area

5.8 Normalized Radionuclide Concentration or Equivalent Multiplicative Factor (EMF)

Comparisons of the results obtained in the current study with those of studies conducted in other parts of the world [49-70,11] are shown in Figures 14 and 15. The results of other studies are compared with those of the current study by defining the equivalent multiplicative factor (EMF) or normalized values of the radionuclide concentrations.

Figure. 14

Equivalent Multiplicative Factor (EMF) of ⁴⁰K for comparison with data in literature

Figure. 15

(Color online) Equivalent Multiplicative Factor (EMF) of ²³²Th and ²²⁶Ra for comparison with data in literature

Equivalent Multiplicative Factor(EMF) = \frac{Current Study Values of^{226} Ra,^{232} Th,^{40} K}{Values of^{226} Ra,^{232} Th,^{40} K in literature}

(18)

These EMF values show that the ²²⁶Ra and ²³²Th concentrations in the sand samples collected from different rivers in Muzaffarabad are somewhat higher than the data available in the literature. Similarly, the radionuclide ⁴⁰K Equivalent Multiplicative Factor for the sand samples collected from different rivers in Muzaffarabad is relatively greater, but in some cases smaller, for some countries [49–70,11].

For the majority of the analyzed sand samples, the presence of ⁴⁰K does not present any possible radiological risks to human health because the concentration of ⁴⁰K is in a homeostatic balance with other isotopes of potassium in the atmosphere.

6. Conclusion

The NORMSs concentration and RER in sand samples with associated radiological hazards were investigated. The ELCR values for outdoor exposure were found to be within the ICRP recommended safety limit of 1, whilst about 43% of the samples were found to have slightly higher values than the recommended limit. On the other hand, 90% of the sand samples were found to have average AGED values, resulting from gamma exposure from the sand samples, that is higher than the world average of 0.298 mSv·y⁻¹ reported for soil. Furthermore, 37% of the samples were found to have a representative level index (I_r) that is slightly higher than recommended value of 1. Hence, the use of sand from the study area, given its elevated concentration of radionuclides, as a construction material requires careful regulation to minimize radiation hazards. The Ra_eq, external, and internal hazard indices, as well as the alpha and gamma index for the majority of the samples were found to be less than 1, while the results for RER shows that most of the sand samples available in area can be utilized as building materials and do not pose any significant radiation hazard to building inhabitants or other persons. It is recommended that the use of sand with high radionuclide concentrations as a construction material should be avoided.

References

NCRP,

Exposures from the Uranium Series with Emphasis on Radon and Its Daughters. (National Council on Radiation Protection and Measurements, 7910 Woodmont Avenue 1Bethesda, Report No. 77, 1984)

Radionuclide concentrations in sand samples from riverbanks of Muzaffarabad, Azad Kashmir

1. INTRODUCTION

2. Experimental Procedures

2.1 Sample Collection

2.2 Sample Preparation

2.3 Gamma spectrometric analysis

3. Determination of Gamma Hazard Indices and Dose Rate

4. Radon exhalation rate (RER) determination

4.1 RER Measurement Methodology

5. Results and Discussion

5.1 Determination of 226Ra, 232Th, and 40K activity levels

5.2 Radium Equivalent Activities (Raeq) and Radiological Indices

5.3 Estimation of Absorbed Dose Rate (ADR) and Total Annual Effective Dose (AED)

5.4 Excess Life Time Cancer Risk (ELCR)

5.5 Annual gonadal equivalent dose (AGED)

5.6 Correlation between radon exhalation rate (RER) and radium content

5.7 Geo-Spatial Analysis

5.8 Normalized Radionuclide Concentration or Equivalent Multiplicative Factor (EMF)

6. Conclusion

5.1 Determination of ²²⁶Ra, ²³²Th, and ⁴⁰K activity levels

5.2 Radium Equivalent Activities (Ra_eq) and Radiological Indices