Home Page Departments Personnel Reports & Data Website Index
Surface Soil Radionuclides and Inorganic Chemicals

Introduction

Results reported herein are from soil samples collected during 1998 from a grid of 16 locations surrounding the WIPP site (the Near Field grid) and a grid of 16 locations approximately 12 miles southeast of the WIPP (the Cactus Flats grid, Fig. 2). Radioanalyses of a subset of these samples for 234U, 235U, 238U, 230Th, 232Th, 228Th and 239,240Pu were reported in the CEMRC 1998 Report using measurements made by Duke Engineering and Services (DES) (Bolton, Massachusetts). Measurements presented herein were made by CEMRC for those same radio-analytes plus 137Cs, 208Tl, 212Pb, 212Bi, 214Pb, 228Ac, 234mPa, 241Am, 40K, 60Co and 7Be. The natural radionuclides 208Tl, 212Bi, 214Pb and 212Pb are measured after allowing for ingrowth and can be used to estimate the concentration of Ra. However, these measurements do not reflect natural levels of those radionuclides in the environment. Results are also presented for 36 non-radiological analytes measured using ICP-MS, AAS and IC.

One finding presented in the CEMRC 1998 Report was that there were significant differences in many analyte concentrations between the Near Field and Cactus Flats grids. Differences in soil texture were postulated as a possible cause for these differences. Aliquots of 73 of the 96 samples collected in 1998 were subsequently analyzed for soil texture in order to test this hypothesis.

Methods

The 16 sampling locations comprising each grid are distributed over approximately 16,580 hectares. At each of the 32 locations, soil was collected at three randomly selected sites within a 50-m radius of the selected reference point. Individual sampling sites were selected on the basis of: relatively flat topography, minimum surface erosion and minimum surface disturbance by human or livestock activity. At each sampling site, approximately 12 L of soil were collected from within two 50-cm x 50-cm areas, to a depth of approximately 2 cm. Soil samples were excavated using a trowel and placed in plastic bags for transport and storage. Sampling equipment was cleaned between samples.

Initial preparation of the samples consisted of passing the soil through a 2-mm sieve to remove rocks, roots and other materials. One half of the 96 samples were originally analyzed for inorganic analytes in 1998. These soil samples were homogenized in the laboratory using a riffler. The method of homogenization was shown previously to yield subsamples that differed from the overall mean count of a radioactive tracer (137Cs) by no more than 7%.

Each sample was homogenized again before aliquots were split for inorganic/metal analyses. New aliquots obtained from samples analyzed and reported in the 1998 CEMRC Report had already been dried at 105° C. Therefore, no analysis for Hg was performed on this set of 48 samples. A 50-g aliquot was removed for radiochemical analyses and ~300 mL aliquots were used for gamma spectroscopy analysis. The samples for gamma analysis were sealed in a ~300 mL paint can and stored for at least 21 days to allow radon progeny to reach equilibrium with parent radionuclides.

Soil samples were analyzed by AAS for As and Se. ICP-MS was used to analyze samples for Ag, Al, As, Ba, Be, Ca, Cd, Co, Cr, Cu, Eu, Fe, Hg, K, La, Li, Mg, Mn, Mo, Na, Ni, Pb, Sb, Sr, Th, Ti, Tl, U, V and Zn. The lower detection limits for both of these systems are in the low parts per billion range (Appendix K). Soil samples were analyzed by IC for chloride, nitrate, nitrite, phosphate and sulfate. A summary of CEMRC QA/QC for inorganic analyses is presented in Appendix K. The mean concentrations of these analytes reported herein for soils include only those values that are above detection levels. Thus, some estimates of the mean may be biased toward larger values.

Gamma spectroscopy analysis was conducted using high purity Ge (HPGe) detector systems for 2-3 d. The systems have HPGe, p-type, coaxial detectors of ~50% nominal efficiency. A set of soil matrix standards was prepared using NIST traceable solutions and used to establish matrix-specific calibration and counting efficiencies. For analyses of alpha-emitting radionuclides, 10-g aliquots of each sample were heated in a muffle furnace to combust organic material and spiked with a radioactive tracer to allow determination of the efficiency of extraction.

