USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes.”

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 289 records · Page 16Linked to original sources

Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 2. Quantitative hydrologic assessment

Groundwater flow modeling was used to quantitatively assess the hydrologic processes affecting ground water and solute movement to drain laterals. Modeling results were used to calculate the depth distribution of groundwater flowing into drain laterals at 1.8 m (drain lateral 1) and 2.7 m (drain lateral 2) below land surface. The simulations indicated that under nonirrigated conditions about 89% of the flow in drain lateral 2 was from groundwater originating from depths greater than 6 m below land surface. The deep groundwater has higher selenium concentrations than shallow groundwater. Simulation of irrigated conditions indicates that as recharge (deep percolation) increases, the proportional contribution of deep groundwater to drain lateral flow decreases. Groundwater flow paths and travel times estimated from the simulation results indicate that groundwater containing high concentrations of selenium (greater than 780 μg L −1 ) probably will continue to enter drain lateral 2 for decades.

California↗

Dating of shallow groundwater: Comparison of the transient tracers 3H/3He, chlorofluorocarbons, and 85Kr

This paper describes a direct comparison of apparent ages derived from 3 H/ 3 He, chlorofluorocarbons (CCl 3 F and CCl 2 F 2 ), and 85 Kr measurements in shallow groundwater. Wells chosen for this study are completed in the unconfined surficial aquifers in late Cenozoic Atlantic Coastal Plain sediments of the Delmarva Peninsula, on the east coast of the United States. Most of the apparent tracer ages agree within 2 years of each other for recharge dates between 1965 and 1990. Discrepancies in apparent tracer ages usually can be explained by hydrological processes such as mixing in a discharge area. Recharge rate calculations based on apparent tracer age gradients at multilevel well locations agree with previous recharge estimates. High recharge rates on the Delmarva Peninsula result in nearly complete dissolved-gas confinement in the groundwater. The remarkable agreement between the different tracer ages indicates negligible mixing of waters of different ages, insignificant dispersion, minimal gas loss to the atmosphere, and insignificant sorption-desorption processes at this location.

Water Resources Research↗

Tritium-helium 3 dating under complex conditions in hydraulically stressed areas of a buried-valley aquifer

The 3 H- 3 He dating method is applied in a buried-valley aquifer near Dayton, Ohio. The study area is large, not all sampling locations lie along well-defined flow paths, and existing wells with variable screen lengths and diameters are used. Reliable use of the method at this site requires addressing several complications: (1) The flow system is disturbed because of high pumping rates and induced infiltration; (2) tritium contamination is present in several areas of the aquifer; and (3) radiogenic helium concentrations are elevated in a significant number of the wells. The 3 H- 3 He ages are examined for self-consistency by comparing the reconstructed tritium evolution to the annual weighted tritium measured in precipitation; deviations result from dispersion, tritium contamination, and mixing. 3 H- 3 He ages are next examined for consistency with chlorofluorocarbon ages; the agreement is poor because of degradation of CFCs. Finally, the 3 H- 3 He ages are examined for consistency with the current understanding of local hydrologic processes; the ages are generally supported by hydrogeologic data and the results of groundwater flow modeling coupled with particle-tracking analyses.

Ohio↗

Metal mobilization from thawing permafrost to aquatic ecosystems is driving rusting of Arctic streams

Climate change in the Arctic is altering watershed hydrologic processes and biogeochemistry. Here, we present an emergent threat to Arctic watersheds based on observations from 75 streams in Alaska’s Brooks Range that recently turned orange, reflecting increased loading of iron and toxic metals. Using remote sensing, we constrain the timing of stream discoloration to the last 10 years, a period of rapid warming and snowfall, suggesting impairment is likely due to permafrost thaw. Thawing permafrost can foster chemical weathering of minerals, microbial reduction of soil iron, and groundwater transport of metals to streams. Compared to clear reference streams, orange streams have lower pH, higher turbidity, and higher sulfate, iron, and trace metal concentrations, supporting sulfide mineral weathering as a primary mobilization process. Stream discoloration was associated with dramatic declines in macroinvertebrate diversity and fish abundance. These findings have considerable implications for drinking water supplies and subsistence fisheries in rural Alaska.

Alaska↗

