Geology topics

S.M. Gorelick

Publications and source records attributed to S.M. Gorelick.

6 recordsLinked to original sources

USGS2011

Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths

This study investigates connectivity in a small portion of the extremely heterogeneous aquifer at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi. A total of 19 fully penetrating soil cores were collected from a rectangular grid of 4 m by 4 m. Detailed grain size analysis was performed on 5 cm segments of each core, yielding 1740 hydraulic conductivity ( K ) estimates. Three different geostatistical simulation methods were used to generate 3‐D conditional realizations of the K field for the sampled block. Particle tracking calculations showed that the fastest particles, as represented by the first 5% to arrive, converge along preferential flow paths and exit the model domain within preferred areas. These 5% fastest flow paths accounted for about 40% of the flow. The distribution of preferential flow paths and particle exit locations is clearly influenced by the occurrence of clusters formed by interconnected cells with K equal to or greater than the 0.9 decile of the data distribution (10% of the volume). The fraction of particle paths within the high‐ K clusters ranges from 43% to 69%. In variogram‐based K fields, some of the fastest paths are through media with lower K values, suggesting that transport connectivity may not require fully connected zones of relatively homogenous K . The high degree of flow and transport connectivity was confirmed by the values of two groups of connectivity indicators. In particular, the ratio between effective and geometric mean K (on average, about 2) and the ratio between the average arrival time and the arrival time of the fastest particles (on average, about 9) are consistent with flow and advective transport behavior characterized by channeling along preferential flow paths.

Water Resources Research↗

USGS2006

Design and cost analysis of rapid aquifer restoration systems using flow simulation and quadratic programming

Detailed two-dimensional flow simulation of a complex ground-water system is combined with quadratic and linear programming to evaluate design alternatives for rapid aquifer restoration. The design model ensures that a contaminant plume is removed and treated within four years at the least possible cost. Rapid restoration is accomplished by maintaining specified velocities around the plume perimeter toward a group of pumping wells located near the plume center. Ground-water velocities are adjusted to include the effects of solute retardation due to sorption. As a simplification, the model does not account for hydrodynamic dispersion. Results show how treatment and pumping costs depend dynamically on the type of treatment process, the capacity of pumping and injection wells, and the number of wells. The design for an inexpensive treatment process minimizes pumping costs, while an expensive process results in the minimization of treatment costs. Substantial reductions in pumping costs occur with increases in injection capacity or in the number of wells. Treatment costs are reduced by expansions in pumping capacity or injection capacity. The analysis identifies maximum pumping and injection capacities.

Groundwater↗

USGS2004

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

New Hampshire↗

USGS2000

Attenuation-difference radar tomography: Results of a multiple-plane experiment at the U.S. Geological Survey Fractured-Rock Research Site, Mirror Lake, New Hampshire

Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride tracer injection tests conducted within the FSE, wellfield at the U.S. Geological Survey Fractured-Rock Hydrology Research Site in Grafton County, New Hampshire, USA. Borehole-radar tomography surveys were conducted using the sequential-scanning and injection method in three boreholes that form a triangular prism of adjoining tomographic image planes. Results indicate that time-lapse tomography methods provide high-resolution images of tracer distribution in permeable zones.

New Hampshire↗

USGS1986

Evaluation of the matrix exponential for use in ground-water-flow and solute-transport simulations; theoretical framework

It is possible to obtain analytic solutions to the groundwater flow and solute transport equations if space variables are discretized but time is left continuous. From these solutions, hydraulic head and concentration fields for any future time can be obtained without ' marching ' through intermediate time steps. This analytical approach involves matrix exponentiation and is referred to as the Matrix Exponential Time Advancement (META) method. Two algorithms are presented for the META method, one for symmetric and the other for non-symmetric exponent matrices. A numerical accuracy indicator, referred to as the matrix condition number, was defined and used to determine the maximum number of significant figures that may be lost in the META method computations. The relative computational and storage requirements of the META method with respect to the time marching method increase with the number of nodes in the discretized problem. The potential greater accuracy of the META method and the associated greater reliability through use of the matrix condition number have to be weighed against this increased relative computational and storage requirements of this approach as the number of nodes becomes large. For a particular number of nodes, the META method may be computationally more efficient than the time-marching method, depending on the size of time steps used in the latter. A numerical example illustrates application of the META method to a sample ground-water-flow problem. (Author 's abstract)

Water-Resources Investigations Report↗

USGS1985

Hydraulic gradient control for groundwater contaminant removal

The Rocky Mountain Arsenal near Denver, Colorado, U.S.A., is used as a realistic setting for a hypothetical test of a procedure that plans the hydraulic stabilization and removal of a groundwater contaminant plume. A two-stage planning procedure successfully selects the best wells and their optimal pumping/recharge schedules to contain the plume while a well or system of wells within the plume removes the contaminated water. In stage I, a combined groundwater flow and solute transport model is used to simulate contaminant removal under an assumed velocity field. The result is the approximated plume boundary location as a function of time. In stage II, a linear program, which includes a groundwater flow model as part of the set of constraints, determines the optimal well selection and their optimal pumping/recharge schedules by minimizing total pumping and recharge. The simulation—management model eliminates wells far from the plume perimeter and activates wells near the perimeter as the plume decreases in size. This successfully stablizes the hydraulic gradient during aquifer cleanup.

Colorado↗