The samples for radiochemical analyses underwent total dissolution followed by NaOH fusion of the insoluble residues. Multiple precipitation, co-precipitation and ion-exchange and/or extraction chromatography procedures were then used to separate and purify the desired elements. The elements of interest were then precipitated with La, deposited onto filters, mounted and counted on an alpha spectroscopy system. A summary of CEMRC QA/QC for radioanalyses is presented in Appendix L.

Aliquots of 73 of the 96 samples collected in 1998 were analyzed by the NMSU SWAT laboratory for soil texture. Analyses were made using the hydrometer method (Soil Conservation Service, 1972, Soil Survey Laboratory Methods and Procedures for Collecting Soil Samples. SCS, USDA; Gee, G. W. and J. W. Bauder, 1986, Particle-size Analysis. In Kline, A. (ed.) Methods of Soil Analysis. Part I. Physical and Mineralogical Methods - Agronomy Monograph no. 9. American Society of Agronomy, Madison, WI) resulting in measurements of the percentages of sand, silt and clay.

Multivariate analysis of variance (MANOVA) was used to test for differences between the Cactus Flats and Near Field grids across analytes. Mean concentrations of all analytes were estimated by grid and by soil type and significant differences between means were identified using t-tests. Correlations of the concentrations of radionuclides to soil texture classes, and to concentrations of non-radioactive analytes, were computed using Pearson correlation coefficients. Student's t-tests were applied to pairwise differences between radionuclide concentrations reported by Duke Engineering (1998 CEMRC Report) and the concentrations determined by CEMRC.

Results and Discussion

Baseline Analyses by CEMRC

Results of the MANOVA showed that there were significant (p<0.05) differences between the two grids, and that the Cactus Flats grid generally had higher concentrations of metals than found on the Near Field grid. Of the inorganic analytes, 20 (Ag, Al, As, Ba, Be, Cd, Co, Cu, Eu, Fe, K, La, Mn, Ni, Pb, Th, Ti, U, V and Zn) showed significantly higher concentrations on the Cactus Flats grid as compared to the Near Field grid. Chloride showed a significantly lower concentration on the Cactus Flats grid.

Many metals, including radionuclides, are known to have an affinity for attaching to small particles in the soil (Muller, R. N. and D. G. Sprugel, 1977, Health Physics 33, 405; Muller, R. N. and G. T. Tisue, 1977, Soil Science 124, 191; Watters, R. L., et al., 1983, Radiochema Acta 32, 89; Little, C. A., 1980, Journal of Environmental Quality 9, 350; Tamura, T., 1975, J.Environ. Qual. 4, 350). Clay minerals are aluminosilicates and hydrated oxides, and usually account for the major adsorptive component of soils (Wild, A., 1994, Soils and the Environment. Cambridge University Press; Whicker, F. W. and V. Schultz, 1982, Radioecology: Nuclear Energy and the Environment. Vol. II. CRC Press). Therefore, the concentration of aluminum in soil can be used as a surrogate for estimating the clay content of soils. A correlation between the concentration of Al and the percent clay in soils from 73 of the samples was highly significant (r=0.51, p<0.001), although the correlation between Al and silt was stronger (r=0.69, p<0.001). Similar correlations can be seen for many of the analytes (Table 12). Thus, comparing the ratios of the metals to Al helps to normalize for the effect of soil texture on the concentrations. The ratios of the metals to Al were similar between the two grids (Table 13). Only the ratios of Mo and Tl to Al were significantly (p<0.01) lower on the Cactus Flats grid.

