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 37 records · Page 2Linked to original sources
USGS2008

Simulating water, solute, and heat transport in the subsurface with the VS2DI software package

The software package VS2DI was developed by the U.S. Geological Survey for simulating water, solute, and heat transport in variably saturated porous media. The package consists of a graphical preprocessor to facilitate construction of a simulation, a postprocessor for visualizing simulation results, and two numerical models that solve for flow and solute transport (VS2DT) and flow and heat transport (VS2DH). The finite-difference method is used to solve the Richards equation for flow and the advection-dispersion equation for solute or heat transport. This study presents a brief description of the VS2DI package, an overview of the various types of problems that have been addressed with the package, and an analysis of the advantages and limitations of the package. A review of other models and modeling approaches for studying water, solute, and heat transport also is provided. ?? Soil Science Society of America. All rights reserved.

Vadose Zone Journal↗
USGS2008

Susceptibility to enhanced chemical migration from depression-focused preferential flow, High Plains aquifer

Aquifer susceptibility to contamination is controlled in part by the inherent hydrogeologic properties of the vadose zone, which includes preferential-flow pathways. The purpose of this study was to investigate the importance of seasonal ponding near leaky irrigation wells as a mechanism for depression-focused preferential flow and enhanced chemical migration through the vadose zone of the High Plains aquifer. Such a mechanism may help explain the widespread presence of agrichemicals in recently recharged groundwater despite estimates of advective chemical transit times through the vadose zone from diffuse recharge that exceed the historical period of agriculture. Using a combination of field observations, vadose zone flow and transport simulations, and probabilistic neural network modeling, we demonstrated that vadose zone transit times near irrigation wells range from 7 to 50 yr, which are one to two orders of magnitude faster than previous estimates based on diffuse recharge. These findings support the concept of fast and slow transport zones and help to explain the previous discordant findings of long vadose zone transit times and the presence of agrichemicals at the water table. Using predictions of aquifer susceptibility from probabilistic neural network models, we delineated approximately 20% of the areal extent of the aquifer to have conditions that may promote advective chemical transit times to the water table of <50 yr if seasonal ponding and depression-focused flow exist. This aquifer-susceptibility map may help managers prioritize areas for groundwater monitoring or implementation of best management practices.

High Plains↗
USGS2009

Rapid measurement of field-saturated hydraulic conductivity for areal characterization

To provide an improved methodology for characterizing the field‐saturated hydraulic conductivity ( K fs ) over broad areas with extreme spatial variability and ordinary limitations of time and resources, we developed and tested a simplified apparatus and procedure, correcting mathematically for the major deficiencies of the simplified implementation. The methodology includes use of a portable, falling‐head, small‐diameter (∼20 cm) single‐ring infiltrometer and an analytical formula for K fs that compensates both for nonconstant falling head and for the subsurface radial spreading that unavoidably occurs with small ring size. We applied this method to alluvial fan deposits varying in degree of pedogenic maturity in the arid Mojave National Preserve, California. The measurements are consistent with a more rigorous and time‐consuming K fs measurement method, produce the expected systematic trends in K fs when compared among soils of contrasting degrees of pedogenic development, and relate in expected ways to results of widely accepted methods.

Vadose Zone Journal↗
USGS2009

Transport of tritium contamination to the atmosphere in an arid environment