Biogeochemical and hydrologic synergy control mercury fate in an arid land river-reservoir system

Reservoirs in arid landscapes provide critical water storage and hydroelectric power but influence the transport and biogeochemical cycling of mercury (Hg). Improved management of reservoirs to mitigate the supply and uptake of bioavailable methylmercury (MeHg) in aquatic food webs will benefit from a mechanistic understanding of inorganic divalent Hg (Hg( II )) and MeHg fate within and downstream of reservoirs. Here, we quantified Hg( II ), MeHg, and other pertinent biogeochemical constituents in water (filtered and associated with particles) at high temporal resolution from 2016–2020. This was done (1) at inflow and outflow locations of three successive hydroelectric reservoirs (Snake River, Idaho, Oregon) and (2) vertically and longitudinally within the first reservoir (Brownlee Reservoir). Under spring high flow, upstream inputs of particulate Hg (Hg( II ) and MeHg) and filter-passing Hg( II ) to Brownlee Reservoir were governed by total suspended solids and dissolved organic matter, respectively. Under redox stratified conditions in summer, net MeHg formation in the meta- and hypolimnion of Brownlee reservoir yielded elevated filter-passing and particulate MeHg concentrations, the latter exceeding 500 ng g −1 on particles. Simultaneously, the organic matter content of particulates increased longitudinally in the reservoir (from 9–29%) and temporally with stratified duration. In late summer and fall, destratification mobilized MeHg from the upgradient metalimnion and the downgradient hypolimnion of Brownlee Reservoir, respectively, resulting in downstream export of elevated filter-passing MeHg and organic-rich particles enriched in MeHg (up to 43% MeHg). We document coupled biogeochemical and hydrologic processes that yield in-reservoir MeHg accumulation and MeHg export in water and particles, which impacts MeHg uptake in aquatic food webs within and downstream of reservoirs.

Environmental Science: Processes & Impacts↗

Climate change and physical disturbance cause similar community shifts in biological soil crusts

Biological soil crusts (biocrusts)—communities of mosses, lichens, cyanobacteria, and heterotrophs living at the soil surface—are fundamental components of drylands worldwide, and destruction of biocrusts dramatically alters biogeochemical processes, hydrology, surface energy balance, and vegetation cover. While there has been long-standing concern over impacts of 5 physical disturbances on biocrusts (e.g., trampling by livestock, damage from vehicles), there is also increasing concern over the potential for climate change to alter biocrust community structure. Using long-term data from the Colorado Plateau, USA, we examined the effects of 10 years of experimental warming and altered precipitation (in full-factorial design) on biocrust communities, and compared the effects of altered climate with those of long-term physical 10 disturbance (>10 years of replicated human trampling). Surprisingly, altered climate and physical disturbance treatments had similar effects on biocrust community structure. Warming, altered precipitation frequency [an increase of small (1.2 mm) summer rainfall events], and physical disturbance from trampling all promoted early successional community states marked by dramatic declines in moss cover and increased cyanobacteria cover, with more variable effects 15 on lichens. While the pace of community change varied significantly among treatments, our results suggest that multiple aspects of climate change will affect biocrusts to the same degree as physical disturbance. This is particularly disconcerting in the context of warming, as temperatures for drylands are projected to increase beyond those imposed by the climate treatments used in our study.

PNAS↗

A remote-sensing, GIS-based approach to identify, characterize, and model spawning habitat for fall-run chum salmon in a sub-arctic, glacially fed river

At northern limits of a species’ distribution, fish habitat requirements are often linked to thermal preferences, and the presence of overwintering habitat. However, logistical challenges and hydrologic processes typical of glacial systems could compromize the identification of these habitats, particularly in large river environments. Our goal was to identify and characterize spawning habitat for fall-run chum salmon Oncorhynchus keta and model habitat selection from spatial distributions of tagged individuals in the Tanana River, Alaska using an approach that combined ground surveys with remote sensing. Models included braiding, sinuosity, ice-free water surface area (indicating groundwater influence), and persistent ice-free water (i.e., consistent presence of ice-free water for a 12-year period according to satellite imagery). Candidate models containing persistent ice-free water were selected as most likely, highlighting the utility of remote sensing for monitoring and identifying salmon habitat in remote areas. A combination of ground and remote surveys revealed spatial and temporal thermal characteristics of these habitats that could have strong biological implications. Persistent ice-free sites identified using synthetic aperture radar appear to serve as core areas for spawning fall chum salmon, and the importance of stability through time suggests a legacy of successful reproductive effort for this homing species. These features would not be captured with a one-visit traditional survey but rather required remote-sensing monitoring of the sites through time.

