Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Vadose Zone Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
USGS2009

UZIG USGS research: Advances through interdisciplinary interaction

B Because vadose zone research relates to diverse disciplines, applications, and modes of research, collaboration across traditional operational and topical divisions is especially likely to yield major advances in understanding. The Unsaturated Zone Interest Group (UZIG) is an informal organization sponsored by the USGS to encourage and support interdisciplinary collaboration in vadose or unsaturated zone hydrologic research across organizational boundaries. It includes both USGS and non-USGS scientists. Formed in 1987, the UZIG operates to promote communication, especially through periodic meetings with presentations, discussions, and field trips. The 10th meeting of the UZIG at Los Alamos, NM, in August 2007 was jointly sponsored by the USGS and Los Alamos National Laboratory. Presentations at this meeting served as the initial basis for selecting papers for this special section of Vadose Zone Journal , the purpose of which is to present noteworthy cutting-edge unsaturated zone research promoted by, facilitated by, or presented in connection with the UZIG.

Vadose Zone Journal↗
USGS2018

UZIG research: Measurement and characterization of unsaturated zone processes under wide-ranging climates and changing conditions

Unsaturated zone properties and processes are central to understanding the interacting effects of land-use change, contamination, and hydroclimate on our ability to grow food, sustain clean water supplies, and minimize loss of life and property. Advances in unsaturated zone science are being achieved through collaborations across traditional boundaries where information from biological, physical, and chemical disciplines is combined for new insights. The Unsaturated Zone Interest Group (UZIG) is an organization that exists principally to promote multidisciplinary collaborations and the sharing of ideas, expertise, and technical assets. Here we summarize key findings from 14 papers, several of which originated from a meeting convened by UZIG in 2017 at the University of Florida in Gainesville titled “Land-Use Change, Climate Change, and Hydrologic Extremes: Unsaturated Zone Responses and Feedbacks.” This special section of Vadose Zone Journal contains multidisciplinary research in three general categories relevant to measuring and understanding unsaturated zone responses to changing land uses and climate: (i) unsaturated zone properties and processes; (ii) soil–plant–atmosphere interactions; and (iii) novel field sampling devices. A strong cross-cutting theme in these papers is the value of continuous monitoring data and ways of utilizing them to discover novel hydrologic, biologic, and pedologic information. As climatic and land-use conditions change and demands for resources and stresses on ecosystems continue to intensify, it is vital to improve our fundamental understanding of the processes at work in the unsaturated zone. Toward that goal, we discuss the need for improved ground-based unsaturated zone monitoring networks.

Vadose Zone Journal↗
USGS2021

Developing a strategy for the national coordinated soil moisture monitoring network

Soil moisture is a critical land surface variable, affecting a wide variety of climatological, agricultural, and hydrological processes. Determining the current soil moisture status is possible via a variety of methods, including in situ monitoring, remote sensing, and numerical modeling. Although all of these approaches are rapidly evolving, there is no cohesive strategy or framework to integrate these diverse information sources to develop and disseminate coordinated national soil moisture products that will improve our ability to understand climate variability. The National Coordinated Soil Moisture Monitoring Network initiative has developed a national strategy for network coordination with NOAA's National Integrated Drought Information System. The strategy is currently in review within NOAA, and work is underway to implement the initial milestones of the strategy. This update reviews the goals and steps being taken to establish this national-scale coordination for soil moisture monitoring in the United States.

Vadose Zone Journal↗
USGS2022

Evapotranspiration covers at uranium mill tailings sites

Waste isolation is a key strategy for mitigating risk from municipal solid waste (MSW) and hazardous waste streams. Conventional covers at MSW facilities are designed for a 30-yr post-closure period where compacted soils and geosynthetics are used to minimize percolation into buried waste. Recently, evapotranspiration (ET) covers have shown beneficial use for MSW management. Evapotranspiration covers encourage infiltration, storage, and transpiration of precipitation to minimize percolation. Such covers may also have beneficial use for long-term waste issues, such as at Uranium Mill Tailings Radiation Control Act (UMTRCA) sites. These sites were covered by a clay radon barrier to create tortuous flow paths that allow radioactive decay and attenuation of short-lived, radon-222 gas. For long-term waste isolation, an ET-radon cover may provide greater resilience by exploiting natural processes instead of resisting them. This update presents a review of the current state-of-the-science regarding ET covers and considerations for long-term applications.

