Geology topics

James J. Butler

Publications and source records attributed to James J. Butler.

12 recordsLinked to original sources

USGS2023

Evaluation of models for estimating hydraulic conductivity in glacial aquifers from NMR logging

Nuclear magnetic resonance (NMR) logging is a promising method for estimating hydraulic conductivity ( K ). During the past ∼60 years, NMR logging has been used for petroleum applications, and different models have been developed for deriving estimates of permeability. These models involve calibration parameters whose values were determined through decades of research on sandstones and carbonates. We assessed the use of five models to derive estimates of K in glacial aquifers from NMR logging data acquired in two wells at each of two field sites in central Wisconsin, USA. Measurements of K , obtained with a direct push permeameter (DPP), K DPP , were used to obtain the calibration parameters in the Schlumberger-Doll Research, Seevers, Timur-Coates, Kozeny-Godefroy, and sum-of-echoes (SOE) models so as to predict K from the NMR data; and were also used to assess the ability of the models to predict K DPP . We obtained four well-scale calibration parameter values for each model using the NMR and DPP measurements in each well; and one study-scale parameter value for each model by using all data. The SOE model achieved an agreement with K DPP that matched or exceeded that of the other models. The Timur-Coates estimates of K were found to be substantially different from K DPP . Although the well-scale parameter values for the Schlumberger-Doll, Seevers, and SOE models were found to vary by less than a factor of 2, more research is needed to confirm their general applicability so that site-specific calibration is not required to obtain accurate estimates of K from NMR logging data.

Wisconsin↗

USGS2021

Assessment of NMR logging for estimating hydraulic conductivity in glacial aquifers

Glacial aquifers are an important source of groundwater in the United States and require accurate characterization to make informed management decisions. One parameter that is crucial for understanding the movement of groundwater is hydraulic conductivity, K . Nuclear magnetic resonance (NMR) logging measures the NMR response associated with the water in geological materials. By utilizing an external magnetic field to manipulate the nuclear spins associated with 1 H, the time‐varying decay of the nuclear magnetization is measured. This logging method could provide an effective way to estimate K at submeter vertical resolution, but the models that relate NMR measurements to K require calibration. At two field sites in a glacial aquifer in central Wisconsin, we collected a total of four NMR logs and obtained measurements of K in their immediate vicinity with a direct‐push permeameter (DPP). Using a bootstrap algorithm to calibrate the Schlumberger‐Doll Research (SDR) NMR‐ K model, we estimated K to within a factor of 5 of the DPP measurements. The lowest levels of accuracy occurred in the lower‐ K ( K < 10 −4 m/s) intervals. We also evaluated the applicability of prior SDR model calibrations. We found the NMR calibration parameters varied with K , suggesting the SDR model does not incorporate all the properties of the pore space that control K . Thus, the expected range of K in an aquifer may need to be considered during calibration of NMR‐ K models. This study is the first step toward establishing NMR logging as an effective method for estimating K in glacial aquifers.

Groundwater↗

USGS2019

Radiometric calibration of a non-imaging airborne spectrometer to measure the Greenland ice sheet surface

