Evaluation of Natural Radioactivity in Marine Sand Deposits from Offshore China

Natural radioactivity is very important for the assessment of the marine sand property and usability. By using gamma spectrometry, the concentration of the natural radionuclides Ra, Th and K have been measured in marine sand deposits from Liaodong Bay (LDB), North Yellow Sea (NYS), Zhoushan area (ZS), Taiwan Shoal (TS) and Pearl River Mouth (PR), offshore China, which are potential marine sand mining areas. The radiation activity equivalent (Raeq), indoor gamma absorbed dose rate (DR), annual effective dose (HR), alpha index (Ia), gamma index (Ig), external radiation hazard index (Hex), internal radiation hazard index (Hin), representative level index (RLI), excess lifetime cancer risk (ELCR) and annual gonadal dose equivalent (AGDE) associated with the natural radionuclides are calculated to assess the radiation hazard of the natural radioactivity in the marine sands offshore China. From the analysis, it is found that these marine sands are safe for the constructions. The Pearson correlation coefficient reveals that the Ra distribution in the marine sands offshore China is controlled by the variation of the K concentration. Principal component analysis (PCA) yields a two-component representation of the entire data from the marine sands, wherein 98.22% of the total variance is explained. Our results provide good baseline data to expand the database of radioactivity of building materials in China and all over the world.


Introduction
With the exhausting on-land sand resources and increasing environmental pres-J.Li et al.
sure, it has become necessary to look for alternative sources.Marine sand resources now stand for a more promising alternative, and contribute significantly to the overall provision of sand material in many countries (e.g., [1] [2] [3]).It was widely utilized as the main construction materials for buildings, road, artificial islands, coastal reclamation and beach nourishment, etc.Today, marine sand mining has become the second most important marine mining activity after offshore oil extraction [1] [3].The annual global production of the aggregate is about 16.5 billion tons, of which approximately 10% is supplied by marine sand mining in coastal waters [4].
Together with the rapid economic development of China, the unprecedented demand for the marine sand mining has increased greatly in the recent years.
After nearly 10 years of extensive exploration of marine sand offshore China, there are five areas were assessed which are suitable for marine sand mining, namely Liaodong Bay (LB), North Yellow Sea (NYS), Zhoushan Area (ZS), Taiwan Shoal (TS) and Pearl River Mouth (PR) (Figure 1), where very large quantity of sandy sediments developed.
Figure 1.Locations of samples (solid yellow dots) and study areas offshore China in this study.A, Liaodong Bay (LDB); B, North Yellow Sea (NYS); C, Zhoushan Area (ZS); D, Taiwan Shoal (TS); E, Pearl River Mouth (PR).The geographic coordinate of the samples are shown in Table 1.

J. Li et al.
However, the quality of marine sand should be scientifically evaluated before the mining and utilization.The quality parameters include the grain-size composition and sorting, heavy metal constituents, the mineral components, and most importantly, the natural radioactive properties ( 226 Ra, 232 Th and 40 K).The study of the concentrations of radionuclides and their distribution in sands enables the assessment of radiological risk due to external human exposure to gamma radiation outdoors and inhalation of airborne radioactivity emanating from building constructions and dwellings (e.g., [5] [6]).Generally, the specific activities of 226 Ra, 232 Th and 40 K in raw building materials and their products depend on their geological and geographical conditions as well as the geochemical features of those materials [7].And the natural radioactivity of marine sand depends on the sediment formation and transport processes that were involved; chemical and biochemical interactions influence the distribution patterns of uranium, thorium and their decay products.However, many studies had been carried out on the radionuclide concentrations in sand beaches around the world, such as India [8] [9], Brazil [10] [11], Thailand [12], Egypt [13], Iran [14], and China (Xiamen, [15]) using the gamma ray spectrometry.In spite of the high number of works carried out around the world on the beach sands, there is a lack of studies about radionuclides of offshore sands.
The aim of this study is to determine natural radioactivity ( 226 Ra, 232 Th, 40 K)   levels in sandy sediments (potential marine sand resources) collected from offshore China.Also, the average radium equivalent activity (Raeq), the total absorbed dose rate (D), the indoor and external hazard index (Hin, Hex), the annual gonadal dose equivalent (AGDE) and the annual effective dose equivalent (AEDE), etc., which will be defined later have been calculated and compared with the results of beach sands in literature all over the world.The results of this study will provide background data on the natural radioactive isotopes and environmental pollution of marine sand deposits offshore China.

