James Robert Brainard

A geostatistically based inverse model for electrical resistivity surveys and its applications to vadose zone hydrology

[1] A sequential, geostatistical inverse approach was developed for electrical resistivity tomography (ERT). Unlike most ERT inverse approaches, this new approach allows inclusion of our prior knowledge of general geological structures of an area and point electrical resistivity measurements to constrain the estimate of the electrical resistivity field. This approach also permits sequential inclusion of different data sets, mimicking the ERT data collection scheme commonly employed in the field survey. Furthermore, using the conditional variance concept, the inverse model quantifies uncertainty of the estimate caused by spatial variability and measurement errors. Using this approach, numerical experiments were conducted to demonstrate the effects of bedding orientation on ERT surveys and to show both the usefulness and uncertainty associated with the inverse approach for delineating the electrical resistivity distribution using down-hole ERT arrays. A statistical analysis was subsequently undertaken to explore the effects of spatial variability of the electrical resistivity-moisture relation on the interpretation of the change in water content in the vadose zone, using the change in electrical resistivity. Core samples were collected from a field site to investigate the spatial variability of the electrical resistivity-moisture relation. Numerical experiments were subsequently conducted to illustrate how the spatially varying relations affect the level of uncertainty in the interpretation of change of moisture content based on the estimated change in electrical resistivity. Other possible complications are also discussed. INDEX TERMS: 0903 Exploration Geophysics: Computational methods, potential fields; 1869 Hydrology: Stochastic processes; 1875 Hydrology: Unsaturated zone; 1866 Hydrology: Soil moisture; 3260 Mathematical Geophysics: Inverse theory; KEYWORDS: geostatistical inverse model, electrical resistivity tomography, vadose zone, resistivitymoisture relation, spatial variability, sequential/successive linear estimator

Infiltration in Unsaturated Layered Fluvial Deposits at Rio Bravo: Macroscopic Anisotropy and Heterogeneous Transport

An infiltration and dye transport experiment was conducted to visualize flow and transport processes in a heterogeneous, layered, sandy-gravelly fluvial deposit adjacent to Rio Bravo Boulevard in Albuquerque, NM. Water containing red dye followed by blue-green dye was ponded in a small horizontal zone (about 0.5 by 0.5 m) above a vertical outcrop (about 4 by 2.5 m). The red dye lagged behind the wetting front due to slight adsorption, thus allowing both the wetting front and dye fronts to be observed in time at the outcrop face. After infiltration, vertical slices were excavated to the midpoint of the infiltrometer, exposing the wetting front and dye distribution in a quasi three-dimensional manner. At small scale, wetting front advancement was influenced by the multitude of local capillary barriers within the deposit. However, at the scale of the experiment, the wetting front appeared smooth with significant lateral spreading, twice that in the vertical, indicating a strong anisotropy due to the pronounced horizontal layering. The dye fronts exhibited appreciably more irregularity than the wetting front, as well as the influence of preferential flow features (a fracture) that moved the dye directly to the front, bypassing the fresh water between. To illustrate the ability of equivalent homogeneous media models to capture the behavior of the wetting front, we performed numerical simulations using equivalent homogeneous media with isotropic, anisotropic, and moisture-dependent anisotropic properties. Those containing anisotropy matched the experimental data best.

A Hydrologic-Geophysical Method for Characterizing Flow and Transport Processes within the Vadose Zone

The objective of this study is to analyze flow within the vadose zone during a mid-scale hydrologic test to and to characterize transport processes in-situ. This project will employ numerical and experimental tools that have been developed under a previously funded EMSP proposal (project number 55332). Geophysical imaging techniques will be employed to image the changes produced by the transport experiments in-situ as they occur. Results will help to better understand flow and transport modes within the vadose zone at DOE sites, including the influence of natural heterogeneities and man-made structures. In addition the data will provide checks against which numerical flow and transport simulations can be compared.

Data collection for cooperative water resources modeling in the Lower Rio Grande Basin, Fort Quitman to the Gulf of Mexico.