The analyte:Al ratios in soils were compared to ratios reported herein for aerosols (Particulate Concentrations and Inorganics in Near Surface Air). Because aerosols are generally smaller than soil particles, the aerosols tend to have higher mass concentrations of metals which adsorb to fine particles than do soils. The ratios of Ag, Ca, Cd, Co, Cr, Cu, La, Mn, Mo, Ni and V to Al were lower by a factor of 2 or more in the soils as compared to the aerosols at both Cactus Flats and Near Field. The ratios of Fe, Hg, K, Mg, Na and Th to Al were higher in the soils than in the aerosols by a factor of 2 or greater.

The average MDC for radionuclides determined by alpha spectroscopy was »0.04 mBq g-1. The average MDC for radionuclides determined by gamma spectrometry was » 0.5 mBq g-1. Radionuclide activities greater than MDC were detected in all soil samples except for 60Co, 234mPa (1 value >MDC), 238Pu (11 of 96 >MDC) and 241Am (72 of 96 >MDC). Only values greater than the MDC were included in the statistical analyses (Table 14). MANOVA confirmed that the Cactus Flats and Near Field grids were significantly (p<0.05) different. Concentrations of the radionuclides were significantly (p<0.05) higher on the Cactus Flats grid than on the Near Field grid for all radionuclides except 238Pu (Fig. 21). The mean 239,240Pu concentration (0.158 mBq g-1) fell within the range reported by Kenny et al., (1995, Radionuclide Baseline in Soil Near Project Gnome and the Waste Isolation Pilot Plant, Environmental Evaluation Group, Carlsbad, New Mexico) at the WIPP (0 - 0.74 mBq g-1) and was lower than background concentrations found at Hueston Woods and Urbana, Ohio (0.7 - 1.0 mBq g-1) (Alberts J. J., et al., 1980, J.Environ.Qual. 9, 592) and at a series of 15 locations between Ft. Collins and Colorado Springs, Colorado (0.6 - 1.7 mBq g-1) (Hodge, V. et al., 1996, Chemosphere 32, 2067).

The soil textures for all of the samples were very similar, with 81.2-95.8% sand, 1.5-13.4% silt and 1.4% to 5.7% clay. Eleven of the 73 samples were classified as loamy sand and the remainder as sand. All of the radionuclides except 238Pu showed significant (p<0.01) correlations with each of the soil texture classes (Table 12). In general, the percentage of sand, or equivalently, the sum of the percentages of silt and clay, was the best correlate for the concentration of radionuclides. However, for 137Cs and 239,240Pu concentrations the percentage of silt gave a slightly better correlation (Table 12).

In general, the mean ratios of concentrations of the radionuclides to Al were similar between the two grids. Only 239,240Pu showed a significant difference in these ratios, with the ratio for Cactus Flats being significantly greater (p<0.05) than that for Near Field. Concentrations of the radionuclides were also positively correlated by location with the concentrations of many of the analytes. All of the radionuclides were significantly (p<0.01) correlated with Al. Pb showed stronger correlations than Al for all radionuclides except 235U.

The correlation between 40K and 39K was only 0.625 (Fig. 22) and the average ratio (± 95% CI) was 0.00165 (± 0.00009) as compared to the expected isotopic ratio of 0.0117 (Browne, E. and R. B. Firestone, 1986, Table of Radioactive Isotopes, John Wiley & Sons, New York). The cause for this significant departure from the expected ratio and for the variability in the ratio is not known, but some of the variability may be due to the difference in the dissolution methods used for radioanalyses and inorganic analyses. The ratios of 40K to 39K are correlated with soil texture (r for percent sand = 0.369, r for percent silt = -0.325 and r for percent clay = -0.350; n = 96, p<0.01). Because the half-life of 40K is about 109 years, it is unlikely that the effect on the ratios is due to input from decay-depleted sources, such as the local Permian-age potash deposits.