Transactions of the American Fisheries Society↗

Efficacy of environmental DNA to detect and quantify Brook Trout populations in headwater streams of the Adirondack Mountains, New York

Environmental DNA (eDNA) analysis is rapidly evolving as a tool for monitoring the distributions of aquatic species. Detection of species’ populations in streams may be challenging because the persistence time for intact DNA fragments is unknown and because eDNA is diluted and dispersed by dynamic hydrological processes. During 2015, the DNA of Brook Trout Salvelinus fontinalis was analyzed from water samples collected at 40 streams across the Adirondack region of upstate New York, where Brook Trout populations were recently quantified. Study objectives were to evaluate different sampling methods and the ability of eDNA to accurately predict the presence and abundance of resident Brook Trout populations. Results from three-pass electrofishing surveys indicated that Brook Trout were absent from 10 sites and were present in low (<100 fish/0.1 ha), moderate (100–300 fish/0.1 ha), and high (>300 fish/0.1 ha) densities at 9, 11, and 10 sites, respectively. The eDNA results correctly predicted the presence and confirmed the absence of Brook Trout at 85.0–92.5% of the study sites; eDNA also explained 44% of the variability in Brook Trout population density and 24% of the variability in biomass. These findings indicate that eDNA surveys will enable researchers to effectively characterize the presence and abundance of Brook Trout and other species’ populations in headwater streams across the Adirondack region and elsewhere.

New York↗

Projecting the land cover change and its environmental impacts in the Cedar River Basin in the Midwestern United States

The physical surface of the Earth is in constant change due to climate forcing and human activities. In the Midwestern United States, urban area, farmland, and dedicated energy crop (e.g., switchgrass) cultivation are predicted to expand in the coming decades, which will lead to changes in hydrological processes. This study is designed to (1) project the land use and land cover (LULC) by mid-century using the FORecasting SCEnarios of future land-use (FORE-SCE) model under the A1B greenhouse gas emission scenario (future condition) and (2) assess its potential impacts on the water cycle and water quality against the 2001 baseline condition in the Cedar River Basin using the physically based soil and water assessment tool (SWAT). We compared the baseline LULC (National Land Cover data 2001) and 2050 projection, indicating substantial expansions of urban area and pastureland (including the cultivation of bioenergy crops) and a decrease in rangeland. We then used the above two LULC maps as the input data to drive the SWAT model, keeping other input data (e.g., climate) unchanged to isolate the LULC change impacts. The modeling results indicate that quick-response surface runoff would increase significantly (about 10.5%) due to the projected urban expansion (i.e., increase in impervious areas), and the baseflow would decrease substantially (about 7.3%) because of the reduced infiltration. Although the net effect may cause an increase in water yield, the increased variability may impede its use for public supply. Additionally, the cultivation of bioenergy crops such as switchgrass in the newly added pasture lands may further reduce the soil water content and lead to an increase in nitrogen loading (about 2.5% increase) due to intensified fertilizer application. These study results will be informative to decision makers for sustainable water resource management when facing LULC change and an increasing demand for biofuel production in this area.

Environmental Research Letters↗

Fire and tree death: Understanding and improving modeling of fire-induced tree mortality