Water resource scarcity around the world is driving the need for the development of simulation models that can assist in water resources management. Transboundary water resources are receiving special attention because of the potential for conflict over scarce shared water resources. The Rio Grande/Rio Bravo along the U.S./Mexican border is an example of a scarce, transboundary water resource over which conflict has already begun. The data collection and modeling effort described in this report aims at developing methods for international collaboration, data collection, data integration and modeling for simulating geographically large and diverse international watersheds, with a special focus on the Rio Grande/Rio Bravo. This report describes the basin, and the data collected. This data collection effort was spatially aggregated across five reaches consisting of Fort Quitman to Presidio, the Rio Conchos, Presidio to Amistad Dam, Amistad Dam to Falcon Dam, and Falcon Dam to the Gulf of Mexico. This report represents a nine-month effort made in FY04, during which time the model was not completed.

Infiltration in Unsaturated Layered Fluvial Deposits at Rio Bravo: Photo Essay and Data Summary

An infiltration and dye transport experiment was conducted to visualize flow and transport processes in a heterogeneous, layered, sandy-gravelly fluvial deposit adjacent to Rio Bravo Boulevard in Albuquerque, NM. Water containing red dye followed by blue-green dye was ponded in a small horizontal zone ({approx}0.5 m x 0.5 m) above a vertical outcrop ({approx}4 m x 2.5 m). The red dye lagged behind the wetting front due to slight adsorption thus allowing both the wetting front and dye fronts to be observed in time at the outcrop face. After infiltration, vertical slices were excavated to the midpoint of the infiltrometer exposing the wetting front and dye distribution in a quasi three-dimensional manner. At small-scale, wetting front advancement was influenced by the multitude of local capillary barriers within the deposit. However at the scale of the experiment, the wetting front appeared smooth with significant lateral spreading {approx} twice that in the vertical, indicating a strong anisotropy due to the pronounced horizontal layering. The dye fronts exhibited appreciably more irregularity than the wetting front, as well as the influence of preferential flow features (a fracture) that moved the dye directly to the front, bypassing the fresh water between.

Vadose Zone Monitoring of Dairy Green Water Lagoons using Soil Solution Samplers.

Over the last decade, dairy farms in New Mexico have become an important component to the economy of many rural ranching and farming communities. Dairy operations are water intensive and use groundwater that otherwise would be used for irrigation purposes. Most dairies reuse their process/green water three times and utilize lined lagoons for temporary storage of green water. Leakage of water from lagoons can pose a risk to groundwater quality. Groundwater resource protection infrastructures at dairies are regulated by the New Mexico Environment Department which currently relies on monitoring wells installed in the saturated zone for detecting leakage of waste water lagoon liners. Here we present a proposal to monitor the unsaturated zone beneath the lagoons with soil water solution samplers to provide early detection of leaking liners. Early detection of leaking liners along with rapid repair can minimize contamination of aquifers and reduce dairy liability for aquifer remediation. Additionally, acceptance of vadose zone monitoring as a NMED requirement over saturated zone monitoring would very likely significantly reduce dairy startup and expansion costs. Acknowledgment Funding for this project was provided by the Sandia National Laboratories Small Business Assistance Program

A Hydrologic-geophysical Method for Characterizing Flow and Transport Processes Within The Vadose Zone