Methods to radiometrically calibrate a non-imaging airborne visible-to-shortwave infrared (VSWIR) spectrometer to measure the Greenland ice sheet surface are presented. Airborne VSWIR measurement performance for bright Greenland ice and dark bare rock/soil targets is compared against the MODerate resolution atmospheric TRANsmission (MODTRAN ® ) radiative transfer code (version 6.0), and a coincident Landsat 8 Operational Land Imager (OLI) acquisition on 29 July 2015 during an in-flight radiometric calibration experiment. Airborne remote sensing flights were carried out in northwestern Greenland in preparation for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimeter mission. A total of nine science flights were conducted over the Greenland ice sheet, sea ice, and open-ocean water. The campaign's primary purpose was to correlate green laser pulse penetration into snow and ice with spectroscopic-derived surface properties. An experimental airborne instrument configuration that included a nadir-viewing (looking downward at the surface) non-imaging Analytical Spectral Devices (ASD) Inc. spectrometer that measured upwelling VSWIR (0.35 to 2.5 µ m) spectral radiance ( W m - 2 sr - 1 µ m - 1 "> W m − 2 sr − 1 µ m − 1 ) in the two-color Slope Imaging Multi-polarization Photon-Counting Lidar's (SIMPL) ground instantaneous field of view, and a zenith-viewing (looking upward at the sky) ASD spectrometer that measured VSWIR spectral irradiance (W m −2 nm −1 ) was flown. National Institute of Standards and Technology (NIST) traceable radiometric calibration procedures for laboratory, in-flight, and field environments are described in detail to achieve a targeted VSWIR measurement requirement of within 5 % to support calibration/validation efforts and remote sensing algorithm development. Our MODTRAN predictions for the 29 July flight line over dark and bright targets indicate that the airborne nadir-viewing spectrometer spectral radiance measurement uncertainty was between 0.6 % and 4.7 % for VSWIR wavelengths (0.4 to 2.0 µ m) with atmospheric transmittance greater than 80 %. MODTRAN predictions for Landsat 8 OLI relative spectral response functions suggest that OLI is measuring 6 % to 16 % more top-of-atmosphere (TOA) spectral radiance from the Greenland ice sheet surface than was predicted using apparent reflectance spectra from the nadir-viewing spectrometer. While more investigation is required to convert airborne VSWIR spectral radiance into atmospherically corrected airborne surface reflectance, it is expected that airborne science flight data products will contribute to spectroscopic determination of Greenland ice sheet surface optical properties to improve understanding of their potential influence on ICESat-2 measurements.

Atmospheric Measurement Techniques↗

USGS2013

A small-diameter NMR logging tool for groundwater investigations

A small-diameter nuclear magnetic resonance (NMR) logging tool has been developed and field tested at various sites in the United States and Australia. A novel design approach has produced relatively inexpensive, small-diameter probes that can be run in open or PVC-cased boreholes as small as 2 inches in diameter. The complete system, including surface electronics and various downhole probes, has been successfully tested in small-diameter monitoring wells in a range of hydrogeological settings. A variant of the probe that can be deployed by a direct-push machine has also been developed and tested in the field. The new NMR logging tool provides reliable, direct, and high-resolution information that is of importance for groundwater studies. Specifically, the technology provides direct measurement of total water content (total porosity in the saturated zone or moisture content in the unsaturated zone), and estimates of relative pore-size distribution (bound vs. mobile water content) and hydraulic conductivity. The NMR measurements show good agreement with ancillary data from lithologic logs, geophysical logs, and hydrogeologic measurements, and provide valuable information for groundwater investigations.

Groundwater↗

USGS2008

Simulation assessment of the direct‐push permeameter for characterizing vertical variations in hydraulic conductivity

The direct‐push permeameter (DPP) is a tool for the in situ characterization of hydraulic conductivity (K) in shallow, unconsolidated formations. This device, which consists of a short screened section with a pair of pressure transducers near the screen, is advanced into the subsurface with direct‐push technology. K is determined through a series of injection tests conducted between advancements. Recent field work by Butler et al. (2007) has shown that the DPP holds great potential for describing vertical variations in K at an unprecedented level of detail, accuracy and speed. In this paper, the fundamental efficacy of the DPP is evaluated through a series of numerical simulations. These simulations demonstrate that the DPP can provide accurate K information under conditions commonly faced in the field. A single DPP test provides an effective K for the domain immediately surrounding the interval between the injection screen and the most distant pressure transducer. Features that are thinner than that interval can be quantified by reducing the vertical distance between successive tests and analyzing the data from all tests simultaneously. A particular advantage of the DPP is that, unlike most other single borehole techniques, a low‐K skin or a clogged screen has a minimal impact on the K estimate. In addition, the requirement that only steady‐shape conditions be attained allows for a dramatic reduction in the time required for each injection test.