Conterminous United States↗
USGS2023

Field evaluation of semi-automated moisture estimation from geophysics using machine learning

Geophysical methods can provide three-dimensional (3D), spatially continuous estimates of soil moisture. However, point-to-point comparisons of geophysical properties to measure soil moisture data are frequently unsatisfactory, resulting in geophysics being used for qualitative purposes only. This is because (1) geophysics requires models that relate geophysical signals to soil moisture, (2) geophysical methods have potential uncertainties resulting from smoothing and artifacts introduced from processing and inversion, and (3) results from multiple geophysical methods are not easily combined within a single soil moisture estimation framework. To investigate these potential limitations, an irrigation experiment was performed wherein soil moisture was monitored through time, and several surface geophysical datasets indirectly sensitive to soil moisture were collected before and after irrigation: ground penetrating radar, electrical resistivity tomography (ERT), and frequency domain electromagnetics (FDEM). Data were exported in both raw and processed form, and then snapped to a common 3D grid to facilitate moisture prediction by standard calibration techniques, multivariate regression, and machine learning. A combination of inverted ERT data, raw FDEM, and inverted FDEM data was most informative for predicting soil moisture using a random regression forest model (one-thousand 60/40 training/test cross-validation folds produced root mean squared errors ranging from 0.025–0.046 cm 3 /cm 3 ). This cross-validated model was further supported by a separate evaluation using a test set from a physically separate portion of the study area. Machine learning was conducive to a semi-automated model-selection process that could be used for other sites and datasets to locally improve accuracy.

Connecticut↗
USGS2025

Identifying preferential flow from soil moisture time series: Review of methodologies

Identifying and quantifying preferential flow (PF) through soil—the rapid movement of water through spatially-distinct pathways in the subsurface—is vital to understanding how the hydrologic cycle responds to climate, land cover, and anthropogenic changes. In recent decades, methods have been developed that use measured soil moisture time series to identify PF. Because they allow for continuous monitoring and are relatively easy to implement, these methods have become an important tool for recognizing when, where, and under what conditions PF occurs. The methods seek to identify a pattern or quantification that indicates the occurrence of PF. Most commonly, the chosen signature is either (1) a nonsequential response to infiltrated water, in which soil moisture responses do not occur in order of shallowest to deepest, or (2) a velocity criterion, in which newly infiltrated water is detected at depth earlier than is possible by nonpreferential flow processes. Alternative signatures have also been developed that have certain advantages but are less commonly utilized. Choosing among these possible signatures requires attention to their pertinent characteristics, including susceptibility to errors, possible bias toward false negatives or false positives, reliance on subjective judgments, and possible requirements for additional types of data. We review 77 studies that have applied such methods, to highlight important information for readers who want to identify PF from soil moisture data, and to inform those who aim to develop new methods or improve existing ones.

Vadose Zone Journal↗
USGS2026

Network-wide assessment of soil water content calibration and sensitivity to biomass proxies using cosmic-ray neutron sensing in the Roaring Fork Basin, Colorado

Soil water content (SWC) is a key state variable of the climate system but is often uncertain in water balance monitoring, especially in alpine environments. SWC measurements can be challenging in alpine environments due to the topography and rocky soils. In 2022, the US Geological Survey's Next Generation Water Observing System Program began research to evaluate water balance monitoring technologies, including cosmic-ray neutron sensors (CRNS). This work evaluated the uncertainty resulting from network-wide calibration of CRNS for SWC monitoring in an alpine watershed and investigated the stability of the calibration parameters across space and time, focusing on potential influence of biomass dynamics. Fifteen stations with moderated and unmoderated (bare) CRNS were deployed and made operational within the Roaring Fork Basin in west-central Colorado. The root mean squared error of the network-wide calibration using the moderated CRNS was 0.042 or 0.047 cm 3 cm −3 , depending on the calibration equation used. Relative SWC dynamics from CRNS were correlated with the in situ probes with a correlation coefficient of 0.91 or 0.87 (depending on calibration equation). We did not find significant relationships between the calibration parameters and stationary site-specific variables. However, the calibration parameters derived from in situ probe SWC dynamics varied over time and were correlated with biomass proxies of cumulative growing degree-day, cumulative growing season index, and bare neutron counts. Future use of the CRNS network can leverage the reliable relative SWC data from network-wide calibration for watershed modeling and continue to research sensitivity of bare neutron measurements to biomass dynamics.