Each year wildland fires kill and injure trees on millions of forested hectares globally, affecting plant and animal biodiversity, carbon storage, hydrologic processes, and ecosystem services. The underlying mechanisms of fire-caused tree mortality remain poorly understood, however, limiting the ability to accurately predict mortality and develop robust modeling applications, especially under novel future climates. Virtually all post-fire tree mortality prediction systems are based on the same underlying empirical model described in Ryan and Reinhardt (1988 Can. J. For. Res. 18 1291–7), which was developed from a limited number of species, stretching model assumptions beyond intended limits. We review the current understanding of the mechanisms of fire-induced tree mortality, provide recommended standardized terminology, describe model applications and limitations, and conclude with key knowledge gaps and future directions for research. We suggest a two-pronged approach to future research: (1) continued improvements and evaluations of empirical models to quantify uncertainty and incorporate new regions and species and (2) acceleration of basic, physiological research on the proximate and ultimate causes of fire-induced tree mortality to incorporate processes of tree death into models. Advances in both empirical and process fire-induced tree modeling will allow creation of hybrid models that could advance understanding of how fire injures and kills trees, while improving prediction accuracy of fire-driven feedbacks on ecosystems and landscapes, particularly under novel future conditions.

Environmental Research Letters↗

Development of perennial thaw zones in boreal hillslopes enhances potential mobilization of permafrost carbon

Permafrost thaw alters subsurface flow in boreal regions that in turn influences the magnitude, seasonality, and chemical composition of streamflow. Prediction of these changes is challenged by incomplete knowledge of timing, flowpath depth, and amount of groundwater discharge to streams in response to thaw. One important phenomenon that may affect flow and transport through boreal hillslopes is development of lateral perennial thaw zones (PTZs), the existence of which is here supported by geophysical observations and cryohydrogeologic modeling. Model results link thaw to enhanced and seasonally-extended baseflow, which have implications for mobilization of soluble constituents. Results demonstrate the sensitivity of PTZ development to organic layer thickness and near-surface factors that mediate heat exchange at the atmosphere/ground-surface interface. Study findings suggest that PTZs serve as a detectable precursor to accelerated permafrost degradation. This study provides important contextual insight on a fundamental thermo-hydrologic process that can enhance terrestrial-to-aquatic transfer of permafrost carbon, nitrogen, and mercury previously sequestered in thawing watersheds.

Alaska↗

Efficient delineation of nested depression hierarchy in digital elevation models for hydrological analysis using level-set method

In terrain analysis and hydrological modeling, surface depressions (or sinks) in a digital elevation model (DEM) are commonly treated as artifacts and thus filled and removed to create a depressionless DEM. Various algorithms have been developed to identify and fill depressions in DEMs during the past decades. However, few studies have attempted to delineate and quantify the nested hierarchy of actual depressions, which can provide crucial information for characterizing surface hydrologic connectivity and simulating the fill‐merge‐spill hydrological process. In this paper, we present an innovative and efficient algorithm for delineating and quantifying nested depressions in DEMs using the level‐set method based on graph theory. The proposed level‐set method emulates water level decreasing from the spill point along the depression boundary to the lowest point at the bottom of a depression. By tracing the dynamic topological changes (i.e., depression splitting/merging) within a compound depression, the level‐set method can construct topological graphs and derive geometric properties of the nested depressions. The experimental results of two fine‐resolution Light Detection and Ranging‐derived DEMs show that the raster‐based level‐set algorithm is much more efficient (~150 times faster) than the vector‐based contour tree method. The proposed level‐set algorithm has great potential for being applied to large‐scale ecohydrological analysis and watershed modeling.

Journal of the American Water Resources Associatio↗

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗

Changes in dissolved organic material determine exposure of stream benthic communities to UV-B radiation and heavy metals: Implications for climate change

Changes in regional climate in the Rocky Mountains over the next 100 years are expected to have significant effects on biogeochemical cycles and hydrological processes. In particular, decreased discharge and lower stream depth during summer when ultraviolet radiation (UVR) is the highest combined with greater photo-oxidation of dissolved organic materials (DOM) will significantly increase exposure of benthic communities to UVR. Communities in many Rocky Mountain streams are simultaneously exposed to elevated metals from abandoned mines, the toxicity and bioavailability of which are also determined by DOM. We integrated field surveys of 19 streams (21 sites) along a gradient of metal contamination with microcosm and field experiments conducted in Colorado, USA, and New Zealand to investigate the influence of DOM on bioavailability of heavy metals and exposure of benthic communities to UVR. Spatial and seasonal variation in DOM were closely related to stream discharge and significantly influenced heavy metal uptake in benthic organisms. Qualitative and quantitative changes in DOM resulting from exposure to sunlight increased UV-B (290-320nm) penetration and toxicity of heavy metals. Results of microcosm experiments showed that benthic communities from a metal-polluted stream were tolerant of metals, but were more sensitive to UV-B than communities from a reference stream. We speculate that the greater sensitivity of these communities to UV-B resulted from costs associated with metal tolerance. Exclusion of UVR from 12 separate Colorado streams and from outdoor stream microcosms in New Zealand increased the abundance of benthic organisms (mayflies, stoneflies, and caddisflies) by 18% and 54%, respectively. Our findings demonstrate the importance of considering changes in regional climate and UV-B exposure when assessing the effects of local anthropogenic stressors. ?? Journal compilation ?? 2008 Blackwell Publishing.