Water Resources Research↗

USGS2007

A field assessment of the value of steady shape hydraulic tomography for characterization of aquifer heterogeneities

Hydraulic tomography is a promising approach for obtaining information on variations in hydraulic conductivity on the scale of relevance for contaminant transport investigations. This approach involves performing a series of pumping tests in a format similar to tomography. We present a field‐scale assessment of hydraulic tomography in a porous aquifer, with an emphasis on the steady shape analysis methodology. The hydraulic conductivity (K) estimates from steady shape and transient analyses of the tomographic data compare well with those from a tracer test and direct‐push permeameter tests, providing a field validation of the method. Zonations based on equal‐thickness layers and cross‐hole radar surveys are used to regularize the inverse problem. The results indicate that the radar surveys provide some useful information regarding the geometry of the K field. The steady shape analysis provides results similar to the transient analysis at a fraction of the computational burden. This study clearly demonstrates the advantages of hydraulic tomography over conventional pumping tests, which provide only large‐scale averages, and small‐scale hydraulic tests (e.g., slug tests), which cannot assess strata connectivity and may fail to sample the most important pathways or barriers to flow.

Water Resources Research↗

USGS2007

A field investigation of phreatophyte‐induced fluctuations in the water table

Hydrographs from shallow wells in vegetated riparian zones frequently display a distinctive pattern of diurnal water table fluctuations produced by variations in plant water use. A multisite investigation assessed the major controls on these fluctuations and the ecohydrologic insights that can be gleaned from them. Spatial and temporal variations in the amplitude of the fluctuations are primarily a function of variations in (1) the meteorological drivers of plant water use, (2) vegetation density, type, and vitality, and (3) the specific yield of sediments in the vicinity of the water table. Past hydrologic conditions experienced by the riparian zone vegetation, either in previous years or earlier within the same growing season, are also an important control. Diurnal water table fluctuations can be considered a diagnostic indicator of groundwater consumption by phreatophytes at most sites, so the information embedded within these fluctuations should be more widely exploited in ecohydrologic studies.

Water Resources Research↗

USGS2005

Estimation of groundwater consumption by phreatophytes using diurnal water table fluctuations: A saturated‐unsaturated flow assessment

Groundwater consumption by phreatophytes is a difficult‐to‐measure but important component of the water budget in many arid and semiarid environments. Over the past 70 years the consumptive use of groundwater by phreatophytes has been estimated using a method that analyzes diurnal trends in hydrographs from wells that are screened across the water table (White, 1932). The reliability of estimates obtained with this approach has never been rigorously evaluated using saturated‐unsaturated flow simulation. We present such an evaluation for common flow geometries and a range of hydraulic properties. Results indicate that the major source of error in the White method is the uncertainty in the estimate of specific yield. Evapotranspirative consumption of groundwater will often be significantly overpredicted with the White method if the effects of drainage time and the depth to the water table on specific yield are ignored. We utilize the concept of readily available specific yield as the basis for estimation of the specific yield value appropriate for use with the White method. Guidelines are defined for estimating readily available specific yield based on sediment texture. Use of these guidelines with the White method should enable the evapotranspirative consumption of groundwater to be more accurately quantified.

Water Resources Research↗

USGS2004

Hydraulic tests in highly permeable aquifers