The primary purpose of this project was to employ two geophysical imaging techniques, electrical resistivity tomography and cross-borehole ground penetrating radar, to image a controlled infiltration of a saline tracer under unsaturated flow conditions. The geophysical techniques have been correlated to other more traditional hydrologic measurements including neutron moisture measurements and induction conductivity logs. Images that resulted during two successive infiltrations indicate the development of what appear to be preferential pathways through the finer grained materials, although the results could also be produced by cationic capture of free ions in clays. In addition the site as well as the developing solute plume exhibits electrical anisotropy which is likely related to flow properties. However the geologic significance of this phenomenon is still under investigation.The primary purpose of this project was to employ two geophysical imaging techniques, electrical resistivity tomography and cross-borehole ground penetrating radar, to image a controlled infiltration of a saline tracer under unsaturated flow conditions. The geophysical techniques have been correlated to other more traditional hydrologic measurements including neutron moisture measurements and induction conductivity logs. Images that resulted during two successive infiltrations indicate the development of what appear to be preferential pathways through the finer grained materials, although the results could also be produced by cationic capture of free ions in clays. In addition the site as well as the developing solute plume exhibits electrical anisotropy which is likely related to flow properties. However the geologic significance of this phenomenon is still under investigation.

Monitoring stream stage, channel profile, and aqueous conductivity with time domain reflectometry (TDR).

Time domain reflectometry (TDR) operates by propagating a radar frequency electromagnetic pulse down a transmission line while monitoring the reflected signal. As the electromagnetic pulse propagates along the transmission line, it is subject to impedance by the dielectric properties of the media along the transmission line (e.g., air, water, sediment), reflection at dielectric discontinuities (e.g., air-water or water-sediment interface), and attenuation by electrically conductive materials (e.g., salts, clays). Taken together, these characteristics provide a basis for integrated stream monitoring; specifically, concurrent measurement of stream stage, channel profile and aqueous conductivity. Here, we make novel application of TDR within the context of stream monitoring. Efforts toward this goal followed three critical phases. First, a means of extracting the desired stream parameters from measured TDR traces was required. Analysis was complicated by the fact that interface location and aqueous conductivity vary concurrently and multiple interfaces may be present at any time. For this reason a physically based multisection model employing the S11 scatter function and Cole-Cole parameters for dielectric dispersion and loss was developed to analyze acquired TDR traces. Second, we explored the capability of this multisection modeling approach for interpreting TDR data acquired from complex environments, such as encountered in stream monitoring. A series of laboratory tank experiments were performed in which the depth of water, depth of sediment, and conductivity were varied systematically. Comparisons between modeled and independently measured data indicate that TDR measurements can be made with an accuracy of {+-}3.4x10{sup -3} m for sensing the location of an air/water or water/sediment interface and {+-}7.4% of actual for the aqueous conductivity. Third, monitoring stations were sited on the Rio Grande and Paria rivers to evaluate performance of the TDR system under normal field conditions. At the Rio Grande site (near Central Bridge in Albuquerque, New Mexico) continuous monitoring of stream stage and aqueous conductivity was performed for 6 months. Additionally, channel profile measurements were acquired at 7 locations across the river. At the Paria site (near Lees Ferry, Arizona) stream stage and aqueous conductivity data were collected over a 4-month period. Comparisons drawn between our TDR measurements and USGS gage data indicate that the stream stage is accurate within {+-}0.88 cm, conductivity is accurate within {+-}11% of actual, and channel profile measurements agree within {+-}1.2 cm.

A RESOLUTION ANALYSIS OF TWO GEOPHYSICAL IMAGING METHODS FOR CHARACTERIZING AND MONITORING HYDROLOGIC CONDITIONS IN THE VADOSE ZONE

This project has been designed to analyze the resolution of two different geophysical imaging techniques (electrical resistivity tomography and cross-borehole ground penetrating radar) for monitoring subsurface flow and transport processes within the vadose zone. This is to be accomplished through a coupled approach involving large scale unsaturated flow modeling, petrophysical conversion of the resulting 2 and 3 Dimensional water content and solute concentration fields to geophysical property models and generation of synthetic geophysical data, followed by the inversion of the synthetic geophysical data. The resolution, strengths, and limitations of the geophysical techniques will then be ascertained through an analysis involving comparisons between the original hydrologic modeling results and inverted geophysical images. Increasing levels of complexity will be added to the models as the project progresses through the addition of heterogeneity in the original hydrologic property model, and by adding uncertainty to the petrophysical relationship that couples the geophysical model to the hydrologic modeling results.