Global Change Biology↗

Effects of linking a soil-water-balance model with a groundwater-flow model

A previously published regional groundwater-flow model in north-central Nebraska was sequentially linked with the recently developed soil-water-balance (SWB) model to analyze effects to groundwater-flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater-level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root-mean-squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB-generated recharge, the RMS difference between simulated and estimated base-flow target values for the groundwater-flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater-level and base-flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.

Nebraska↗

Neogene geomorphic and climatic evolution of the central San Juan Mountains, Colorado: K/Ar age and stable isotope data on supergene alunite and jarosite from the Creede mining district

K/Ar age determinations or supergene alunite and jarosite, formed during Neogene weathering of the epithermal silver and base-metal ores of the Creede mining district, have been combined with geologic evidence to estimate the timing of regional uplift of the southern Rocky Mountains and related canyon cutting. In addition, oxygen and hydrogen isotopic studies suggest climate changes in the central San Juan Mountains during the past 5 m.y. Alunite [ideally (K,Na)Al 3 (SO 4 ) 2 (OH) 6 ] and jarosite [ideally KFe 3 (SO 4 ) 2 (OH) 6 ] can be dated by K/Ar or 40 Ar/ 39 Ar techniques and both contain OH and SO 4 sites that enable four stable isotope analyses (&delta;D, &delta; 18 O OH , and &delta; 34 S) to be made. This supergene alunite and jarosite formed by weathering of sulfide-rich ore bodies may record the evolution of the chemical and hydrologic processes affecting ancient oxidized acid ground water, as well as details of climate history and geomorphic evolution. Fine-grained (1-10 &mu;m) supergene alunite and jarosite occur in minor fractures in the upper, oxidized parts of the 25 Ma sulfide-bearing veins of the Creede mining district, and jarosite also occurs in adjacent oxidized Ag-bearing clastic sediments. K/Ar ages for alunite range from 4.8 to 3.1 Ma, and for jarosite range from 2.6 to 0.9 Ma. The &delta;D values for alunite and jarosite show opposite correlations with elevation, and values for jarosite correlate with age. Calculated &delta;D H2O values of alunite fluids approach but are larger than those of present-day meteoric water. Calculated &delta;D H2O values for jarosite fluids are more variable; the values of the youngest jarosites are lowest and are similar to those of present-day meteoric water in the district. The narrow &delta;D-&delta; 18 O SO4 values of alunites reflects oxidation of sulfide below the water table. The greater range in these values for jarosites reflects oxidation of sulfide under vadose conditions. The ages of alunite mark the position of the paleo-water table at the end of a period of moderate erosion from ca. 25 to 5 Ma that exposed the tops of the ore bodies to oxidation. The younger jarosite formed in the vadose zone during or following subsequent canyon cutting related to regional uplift of the southern Rocky Mountains, The &delta;D values suggest that climates in the area were similar to those of the present day prior to regional uplift but went through a warm period before returning to present conditions during or after regional uplift. The results of this study indicate that the combined stable and radiogenic isotope analysis of supergene alunite and jarosite has broad application in understanding climate and geomorphic evolution of selected areas.

Colorado↗

Annual cycle of magmatic CO2 in a tree-kill soil at Mammoth Mountain, California: Implications for soil acidification