Soil–plant–atmosphere interactions strongly influence water movement in desert unsaturated zones, but little is known about how such interactions affect atmospheric release of subsurface water-borne contaminants. This 2-yr study, performed at the U.S. Geological Survey's Amargosa Desert Research Site in southern Nevada, quantified the magnitude and spatiotemporal variability of tritium (3H) transport from the shallow unsaturated zone to the atmosphere adjacent to a low-level radioactive waste (LLRW) facility. Tritium fluxes were calculated as the product of 3H concentrations in water vapor and respective evaporation and transpiration water-vapor fluxes. Quarterly measured 3H concentrations in soil water vapor and in leaf water of the dominant creosote-bush [ Larrea tridentat a (DC.) Coville] were spatially extrapolated and temporally interpolated to develop daily maps of contamination across the 0.76-km2 study area. Maximum plant and root-zone soil concentrations (4200 and 8700 Bq L−1, respectively) were measured 25 m from the LLRW facility boundary. Continuous evaporation was estimated using a Priestley–Taylor model and transpiration was computed as the difference between measured eddy-covariance evapotranspiration and estimated evaporation. The mean evaporation/transpiration ratio was 3:1. Tritium released from the study area ranged from 0.12 to 12 μg d−1 and totaled 1.5 mg (8.2 × 1010 Bq) over 2 yr. Tritium flux variability was driven spatially by proximity to 3H source areas and temporally by changes in 3H concentrations and in the partitioning between evaporation and transpiration. Evapotranspiration removed and limited penetration of precipitation beneath native vegetation and fostered upward movement and release of 3H from below the root zone.

Vadose Zone Journal↗
USGS2009

Hydrologic characterization of desert soils with varying degrees of pedogenesis: 1. field experiments evaluating plant-relevant soil water behavior

To assess the eff ect of pedogenesis on the soil moisture dynamics infl uencing the character and quality of ecological habitat, we conducted infi ltration and redistribution experiments on three alluvial deposits in the Mojave National Preserve: (i) recently deposited active wash sediments, (ii) a soil of early Holocene age, and (iii) a highly developed soil of late Pleistocene age. At each, we ponded water in a 1-m-diameter infi ltration ring for 2.3 h and monitored soil water content and matric pressure during and atier infi ltration, using probes and electrical resistivity imaging (ERI). Infi ltration and downward fl ow rates were greater in younger material, favoring deep-rooted species. Deep-rooted species tend to colonize the margins of washes, where they are unaff ected by sediment transport that inhibits colonization. The ERI results support important generalizations, for example that shallower than 0.5 m, infi ltrated water persists longer in highly developed soil, favoring shallow-rooted species. Soil moisture data for the two youngest soils suggested that saturation overshoot, which may have signifi cant but unexplored hydroecologic and pedogenic eff ects, occurred at the horizontally advancing weting front. Spatial heterogeneity of soil properties generally increased with pedogenic development. Evidence suggested that some early-stage developmental processes may promote uniformity; the intermediate- age soil appeared to have the least heterogeneity in terms of textural variation with depth, and also the least anisotropy. Lateral heterogeneity was pronounced in older soil, having a multitude of eff ects on the distribution and retention of soil water, and may facilitate certain water-conserving strategies of plants over what would be possible in a laterally homogeneous soil. ?? Soil Science Society of America.

California↗
USGS2009

Effect of grain-coating mineralogy on nitrate and sulfate storage in the unsaturated zone

Unsaturated-zone sediments and the chemistry of shallow groundwater underlying a small (∼8-km 2 ) watershed were studied to identify the mechanisms responsible for anion storage within the Miocene Bridgeton Formation and weathered Coastal Plain deposits in southern New Jersey. Lower unsaturated-zone sediments and shallow groundwater samples were collected and concentrations of selected ions (including NO 3 − and SO 4 2− ) from 11 locations were determined. Grain size, sorting, and color of the lower unsaturated-zone sediments were determined and the mineralogy of these grains and the composition of coatings were analyzed by petrographic examination, scanning electron microscopy and energy dispersive analysis of x-rays, and quantitative whole-rock x-ray diffraction. The sediment grains, largely quartz and chert (80–94% w/w), are coated with a very fine-grained (<20 μm), complex mixture of kaolinite, halloysite, goethite, and possibly gibbsite and lepidocrocite. The mineral coatings are present as an open fabric, resulting in a large surface area in contact with pore water. Significant correlations between the amount of goethite in the grain coatings and the concentration of sediment-bound SO 4 2− were observed, indicative of anion sorption. Other mineral–chemical relations indicate that negatively charged surfaces and competition with SO 4 2− results in exclusion of NO 3 − from inner sphere exchange sites. The observed NO 3 − storage may be a result of matrix forces within the grain coatings and outer sphere complexation. The results of this study indicate that the mineralogy of grain coatings can have demonstrable effects on the storage of NO 3 − and SO 4 2− in the unsaturated zone.