Colorado↗
USGS2002

Kilometer-scale rapid transport of naphthalene sulfonate tracer in the unsaturated zone at the Idaho National Engineering and Environmental Laboratory

To investigate possible long-range flow paths through the interbedded basalts and sediments of a 200-m-thick unsaturated zone, we applied a chemical tracer to seasonally filled infiltration ponds on the Snake River Plain in Idaho. This site is near the Subsurface Disposal Area for radioactive and other hazardous waste at the Idaho National Engineering and Environmental Laboratory. Within 4 mo, we detected tracer in one of 13 sampled aquifer wells, and in eight of 11 sampled perched-water wells as far as 1.3 km away. These detections show that (i) low-permeability layers in the unsaturated zone divert some flow horizontally, but do not prevent rapid transport to the aquifer; (ii) horizontal convective transport rates within the unsaturated zone may exceed 14 m d −1 , perhaps through essentially saturated basalt fractures, tension cracks, lava tubes, or rubble zones; and (iii) some perched water beneath the Subsurface Disposal Area derives from episodic surface water more than 1 km away. Such rapid and far-reaching flow may be common throughout the Snake River Plain, and possibly occurs in other locations that have a geologically complex unsaturated zone and comparable sources of infiltrating water.

Idaho↗
USGS2003

Temperature-profile methods for estimating percolation rates in arid environments

Percolation rates are estimated using vertical temperature profiles from sequentially deeper vadose environments, progressing from sediments beneath stream channels, to expansive basin-fill materials, and finally to deep fractured bedrock underlying mountainous terrain. Beneath stream channels, vertical temperature profiles vary over time in response to downward heat transport, which is generally controlled by conductive heat transport during dry periods, or by advective transport during channel infiltration. During periods of stream-channel infiltration, two relatively simple approaches are possible: a heat-pulse technique, or a heat and liquid-water transport simulation code. Focused percolation rates beneath stream channels are examined for perennial, seasonal, and ephemeral channels in central New Mexico, with estimated percolation rates ranging from 100 to 2100 mm d −1 Deep within basin-fill and underlying mountainous terrain, vertical temperature gradients are dominated by the local geothermal gradient, which creates a profile with decreasing temperatures toward the surface. If simplifying assumptions are employed regarding stratigraphy and vapor fluxes, an analytical solution to the heat transport problem can be used to generate temperature profiles at specified percolation rates for comparison to the observed geothermal gradient. Comparisons to an observed temperature profile in the basin-fill sediments beneath Frenchman Flat, Nevada, yielded water fluxes near zero, with absolute values <10 mm yr −1 For the deep vadose environment beneath Yucca Mountain, Nevada, the complexities of stratigraphy and vapor movement are incorporated into a more elaborate heat and water transport model to compare simulated and observed temperature profiles for a pair of deep boreholes. Best matches resulted in a percolation rate near zero for one borehole and 11 mm yr −1 for the second borehole.

Vadose Zone Journal↗
USGS2003

Spatial variability of groundwater recharge and its effect on shallow groundwater quality in southern New Jersey

Point estimates of groundwater recharge at 48 sediment-coring locations vary substantially (−18.5–1840 cm yr −1 ) in a 930-km 2 area of southern New Jersey. Darcian estimates of steady, long-term recharge made at depth in the unsaturated zone were estimated using pedotransfer functions of soil texture and interpolated (mapped) with nonparametric methods to assess aquifer vulnerability in the area. The probability of exceeding the median recharge (29.1 cm yr −1 ) is low in the southwestern and northeastern portions of the study area and high in the eastern and southeastern portions. Estimated recharge is inversely related to measured percentage clay and positively related to the percentage of well-drained soils near wells. Spatial patterns of recharge estimates, exceedance probabilities, and clay content indicate that sediment texture controls recharge in the study area. Relations with land elevation and a topographic wetness index were statistically insignificant. Nitrate concentration and atrazine (6-chloro- N 2 -ethyl- N 4 -isopropyl-1,3,5-triazine-2,4-diamine) percentage detection in samples of shallow groundwater (typically <10 m) are higher for low recharge sites (≤29.1 cm yr −1 ) than for high recharge sites (>29.1 cm yr −1 ) in agricultural and urban areas. Differences between high and low recharge sites in these areas are highly significant for NO 3 concentration, but not for atrazine concentration.