These results demonstrate that significant levels of variability in background levels of soil contaminants and constituents can occur in areas having relatively low variability in soil texture. The high correlation of the radionuclides and many of the non-radioactive metals to the percentages of sand, silt and clay in the soil explains much of the between-sample variability. Actinides can form strong complexes with oxygen-containing ligands (Beal, G. W. and B. Allard, 1981, In Tewari, P. H. (ed.), Adsorption from Aqueous Solutions, Plenum Press, New York; Allard, B., 1982, In Edlestein, N. M. (ed.), Actinides in Perspective, Pergamon Press, Oxford). The fact that 239,240Pu is more strongly correlated with the percentage of silt than with the percentage of clay suggests that many of the clay size particles provide less suitable binding sites for the adsorption of metals than those of typical clays. Actinides also form complexes with humic molecules and these can be significantly more stable than their complexes with simple inorganic ligands (Livens, F. R. and D. L. Singleton, 1991, J. Environ. Rad.13, 323), although subject to a significant concentration effect (Hummel, W. et al., 1999, Radiochim. Acta 84, 111). Therefore, the affinity of 239,240Pu with silt could also be due to a larger concentration of organic material in the silt than in the clay fractions. An investigation of soil particle mineralogy may help explain why the silt fraction is a better correlate of metal concentrations than is the clay fraction.

These data also suggest that the variability in concentrations across locations may arise from a redistribution of contaminated fine soil particles or from a greater degree of entrapment of the contaminants in the upper layer of the soil. Further elucidation of these complex relationships may be possible through selective soil profile analyses.

Inter- and Intra-Laboratory Comparisons

Concentrations of two radioactive elements, uranium (234U, 235U and 238U) and thorium (232Th), were measured independently using ICP-MS and alpha counting methods. A comparison of the estimated concentrations for these analytes showed that the average ratio of the alpha estimate to the ICP-MS estimate for Th was 2.22. The average ratio of the alpha estimate of uranium to that of the ICP-MS method was 8.20. The larger estimates from the actinide methodology may reflect that the samples prepared for actinide analysis undergo more complete dissolution, including fusion of any residues, whereas the samples prepared for ICP-MS analysis undergo only acid leaching. In drinking water analyses reported herein (Radiological and Non-radiological Constituents in Selected Drinking Water Sources), where dissolved forms of the elements are expected to dominate, the average ratio (± SE) for U by the two methods was 1.02 (± 0.09), but 232Th was not detected in drinking water samples. A correlation analysis failed to show any relation between the alpha counting to ICP-MS ratios and soil texture. If the difference in the estimates is due to the difference in dissolution, it appears that the refractory components are not strongly associated with a particular particle size class.

CEMRC measurements of radionuclides in soils were compared to previous measurements presented in the CEMRC 1997 Report performed by Accu-Labs Research, Inc. (Golden, Colorado) for samples collected during 1997 at the Near Field. Comparison of the means of the Accu-Lab and CEMRC results is somewhat difficult because of the limited precision of the Accu-Lab data. However, there were no significant (p<0.05) differences between the estimated activity concentrations except for 230Th, for which the difference between the Accu-Lab mean (14.8 mBq g-1) and the CEMRC mean (10.4 mBq g-1) was significant at p<0.001. The mean 239,240Pu activity concentration (± 95% CI) of the 48 samples measured by CEMRC for the Near Field grid was 0.101 (± 0.014)mBq g-1. This mean is not significantly different from the mean of the Accu-Lab 1997 results [0.14 (± 0.041) mBq g-1].

CEMRC measurements of radionuclides in soils were also compared to measurements presented in the CEMRC 1998 Report performed by DES. Comparisons used (1) a paired t-test on the difference of the CEMRC and DES results paired by sample and (2) correlations between the DES and CEMRC results. The paired t-test showed a significant difference only between the estimates of 232Th concentrations (p<0.05, n=36) (Table 15). The DES results for each radionuclide were found to be uncorrelated with the Accu-Lab results when paired by location (CEMRC 1998 Report) but were well correlated with the CEMRC results when paired by sample. It should be noted that the DES/CEMRC comparison is based on pairing aliquots from the same samples, not on independent samples from the same grid locations, as was the case with the Accu-Labs/DES comparison.

