Advances in Materials Science and Engineering Assessment of Gamma Radiation Hazards Related to Geologic Materials: Comparison of Results by Field Gamma Spectrometry and Laboratory Methods

Earth Sciences Department, School of Sciences, University of Minho, 4710-057 Braga, Portugal Instituto Universitario de Xeoloxı́a “Isidro Parga Pondal”, University of A Coruña, 915001 A Coruña, Spain Earth Sciences Department, LandS/Lab2PT-Landscape, Heritage and Territory Laboratory (FCT UID/AUR/04509/2013, FEDER COMPETE POCI-01-0145-FEDER-007528), School of Sciences, University of Minho, 4710-057 Braga, Portugal


Introduction
Natural ionizing radiation includes external and internal sources.While internal sources comprise naturally occurring radioisotopes that are taken into the human body, external ones are mainly cosmic rays and gamma radiation emissions from both terrains and building materials [1].Isotopes present different activity values (radioactive decay rate, which is expressed in SI units as Becquerel, Bq, with 1 Bq � 1 s −1 ), and the ratio of the activity of the isotope to the total mass can be considered as an expression of its concentration, referred as activity concentration or specific activity (Bq/kg).e external gamma radiation related to terrestrial sources is mainly due to primordial radionuclides with radioactive decay half-lives that are approximately Earth's age or older [2].ey are mainly 40 K and radioisotopes of the radioactive decay series of 238 U and 232 (daughter isotopes of these radioisotopes) present in minerals and rocks used in building materials, being the main cause of indoor radiation exposure, which is about 40% greater than outdoor exposures in middle and high latitudes [3].
Different approaches have been considered to assess the dose that a person will receive from building materials in a given indoor space.e 2013/59/EURATOM directive of the Council of Europe (CEU [4]), which is directed to setting "basic safety standards for protection against the dangers arising from exposure to ionising radiation," indicates the use of an activity concentration index (1) which is based on the sum of factors obtained by dividing the activity concentrations (or specific activities) of isotopes by constants that correspond to certain reference conditions in terms of room characteristics, materials application, and bulk density and exposition time.
where C Ra-226 is the activity concentration of 226 Ra in Bq/kg, C -232 is the activity concentration of 232 in Bq/kg, and C K-40 is the activity concentration of 40 K in Bq/kg.A value of one for this activity concentration index is indicated as a conservative screening tool with the main goal of distinguishing building materials with radioisotope levels that may have a contribution to external gamma radiation in excess of typical outdoor exposure above the reference level.
is contribution of the materials will be different for different conditions of those considered to obtain the index, and the question of the typical outdoor exposure might also be a subject of discussion as it can show marked differences depending on the terrain characteristics [5].
ere is hence an interest (that might be increased in the future as this directive is implemented) in assessing the concentrations of radioisotopes in building materials.However, direct measurements of isotopes activities are costly and time-consuming and can be considered destructive in the sense that they require sampling, which limits its application to building materials applied in existing structures.is limitation of the extraction of a sample also applies to other, eventually faster, laboratory techniques such as inductively coupled mass spectrometry (ICP-MS) and X-ray fluorescence (XRF) that allow the measuring of the elements amounts (from which specific activities of isotopes can be calculated, e.g., for U, , and K by the factors indicated in IAEA [6]).
Portable gamma-ray spectrometry (which is referred here as field gamma spectrometry-FGS) is a classical geophysical technique that has been used for geological exploration as well as for assessing the dose rate in terrains.It can be a fast procedure, and there are variants (such as the one that is considered here) that can be used by a one-person team.When FGS measurements are carried out outdoor on rock outcrops, a 2πr geometry (plain areas) is needed to avoid over or underestimation of the radioisotopes (K, U, and ) and gamma radiation [6].ere is, consequently, interest in studying the usability of FGS at several perspectives related to building materials and their contribution to external gamma radiation: (i) To assess natural variability in quarries (which can be used to avoid portions of the rock mass with undesirable amounts of the radioisotopes) (ii) To control the products which are going out (again to avoid situations above certain values) (iii) To assess the potential hazards levels related to a given material applied in existing structures ere are several studies comparing field and laboratory analyses of radioisotopes, some of them regarding building materials, such as dose measurements in spaces and in relation to the applications of materials [7,8].In the studies of Mueller et al. [9] and Nuccetelli et al. [10], one can find detailed discussions regarding the issues concerning the comparison of laboratory and field measurements of radioisotopes.However, studies with the perspective that is presented here (assessment in relation to the activity concentration index as a screening tool of building materials) are much scarcer.As far as we know, Stals et al. [11], concerning different types of building materials, and Döse [12], concerning concrete slabs with different types of rock aggregates, are the only works oriented to the assessment of activity concentration index of building materials with portable gamma spectrometry.ese two publications report excellent positive linear correlations between portable and laboratory results.Stals et al. [11] found that 226 Ra and (especially) 40 K results by portable equipment were lower than results from laboratory analyses while 232 results were higher in the portable equipment results.e final result in terms of activity concentration index (I) was slightly lower for portable measurements (between <0.1 and 3.65) than laboratory analyses (<0.1 to 3.91), getting activity concentration index values between 0.15 and 2.04 for laboratory measurements and between 0.28 and 1.6 for the portable measurements.e results of laboratory analyses were lower for the lower I values and higher by 15%-20% for the higher I range (above 1 in the laboratory results).ese studies [11,12] were both made with laboratory well-defined geometrical specimens, while our study considers outcrops with variable geometrical conditions (hence harder conditions for correlations), which are nearer the field conditions that one can find in quarry mapping.In a field study with a perspective near to this paper (but not in terms of the activity index as screening tool of building materials), Bochiolo et al. [13] found that field measurements gave higher measurements in all the considered radioisotopes (which means, consequently, higher activity concentration index values), in a study performed in an old mine tunnel (therefore in geometrical conditions that are more homogeneous than those considered here).
Our main goal in this paper was to assess the performance of field measurements by portable gamma spectrometry in rock outcrops (including their weathering products) by comparison with laboratory analyses of samples from the same places where the measurements were made.

