Estimation of external dose by car-borne survey in Kerala, India.

A car-borne survey was carried out in Kerala, India to estimate external dose. Measurements were made with a 3-in × 3-in NaI(Tl) scintillation spectrometer from September 23 to 27, 2013. The routes were selected from 12 Panchayats in Karunagappally Taluk which were classified into high level, mid-level and low level high background radiation (HBR) areas. A heterogeneous distribution of air kerma rates was seen in the dose rate distribution map. The maximum air kerma rate, 2.1 μGy/h, was observed on a beach sand surface. 232Th activity concentration for the beach sand was higher than that for soil and grass surfaces, and the range of activity concentration was estimated to be 0.7-2.3 kBq/kg. The contribution of 232Th to air kerma rate was over 70% at the measurement points with values larger than 0.34 μGy/h. The maximum value of the annual effective dose in Karunagappally Taluk was observed around coastal areas, and it was estimated to be 13 mSv/y. More than 30% of all the annual effective doses obtained in this survey exceeded 1 mSv/y.

Introduction

Large amounts of artificial radionuclides were released from the reactor buildings into the environment by the 2011 Fukushima Dai-ichi Nuclear Power Plant accident [1]. External and internal dose estimations were carried out by government and university researchers immediately after the accident [2-4]. In a 2013 report, the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) reported that no apparent increase of cancer incidence due to this accident is expected in the future [5]. However, the biological effects due to radiation exposure have become a large concern among Japanese residents. Recently, a review paper was published on the effects of natural variation in natural radionuclides on humans, animals and other organisms [6]. Since the available studies are few in number however, it is difficult to conclude what the effect is of chronic exposure by low level radiation. The authors proposed a three-year project (from July 2012 through March 2015) to look at the health risk associated with exposure to low-dose-rate ionizing radiation. The high background radiation (HBR) area of Kerala state was selected from which prediction data could be obtained for the project. Kerala state is located in southern India and is well-documented as a high background radiation area [7]. The annual absorbed dose for the coastal areas of Karunagappally Taluk has been observed to have especially high values due to thorium-containing monazite sand found there and external dose estimation was carried out for this area as a part of an epidemiological study [8]. However, it is necessary to make a dose distribution map for a detailed evaluation of external dose and no such map has been made yet for this area. Thus, a car-borne survey was carried out in Karunagappally Taluk for the purposes of making a dose distribution map and estimating the annual external dose there.

Materials and Methods

Ethics Statement

No specific permissions were required for these locations/activities for the following reasons; a) The field studies did not involve endangered or protected species; b) We provide the specific locations of our study (GPS coordinates).

Survey area

Karunagappally Taluk, in Kerala state, India consists of 12 Panchayats which are administrative sub units [8]. The survey route which covered all the Panchayats of Karunagappally Taluk is shown in Fig 1. This route map was drawn using the Generic Mapping Tools (GMT) created by Wessel and Smith [9]. The 12 Panchayats were classified as low level, mid-level, and high level HBR areas according to the annual absorbed doses. Oachira, Thazhava, Thevalakkara, and Thodiyoor Panchayats were low level HBR areas; Clappana, Karunagappally, K. S. Puram, and Thekkumbhagam Panchayats were mid-level HBR areas; and Alappad, Chavara, Neendakara, and Panmana Panchayats were high level HBR areas.

Fig 1

The route map in Karunagappally Taluk.

The abbreviations of each Panchayat are shown in the map and are as follows. AL: Alappad, CH: Chavara, CL: Clappana, KA: Karunagappally, KS: K. S. Puram, NE: Neendakara, OA: Oachira, PA: Panmana, TD: Thodiyoor, TK: Thekkumbhagam, TV: Thevalakkara, TZ: Thazhava.

Car-borne survey

A car-born survey technique is a convenient method for the evaluation of radiation dose in a wide area in a short period [10]. Thus, a car-borne survey which used a 3-in × 3-in NaI(Tl) scintillation spectrometer (EMF-211, EMF Japan Co., Japan) was carried out in Karunagappally Taluk from September 23 to 27 in 2013. Measurements of the counts inside the car were carried out every 30 s along the route. Simultaneously with the gamma-ray count measurements, the latitude and longitude at each measurement point were measured with a global positioning system. Car speed was kept around 40 km/h. Since count rate was measured inside the car, it was corrected by multiplying with a shielding factor in order to represent the unshielded external dose rate. The shielding factor of the car body was estimated by making measurements inside and outside the car at 34 points and correcting them with the inside count rates.