Vadose Zone Journal↗
USGS2009

Comment on "Evaluating interactions between groundwater and vadose zone using the HYDRUS-based flow package for MODFLOW" by Navin Kumar C. Twarakavi, Jirka Šimůnek and Sophia Seo

Twarakavi et al (2008) compared four packages that can be used to estimate recharge for regional-scale groundwater flow simulations using MODFLOW (Harbaugh, 2005). This comment is focused on the comparisons made between two of these packages, namely, UZF1 (Niswonger et al., 2006) and a derivative of HYDRUS referred to herein as HYDRUS (Seo et al., 2007). In their paper, Twarakavi et al. (2008) stated that HYDRUS more accurately simulates unsaturated flow processes and groundwater recharge as compared to UZF1. However, Twarakavi et al. (2008) did not address several important differences between these models that undermine the advantages of HYDRUS as compared to UZF1 for simulating recharge. These differences were not revealed by the comparisons presented by Twarakavi et al. because the test simulations used to compare the models were too simple

Vadose Zone Journal↗
USGS2011

Interacting vegetative and thermal contributions to water movement in desert soil

Thermally driven water-vapor flow can be an important component of total water movement in bare soil and in deep unsaturated zones, but this process is often neglected when considering the effects of soil–plant–atmosphere interactions on shallow water movement. The objectives of this study were to evaluate the coupled and separate effects of vegetative and thermal-gradient contributions to soil water movement in desert environments. The evaluation was done by comparing a series of simulations with and without vegetation and thermal forcing during a 4.7-yr period (May 2001–December 2005). For vegetated soil, evapotranspiration alone reduced root-zone (upper 1 m) moisture to a minimum value (25 mm) each year under both isothermal and nonisothermal conditions. Variations in the leaf area index altered the minimum storage values by up to 10 mm. For unvegetated isothermal and nonisothermal simulations, root-zone water storage nearly doubled during the simulation period and created a persistent driving force for downward liquid fluxes below the root zone (total net flux ~1 mm). Total soil water movement during the study period was dominated by thermally driven vapor fluxes. Thermally driven vapor flow and condensation supplemented moisture supplies to plant roots during the driest times of each year. The results show how nonisothermal flow is coupled with plant water uptake, potentially influencing ecohydrologic relations in desert environments.

Vadose Zone Journal↗
USGS2012

Modeling a thick unsaturated zone at San Gorgonio Pass, California: lessons learned after five years of artificial recharge

The information flow among the tasks of framework assessment, numerical modeling, model forecasting and hind casting, and system-performance monitoring is illustrated. Results provide an understanding of artificial recharge in high-altitude desert settings where large vertical distances may separate application ponds from their target aquifers. Approximately 3.8 million cubic meters of surface water was applied to spreading ponds from 2003–2007 to artificially recharge the underlying aquifer through a 200-meter thick unsaturated zone in the San Gorgonio Pass area in southern California. A study was conducted between 1997 and 2003, and a numerical model was developed to help determine the suitability of the site for artificial recharge. Ongoing monitoring results indicated that the existing model needed to be modified and recalibrated to more accurately predict artificial recharge at the site. The objective of this work was to recalibrate the model by using observation of the application rates, the rise and fall of the water level above a perching layer, and the approximate arrival time to the water table during the 5-yr monitoring period following initiation of long-term artificial recharge. Continuous monitoring of soil-matric potential, temperature, and water levels beneath the site indicated that artificial recharge reached the underlying water table between 3.75 and 4.5 yr after the initial application of the recharge water. The model was modified to allow the simulation to more adequately match the perching layer dynamics and the time of arrival at the water table. The instrumentation also showed that the lag time between changes in application of water at the surface and the response at the perching layer decreased from about 4 mo to less than 1 mo due to the wet-up of the unsaturated zone and the increase in relative permeability. The results of this study demonstrate the importance of iteratively monitoring and modeling the unsaturated zone in layered alluvial systems in the context of artificial recharge. They show that adequate geologic and hydraulic-property data on perching layers are critical to success. Continuous monitoring in the unsaturated and saturated zones beneath a site provides data to develop and constrain numerical models, better understand local unsaturated zone process, manage artificial recharge operations, and to determine the timing and volume of recoverable water for consumptive use.