Studied Geologic Objects.
e present study considered samples from igneous and metamorphic rocks and their weathering products in outcrops near the towns of A Coruña (NW Spain) and Braga (NW Portugal), which were based on the available geological mapping (IGME [14,15]; Ferreira et al., [16]), and our field studies can be described as gabbro (3 samples in A Coruña), leucocratic body in gabbro (1 sample in A Coruña), biotite-rich granites (4 samples in Braga), leucocratic granite (3 samples in A Coruña and 5 samples in Braga), and Silurian siliciclastic metamorphic rocks (7 samples in Braga).

Analytical Methods.
e present work considers results from different methods to assess the concentration of radionuclides in such rocks and weathering products to assess their activity index and possible limitations in their use as building materials under the present European regulations established by the directive CEU [4].e main aim of the work is to compare the assessments of radioisotopes concentrations in these geological materials, namely, in the perspective of the European directive trough the results of 2 Advances in Materials Science and Engineering the concentration activity index (which can be considered a synthesis parameter of the contribution of the different isotopes in relation to a certain reference state).

Field Gamma Spectrometry (FGS).
Gamma spectra were acquired in situ in the field with a portable gamma spectrometer GF Instruments GRS-2000, equipped with a probe BGO (Bi 4 Ge 3 O 12 ) having a 51 mm × 51 mm detector and a shielded photomultiplier (Figure 1).e equipment allows measuring energies up to 3 MeV.Measurements were made by direct probe contact at outcrops with variable geometry.e recommendations given by the IAEA [6] for static measurements (recording a spectrum in a given point during a certain amount of time) were followed.Each acquired spectrum corresponds to a counting period of 3 minutes.
is time period is recommended by the manufacturer to optimize time and precision.According to the manufacturer's specifications for low gamma-ray emissions, statistical uncertainties are estimated to be 6% for K, 30% for U, and 16% for by measuring for 180 s. is error is due to background radiation that is basically caused by three main sources: atmospheric radon and its daughters, cosmic rays, and instrument background.Moreover, the geometry of measured surface must be planar (2πr geometry) being both possible surface roughness of the measured terrains and any nearby topographic feature possible causes of inaccuracies in the measurements.
From the obtained spectra (an example is presented in Figure 2), we have used some peaks to assess the contents of K (potassium, mass percent), eU (uranium equivalent, in mass parts per million or ppm), and e (thorium equivalent, also in ppm).While potassium contents are estimated from the peak of 40 K (energy range 1.366-1.564keV), there are several radioactive isotopes in the uranium and thorium decay series, and hence, "equivalents" of uranium and thorium (eU and e , resp.) can be used, assuming secular equilibrium.Estimations of eU and e were made from peaks of the spectra corresponding to 214 Bi (energy range 1.57-1.959keV) and 208 Tl (energy range 2.42-2.81),respectively, on the decay chain of uranium and thorium.e specific activities of 226 Ra (or 238 U), 232 , and 40 K were calculated from eU, e , and K with the conversion factors given in IAEA [6].