In this study, a response matrix method was used to convert the gamma-ray pulse height distribution measured with a NaI(Tl) scintillation spectrometer to the energy spectrum of incident gamma-rays [11, 12]. The details of the response matrix method are described as follows. The photon peaks of 40K (E γ = 1.464 MeV) and 208Tl (E γ = 2.615 MeV) were used for calibration from the channel number to gamma-ray energy. The photon peak positions of 40K and 208Tl were determined accurately by smoothing the gamma-ray pulse height distribution in Fourier expansion [13].

Since it is difficult to obtain the photon peak for each gamma-ray in a 30-s measurement in the natural environment, the accuracy of the energy calibration will be low. Furthermore, photon peaks obtained with the NaI(Tl) scintillation detector are affected by temperature, and photon peak positions of gamma-rays from 40K and 208Tl will drift from side to side [13]. Thus, accuracy of the air kerma rate which is estimated using a gamma-ray pulse height distribution obtained by 30-s measurements will be low. Accordingly, the relationship between the total counts of a gamma-ray pulse height distribution and air kerma rates was examined for the estimation of dose rate conversion factor. Measurements of gamma-ray pulse height distribution were carried out 1 m above the ground surface at 35 points. At one measurement site, the gamma-ray pulse height distribution was obtained at two different points. So the number of measurements for the estimation of dose rate conversion factor differed from the number of measurements for the estimation of shielding factor. Counting times at high level areas, mid-level areas and low level areas were set as 300 s, 600 s and 900 s, respectively. The weather condition was sunny throughout the entire measurement period. Thus, the estimated air kerma rates were not affected by rainfall. Air kerma rate due to the gamma-ray energy of natural radiation can be regarded practically as equivalent with the absorbed dose rate in air.

Activity concentration and the contribution of 40K, 238U and 232Th to air kerma rate

Three measurement points for asphalt-covered surfaces were excluded and the gamma-ray pulse height distributions for the remaining 32 points on soil, grass and beach sand surfaces were unfolded using a 22 × 22 response matrix for the estimations of activity concentrations and the contributions of 40K, 238U and 232Th to air kerma rate [12, 13]. By following the method presented in reference [12, 13], a gamma-ray pulse height distribution obtained by measurements was converted to the energy bin spectrum of incident gamma-rays which is a distribution of gamma-ray flux density to each energy bin. In this method, each channel of a multi-channel analyzer was divided into unequal intervals as energy bins in the gamma-ray energy range from 0 to 3.2 MeV. The gamma-ray energies over 3.2 MeV were not included for evaluation, since the maximum value of the gamma-ray energy from natural radionuclides is 2.615 MeV which is emitted from 208Tl. The energy ranges were set up so that the gamma-rays of 1.464 MeV from 40K, 1.765 MeV and 2.205 MeV from 214Bi (238U-series) and 2.615 MeV from 208Tl (232Th-series) could be stored in each bin, respectively. A 22 × 22 matrix for the 3-in × 3-in NaI(Tl) scintillator for an isotropic field was calculated using the Monte Carlo code, SPHERIX [14, 15]. The gamma-ray flux density and dose rate per unit solid angle are considered almost isotropic in the natural environment [16]. The contribution of cosmic-rays to the gamma-ray pulse height distribution was subtracted using the energy bin number 22 which was stored as the gamma-ray energy range from 3.0 to 3.2 MeV. Moreover, the contribution of 40K contamination due to photo-multiplier tube was also subtracted from the gamma-ray pulse height distribution. Clear peaks from 40K (bin number 14, energy range: 1.39–1.54 MeV), 214Bi (bin numbers 16 and 18; energy range, 1.69–1.84 MeV and 2.10–2.31 MeV) and 208Tl (bin number 20; energy range, 2.51–2.72 MeV) were observed in the unfolded spectrum of the gamma-ray pulse height distribution. The activity concentrations of the 238U-series, 232Th-series, and 40K could be estimated by comparing them with the theoretically evaluated gamma-ray flux density spectrum due to these nuclides. In order to evaluate each activity concentration of the natural radionuclides from an energy bin spectrum, it is necessary to calculate the gamma-ray flux densities per unit activity concentrations of the 238U-series, 232Th-series, and 40K. This calculation assumed that a semi-infinite volume source was formed in the ground [13]. The primary and scattered gamma-ray flux density per unit activity concentrations could be calculated using the one-dimensional Monte Carlo gamma transport code, MONARIZA/G2 [17, 18]. A total of a million histories were traced for each natural radionuclide. The nuclear data of gamma-ray energies and disintegration rates used the reported values by Beck [19] and Beck et al. [20] for this Monte Carlo simulation. The activity concentration of each natural radionuclide was evaluated by a successive approximation which used a 3 × 3 matrix determined by Minato [13] to the values of energy bins for 40K, 238U-series, and 232Th-series. Incidentally, the energy resolutions of the photon peak from each natural radionuclide have wide variation for commercially available scintillators. So a diagonal elements fitting (DEF) technique which was developed by Minato [13] was used in this calculation. A new matrix was reconstructed easily and quickly from the standard response matrix using the technique of DEF [13]. The statistical errors for air kerma rate and activity concentrations for 40K, 238U-series and 232Th-series obtained using this software depend on the integral air kerma (nGy) at each measurement point [21], and these were evaluated in this study as 2%, 2%, 6–8% and 4–5%, respectively.

