3-D Transient Hydraulic Tomography in Unconfined Aquifers with Fast Drainage Response

We investigate, through numerical experiments, the viability of three-dimensional transient hydraulic tomography (3DTHT) for identifying the spatial distribution of groundwater flow parameters (primarily, hydraulic conductivity K) in permeable, unconfined aquifers. To invert the large amount of transient data collected from 3DTHT surveys, we utilize an iterative geostatistical inversion strategy in which outer iterations progressively increase the number of data points fitted and inner iterations solve the quasi-linear geostatistical formulas of Kitanidis. In order to base our numerical experiments around realistic scenarios, we utilize pumping rates, geometries, and test lengths similar to those attainable during 3DTHT field campaigns performed at the Boise Hydrogeophysical Research Site (BHRS). We also utilize hydrologic parameters that are similar to those observed at the BHRS and in other unconsolidated, unconfined fluvial aquifers. In addition to estimating K, we test the ability of 3DTHT to estimate both average storage values (specific storage Ss and specific yield Sy) as well as spatial variability in storage coefficients. The effects of model conceptualization errors during unconfined 3DTHT are investigated including: (1) assuming constant storage coefficients during inversion and (2) assuming stationary geostatistical parameter variability. Overall, our findings indicate that estimation of K is slightly degraded if storage parameters must be jointly estimated, but that this effect is quite small compared with the degradation of estimates due to violation of “structural” geostatistical assumptions. Practically, we find for our scenarios that assuming constant storage values during inversion does not appear to have a significant effect on K estimates or uncertainty bounds

A View Toward the Future of Subsurface Characterization: CAT Scanning Groundwater Basins

In this opinion paper we contend that high-resolution characterization, monitoring, and prediction are the key elements to advancing and reducing uncertainty in our understanding and prediction of subsurface processes at basin scales. First, we advocate that recently developed tomographic surveying is an effective and high-resolution approach for characterizing the field-scale subsurface. Fusion of different types of tomographic surveys further enhances the characterization. A basin is an appropriate scale for many water resources management purposes. We thereby propose the expansion of the tomographic surveying and data fusion concept to basin-scale characterization. In order to facilitate basin-scale tomographic surveys, different types of passive, basin-scale, CAT scan technologies are suggested that exploit recurrent natural stimuli (e.g., lightning, earthquakes, storm events, barometric variations, river-stage variations, etc.) as sources of excitations, along with implementation of sensor networks that provide long-term and spatially distributed monitoring of excitation as well as response signals on the land surface and in the subsurface. This vision for basin-scale subsurface characterization faces many significant technological challenges and requires interdisciplinary collaborations (e.g., surface and subsurface hydrology, geophysics, geology, geochemistry, information and sensor technology, applied mathematics, atmospheric science, etc.). We nevertheless contend that this should be a future direction for subsurface science research

Hierarchical Geostatistics and Multifacies Systems: Boise Hydrogeophysical Research Site, Boise, Idaho

The geostatistical structure of a heterogeneous coarse fluvial aquifer is investigated with porosity data derived from neutron logs at a research well field (Boise Hydrogeophysical Research Site, or BHRS) that was designed, in part, to support three-dimensional geostatistical analysis of hydrologic and geophysical parameters. Recognizing that the coarse fluvial deposits include subdivisions (units between bounding surfaces), we adopt a hierarchical approach and examine the porosity geostatistics of the aquifer at three scales. At the BHRS, the saturated fluvial deposits as a whole (maximum interwell spacing ~80 m, thickness ~16–18 m) are at hierarchical level 1; five subhorizontal units within these deposits (four cobble-dominated units and a channel sand) can be traced across the central area of the BHRS and are at hierarchical level 2; and subunits (patches or lenses) in one of the level 2 units (Unit 4), are at hierarchical level 3. We use variography and porosity statistics to recognize nonstationarity at hierarchical level 1 and in one of the level 2 units (Unit 4) where the means and variances of porosity differences as a function of lag are not constant between distinct units and subunits, respectively. The geostatistical structure at level 1 is modeled with different horizontal and vertical structures that have different sills (vertical sill greater than horizontal sill). The difference in sills can be explained quantitatively by the summing of weighted sills from all individual units and combined units (i.e., a given pair of different units), where the weights are the proportions of data pairs contributing to the sills at each lag from the individual and combined units. Extension of this analysis leads to a weighted, multistructure form of the variogram function whereby a global experimental variogram in a hierarchical system can be decomposed quantitatively into weighted component individual- and combined-unit (or facies) structures for any number of units or hierarchical levels. Such decomposition of the global horizontal variogram from the BHRS indicates that short-range periodicity in that structure is due to both (1) combined-unit structures associated with patches or lenses at hierarchical level 3 in Unit 4 and (2) variations in thickness of Unit 2. For hierarchical multifacies systems, structure models fit to global horizontal and vertical experimental variograms may not be useful for subsequent stochastic modeling if the system on which the structure models are based is nonstationary