Materials and Methods
Totally 141 sandy sediments were collected in five potential marine sand deposits offshore China, in which 8 samples are from Liaodong Bay, 68 samples from North Yellow Sea, 12 samples from Zhoushan Area, 41 samples from Taiwan Shoal and 12 samples from Pearl River Mouth, respectively.The locations of each sample were shown in Figure 1 and Table 1.The water depths of the samples range from 10 m to 50 m.After all samples were dried at room temperature, samples were pulverized by Retsch mill and sieved through a 100 mesh to be homogenized, then weighed and transferred to Marinelli beakers of 1000 ml volume.Each sample was sealed for 30 days to reach radioactive equilibrium where the decay rate of the daughters becomes equal to that of the parent [16].Sample preparation and all radioactivity measurements were carried out by using a Gamma-ray spectrometer (BE3830, Canbarra Industries, Inc.) in Qingdao Institute of Marine Geology.
Gamma-ray spectrometer was used to determine the activities of 226

Specific Radioactivity
The activity concentrations of the detected radionuclide 226 Ra, 232 Th and 40 K in the five offshore sand deposits are presented in Table 1.
In the Liaodong Bay, the concentrations of 226    2).The 226 Ra concentrations from the five offshore sandy deposits in this study are all lower than that of world average value.However, the 40 K concentrations in Liaodong Bay, North Yellow Sea, Zhoushan Area and Taiwan Shoal are all much higher than that of the world mean value, and the 232 Th concentrations from Taiwan Shoal and Pearl River Mouth are lower than that of the world average (Figure 2).

Radium Equivalent Activity (Raeq)
Since the distribution of natural radionuclides in the samples is not uniform, a common radiological index has been introduced to evaluate the actual activity level of 226 Ra, 232 Th and 40 K in the samples and the radiation hazards where  The world average values were also shown for comparisons [17], shown as the dashed line for 226 Ra (green), 232 Th (red) and 40 K (blue), respectively.associated with these radionuclides, the radium equivalent activity (Raeq), which can be calculated from the relation [22] [31]: CRa, CTh, and CK are the specific activities of 226 Ra, 232 Th and 40 K, respectively, in units of Bq•kg −1 .In the definition of the radium equivalent, it is assumed that 10 Bq/kg of 226 Ra, 7 Bq/kg of 232 Th and 130 Bq/kg of 40 K each produce an equal gamma-ray dose rate [32] [33].
The calculated values of Raeq for the five sand deposits in investigation are shown in Figure 3.The calculated values of Raeq range from 52.0 (Pearl River Mouth) to 139.8 (North Yellow Sea), with a trend that the Raeq is much higher in North Part of the Chinese Seas than that in the South.All values of Raeq in the studied samples are found to be lower than the criterion limit of 370 Bq/kg [34], and therefore, do not pose any radiological hazard when used for construction of buildings.

Representative Level Index (RLI)
In order to estimate the level of gamma radioactivity associated with different concentrations of certain specific radionuclides, known as the representative lev- (Pearl River Mouth) to 1.08 (North Yeloow Sea) (Figure 4).It is clear that these values do not exceed the upper limit for RLI, which is unity [5].
J. Li et al.

Absorbed Gamma Dose Rate (DR), and Annual Effective Dose Rate (HR)
The absorbed dose rates in indoor air (DR) and the corresponding annual effective doses (HR) attributed to gamma-ray emission from the radionuclides ( 226 Ra, 232 Th and 40 K) in building materials were evaluated using data and formula provided by UNSCEAR (2000) [17] and the EC (1999) [36].In the UNSCEAR and EC reports, the dose conversion coefficients were calculated for the center of a (concrete), respectively.The resulting dose coefficients were found to be 0.92 To estimate the annual effective dose rates, it is necessary to use the conversion coefficient from the absorbed dose in air to the effective dose (0.7 Sv/Gy) and the outdoor occupancy factor (0.2) proposed by UNSCEAR (2000).Therefore, the effective dose rate is determined as follows: Outdoor where DR is given by Equation ( 4).
The estimated results of DR and HR for all the studied marine sands range from 31.4 nGy•h −1 (Taiwan Shoal) to 327.1 nGy•h −1 (North Yellow Sea) and from 0.04 mSv•y −1 (Taiwan Shoal) to 0.20 mSv•y −1 (North Yellow Sea), respectively (Figure 5).And the estimated mean value of DR in all of studied samples is 107.9 nGy•h −1 , which is little bit higher than world average indoor absorbed gamma dose rate of 84 nGy•h −1 [17].Additionally, the estimated mean value of the annual effective dose rate of 0.13 mSv•y −1 is also higher than the world average value (0.07 mSv•y −1 , [17]).