Time-series sensor data reveal significant short-term and seasonal variations of magmatic CO 2 in soil over a 12 month period in 1995–1996 at the largest tree-kill site on Mammoth Mountain, central-eastern California. Short-term variations leading to ground-level soil CO 2 concentrations hazardous and lethal to humans were triggered by shallow faulting in the absence of increased seismicity or intrusion, consistent with tapping a reservoir of accumulated CO 2 , rather than direct magma degassing. Hydrologic processes closely modulated seasonal variations in CO 2 concentrations, which rose to 65%–100% in soil gas under winter snowpack and plunged more than 25% in just days as the CO 2 dissolved in spring snowmelt. The high efflux of CO 2 through the tree-kill soils acts as an open-system CO 2 buffer causing infiltration of waters with pH values commonly of <4.2, acid loading of up to 7 keq H+ ṁha −1 ṁyr −1 , mobilization of toxic Al 3+ , and long-term decline of soil fertility.

Geology↗

Hydrology of the unsaturated zone, Yucca Mountain, Nevada

The unsaturated zone at Yucca Mountain was investigated as a possible site for the nation's first high-level nuclear waste repository. Scientific investigations included infiltration studies, matrix properties testing, borehole testing and monitoring, underground excavation and testing, and the development of conceptual and numerical models of the hydrologic processes at Yucca Mountain. Infiltration estimates by empirical and geochemical methods range from 0.2 to 1.4 mm/yr and 0.2–6.0 mm/yr, respectively. Infiltration estimates from numerical models range from 4.5 mm/yr to 17.6 mm/yr. Rock matrix properties vary vertically and laterally as the result of depositional processes and subsequent postdepositional alteration. Laboratory tests indicate that the average matrix porosity and hydraulic conductivity values for the main level of the proposed repository (Topopah Spring Tuff middle nonlithophysal zone) are 0.08 and 4.7 × 10 −12 m/s, respectively. In situ fracture hydraulic conductivity values are 3–6 orders of magnitude greater. The permeability of fault zones is approximately an order of magnitude greater than that of the surrounding rock unit. Water samples from the fault zones have tritium concentrations that indicate some component of postnuclear testing. Gas and water vapor movement through the unsaturated zone is driven by changes in barometric pressure, temperature-induced density differences, and wind effects. The subsurface pressure response to surface barometric changes is controlled by the distribution and interconnectedness of fractures, the presence of faults and their ability to conduct gas and vapor, and the moisture content and matrix permeability of the rock units. In situ water potential values are generally less than −0.2 MPa (−2 bar), and the water potential gradients in the Topopah Spring Tuff units are very small. Perched-water zones at Yucca Mountain are associated with the basal vitrophyre of the Topopah Spring Tuff or the Calico Hills bedded tuff. Thermal gradients in the unsaturated zone vary with location, and range from ~2.0 °C to 6.0 °C per 100 m; the variability appears to be associated with topography. Large-scale heater testing identified a heat-pipe signature at ~97 °C, and identified thermally induced and excavation-induced changes in the stress field. Elevated gas-phase CO 2 concentrations and a decrease in the pH of water from the condensation zone also were identified. Conceptual and numerical flow and transport models of Yucca Mountain indicate that infiltration is highly variable, both spatially and temporally. Flow in the unsaturated zone is predominately through fractures in the welded units of the Tiva Canyon and Topopah Spring Tuffs and predominately through the matrix in the Paintbrush Tuff nonwelded units and Calico Hills Formation. Isolated, transient, fast-flow paths, such as faults, do exist but probably carry only a small portion of the total liquid-water flux at Yucca Mountain. The Paintbrush Tuff nonwelded units act as a storage buffer for transient infiltration pulses. Faults may act as flow boundaries and/or fast pathways. Below the proposed repository horizon, low-permeability lithostratigraphic units of the Topopah Spring Tuff and/or the Calico Hills Formation may divert flow laterally to faults that act as conduits to the water table. Advective transport pathways are consistent with flow pathways. Matrix diffusion is the major mechanism for mass transfer between fractures and the matrix and may contribute to retardation of radionuclide transport when fracture flow is dominant. Sorption may retard the movement of radionuclides in the unsaturated zone; however, sorption on mobile colloids may enhance radionuclide transport. Dispersion is not expected to be a major transport mechanism in the unsaturated zone at Yucca Mountain. Natural analogue studies support the concepts that percolating water may be diverted around underground openings and that the percentage of infiltration that becomes seepage decreases as infiltration decreases.

Nevada↗