New Jersey↗
USGS2005

Modeling tritium transport through a deep unsaturated zone in an arid environment

Understanding transport of tritium ( 3 H) in unsaturated zones is critical to evaluating options for waste isolation. Tritium typically is a large component of low-level radioactive waste (LLRW). Studies at the U.S. Geological Survey's Amargosa Desert Research Site (ADRS) in Nevada investigate 3 H transport from a closed LLRW facility. Two boreholes are 100 and 160 m from the nearest waste trench and extend to the water table at 110 m. Soil-water vapor samples from the deep boreholes show elevated levels of 3 H at all depths. The objectives of this study were to (i) test source thermal and gas-advection mechanisms driving 3 H transport and (ii) evaluate model sensitivity to these mechanisms and to selected physical and hydraulic properties including porosity, tortuosity, and anisotropy. A two-dimensional numerical model incorporated a non-isothermal, heterogeneous domain of the unsaturated zone and instantaneous isotopic equilibrium. The TOUGH2 code was used; however, it required modification to account for temperature dependence of both the Henry's law equilibrium constant and isotopic fractionation with respect to tritiated water. Increases in source temperature, pressure, and porosity enhanced 3 H migration, but failed to match measured 3 H distributions. All anisotropic simulations with a source pressure component resembled, in shape, the upper portion of the 3 H distribution of the nearest borehole. Isotopic equilibrium limited migration of 3 H, while effects of radioactive decay were negligible. A 500 Pa pressure increase above ambient pressure in conjunction with a high degree of anisotropy (1:100) was necessary for simulated 3 H transport to reach the nearest borehole.

Vadose Zone Journal↗
USGS2006

Characterization of unsaturated zone hydrogeologic units using matrix properties and depositional history in a complex volcanic environment

Characterization of the physical and unsaturated hydrologic properties of subsurface materials is necessary to calculate flow and transport for land use practices and to evaluate subsurface processes such as perched water or lateral diversion of water, which are influenced by features such as faults, fractures, and abrupt changes in lithology. Input for numerical flow models typically includes parameters that describe hydrologic properties and the initial and boundary conditions for all materials in the unsaturated zone, such as bulk density, porosity, and particle density, saturated hydraulic conductivity, moisture-retention characteristics, and field water content. We describe an approach for systematically evaluating the site features that contribute to water flow, using physical and hydraulic data collected at the laboratory scale, to provide a representative set of physical and hydraulic parameters for numerically calculating flow of water through the materials at a site. An example case study from analyses done for the heterogeneous, layered, volcanic rocks at Yucca Mountain is presented, but the general approach for parameterization could be applied at any site where depositional processes follow deterministic patterns. Hydrogeologic units at this site were defined using (i) a database developed from 5320 rock samples collected from the coring of 23 shallow (<100 m) and 10 deep (500–1000 m) boreholes, (ii) lithostratigraphic boundaries and corresponding relations to porosity, (iii) transition zones with pronounced changes in properties over short vertical distances, (iv) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (v) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. Model parameters developed in this study, and the relation of flow properties to porosity, can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Vadose Zone Journal↗
USGS2004

Constraining the inferred paleohydrologic evolution of a deep unsaturated zone in the Amargosa Desert

Natural flow regimes in deep unsaturated zones of arid interfluvial environments are rarely in hydraulic equilibrium with near-surface boundary conditions imposed by present-day plant–soil–atmosphere dynamics. Nevertheless, assessments of water resources and contaminant transport require realistic estimates of gas, water, and solute fluxes under past, present, and projected conditions. Multimillennial transients that are captured in current hydraulic, chemical, and isotopic profiles can be interpreted to constrain alternative scenarios of paleohydrologic evolution following climatic and vegetational shifts from pluvial to arid conditions. However, interpreting profile data with numerical models presents formidable challenges in that boundary conditions must be prescribed throughout the entire Holocene, when we have at most a few decades of actual records. Models of profile development at the Amargosa Desert Research Site include substantial uncertainties from imperfectly known initial and boundary conditions when simulating flow and solute transport over millennial timescales. We show how multiple types of profile data, including matric potentials and porewater concentrations of Cl − , δD, δ 18 O, can be used in multiphase heat, flow, and transport models to expose and reduce uncertainty in paleohydrologic reconstructions. Results indicate that a dramatic shift in the near-surface water balance occurred approximately 16000 yr ago, but that transitions in precipitation, temperature, and vegetation were not necessarily synchronous. The timing of the hydraulic transition imparts the largest uncertainty to model-predicted contemporary fluxes. In contrast, the uncertainties associated with initial (late Pleistocene) conditions and boundary conditions during the Holocene impart only small uncertainties to model-predicted contemporaneous fluxes.