Alpha Index (Ia)
The alpha index was developed as an assessment of the excess alpha radiation exposure caused by inhalation originating from building materials.The alpha index (Ia) is determined by the following formula [37]: where CRa is the 226 Ra activity concentration (Bq•kg −1 ) in the building materials.
The recommended exemption level and recommended upper level for the 226 Ra activity concentration in building materials as are 100 Bq•kg −1 and 200 Bq•kg −1 , respectively, as suggested by the Radiation Protection Authorities in Denmark, Finland, Iceland, Norway and Sweden and the upper level is in agreement with the action level given by the ICRP in Publication 65 (1994) and by the European Commission (1990) [39].It was observed that the values of the alpha index in the studied marine sand samples are below the recommended unity (Figure 6(a)).

Gamma Index (Ig)
Another radiation hazard index, the gamma activity concentration index, Ig, has been defined by the European Commission (1990) [38] and Righi and Bruzzi (2006) [37], which is given as: The index Ig is corrected with the annual dose rate attributed to excess external gamma radiation caused by superficial material.Values of Ig ≤ 2 correspond to the dose rate criterion of 0.3 mSV•y −1 , whereas 2 < Ig ≤ 6 correspond to a criterion of 1 mSv•y −1 ( [36]).Therefore, the activity concentration index should be used only as a screening tool for identifying materials that might be of concern when used as construction materials; although material with Ig > 6 should be avoided, in that these values correspond to dose rates higher than 1 mSv•y −1 [36], which is the highest dose rate value recommended for the population [17].
The gamma index Ig for the marine sands varies between 0.13 (Taiwan Shoal) and 0.62 (Taiwan Shoal) with an average of 0.45 (Figure 6(b)).Therefore, these marine sands can be exempted from all restrictions concerning radioactivity.J. Li et al.

Radiation Hazard Indices
The external radiation hazard (Hex) and the internal radiation hazard, (Hin) was developed for the additional criteria to assess the radiological suitability of a building material [22].And they are defined as follows: and Hin = CRa/185 + CTh/259 + CK/4810 (9) where CRa, CTh and CK are the activities of 226 Ra, 232 Th and 40 K, respectively, in units of Bq•kg −1 .
The determined values of Hex vary from 0.09 (Taiwan Shoal) to 1.05 (North Yellow Sea) with an average of 0.32 (Figure 7).The Hin values range between 0.11 (Taiwan Shoal) and 1.23 (North Yellow Sea) with an average of 0.37 (Figure 7).However, the highest values of Hex and Hin indices are both found in the same sample from North Yellow Sea, with highest value of 226 Ra activity in all samples, which make them larger than 1 (the criterion, [31]).

Excess Lifetime Cancer Risk (ELCR)
The excess lifetime cancer risk (ELCR) was determined using the following equation [39]: where HR, DL and RF are the annual effective dose equivalent, duration of life (70 years) and risk factor (0.05 Sv −1 ), respectively.The risk factor is defined as the fatal cancer risk per Sievert.For stochastic effects, the ICRP 60 uses a value of 0.05 for the public [39].
The calculated range of ELCR is from 0.14 × 10 −3 (Taiwan Shoal) to 1.41 × 10 −3 (North Yellow Sea).The average ELCR values for the five marine sand Mouth), respectively (Figure 8).And it is very clear that most of the marine sands offshore China is slightly higher than the world average (0.29 × 10 −3 ) [17].

Annual Gonadal Dose Equivalent (AGDE)
The activity of bone marrow and bone surface cells are considered to be origins   (North Yellow Sea) and the average value was found to be 409.0mSv•y −1 .The average values do not generally exceed the permissible recommended limits, indicating that the hazardous effects of the radiation are negligible.In the literature, the average AGDE value for the Eastern Desert of Egypt was found to be 2398 mSv•y −1 [41], for the Tamilnadu of 350.63 mSv•y −1 [42], for the Fırtına Valley (Turkey) of 550.5 mSv•y −1 [43].