Crosshole Radar Tomography in a Fluvial Aquifer Near Boise, Idaho

To determine the distribution of heterogeneities in the saturated zone of an unconfined aquifer in Boise, ID, we compute tomograms for three adjacent well pairs. The fluvial deposits consist of unconsolidated cobbles and sands. We used a curved-ray, finite-difference approximation to the eikonal equation to generate the forward model. The inversion uses a linearized, iterative scheme to determine the slowness distribution from the first arrival traveltimes. The tomograms consist of a layered zone representing the saturated aquifer. The velocities in this saturated zone range between 0.06 to 0.10 m/ns. We use a variety of methods to assess the reliability of our velocity models. Finally, we compare our results to neutron-derived porosity logs in the wells used for the tomograms. The comparison shows that the trends in porosity derived from the tomograms match the trends in porosity measured with the neutron probe

Hydraulic Conductivity Imaging from 3-D Transient Hydraulic Tomography at Several Pumping/Observation Densities

[1] 3-D Hydraulic tomography (3-D HT) is a method for aquifer characterization whereby the 3-D spatial distribution of aquifer flow parameters (primarily hydraulic conductivity, K) is estimated by joint inversion of head change data from multiple partially penetrating pumping tests. While performance of 3-D HT has been studied extensively in numerical experiments, few field studies have demonstrated the real-world performance of 3-D HT. Here we report on a 3-D transient hydraulic tomography (3-D THT) field experiment at the Boise Hydrogeophysical Research Site which is different from prior approaches in that it represents a “baseline” analysis of 3-D THT performance using only a single arrangement of a central pumping well and five observation wells with nearly complete pumping and observation coverage at 1 m intervals. We jointly analyze all pumping tests using a geostatistical approach based on the quasi-linear estimator of Kitanidis (1995). We reanalyze the system after progressively removing pumping and/or observation intervals; significant progressive loss of information about heterogeneity is quantified as reduced variance of the K field overall, reduced correlation with slug test K estimates at wells, and reduced ability to accurately predict independent pumping tests. We verify that imaging accuracy is strongly improved by pumping and observational densities comparable to the aquifer heterogeneity geostatistical correlation lengths. Discrepancies between K profiles at wells, as obtained from HT and slug tests, are greatest at the tops and bottoms of wells where HT observation coverage was lacking

A Field Proof-of-Concept of Aquifer Imaging Using 3-D Transient Hydraulic Tomography with Modular, Temporarily-Emplaced Equipment

Hydraulic tomography is a field scale aquifer characterization method capable of estimating 3-D heterogeneous parameter distributions, and is directly sensitive to hydraulic conductivity (K), thus providing a useful data source for improving flow and transport models. We present results from a proof-of-concept field and modeling study in which we apply 3-D transient hydraulic tomography (3DTHT) to the relatively high-K and moderately heterogeneous unconfined aquifer at the Boise Hydrogeophysical Research Site. Short-duration (20 min) partially penetrating pumping tests, for which observed responses do not reach steady state, are used as the aquifer stimulation. To collect field data, we utilize a system of temporarily emplaced packer equipment to isolate multiple discrete intervals in boreholes. To analyze the data, we utilize MODFLOW combined with geostatistical inversion code based on the quasilinear approach of Kitanidis (1995). This combination of practical software allows inversion of large datasets (\u3e250 drawdown curves, and almost 1000 individual data points) and estimation of K at \u3e100,000 locations; reasonable runtimes are obtained using a single multicore computer with 12 GB of RAM. The K heterogeneity results from 3DTHT are cross-validated against K characterization from a large set of partially penetrating slug tests, and found to be quite consistent. The use of portable, modular equipment for field implementation means that 3DTHT data collection can be performed (including mobilization/demobilization) within a matter of days. Likewise, use of a practical, efficient and scalable numerical modeling and inversion strategy means that computational effort is drastically reduced, such that 3-D aquifer property distributions can be estimated quickly

Estimating Porosity with Ground-Penetrating Radar Reflection Tomography: A Controlled 3-D Experiment at the Boise Hydrogeophysical Research Site