Results and Discussion

Shielding factor and dose rate conversion factor

Fig 2 shows the relationship between the outside and inside count rates, yielding a shielding factor of 1.47. This value was similar to the previously reported value [11]. Fig 3 shows the relationship between air kerma rate (nGy/h) evaluated using the 22 × 22 response matrix method and total count rate (cpm). The conversion factor was evaluated as 0.00244 (nGy/h/cpm). Thus, K out the air kerma rate 1 m above the ground surface at each measurement point can be estimated using Eq (1).

Fig 2

Correlation between count rates outside and inside the car.

This regression formula was used as the shielding factor of the car body.

Fig 3

Correlation between air kerma rate which was calculated by software using the response matrix method and total count rate outside the car.

This regression formula was used as the dose rate conversion factor.

In this study, the counts (N in) inside the car were obtained by the measurements for 30 s. Since the dose rate conversion factor was given as dose rate (nGy/h) for counts per minute (cpm) it is necessary to double N in in order to convert into the counts per minute.

Air kerma rate distribution of Karunagappally Taluk

A dose rate distribution map of Karunagappally Taluk is shown in Fig 4. This map was drawn using 2053 data which included the spot measurements. A heterogeneous distribution of air kerma rate was seen. Air kerma rates of over 0.5 μGy/h were found along the coast. However, values of over 0.3 μGy/h were observed in the interior of Karunagappally Taluk. The maximum air kerma rate, 2.1 μGy/h, was observed near the rare earth mining and separation facility in Chavara Panchayat. Moreover, the air kerma rate of 1.9 μGy/h was observed for the coastal area of Neendakara Panchayat. Both these values were observed on sand which likely contained monazite. According to the reports by Nair et al. [8] and Christa et al. [22], an air kerma rate of about 9 μGy/h was observed for the coastal area of Neendakara Panchayat. Furthermore, Derin, et al. [22] observed the maximum value of 28 μGy/h, with an average value of 9.8 μGy/h for the Chavara and Neendakara coastal areas. Since the air kerma rates on the beach sand could not be obtained by the car-borne survey, the maximum values of the present study might be much lower than the value obtained by Derin, et al. Thus, the present result suggests that the highest dose rate might be observed in the immediate vicinity of the Lakshadweep Sea. The air kerma rates measured by the car-borne survey for each Panchayat are shown as a box plot in Fig 5. Air kerma rate for high radiation level Panchayats had a large variance. The range of the coefficient of variation for these Panchayats was evaluated to be from 39% to 95%, and these values were relatively higher than those at mid-level and low level Panchayats which ranged from 12% to 32% and from 19% to 58%, respectively. It seems that the higher 232Th activity concentration for the coastal areas was a spatially localized distribution.

Fig 4

The distribution map of air kerma rate in Karunagappally Taluk.

This map was also drawn using GMT [9].

Fig 5

The air kerma rates for each Panchayat as measured by the car-borne survey.