A semianalytical solution is presented for a mathematical model describing the flow of groundwater in response to a slug or pumping test in a highly permeable, confined aquifer. This solution, which is appropriate for wells of any degree of penetration and incorporates inertial mechanisms at both the test and observation wells, can be used to gain new insights into hydraulic tests in highly permeable settings. The oscillatory character of slug‐ and pumping‐induced responses will vary considerably across a site, even in an essentially homogeneous formation, when wells of different radii, depths, and screen lengths are used. Thus variations in the oscillatory character of responses do not necessarily indicate variations in hydraulic conductivity (K). Existing models for slug tests in partially penetrating wells in high‐K aquifers neglect the storage properties of the media. That assumption, however, appears reasonable for a wide range of common conditions. Unlike in less permeable formations, drawdown at an observation well in a high‐K aquifer will be affected by head losses in the pumping well. Those losses, which affect the form of the pumping‐induced oscillations, can be difficult to characterize. Thus analyses of observation‐well drawdown should utilize data from the period after the oscillations have dissipated whenever possible. Although inertial mechanisms can have a large impact on early‐time drawdown, that impact decreases rapidly with duration of pumping and distance to the observation well. Conventional methods that do not consider inertial mechanisms should therefore be viable options for the analysis of drawdown data at moderate to large times.

Water Resources Research↗

USGS2003

Pumping‐induced leakage in a bounded aquifer: An example of a scale‐invariant phenomenon

A new approach is presented for calculation of the volume of pumping‐induced leakage entering an aquifer as a function of time. This approach simplifies the total leakage calculation by extending analytical‐based methods developed for infinite systems to bounded aquifers of any size. The simplification is possible because of the relationship between drawdown and leakage in aquifers laterally bounded by impermeable formations. This relationship produces a scale‐invariant total leakage; i.e., the volume of leakage as a function of time does not change with the size of the aquifer or with the location of the pumping well. Two examples and image well theory are used to demonstrate and prove, respectively, the generality of this interesting phenomenon.

Water Resources Research↗

USGS2002

Steady shape analysis of tomographic pumping tests for characterization of aquifer heterogeneities

Hydraulic tomography, a procedure involving the performance of a suite of pumping tests in a tomographic format, provides information about variations in hydraulic conductivity at a level of detail not obtainable with traditional well tests. However, analysis of transient data from such a suite of pumping tests represents a substantial computational burden. Although steady state responses can be analyzed to reduce this computational burden significantly, the time required to reach steady state will often be too long for practical applications of the tomography concept. In addition, uncertainty regarding the mechanisms driving the system to steady state can propagate to adversely impact the resulting hydraulic conductivity estimates. These disadvantages of a steady state analysis can be overcome by exploiting the simplifications possible under the steady shape flow regime. At steady shape conditions, drawdown varies with time but the hydraulic gradient does not. Thus transient data can be analyzed with the computational efficiency of a steady state model. In this study, we demonstrate the value of the steady shape concept for inversion of hydraulic tomography data and investigate its robustness with respect to improperly specified boundary conditions.

Water Resources Research↗

USGS2002

Dynamic interpretation of slug tests in highly permeable aquifers

Considerable progress has been made in developing a theoretical framework for modeling slug test responses in formations with high hydraulic conductivity K . However, several questions of practical significance remain unresolved. Given the rapid and often oscillatory nature of test responses, the traditional hydrostatic relationship between the water level and the transducer‐measured head in the water column may not be appropriate. A general dynamic interpretation is proposed that describes the relationship between water level response and transducer‐measured head. This theory is utilized to develop a procedure for transforming model‐generated water level responses to transducer readings. The magnitude of the difference between the actual water level position and the apparent position based on the transducer measurement is a function of the acceleration and velocity of the water column, test geometry, and depth of the transducer. The dynamic approach explains the entire slug test response, including the often‐noted discrepancy between the actual initial water level displacement and that measured by a transducer in the water column. Failure to use this approach can lead to a significant underestimation of K when the transducer is a considerable distance below the static water level. Previous investigators have noted a dependence of test responses on the magnitude of the initial water level displacement and have developed various approximate methods for analyzing such data. These methods are re‐examined and their limitations clarified. Practical field guidelines are proposed on the basis of findings of this work. The soundness of the dynamic approach is demonstrated through a comparison of K profiles from a series of multilevel slug tests with those from dipole‐flow tests performed in the same wells.

Water Resources Research↗