Multivariate Statistical Analysis
Correlation analysis was carried out in terms of bivariate statistics to determine the mutual relations and strengths of association between pairs of variables through the calculation of the linear Pearson correlation coefficients.The results for Pearson correlation coefficients between all the studied radioactive variables of the marine sand deposits offshore China are shown in Table 3.
A high positive correlation coefficient is observed between 232 Th and 226 Ra (Figure 10(a)), because the radium and thorium decay series occur together in nature [42].In contrast, a very weak negative correlation coefficient was observed between these two nuclides and 40 K (Figure 10(b) and Figure 10(c)), because 40 K is from different origin [42].In addition, all radioactive variables have strong positive correlation coefficients with 226 Ra and 232 Th, while they are weakly negatively correlated with 40 K.All radioactive variables calculated are positively correlated with one another (Table 3).
Principal component analysis was performed on the whole data set (13 variables) to assess the relations between them.The rotated factor analysis was carried out via varimax rotation with Kaiser normalization.The rotated factor 1 and factor 2 values are shown in Table 4. Two principal components were yielded with eigenvalues > 1, explaining 98.26% of the total variance.From the

Conclusion
The natural radionuclide content, radium equivalent activity (Raeq), indoor gamma absorbed dose rate (DR), annual effective dose (HR), alpha index (Ia), gamma index (Ig), external radiation hazard level index (RLI), excess lifetime cancer risk (ELCR) and annual gonadal dose equivalent (AGDE) of five marine sand deposits offshore China were determined.The values obtained in the present study are mostly within the recommended safety limits in spite of just one sample from North Yellow Sea, demonstrating that these marine sands will not pose any significant radiation hazard; thus, the use of these marine sands even in the construction of buildings can be considered safe for the human being.From the statistical analysis, the marine sands in the northern part of offshore China (Liaodong Bay, North Yellow Sea) have 226 Ra, 232 Th and 40 K concentrations that are higher than those from southern part of offshore China (Zhoushan, Taiwan Shoal and Pearl River Mouth).In addition, 40 K and 232 Th are primarily responsible for radioactivity levels of marine sands offshore China.
J. Li et al.

Figure 6 .
Figure 6.Mean (bar) and Maximum (line/data) values of alpha index (Ia) and gamma index (Ig) for different marine sand deposits.

Figure 7 .
Figure 7. Mean (bar) and maximum (line) values of Hex and Hin of different marine sand deposits.

rotation space of component 1 and the component 2 (
Figure 11), the first component accounts for 88.22% of the total variance and is mainly characterized by high positive loadings of concentration of40 K, 232 Th and most of the radioactive variables.While the second component accounts for 10.41% of the total variance and is mainly corresponds to positive loading of 226 Ra and Ia.From the overall component analysis, it can be deduced that40 K and 232 Th dominantly increase the radioactivity in the entire marine sand deposits offshore China.

Figure 11 .
Figure 11.Graphical representations of component 1 and Component 2 (a) and factor loadings for each sample from different marine sands (b).

Table 1 .
Ra, 232 Th and J. Li et al.Location and Activity concentration (Bq•kg −1 ) of 226 Ra, 232 Th, and 40 K in marine sands offshore China.The gamma ray transitions of energies 186.3 keV and 1460 keV gamma-ray transition were used to determine the concentration of 226 Ra and 40 K, respectively, while the gamma-ray lines at 911.0 keV ( 228 Ac) and 583.3 keV ( 208 Tl) were used to determine the concentration of the 232 Th series.The activity levels of the samples obtained for 226 Ra, 232 Th and 40 K are expressed in Bq•kg −1 .
J. Li et al.J. Li et al. 40 K.

Table 2 .
Comparison of activity concentrations of 226 Ra, 232 Th, and 40 K dose rates of sandy sediments found in present study and different areas around the world.
J. Li et al. nGy/h per Bq/kg for 226 Ra, 1.1 nGy/h per Bq/kg for 232 Th, and 0.080 nGy/h per CTh and CK are the activity concentrations of 226 Ra, 232 Th and 40 K, respectively.

Table 3 .
Pearson correlation matrix among the variables.

Table 4 .
Rotated factor loading of the variables.