Gray shows measurements in high level background radiation areas. White shows measurements for mid-level and low level background radiation areas.

Activity concentration

The spot measurement results are shown in Table 1. 232Th activity concentration on beach sand was higher than that on soil and grass surfaces, and the range of the activity concentration was estimated to be 0.7–2.3 kBq/kg. According to Derin, et al. [23], the highest 232Th activity concentration of 41 kBq/kg in beach sand was obtained from coastal areas in Chavara and Neendakara Panchayats and the geometric mean in these areas was estimated to be 14 kBq/kg. In the present study, the contribution of 232Th to air kerma rate ranged from 52% to 79%. Especially, the contribution of 232Th to air kerma rate was over 70% at measured points with values over 0.34 μGy/h.

Table 1

The measured activity concentrations and the contributions of 40K, 238U and 232Th to air kerma rates in Karunagappally Taluk.

Panchayat	Location		Air kerma rate (nGy/h)				Contribution to air kerma rate (%)			Activity concentration (Bq/kg)			Surface
	Latitude (°)	Longitude (°)	40K	238U	232Th	Total	40K	238U	232Th	40K	238U	232Th	Classification (soil, grass, sand)
Alappad	9.0369	76.5087	4	159	471	634	1	25	74	99 ± 2	359 ± 22	716 ± 29	Soil
	9.0524	76.5068	10	25	69	104	10	24	66	240 ± 5	57 ± 3	105 ± 4	Soil
	9.1352	76.4639	28	16	47	91	31	17	52	672 ± 13	36 ± 2	71 ± 3	Soil
Chavara	8.9696	76.5361	1	55	198	254	1	22	78	32 ± 1	124 ± 7	300 ± 12	Soil
	8.9704	76.5450	2	95	334	431	0	22	77	44 ± 1	215 ± 13	507 ± 20	Soil
	8.9877	76.5616	9	24	58	91	10	26	64	214 ± 4	54 ± 3	88 ± 4	Soil
	8.9848	76.5527	4	41	92	137	3	30	67	92 ± 2	93 ± 6	142 ± 6	Grass
	8.9889	76.5314	2	46	166	214	1	21	78	46 ± 1	103 ± 6	253 ± 10	Soil
	8.9698	76.5294	17	163	461	641	3	25	72	407 ± 8	370 ± 22	709 ± 28	Soil
	8.9928	76.5242	40	558	1543	2141	2	26	72	955 ± 19	1269 ± 76	2374 ± 95	Sand
Neendakara	8.9680	76.5316	6	39	113	158	4	25	71	145 ± 3	88 ± 5	172 ± 7	Soil
	8.9634	76.5317	6	98	220	324	2	30	68	151 ± 3	222 ± 13	334 ± 13	Soil
	8.9557	76.5328	25	391	967	1383	2	28	70	594 ± 12	889 ± 53	1488 ± 60	Sand
	8.9557	76.5328	40	521	1460	2021	2	26	72	964 ± 19	1185 ± 71	2247 ± 90	Sand
	8.9534	76.5418	4	76	277	357	1	21	77	102 ± 2	172 ± 10	420 ± 17	Grass
	8.9420	76.5360	5	125	430	560	1	22	77	113 ± 2	283 ± 17	654 ± 26	Sand
Oachira	9.1352	76.5281	6	28	70	104	6	27	67	139 ± 3	64 ± 4	107 ± 4	Soil
	9.1329	76.5338	3	11	28	42	8	26	66	81 ± 2	25 ± 2	42 ± 2	Soil
	9.1167	76.5293	3	18	52	73	4	25	71	66 ± 1	41 ± 2	79 ± 3	Soil
	9.1170	76.5180	2	22	66	90	2	25	73	37 ± 1	51 ± 3	100 ± 4	Soil
	9.1351	76.5088	2	61	215	278	1	22	77	56 ± 1	138 ± 8	327 ± 13	Soil
Panmana	9.0309	76.5501	10	105	322	437	2	24	74	232 ± 5	238 ± 14	496 ± 20	Soil
	9.0019	76.5429	2	69	269	340	1	20	79	53 ± 1	156 ± 9	408 ± 16	Soil
	9.0041	76.5365	4	28	74	106	3	26	70	84 ± 2	64 ± 4	114 ± 5	Soil
	9.0008	76.5254	9	93	244	346	3	27	70	218 ± 4	212 ± 13	375 ± 15	Soil
	9.0162	76.5340	6	37	117	160	4	23	73	135 ± 3	84 ± 5	177 ± 7	Soil
	9.0250	76.5378	1	54	185	240	0	23	77	22 ± 0	123 ± 7	281 ± 11	Soil
