|
|
||||||||
a U.S. Geological Survey, 333 West Nye Lane, Suite 203, Carson City, NV 89706
b U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025
* Corresponding author (andraski{at}usgs.gov)
1 The use of trade names is for identification purposes only and does not imply endorsement by the U.S. Geological Survey. ![]()
Received 30 March 2005.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: ADRS, Amargosa Desert Research Site LLRW, low-level radioactive waste
| INTRODUCTION |
|---|
|
|
|---|
Subsurface monitoring helps us quantify and manage environmental risks. Traditional monitoring methods include sampling of bulk soil, soil water, soil gas, and groundwater (Faybishenko, 2000; Lindgren, 2000). While these methods provide detailed information on the presence of contaminants, the associated site disturbance and costs of equipment, installation, and maintenance usually limit the number of samples. Because traditional methods provide only point measurements, limited sampling makes it difficult to adequately assess potentially contaminated areas. Novel approaches are thus needed for detection and monitoring. Such approaches must be cost effective, sensitive enough to provide early warning of contaminant releases, and accurate and robust enough for routine use in assessments of containment and remediation strategies.
In the late 1980s, trees were found to be useful indicators of subsurface tritium movement from the Maxey Flats, KY LLRW disposal site (Rickard and Kirby, 1987; Kalisz et al., 1988). Water taken up by tree roots showed tritium to be migrating along the regolithbedrock interface. Andraski et al. (2003) recently introduced an efficient method of collecting plant water for tritium analysis. That study also documented a high correlation between tritium levels in plant water and in root-zone soil-water vapor.
The present study builds on previous work to evaluate the applicability of plant-based methods for large-scale contaminant-plume mapping in an arid environment. Testing was done adjacent to a closed LLRW facility in the Amargosa Desert, Nevada. Specific objectives were to (i) characterize and map the spatial variability of tritium in plant water (creosote bush), (ii) develop empirical relations to predict and map subsurface tritium contamination from plant-water concentrations, and (iii) gain insight into tritium migration pathways and processes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The ADRS is in the Mojave Desert ecosystem, one of the driest regions of the United States. Annual average precipitation (19812001) is 108 mm yr1. The volume-weighted tritium concentration in precipitation (November 1997April 2001) averaged 1.1 Bq L1. Surface soils belong to the Yermo (loamy-skeletal, mixed, superactive, calcareous, thermic Typic Torriorthent)Arizo (sandy-skeletal, mixed, thermic Typic Torriorthent) association. Underlying sediments are primarily fluvial deposits. Detailed information on soil properties is given in Andraski (1996). Depth to the water table ranges from 85 to 115 m below land surface.
Vegetation is sparse. Creosote bush, the most abundant shrub of the North American warm deserts, is the dominant species. Creosote bush forms nearly pure stands in many desert areas (Smith et al., 1997). At the ADRS, total shrub-species ground cover is 6.4% [based on April 2001 data collected using the line-transect method (Smith, 1974)]; ground cover by species is 4.0% creosote bush, 1.2% shadscale [or spiny saltbush, Atriplex confertifolia (Torr. & Frém.) S. Wats.], 0.8% burrobush [Ambrosia dumosa (Gray) Payne], and 0.4% wolfberry [Lycium pallidum Miers]. Observations indicate creosote-bush roots are horizontally extensive; literature values for the radial extent of creosote-bush roots exceed 4 m (Cannon, 1911; Gile et al., 1998). A 4-m rooting radius is about 10 times the measured creosote-bush canopy radius (0.43 m). The rooting depth of creosote bush generally corresponds to the maximum-annual penetration depth of precipitation, which is about 0.75 to 1 m at the ADRS (Andraski, 1997). The penetration depth of precipitation is limited by plant-water uptake and a laterally extensive gravelly layer (very gravelly sand, about 1 m thick) that underlies the surface soil and impedes downward percolation of liquid water under unsaturated conditions. Multiyear measurements have shown that, in vegetated parts of the study area, root-zone soil-water content varies temporally between 0.02 and 0.14 m3 m3 (Fischer, 1992; Andraski, 1997; Johnson et al., 2002). The gravelly sub-root-zone layer has a temporally invariant soil water content of approximately 0.05 m3 m3. In devegetated areas, soil moisture accumulates and penetrates below the root zone.
Sampling and Tritium Analysis
Creosote bush samples (n = 103) were collected within a 63-ha area adjacent to the waste facility. Creosote bush was selected for this study because of its dominance and drought-resistant evergreen character. Most plant samples were collected along transects that ran parallel to the perimeter of the waste facility. Spacing between adjacent samples usually ranged between about 50 and 100 m. Closer spacing was used nearer the facility to better allow for identification of potentially localized areas of contamination. The ultimate combination of sample number and locations provided reasonable coverage of the study area and allowed all plant samples to be collected in 1 d. Soil water vapor samples were collected at selected locations to develop empirical relations between plant and subsurface tritium concentrations. Root-zone soil-gas tubes (n = 10) were installed to a nominal depth of 0.5 m. Sub-root-zone soil-gas tubes (n = 17) were installed into the underlying gravelly layer, to a nominal depth of 1.5 m. Each soil-gas tube was located 0.5 to 5 m from the center of a sampled plant. Canopy air samples were collected 0.5 m above land surface along a transect perpendicular to the LLRW facility. Background tritium concentrations in plant water, soil water vapor, and air were determined at the control site.
Andraski et al. (2003) detailed field and laboratory procedures used in the collection, preparation, and analysis of samples. Briefly, the plant-water method involved hand-stripping of foliage (170 g on average) from randomly selected branches, solar distillation in plastic bags, collection of distillate (23 mL on average) by pipet, and extraction of scintillation-interfering constituents by filtration and adsorption on a graphite-based medium before tritium analysis. Foliage samples were weighed before distillation, pipet volumes were recorded, and foliage samples were reweighed after pipeting. These data indicated that, on average, about 2 to 3 g of liquid was lost during the field distillation and sample transfer process. Soil water vapor and canopy water vapor were collected by pulling air through a dry-ice freeze trap using a downstream pump. Water vapor collected as ice in the trap. Tritium analyses were done according to the direct liquid-scintillation counting method described by Thatcher et al. (1977) at the USGS tritium laboratory in Menlo Park, CA. Scintillation counting uncertainty (±one sigma) for canopy air, plant water, and soil water vapor samples averaged 1, 7, and 10%, respectively.
Data and Spatial Analyses
Univariate statistics on tritium concentrations in plant water were used to summarize the data and determine whether transformation was needed before further statistical analysis. Geostatistical techniques (Journel and Huijbregts, 1978; Isaaks and Srivastava, 1989) were used to evaluate spatial correlation in the data and to generate a complete map of tritium concentrations in plant water. The semivariogram was fitted with a Gaussian model:
![]() | [1] |
is the semivariogram of plant-water tritium, h is separation distance, Cn is nugget variance, C0 is structural variance, Cn + C0 is sill variance, and a is correlation length. Equation [1] describes the asymptotic approach to a sill; the effective range (re) is the distance at which the semivariogram reaches 95% of the sill (re
31/2a). Parameters for the model semivariogram were optimized using nonlinear regression procedures. The modeled semivariogram was used for ordinary-kriging estimation and mapping of plant-water tritium concentrations throughout the study area. Statistical analyses were done using SAS software (SAS Institute, 2001).1 Data from collocated plant-water and soil water vapor sampling sites were analyzed using linear regression procedures. The resulting empirical relations were used to predict subsurface tritium concentrations from the more abundant (and more easily determined) plant-water concentrations. Root-zone and sub-root-zone tritium were mapped by combining the kriged plant-water concentrations with the appropriate regression equation.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
75100 m between adjacent samples) would give results similar to those shown in Fig. 2. For example, model-semivariogram parameters for the subset-data analysis were Cn = 0.0 and C0 = 1.29 in units of (log10 Bq L1)2; a = 225 m, giving an effective range of 390 m.
|
|
Andraski et al. (2003) found a good correlation between tritium in plant water and root-zone soil water vapor. Thus, we developed empirical relations to predict root-zone and sub-root-zone tritium concentrations from the more easily determined plant-water concentrations. These results are shown in Fig. 4a and 4b , respectively. In both cases, the regression relation between plant- and subsurface-tritium concentrations is highly significant (p < 0.01). Correlation coefficients (r2) showed that the empirical regressions explained 96% of the variation in the root-zone data and 90% of the variation in the subroot-zone data. The decrease in the r2 value for the sub-root zone prediction reflects the lack of direct contact between plant roots and soil water in that zone. Nevertheless, these results show that plant-water data can be used to predict specific tritium concentrations in both root-zone and sub-root-zone soil water vapor.
|
|
527335 m), along with background data collected 3 km south of the waste facility. At distances of 0, 25, 100, and 200 m from the perimeter of the facility, the sub-root-zone tritium concentration exceeded the root-zone concentration by a factor of two to six. At a distance of 300 m, all measured concentrations were about twice background levels, but differences among soil and plant samples were not readily apparent (compared with the one-sigma scintillation counting errors of ± 0.6 Bq L1 for this location). At a distance of 400 m, all tritium concentrations were statistically equivalent to those measured at the background site.
|
Soil water vapor sampling was limited in the western hot spot, but nevertheless showed a vertical concentration gradient like that observed in the southern hot spot. For example, concentrations approximately 25 m west of the LLRW area (easting 527208 m, northing 4068945 m) decreased as follows: sub-root zone (3623 Bq L1), root zone (3505 Bq L1), and plant water (2960 Bq L1).
The results show that plant sampling in combination with a limited number of in situ soil water vapor samples is useful for plume-scale mapping of subsurface tritium contamination. The plant sampling method is cost effective. Collection and preparation of plant samples for tritium analysis required only one-fifth the time of soil water vapor samples. Plant sampling is relatively noninvasive with respect to plants and completely noninvasive with respect to the soil. Plant water provides a volume-integrated (versus point) sample that reflects the relatively large soil volume exploited by the active root system. A survey of more than 200 species indicates many plants can develop far-reaching root systems in the absence of restrictive subsoil characteristics (Stone and Kalisz, 1991). This underscores the potential for other plants to act as volume-integrated samplers of subsurface contamination.
The greatest source of uncertainty in plant-based plume-scale mapping stems from the empirical relations between plant- and subsurface-tritium concentrations. Other, but relatively minor, sources of uncertainty are those associated with the solar-distillation and solid-phase extraction of plant water, and scintillation counting (Andraski et al., 2003). Future work must establish the presumably species- and environment-specific relation between plant- and soil water contamination. Nevertheless, our approach is likely transferable to additional species and environments. Note that site-specific testing will be needed to establish solid-phase-extraction requirements for the plant water of interest and to establish empirical relations between plant and subsurface contaminant concentrations. Depending on the objectives and results of such testing, the end-user may have the option to use a plant-based method as a general indicator or a quantitative predictor of subsurface tritium contamination. In addition, the resultant tritium data may serve as a proxy for constraining the migration patterns of other liquid-phase or gas-phase (e.g., mercury) contaminants that can be associated with buried LLRW.
Transport Pathways and Processes
The distribution of tritium in relation to site features provides insight into field-scale transport. Data collected as part of this study, together with data from deep unsaturated-zone boreholes (Mayers, 2003; Stonestrom et al., 2004), show that tritium moves preferentially away from the waste source through dry, coarse-textured layers beneath the root zone, with subsequent upward release to the surface through the root zone. Near the southern hot spot at Borehole UZB-3 (
30 m south of LLRW area; easting 527336 m, northing 4068713 m), the tritium concentration in the shallow sub-root-zone gravel layer (2248 Bq L1 at 1.5-m depth) exceeded values at the 0.50-m depth (710 Bq L1) and 4.95-m depth (335 Bq L1) by factors of 3 and 7, respectively. Likewise, at Borehole UZB-2 (
100 m south of LLRW area; easting 527257 m, northing 4068645 m), the tritium concentration in the shallow sub-root-zone gravel layer (1429 Bq L1 at 1.5-m depth) exceeded values at the 0.50-m depth (375 Bq L1) and 5.5-m depth (85 Bq L1) by factors of 4 and 17, respectively. Earlier work postulated subsurface liquid transport along preferential flow paths as a mechanism for tritium contamination observed at the UZB-2 borehole (Striegl et al., 1996). However, multiple-year soil moisture measurements (water content and water potentialFischer, 1992; Andraski, 1997; Johnson et al., 2002) and electromagnetic surveys (B.R. Scanlon, B.J. Andraski, and D.A. Stonestrom, unpublished data, 1999) showed no evidence of elevated soil moisture in the vicinity of the southern hot spot. Thus, the mapped contamination (Fig. 3a and 5) appears to be primarily the result of subsurface vapor-phase transport.
Liquid and vapor transport are coupled processes, requiring integrated analysis for accurate work (Thorstenson and Pollock, 1989). For purposes of initial data evaluation, a single-phase independent-process approach was used to estimate the relative magnitudes of vertical, diffusive and advective fluxes of tritiated water in the liquid and vapor phases. Flux estimates were based on measured root-zone and sub-root-zone soil water vapor tritium concentrations at distances of 0, 25, 100, and 200 m from the perimeter of the LLRW area (Fig. 6). Tritium concentrations in the liquid phase were calculated assuming isotopic equilibrium between water vapor and liquid water and using a fractionation factor of 1.108 (Ferronsky and Polyakov, 1982).
Diffusive fluxes of tritiated water were estimated using (Pruess et al., 1999):
![]() | [2] |
is porosity,
o
ß is a tortuosity factor (
o
ß =
1/3 Sß10/3), Sß is phase saturation,
ß is phase density, Ddß is diffusion coefficient of tritiated water in phase ß, and Xß is mass fraction in phase ß. We used
= 0.25, the average total porosity for root-zone and sub-root-zone soil measurements. For the 25-, 100-, and 200-m locations, calculated Sl and Sg values were based on an average measured soil water content of 0.05 m3 m3. For the 0-m (bare soil) location, Sl and Sg values were based on a soil water content of 0.08 m3 m3 reported for devegetated soil (Johnson et al., 2002). Values used for the other parameters were based on published information:
l = 997.9 kg m3,
g = 1.077 kg m3, Dl = 2.05 x 109 m2 s1, and Dg = 2.60 x 105 m2 s1. Kinetic fractionation effects during sample collection were assumed to produce small errors relative to the range of absolute tritium concentrations.
The advective liquid-phase flux of tritiated water, Jal, was estimated using
![]() | [3] |
l is depth-averaged mass fraction in the liquid phase. For the 25-, 100-, and 200-m locations, Jw was based on the average value beneath native vegetation (1.2 x 1013 kg m2 s1) reported by Andraski (1997). For the 0-m location, Jw (4.3 x 105 kg m2 s1) was based on Darcy's Law using average measured water-potential gradients beneath devegetated soil (Andraski, 1997) and the estimated unsaturated hydraulic conductivity for sub-root-zone-gravel at a soil-water content of 0.08 m3 m3 (Andraski, 1996).
The advective vapor-phase flux of tritiated water due to barometric "pumping" (a dispersive flux) was estimated using (Parker, 2003)
![]() | [4] |
a
Pdrv)/(Pot) where
a is air-filled porosity,
P is average diurnal barometric forcing, d is depth of barometric pressure propagation, r is dispersivity to travel distance ratio, v is vertical travel distance to top of contaminant zone, Po is mean barometric pressure, and t is period We used v = 1 m (distance to the top of the sub-root-zone gravel); d = 108 m (depth of capillary fringe); and t = 1 d,
P = 0.4 kPa, and Po = 92 kPa, based on data from Smith et al. (1999). Following Parker (2003) we assumed r = 0.02. Estimated diffusive and advective fluxes of tritiated water calculated using Eq. [2] through [4] are given in Table 2. For the 0-m location, water accumulation beneath devegetated soil produces a downward-directed, advective-liquid flux more than 103 times greater than the other fluxes. In contrast, fluxes beneath vegetation at distances of 25, 100, and 200 m are upward and dominated by vapor-phase diffusion (Table 2). The diffusive vapor-phase fluxes are one to five orders of magnitude greater than other fluxes in vegetated areas. The diffusion coefficient of tritiated water in vapor is about 104 times greater than that for tritiated water in liquid; this more than offsets differences in tritiated water mass fraction between soil gas and liquid soil water. The diffusive/advective vapor flux ratio is approximately 67 at the 0-m location and about 100 elsewhere. Consistent with results of Parker (2003), relative to advective vapor transport, diffusive vapor transport is more important at the shallow depths considered here and at sites with lower soil water contents. Under present climatic conditions, the efficiency of root-zone water depletion by native plants (Andraski, 1997) and the associated transpiration of tritiated-liquid water will likely support relatively continuous expulsion of tritiated water vapor from the sub-root zone.
|
Evaluation of the western hot spot is speculative, as no deep data are available. As shown, near-surface vertical tritium-concentration gradients are similar to those in the southern hot spot. The drainage diversion ditch (Fig. 1), excavated during 19871988, cut through the near-surface gravel layer. The fact that the plant and shallow-subsurface tritium distributions do not reflect the discontinuity (Fig. 1, 3a, and 5) suggests that (i) the leading edge of a shallow plume migrated beyond the ditch before 1987-88 or (ii) the primary path of lateral migration is through gravel layers at depth. Recognizing the upward migration of tritium from sub-root-zone to root-zone soil layers documented in this study, as well as the importance of upward flow processes in deep unsaturated zones, it is plausible that tritium contamination from a subsurface source well below the root zone could move upward and be detected in near-surface soils and in plants. Insight into the historic record and rate of plume advancement might be gained by tritium analysis of plant-growth rings (Love et al., 2002).
Although vapor-phase transport appears to be the primary mechanism of near-surface tritium migration from the waste area, the correspondence between the finger-like projection in the 10 Bq L1 contour line and the trail south of the chemical waste area (Fig. 1) suggests an off-site surface spill may have contributed to elevated concentrations in this area. Increased subsurface concentrations in the finger area (Fig. 5) are reasonable if a surface spill caused focused infiltration and percolation of tritiated water.
To our knowledge, this study is the first to document extensive lateral, subsurface vapor-phase transport of tritium from a LLRW area. The observed extent of tritium transport, to distances >300 m, is inconsistent with predictions from process-based numerical models (e.g., Striegl et al., 1996; Mayers et al., 2005). Results of this study will help test and refine models of tritium transport. Better understanding of the relevant processes will aid containment and monitoring efforts at contaminated sites.
The strong relation between tritium in plant water and tritium in sub-root-zone soil water vapor demonstrates the role of desert vegetation in the upward migration of waterborne contaminants from beneath the root zone. This has implications beyond monitoring. Phytoremediation is often assumed to be limited to material that is in direct contact with plant roots (Rock, 1997). The present study shows that the remedial effects of desert vegetation on tritium contamination can extend beneath the root zone.
| CONCLUSIONS |
|---|
|
|
|---|
Plant data revealed hot spots and other trends in the distribution of tritium contamination. Geostatistical analysis showed that plant tritium was spatially correlated to a separation distance of 380 m. Measurement uncertainty accounted for <0.1% of the total variability in the data. Kriging produced full-coverage maps of estimated plant contamination. Kriging standard errors were small, typically <0.2 log10 Bq L1. Results provided a synoptic view of the plant tritium distributions. The method is simple and robust enough for repeat sampling to determine evolution of a plume.
Tritium concentrations in plant water accurately reflected those in root-zone and sub-root-zone soil water vapor. The linear regression relations between measured plant and subsurface tritium concentrations were highly significant (p < 0.01). Regression equations explained 96 and 90% of the variation in root-zone and sub-root-zone data, respectively. The empirical extrapolation from plant tritium to subsurface tritium represents the greatest source of uncertainty in the plant-based approach to subsurface contaminant mapping. Other relatively minor sources of uncertainty include those associated with collection and preparation of plant water for tritium analysis.
Preliminary analysis of results in relation to site features provided insight into tritium transport pathways and processes. Tritiated water vapor is inferred to be moving preferentially from the waste source through a coarse-textured layer beneath the root zone. To our knowledge, this study is the first to document large-scale subsurface vapor-phase tritium migration from a LLRW disposal area. Transport beneath native vegetation is upward from the sub-root-zone gravel to the root zone, where diffusive vapor-phase fluxes exceed the advective vapor-phase fluxes by a factor of about 100 and liquid-phase (diffusive and advective) fluxes by a factor of about 20 to 105. Plants can play an important role in the detection, movement, and potential remediation of tritium contamination in desert areas.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. A. Garcia, B. J. Andraski, D. A. Stonestrom, C. A. Cooper, M. J. Johnson, R. L. Michel, and S. W. Wheatcraft Transport of Tritium Contamination to the Atmosphere in an Arid Environment Vadose Zone J., May 21, 2009; 8(2): 450 - 461. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Mayers, B. J. Andraski, C. A. Cooper, S. W. Wheatcraft, D. A. Stonestrom, and R. L. Michel Modeling Tritium Transport Through a Deep Unsaturated Zone in an Arid Environment Vadose Zone J., October 10, 2005; 4(4): 967 - 976. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||