Sample Preparation for Laboratory Analyses.
For laboratory analyses of the samples, the obtained samples were powdered at a particle size less than 0.5 mm and dried at 100 °C to constant weight.Once dried, the powder was fired at 500 °C for 20 hours to remove organic matter.A part of the powder was separated for elemental analyses and most of them for low background gamma spectrometry.

Laboratory Gamma Spectrometry (LGS).
For laboratory gamma spectrometry (which is referred here as LGS), a mass of 250 g of each sample was placed in Marinelli polystyrene containers and sealed with a sellotape and left for more than 4 weeks in order to allow for Ra and its short-lived progeny to reach radioactive equilibrium.e radioactivity concentrations were measured in coaxial HPGe Canberra GR6022 of the "Servicio de Radioisotopos (CITIUS)" of the University of Seville (Spain).e detector was shielded with 10 cm thick lead on all sides including the top, to reduce the background contribution from the surroundings.e shielded volume is 180 cm 3 ; the relative efficiency of the detector is 60%, and it has a resolution of 1.4 keV at 2.3 keV for energies 122 keV and 1332 keV, respectively.

Element Analyses.
e determination of both major and minor elements in the samples was obtained by X-ray fluorescence spectrometry (XRF).Powdered samples were measured in a Bruker-Nonius S4 Pioneer wavelength dispersive fluorescence spectrometer under helium purge at the University of A Coruña.To assess trace elements, inductively coupled plasma-mass spectrometry (ICP-MS) was used.Advances in Materials Science and Engineering 3 ICP-MS analyses were carried out in a ermo Scientific ™ ELEMENT XR ™ ICP-MS which combines a dual mode secondary electron multiplier with a Faraday detector.Lithium metaborate fusion was used for sample preparation.e powdered samples were mixed with an equal amount of lithium tetraborate flux, placed in a carbon crucible, and fused at 1000 °C in a furnace for 30 minutes.Once the melt was cooled, the resultant fusion bead was ground and dissolved in 100 mL of 4% HNO 3 /2% HCl 3 solution which was then analysed by ICP-MS.e same factors (given in IAEA [6]) were used to convert U, , and K measurements to specific activities of 226 Ra (or 238 U), 232 , and 40 K.

Data
Treatment.Some notes on the statistical treatment are presented here (additional information can be found in statistical texts such as Dodge [17]).e statistical summary includes diverse measures of points of the value distribution: arithmetic mean (referred simply as mean from now on), harmonic mean, median, geometric mean, minimum, first (lower) quartile, third (upper) quartile, and maximum.While the minimum, mean, and maximum are certainly widely known, some notes are presented in relation to the other measures.e geometric mean is obtained by the n root of the product of the values (or by the arithmetic mean of their logarithms).e harmonic mean is obtained by the quotient between the number of elements of the data set and the sum of the inverse of each element of the data set.e quartiles are the points that divide a set of values into 4 equal parts.e middle quartile divides the data set into two equal parts and is equal to the median.
e lower and upper quartiles correspond to the upper limit of the lower 25% and lower limit of the upper 25%.While it is not our goal to characterize the distribution of the sets of values from the studied rocks, dispersion statistics (standard deviation, standard error of the mean, and interquartile range) are presented for the reader to assess the dispersion of the data sets under study, which could be useful for comparison.While we will not discuss the question of the normal distribution of such a geologically heterogeneous set of measurements, one can highlight that, for a normal distribution, the 95% confidence limits of the mean are obtained by the mean more or less the product of the standard deviation by 1.96 (for a higher confidence limit, the error interval is, of course, higher).
As for there are groups of samples from some rock types and since our goal here is to compare techniques and not the geochemical characterization of the geological materials, data treatment is based on robust statistical techniques based on nonparametric approaches (that do not assume any premise concerning the shape of the data distribution or its parameters).e correlation between different measures is assessed by the rank-based Spearman correlation coefficient, which does not require each studied value set that presents a normal distribution, which is not affected by the specific values nor is it affected by arithmetic operations with constants (e.g., dividing by a constant) since it does not depend on the variance of the value sets.One will also present results of the coefficient of determination that assess the fit of the data to a proposed regression equation (it is a measure of the contribution of a given factor to the variance of another under a given statistical model such as a regression straight line).e distribution of the values of the concentration activity index (I) obtained from results of each analytical technique is represented by a boxplot, representing the lower, middle, and upper quartile and raw data.Since the realization of several statistical tests with the same data set can cause an error I type inflation, we limit the use of other statistical tests to the comparison of the calculated values of the concentration activity index.Following the advice of the American Statistical Association (Wasserstein and Lazar [18]), the obtained p values are reported.Comparisons of I values from the different analytical techniques are based on two nonparametric statistical tests for paired statistical samples: the sign test and the Wilcoxon test.e sign test is based on the proportion of the signs of the differences between the paired results.e Wilcoxon test is based on the ranking of the differences between the compared statistical samples.
All statistic tests, the boxplot and the mean and standard error plots were prepared with Statistica 11 (Statsoft).