Tables presenting soil data summarized herein are available on the CEMRC web site at http://www.cemrc.org.





Table 12.  Correlations between Analytes and Soil Texture Components


Analyte a,bCorrelation Coefficients Analyte Correlation Coefficients
%
Sand
%
Silt
%
Clay
%
Sand
%
Silt
%
Clay
137Cs -.87 .88 .49 Chloride .19 -.19 -.11
208Tl -.91 .89 .63 Co -.77 .69 .70
212Bi -.88 .83 .66 Cr -.37 .29 .45
212Pb -.90 .86 .66 Cu -.80 .74 .67
214Bi -.85 .83 .59 Eu -.66 .59 .62
214Pb -.87 .85 .59 Fe -.73 .68 .57
228Ac -.92 .89 .64 K -.73 .71 .48
228Th -.93 .90 .66 La -.71 .63 .65
230Th -.90 .88 .60 Li -.39 .32 .43
232Th -.93 .89 .67 Mg -.58 .56 .39
234U -.74 .71 .54 Mn -.74 .72 .52
235U -.49 .51 .24 Ni -.69 .63 .60
238U -.74 .70 .54 Pb -.81 .73 .69
239,240Pu -.84 .85 .49 Sb .02 -.07 .11
40K -.78 .77 .48 Sr -.23 .18 .29
Ag -.28 .28 .16 Th -.71 .63 .68
Al -.71 .69 .51 Ti -.49 .43 .48
As -.67 .59 .63 U -.58 .49 .63
Ba -.49 .43 .45 V -.45 .32 .63
Ca -.25 .25 .16 Zn -.73 .68 .57
Cd -.63 .59 .48        
aPearson product moment correlation coefficients; all correlations are significant at p<0.05 except those in bold
bCorrelations for radionuclides based on 73 samples; correlations for other constituents based on 66 samples




Table 13.  Mean Ratios of Analytes to Al in Soil Samples
from Cactus Flats and Near Field Grids


Analyte Near
Field
aN Cactus
Flats
N Analyte Near
Field
N Cactus
Flats
N
137Cs 2.1E-03 48 3.0E-03 48 Cu 1.2E-03 48 1.2E-03 48
208Tl 1.9E-03 48 2.0E-03 48 Eu 5.4E-05 48 5.2E-05 48
212Bi 6.0E-03 48 6.7E-03 48 Fe 1.5E+00 48 1.6E+00 48
212Pb 5.6E-03 48 5.9E-03 48 Hg 2.9E-03 22 3.9E-03 13
214Bi 5.7E-03 48 5.7E-03 48 K 2.9E-01 48 3.3E-01 48
214Pb 5.9E-03 48 6.0E-03 48 La 1.8E-03 48 1.9E-03 48
228Ac 5.6E-03 48 6.0E-03 48 Li 1.2E-03 48 1.0E-03 48
228Th 5.9E-03 48 6.2E-03 48 Mg 2.1E-01 48 2.1E-01 48
230Th 6.2E-03 48 6.4E-03 48 Mn 2.1E-02 48 2.4E-02 48
232Th 5.6E-03 48 6.0E-03 48 Mo 8.7E-05 29 5.1E-05 38
234U 5.2E-03 48 4.8E-03 48 Na 3.6E-02 7 8.2E-02 9
235U 2.7E-04 48 2.6E-04 48 Ni 1.1E-03 48 1.2E-03 48
238Pu 3.3E-05 4 2.0E-05 7 Nitrate 5.4E-03 48 5.7E-03 48
238U 5.4E-03 48 5.0E-03 48 Nitrite 8.9E-05 4 2.0E-04 5
239,240Pu 6.9E-05 48 1.1E-04 48 Pb 1.6E-03 48 1.9E-03 48
241Am 3.5E-05 31 4.3E-05 40 Phosphate 3.3E-03 36 3.5E-03 41