California, Nevada↗
USGS2004

The role of the unsaturated zone in artificial recharge at San Gorgonio Pass, California

The hydrogeology of the unsaturated zone plays a critical role in determining the suitability of a site for artificial recharge. Optimally, a suitable site has highly permeable soils, a capacity for horizontal flow at the aquifer boundary, a lack of impeding layers, and a thick unsaturated zone. The suitability of a site is often determined by field and laboratory measurements of soil properties, field experiments, and numerical modeling. An artificial recharge site in the San Gorgonio Pass area in southern California, USA was studied to better understand the role of the unsaturated zone in artificial recharge by surface spreading. Field measurements and observations were used to characterize the site and to develop a conceptual model of the unsaturated zone. A numerical model was developed based on the conceptual model and calibrated using data from a 50-d artificial recharge experiment conducted in 1991 and borehole data collected between 1997 and 2002. Results indicate that an impeding layer exists 70 m below land surface that will cause lateral diversion of artificially recharged water, which would spread out and delay recharge to the water table 185 m below land surface.

California↗
USGS2004

Transport and time lag of chlorofluorocarbon gases in the unsaturated zone, Rabis Creek, Denmark

Transport of chlorofluorocarbon (CFC) gases through the unsaturated zone to the water table is affected by gas diffusion, air–water exchange (solubility), sorption to the soil matrix, advective–dispersive transport in the water phase, and, in some cases, anaerobic degradation. In deep unsaturated zones, this may lead to a time lag between entry of gases at the land surface and recharge to groundwater. Data from a Danish field site were used to investigate how time lag is affected by variations in water content and to explore the use of simple analytical solutions to calculate time lag. Numerical simulations demonstrate that either degradation or sorption of CFC-11 takes place, whereas CFC-12 and CFC-113 are nonreactive. Water flow did not appreciably affect transport. An analytical solution for the period with a linear increase in atmospheric CFC concentrations (approximately early 1970s to early 1990s) was used to calculate CFC profiles and time lags. We compared the analytical results with numerical simulations. The time lags in the 15-m-deep unsaturated zone increase from 4.2 to between 5.2 and 6.1 yr and from 3.4 to 3.9 yr for CFC-11 and CFC-12, respectively, when simulations change from use of an exponential to a linear increase in atmospheric concentrations. The CFC concentrations at the water table before the early 1990s can be estimated by displacing the atmospheric input function by these fixed time lags. A sensitivity study demonstrates conditions under which a time lag in the unsaturated zone becomes important. The most critical parameter is the tortuosity coefficient. The analytical approach is valid for the low range of tortuosity coefficients (τ = 0.1–0.4) and unsaturated zones greater than approximately 20 m in thickness. In these cases the CFC distribution may still be from either the exponential or linear phase. In other cases, the use of numerical models, as described in our work and elsewhere, is an option.

Rabis Creek↗
USGS2004

Hydraulic and geochemical framework of the Idaho National Engineering and Environmental Laboratory vadose zone