Results and Discussion
Figure 3 presents a combined boxplot (quartiles) and raw data plot of the logarithms of ratios, in relation to estimations from field measurements, of laboratory estimations of 226 Ra, 232 , and 40 K specific activities by laboratory gamma spectrometry and those calculated from uranium, thorium, and potassium quantifications (the first two by MS-ICP and the later by XRF).It can be seen that, in almost all the cases, the ratios of laboratory estimations to field estimations are below one (negative logarithms), indicating that, for the different radioisotopes, estimations from field measurements are higher than laboratory ones.
Bivariate plots of laboratory measurements against the field ones as well as the values of the activity concentration index calculated from the different sets of isotopes estimations (laboratory and field) are plotted in Figure 4.It can be seen that there are clear trends toward correlations between the different sets of measurements, but there is some wider scatter in the case of 40 K estimations as shown by the coefficient of determination values, which are lower than those obtained by Stals et al. [11].e plots also include the line of equality between plotted variables, showing that, in general (with few exceptions), values of radioisotopes estimations from laboratory measurements are less than those from field measurements and that this happens in the comparison of calculated values of the activity concentration index also (in this case, only one value obtained from field estimations is slightly lower).
e p values for the Spearman coefficients calculated for I value estimations from the three sets of analyses are presented in Table 1, showing that they markedly below 0.05 (the highest p value equals to 0.000001).It can also be seen from the Spearman coefficients in Table 1 that the correlation between the estimations from the laboratory analyses is better than the correlations of the laboratory analyses with the field measurements.
Table 2 presents the univariate statistics of the activity concentration index (I) values calculated from the three procedures.It can be seen that all location measures are higher in the field set than in the ICP/XRF set that, in turn, are higher than those from the laboratory gamma spectrometry measurements.It can also be seen that I values based on laboratory analyses are always nearer each other than in relation to I values based on field measurements.One can also highlight the relation of the statistics presented in Table 2 in relation to the screening criteria indicated in the European directive (I � 1).Most of the average estimations (mean, geometric mean, and median) are above 1.0 for the set of I values from field measurements but does not happen in any of the averages for I values from laboratory measurements.
e lower quartile of I values based on field measurements is very near this screening value.ese relations are also illustrated in the plot of Figure 5(a) that combines a boxplot and the plotting of raw data. is plot also shows that the number of results that exceeds the I � 1 condition is clearly much higher for the field estimations.
e nonparametric statistical tests in the comparison of these paired values (Table 3) show very low p values (the higher one being 0.0005) and being noticeable that the comparison of field values with laboratory values produces lower p values.However, the statistical assessment of higher mean in the field set deserves more caution since for p � 0.05, the intervals of mean more or less 1.96 times the standard error of the mean of the different data sets overlap, as can be seen in the plot of Figure 5(b).Statistical nonparametric tests of paired samples also support these relations, with low p values, namely, in the comparison of I values from field and laboratory measurements (with the highest being p � 0.00003).
To explain the scatter observed among 40 K estimates obtained by both XRF and LGS, two causes can be considered.On the one hand, the 40 K and 214 Bi peaks (i.e., used to assess the U concentration) are very close in the gammaray spectrum.e FGS equipment has low resolution when compared with HPGe used in LGS.Hence, when one of the 214 Bi or 40 K peaks is high, it can overlap to part of the nearby peak providing an overestimation of the other radioisotope [6].To check if this effect is a cause of scatter in our results, we have compared eU and e estimates from FGS and the ratio of K assessed from FGS to K measured by XRF and LGS (Figure 5).e ratios K (FGS)/K (XRF) and K (FGS)/K (LGS) are near 1 when the estimates of K fit, being higher when the K (FGS) overestimates the K (XRF) or K (LGS) estimates.If the overlap of a high 214 Bi peak with the 40 K peak is the cause of K (FGS) overestimation, the ratios should be higher and above 1 when higher eU is measured.However, this is not observed when the K/K ratios are plotted versus eU (Figure 6).In fact, when the ratios are plotted versus the 208 Tl peak estimates (e ), the results are similar to those in the former plot, being negligible any peak overlaps on the K and eU estimates by FGS.Anyway, it must not be forgotten that the large uncertainty of the FGS (6%, 30%, and 16% for K, U, and , resp.) hinders the assessment of this cause of scatter, and a possible larger data set could be required to get concluding results.Anyway, the disequilibrium in the 238 U decay chain should be expected in weathered rocks but not in nonweathered of slightly weathered rocks.In such rocks, secular equilibrium is observed [19].Moreover, the weathering processes and derived disequilibrium are not necessarily similar in different rocks  Advances in Materials Science and Engineering (e.g., gabbro and granite), being necessary to consider any rock type to assess this factor as a possible cause of the scatter results.Other sources of inaccuracy must be taken into account in the estimates, as it is not possible to assess its effect on the results.A flat geometry is required for FGS measurements (2πr).Although planar surfaces have been chosen for the measurements, the work areas are not completely flat, being surface roughness and topographic features a possible cause of such inaccuracies and scattering.