California↗
USGS2012

Response to "Comments on 'Theory for source-responsive and free-surface film modeling of unsaturated flow'"

I am grateful to Masciopinto (2012) for raising several issues from my study (Nimmo, 2010) that deserve elaboration or clarification. In this reply, I address these in what I judge to be the order of importance, the main ones being (i) the discrepancy of scales between the two domains that treat preferential and diffuse unsaturated flow, and (ii) the properties that must be evaluated to characterize the medium.

Vadose Zone Journal↗
USGS2013

Wildfire and aspect effects on hydrologic states after the 2010 Fourmile Canyon Fire

Wildfire can change how soils take in, store, and release water. This study examined differences in how burned and unburned plots on north versus south-facing slope aspects respond to rainfall. The largest wildfire impacts were litter/duff combustion on burned north-facing slopes versus soil-water retention reduction on burned south-facing slopes. Wildfire is one of the most significant disturbances in mountainous landscapes, affecting water supply and ecologic function and setting the stage for natural hazards such as flash floods. The impacts of wildfire can affect the entire hydrologic cycle. Measurements of soil-water content and matric potential in the near surface (top 30 cm) captured the hydrologic state in both burned and unburned hillslopes during the first spring through fall period (1 June–1 Oct. 2011) after the 2010 Fourmile Canyon Fire near Boulder, CO. This time span included different hydrologic periods characterized by cyclonic frontal storms (low-intensity, long duration), convective storms (high-intensity, short duration), and dry periods. In mountainous environments, aspect can also control hydrologic states, so north- vs. south-facing slopes were compared. Wildfire tended to homogenize soil-water contents across aspects and with depth in the soil, yet it also may have introduced an aspect control on matric potential that was not observed in unburned soils. Post-wildfire changes in hydrologic state were observed in south-facing soils, probably reflecting decreased soil-water retention after wildfire. North-facing soils were impacted the most, in terms of hydrologic state, by the loss of water storage in the combusted litter–duff layer and forest canopy, which had provided a large “hydrologic buffering” capacity when unburned. Unsaturated zone measurements showed increased variability in hydrologic states and more rapid state transitions in wildfire-impacted soils. A simple, qualitative analysis suggested that the range of unsaturated-zone processes along the gravity–capillarity–adsorption continuum was expanded by wildfire for a given soil. The small number of experimental plots in this study suggests that further work is needed before these conclusions can be generalized to other geographic areas.

Vadose Zone Journal↗
USGS2013

Predicting vertically-nonsequential wetting patterns with a source-responsive model

Water infiltrating into soil of natural structure often causes wetting patterns that do not develop in an orderly sequence. Because traditional unsaturated flow models represent a water advance that proceeds sequentially, they fail to predict irregular development of water distribution. In the source-responsive model, a diffuse domain (D) represents flow within soil matrix material following traditional formulations, and a source-responsive domain (S), characterized in terms of the capacity for preferential flow and its degree of activation, represents preferential flow as it responds to changing water-source conditions. In this paper we assume water undergoing rapid source-responsive transport at any particular time is of negligibly small volume; it becomes sensible at the time and depth where domain transfer occurs. A first-order transfer term represents abstraction from the S to the D domain which renders the water sensible. In tests with lab and field data, for some cases the model shows good quantitative agreement, and in all cases it captures the characteristic patterns of wetting that proceed nonsequentially in the vertical direction. In these tests we determined the values of the essential characterizing functions by inverse modeling. These functions relate directly to observable soil characteristics, rendering them amenable to evaluation and improvement through hydropedologic development.

Vadose Zone Journal↗
USGS2013

Tritium plume dynamics in the shallow unsaturated zone in an arid environment