Thevalakkara	9.0177	76.5617	5	20	57	82	6	24	70	115 ± 2	45 ± 3	86 ± 3	Soil
	9.0135	76.5881	5	30	93	128	4	24	73	108 ± 2	68 ± 4	142 ± 6	Soil
	9.0106	76.5934	5	33	94	132	4	25	72	112 ± 2	74 ± 4	143 ± 6	Soil
	9.0061	76.5909	2	27	92	121	1	23	76	36 ± 1	62 ± 4	139 ± 6	Soil
	8.9920	76.5746	4	24	72	100	4	24	72	100 ± 2	54 ± 3	110 ± 4	Soil

External dose estimation

Mean, median, minimum and maximum values of the annual effective dose of each Panchayat are shown in Table 2. It is difficult to make individual dose estimation without considering occupancy factor that depends on the age and sex of the residents [8]. Moreover, it is necessary to consider the shielding factor by the building materials that is used to make the individual dose estimation based on an outdoor dose rate. However, the reported individual dose rate ratios of indoors to outdoors were to be between 0.08 and 1.20 among 197 houses investigated in Kerala state, India [24]. This result suggests that the ratio of 1.20 ± 0.24, which was reported as a typical value for normal background radiation areas in India [25], cannot be used for dose estimation. Individual dose estimation which considers each occupancy factor indoors and outdoors is an issue for future work. Thus, only a conservative external dose estimation was made excluding the occupancy factor and the shielding factor by building materials. The result by this simple approach might be described as an overestimated value. Here, annual effective dose H (mSv/y) was estimated using Eq (2): where, K out is air kerma rate outdoors (nGy/h), DCF is the dose conversion factor (0.7 Sv/Gy) from the air kerma to the external effective dose to adults [5], and T is 8760 h (24 h × 365 d). The arithmetic mean of each Panchayat in Karunagappally Taluk ranged from 0.6 mSv/y (Oachira and Thazhava) to 2.3 mSv/y (Neendakara). The highest value (13 mSv/y) was observed in Chavara Panchayat, and a similar value (12 mSv/y) was observed in Neendakara Panchayat. Both these values were observed in beach sand. The annual effective dose less than or equal to 1 mSv/y in Karunagappally Taluk was observed for 69% of all data.

Table 2

The estimated annual effective doses of each Panchayat in Karunagappally Taluk.

Panchayat	No. of data	Average (mSv/y)	Maximum (mSv/y)	Minimum (mSv/y)	Median (mSv/y)	Radiation level
Alappad	136	1.9	5.1	0.4	1.9	High
Chavara	214	1.2	13	0.5	0.9	High
Clappana	50	0.8	1.1	0.6	0.8	Middle
Karunagappally	112	0.9	1.6	0.5	0.9	Middle
K. S. Puram	74	0.8	1.1	0.5	0.9	Middle
Neendakara	125	2.3	12	0.6	1.5	High
Oachira	115	0.6	1.5	0.2	0.5	Low
Panmana	273	1.1	2.9	0.6	1	High
Thazhava	138	0.6	0.8	0.5	0.6	Low
Thekkumbhagam	112	1.1	4.3	0.4	0.9	Middle
Thevalakkara	215	0.7	1	0.5	0.7	Low
Thodiyoor	109	0.7	0.9	0.5	0.7	Low

Conclusions

The car-borne survey with a NaI(Tl) scintillation spectrometer was carried out in Karunagappally Taluk, India, for the purposes of making a dose distribution map and estimating the external dose. The dose rate distribution map showed a heterogeneous distribution of the air kerma rates. Especially, air kerma rates of over 0.5 μGy/h were distributed along coastal areas. The contribution of 232Th to air kerma rate at measurement points with values exceeding 0.34 μGy/h was over 70%. The highest annual effective dose of 13 mSv/y was found near the rare earth mining and separation facility in Chavara Panchayat. More than 30% of all the annual effective dose measurements obtained in this study exceeded 1 mSv/y.