Conclusions
Field portable gamma spectrometry measurements in igneous and metamorphic rocks and their weathering products showed good ordinal correlations with two laboratory techniques of assessment of radioisotopes specific activity, suggesting that it can be used to rank materials in terms of radiological hazard concerning the contribution of materials to the external gamma radiation dose.e results tend to be  LGS: laboratory gamma spectrometry. 6 Advances in Materials Science and Engineering   Advances in Materials Science and Engineering higher in the case of field measurements, which is reflected in a higher set of results above the guideline criteria I � 1.
While one should read carefully in this matter, namely, given the limited sample set considered, these results strongly suggest that portable gamma spectrometry could be a conservative tool in assessing external gamma radiation posed by materials and, of this is confirmed, be used as a screening tool to separated samples to be sent for laboratory analyses.It is possible that these outcrop measurements could work in a conservative way given the great amount of the geological substance that is present in the measuring point.But it is admitted that portable measurements might not be as efficient in the case of measuring specimens of materials as the effect of the surroundings could behave a higher impact (due to the smaller amount of the measured object).e Lab2PT-Landscapes, Heritage and Territory laboratory-AUR/04509 is supported by the Portuguese "Fundação para a Ciência e a Tecnologia" (Portuguese funds and where applicable the FEDER co-financing, in the aim of the new partnership agreement PT2020 and COMPETE2020-POCI 01 0145 FEDER 007528).e University Institute of Geology of the University of A Coruña (Spain) received support from Xunta de Galicia with funds from "Consolidación y estructuración de unidades de investigación competitivas-Grupo de potencial de crecimiento" (GPC2015/024).e authors would like to thank A. Pinto, F. Dias, M. Cunha, and P. Pereira for help during the field work.

Figure 1 :Figure 2 :
Figure 1: Picture of the portable gamma-ray spectrometer GF Instruments model GRS-2000 used for field measurements in this work.On the left side, the probe with the detector can be seen, and on the right side, the computer system with electric supply can be seen.

Figure 3 :
Figure 3: Combined boxplots (quartiles) and raw data plots of logarithms of radioisotopes estimations from laboratory analyses to those from field measurements: (a) comparison of laboratory gamma spectrometry (LGS) to field gamma spectrometry (FGS); (b) comparison of estimations from ICP and XRF to FGS.

Figure 4 :
Figure 4: Plots of the results from laboratory analyses (ICP, XRF, and LGS: laboratory gamma spectrometry) against field estimations: (a) 226 Ra activity concentration; (b) 232 activity concentration; (c) 40 K activity concentration; (d) activity concentration index calculated from the isotopes estimations.

Figure 5 :
Figure 5: Plots for comparison of the activity concentration index values (I) calculated from the laboratory analyses (ICP, XRF, and LGS: laboratory gamma spectrometry) and field estimations: (a) boxplot (quartiles) and raw data plot; (b) mean and standard error plot for 95% confidence; (c) mean and standard error plot for 99.5% confidence.

Figure
Figure Plots for the comparison of estimates of eU (a) and e (b) versus K (%) ratios obtained by FGS, XRF, and LGS. e K ratios correspond to K (FGS)/K (XRF) and K (FGS)/K (LGS).e insets show a detail of the whole plots in the range 0-3 of the ratio and 0 to 35 ppm for eU and e .No correlation is observed between the data.

Table 1 :
Spearman correlation coefficients between the calculated values of the concentration activity index (I) from different procedures for estimation of isotopes specific activity and related p values.

Table 2 :
Statistic measures of the activity concentration index (I) sets from different procedures for estimation of isotopes specific activity and related p values.

Table 3 :
Results (p values) of statistical tests for paired comparisons of estimations of the activity concentration index (I) values.