The spatiotemporal variability of a tritium plume in the shallow unsaturated zone and the mechanisms controlling its transport were evaluated during a 10-yr study. Plume movement was minimal and its mass declined by 68%. Upward-directed diffusive-vapor tritium fluxes and radioactive decay accounted for most of the observed plume-mass declines. Effective isolation of tritium ( 3 H) and other contaminants at waste-burial facilities requires improved understanding of transport processes and pathways. Previous studies documented an anomalously widespread (i.e., theoretically unexpected) distribution of 3 H (>400 m from burial trenches) in a dry, sub-root-zone gravelly layer (1–2-m depth) adjacent to a low-level radioactive waste (LLRW) burial facility in the Amargosa Desert, Nevada, that closed in 1992. The objectives of this study were to: (i) characterize long-term, spatiotemporal variability of 3 H plumes; and (ii) quantify the processes controlling 3 H behavior in the sub-root-zone gravelly layer beneath native vegetation adjacent to the facility. Geostatistical methods, spatial moment analyses, and mass flux calculations were applied to a spatiotemporally comprehensive, 10-yr data set (2001–2011). Results showed minimal bulk-plume advancement during the study period and limited Fickian spreading of mass. Observed spreading rates were generally consistent with theoretical vapor-phase dispersion. The plume mass diminished more rapidly than would be expected from radioactive decay alone, indicating net efflux from the plume. Estimates of upward 3 H efflux via diffusive-vapor movement were >10× greater than by dispersive-vapor or total-liquid movement. Total vertical fluxes were >20× greater than lateral diffusive-vapor fluxes, highlighting the importance of upward migration toward the land surface. Mass-balance calculations showed that radioactive decay and upward diffusive-vapor fluxes contributed the majority of plume loss. Results indicate that plume losses substantially exceeded any continuing 3 H contribution to the plume from the LLRW facility during 2001 to 2011 and suggest that the widespread 3 H distribution resulted from transport before 2001.

Nevada↗
USGS2014

A screening tool for delineating subregions of steady recharge within groundwater models

We have developed a screening method for simplifying groundwater models by delineating areas within the domain that can be represented using steady-state groundwater recharge. The screening method is based on an analytical solution for the damping of sinusoidal infiltration variations in homogeneous soils in the vadose zone. The damping depth is defined as the depth at which the flux variation damps to 5% of the variation at the land surface. Groundwater recharge may be considered steady where the damping depth is above the depth of the water table. The analytical solution approximates the vadose zone diffusivity as constant, and we evaluated when this approximation is reasonable. We evaluated the analytical solution through comparison of the damping depth computed by the analytic solution with the damping depth simulated by a numerical model that allows variable diffusivity. This comparison showed that the screening method conservatively identifies areas of steady recharge and is more accurate when water content and diffusivity are nearly constant. Nomograms of the damping factor (the ratio of the flux amplitude at any depth to the amplitude at the land surface) and the damping depth were constructed for clay and sand for periodic variations between 1 and 365 d and flux means and amplitudes from nearly 0 to 1 × 10 −3 m d −1 . We applied the screening tool to Central Valley, California, to identify areas of steady recharge. A MATLAB script was developed to compute the damping factor for any soil and any sinusoidal flux variation.

California↗
USGS2014

Field-scale sulfur hexafluoride tracer experiment to understand long distance gas transport in the deep unsaturated zone