To evaluate the uncertainty of water-saturated sediment velocity and porosity estimates derived from surface-based, ground-penetrating radar reflection tomography, we conducted a controlled field experiment at the Boise Hydrogeophysical Research Site (BHRS). The BHRS is an experimental well field located near Boise, Idaho. The experimental data set consisted of 3-D multioffset radar acquired on an orthogonal 20 × 30 m surface grid that encompassed a set of 13 boreholes. Experimental control included (1) 1-D vertical velocity functions determined from traveltime inversion of vertical radar profiles (VRP) and (2) neutron porosity logs. We estimated the porosity distribution in the saturated zone using both the Topp and Complex Refractive Index Method (CRIM) equations and found the CRIM estimates in better agreement with the neutron logs. We found that when averaged over the length of the borehole, surface-derived velocity measurements were within 5% of the VRP velocities and that the porosity differed from the neutron log by less than 0.05. The uncertainty, however, is scale dependent. We found that the standard deviation of differences between ground-penetrating-radar-derived and neutron-log-derived porosity values was as high as 0.06 at an averaging length of 0.25 m but decreased to less than 0.02 at length scale of 11 m. Additionally, we used the 3-D porosity distribution to identify a relatively high-porosity anomaly (i.e., local sedimentary body) within a lower-porosity unit and verified the presence of the anomaly using the neutron porosity logs. Since the reflection tomography approach requires only surface data, it can provide rapid assessment of bulk hydrologic properties, identify meterscale anomalies of hydrologic significance, and may provide input for other higherresolution measurement methods

Multivariate Analysis of Cross-Hole Georadar Velocity and Attenuation Tomograms for Aquifer Zonation

We have investigated the potential of combining cross-hole georadar velocity and attenuation tomography as a method for characterizing heterogeneous alluvial aquifers. A multivariate statistical technique, known as k-means cluster analysis, is used to correlate and integrate information contained in velocity and attenuation tomograms. Cluster analysis allows us to identify objectively the major common trends in the tomographic data and thus to ‘‘reduce’’ the information to a limited number of characteristic parameter combinations. The application of this procedure to two synthetic data sets indicates that it is a powerful tool for converting the complex relationships between the tomographically derived velocity and attenuation structures into a lithologically and hydrologically meaningful zonation of the probed region. In addition, these synthetic examples allow us to evaluate the reliability of further petrophysical parameter estimates. We find that although absolute values of the tomographically inferred petrophysical parameters often differ significantly from the actual parameters, the clustering approach enables us to reliably identify the major trends in the petrophysical properties. Finally, we have applied the approach to a cross-hole georadar data set collected in a well-studied alluvial aquifer. A comparison of the clustered tomographic section with well-log data demonstrates that our approach delineates the hydrostratigraphic zonation

Reflectivity Modeling of a Ground-Penetrating-Radar Profile of a Saturated Fluvial Formation

Major horizons in radar reflection profiles may correlate with contacts between stratigraphic units or with structural breaks such as fault surfaces. Minor reflections may be caused by clutter or, in some cases, may indicate material properties or structure within stratigraphic units. In this study, we examine the physical basis for major and minor reflections observed in a shallow, unconfined, fluvial aquifer near Boise, Idaho, U. S. A. We compare a 2D profile from a surface ground-penetrating-radar reflection transect with the 1D modeled reflection profiles at three wells adjacent to the surface-reflection profile. The 1D models are based on dielectric constant and electrical conductivity values from borehole logs and vertical radar profile data. Reflections at the water table/capillary fringe, at the base of a sand-filled channel, and at the base of two sand-rich lenses in a cobble-dominated unit are recognizable in the surfacereflection profile and in all 1D reflectivity models. Less prominent reflections in stratigraphic units occur in both the surface-profile model and the reflectivity model. Although such minor reflections are not correlated easily, general similarities in their presence and location indicate that sometimes the reflections may be useful for recognizing internal facies structure or character

Hydrological Parameter Estimations from a Conservative Tracer Test with Variable-Density Effects at the Boise Hydrogeophysical Research Site

Reliable predictions of groundwater flow and solute transport require an estimation of the detailed distribution of the parameters (e.g., hydraulic conductivity, effective porosity) controlling these processes. However, such parameters are difficult to estimate because of the inaccessibility and complexity of the subsurface. In this regard, developments in parameter estimation techniques and investigations of field experiments are still challenging and necessary to improve our understanding and the prediction of hydrological processes. Here we analyze a conservative tracer test conducted at the Boise Hydrogeophysical Research Site in 2001 in a heterogeneous unconfined fluvial aquifer. Some relevant characteristics of this test include: variable-density (sinking) effects because of the injection concentration of the bromide tracer, the relatively small size of the experiment, and the availability of various sources of geophysical and hydrological information. The information contained in this experiment is evaluated through several parameter estimation approaches, including a grid-search-based strategy, stochastic simulation of hydrological property distributions, and deterministic inversion using regularization and pilot-point techniques. Doing this allows us to investigate hydraulic conductivity and effective porosity distributions and to compare the effects of assumptions from several methods and parameterizations. Our results provide new insights into the understanding of variable-density transport processes and the hydrological relevance of incorporating various sources of information in parameter estimation approaches. Among others, the variable-density effect and the effective porosity distribution, as well as their coupling with the hydraulic conductivity structure, are seen to be significant in the transport process. The results also show that assumed prior information can strongly influence the estimated distributions of hydrological properties