Questions of major importance for subsurface contaminant transport at the Idaho National Engineering and Environmental Laboratory (INEEL) include (i) travel times to the aquifer, both average or typical values and the range of values to be expected, and (ii) modes of contaminant transport, especially sorption processes. The hydraulic and geochemical framework within which these questions are addressed is dominated by extreme heterogeneity in a vadose zone and aquifer consisting of interbedded basalts and sediments. Hydraulically, major issues include diverse possible types of flow pathways, extreme anisotropy, preferential flow, combined vertical and horizontal flow, and temporary saturation or perching. Geochemically, major issues include contaminant mobility as influenced by redox conditions, the concentration of organic and inorganic complexing solutes and other local variables, the interaction with infiltrating waters and with the contaminant source environment, and the aqueous speciation of contaminants such as actinides. Another major issue is the possibility of colloid transport, which inverts some of the traditional concepts of mobility, as sorbed contaminants on mobile colloids may be transported with ease compared with contaminants that are not sorbed. With respect to the goal of minimizing aquifer concentrations of contaminants, some characteristics of the vadose zone are essentially completely favorable. Examples include the great thickness (200 m) of the vadose zone, and the presence of substantial quantities of fine sediments that can retard contaminant transport both hydraulically and chemically. Most characteristics, however, have both favorable and unfavorable aspects. For example, preferential flow, as promoted by several notable features of the vadose zone at the INEEL, can provide fast, minimally sorbing pathways for contaminants to reach the aquifer easily, but it also leads to a wide dispersal of contaminants in a large volume of subsurface material, thus increasing the opportunity for dilution and sorption.

Vadose Zone Journal↗
USGS2006

Borehole environmental tracers for evaluating net infiltration and recharge through desert bedrock

Permeable bedrock aquifers in arid regions are being increasingly developed as water supplies, yet little is generally known about recharge processes and spatial and temporal variability. Environmental tracers from boreholes were used in this study to investigate net infiltration and recharge to the fractured Navajo Sandstone aquifer. Vadose zone tracer profiles at the Sand Hollow study site in southwestern Utah look similar to those of desert soils at other sites, indicating the predominance of matrix flow. However, recharge rates are generally higher in the Navajo Sandstone than in unconsolidated soils in similar climates because the sandstone matrix allows water movement but not root penetration. Water enters the vadose zone either as direct infiltration of precipitation through exposed sandstone and sandy soils or as focused infiltration of runoff. Net infiltration and recharge exhibit extreme spatial variability. High-recharge borehole sites generally have large amounts of vadose zone tritium, low chloride concentrations, and small vadose zone oxygen-18 evaporative shifts. Annual net-infiltration and recharge rates at different locations range from about 1 to 60 mm as determined using vadose zone tritium, 0 to 15 mm using vadose zone chloride, and 3 to 60 mm using groundwater chloride. Environmental tracers indicate a cyclical net-infiltration and recharge pattern, with higher rates earlier in the Holocene and lower rates during the late Holocene, and a return to higher rates during recent decades associated with anomalously high precipitation during the latter part of the 20th century. The slightly enriched stable isotopic composition of modern groundwater indicates this recent increase in precipitation may be caused by a stronger summer monsoon or winter southern Pacific El Nin??o storm track. ?? Soil Science Society of America.

Vadose Zone Journal↗
USGS2005

Plant-based plume-scale mapping of tritium contamination in desert soils

Plant-based techniques were tested for field-scale evaluation of tritium contamination adjacent to a low-level radioactive waste (LLRW) facility in the Amargosa Desert, Nevada. Objectives were to (i) characterize and map the spatial variability of tritium in plant water, (ii) develop empirical relations to predict and map subsurface contamination from plant-water concentrations, and (iii) gain insight into tritium migration pathways and processes. Plant sampling [creosote bush, Larrea tridentata (Sessé & Moc. ex DC.) Coville] required one-fifth the time of soil water vapor sampling. Plant concentrations were spatially correlated to a separation distance of 380 m; measurement uncertainty accounted for <0.1% of the total variability in the data. Regression equations based on plant tritium explained 96 and 90% of the variation in root-zone and sub-root-zone soil water vapor concentrations, respectively. The equations were combined with kriged plant-water concentrations to map subsurface contamination. Mapping showed preferential lateral movement of tritium through a dry, coarse-textured layer beneath the root zone, with concurrent upward movement through the root zone. Analysis of subsurface fluxes along a transect perpendicular to the LLRW facility showed that upward diffusive-vapor transport dominates other transport modes beneath native vegetation. Downward advective-liquid transport dominates at one endpoint of the transect, beneath a devegetated road immediately adjacent to the facility. To our knowledge, this study is the first to document large-scale subsurface vapor-phase tritium migration from a LLRW facility. Plant-based methods provide a noninvasive, cost-effective approach to mapping subsurface tritium migration in desert areas.

Vadose Zone Journal↗