A gas-tracer test in a deep arid unsaturated zone demonstrates that standard estimates of effective diffusivity from sediment properties allow a reasonable first-cut assessment of gas contaminant transport. Apparent anomalies in historic transport behavior at this and other waste disposal sites may result from factors other than nonreactive gas transport properties. A natural gradient SF 6 tracer experiment provided an unprecedented evaluation of long distance gas transport in the deep unsaturated zone (UZ) under controlled (known) conditions. The field-scale gas tracer test in the 110-m-thick UZ was conducted at the U.S. Geological Survey&rsquo;s Amargosa Desert Research Site (ADRS) in southwestern Nevada. A history of anomalous (theoretically unexpected) contaminant gas transport observed at the ADRS, next to the first commercial low-level radioactive waste disposal facility in the United States, provided motivation for the SF 6 tracer study. Tracer was injected into a deep UZ borehole at depths of 15 and 48 m, and plume migration was observed in a monitoring borehole 9 m away at various depths (0.5&ndash;109 m) over the course of 1 yr. Tracer results yielded useful information about gas transport as applicable to the spatial scales of interest for off-site contaminant transport in arid unsaturated zones. Modeling gas diffusion with standard empirical expressions reasonably explained SF 6 plume migration, but tended to underpredict peak concentrations for the field-scale experiment given previously determined porosity information. Despite some discrepancies between observations and model results, rapid SF 6 gas transport commensurate with previous contaminant migration was not observed. The results provide ancillary support for the concept that apparent anomalies in historic transport behavior at the ADRS are the result of factors other than nonreactive gas transport properties or processes currently in effect in the undisturbed UZ.

California, Nevada↗
USGS2016

Quantitative framework for preferential flow initiation and partitioning

A model for preferential flow in macropores is based on the short-range spatial distribution of soil matrix infiltrability. It uses elementary areas at two different scales. One is the traditional representative elementary area (REA), which includes a sufficient heterogeneity to typify larger areas, as for measuring field-scale infiltrability. The other, called an elementary matrix area (EMA), is smaller, but large enough to represent the local infiltrability of soil matrix material, between macropores. When water is applied to the land surface, each EMA absorbs water up to the rate of its matrix infiltrability. Excess water flows into a macropore, becoming preferential flow. The land surface then can be represented by a mesoscale (EMA-scale) distribution of matrix infiltrabilities. Total preferential flow at a given depth is the sum of contributions from all EMAs. Applying the model, one case study with multi-year field measurements of both preferential and diffuse fluxes at a specific depth was used to obtain parameter values by inverse calculation. The results quantify the preferential&ndash;diffuse partition of flow from individual storms that differed in rainfall amount, intensity, antecedent soil water, and other factors. Another case study provided measured values of matrix infiltrability to estimate parameter values for comparison and illustrative predictions. These examples give a self-consistent picture from the combination of parameter values, directions of sensitivities, and magnitudes of differences caused by different variables. One major practical use of this model is to calculate the dependence of preferential flow on climate-related factors, such as varying soil wetness and rainfall intensity.

Vadose Zone Journal↗
USGS2016

Experimental studies and model analysis of noble gas fractionation in porous media

The noble gases, which are chemically inert under normal terrestrial conditions but vary systematically across a wide range of atomic mass and diffusivity, offer a multicomponent approach to investigating gas dynamics in unsaturated soil horizons, including transfer of gas between saturated zones, unsaturated zones, and the atmosphere. To evaluate the degree to which fractionation of noble gases in the presence of an advective&ndash;diffusive flux agrees with existing theory, a simple laboratory sand column experiment was conducted. Pure CO2 was injected at the base of the column, providing a series of constant CO2 fluxes through the column. At five fixed sampling depths within the system, samples were collected for CO2 and noble gas analyses, and ambient pressures were measured. Both the advection&ndash;diffusion and dusty gas models were used to simulate the behavior of CO2 and noble gases under the experimental conditions, and the simulations were compared with the measured depth-dependent concentration profiles of the gases. Given the relatively high permeability of the sand column (5 &acute; 10&minus;11 m2), Knudsen diffusion terms were small, and both the dusty gas model and the advection&ndash;diffusion model accurately predicted the concentration profiles of the CO2 and atmospheric noble gases across a range of CO2 flux from ?700 to 10,000 g m&minus;2 d&minus;1. The agreement between predicted and measured gas concentrations demonstrated that, when applied to natural systems, the multi-component capability provided by the noble gases can be exploited to constrain component and total gas fluxes of non-conserved (CO2) and conserved (noble gas) species or attributes of the soil column relevant to gas transport, such as porosity, tortuosity, and gas saturation.

Vadose Zone Journal↗