Geophysics and Seismology Commons^™

Crosshole Radar Tomography In A Fluvial Aquifer Near Boise, Idaho, William P. Clement, Warren Barrash

Geosciences Faculty Publications and Presentations

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 …

Ground-Penetrating Radar Theory And Application Of Thin-Bed Offset-Dependent Reflectivity, John H. Bradford, Jacob C. Deeds

CGISS Publications and Presentations

Offset-dependent reflectivity or amplitude-variationwith- offset (AVO) analysis of ground-penetrating radar (GPR) data may improve the resolution of subsurface dielectric permittivity estimates. A horizontally stratified medium has a limiting layer thickness below which thin-bed AVO analysis is necessary. For a typical GPR signal, this limit is approximately 0.75 of the characteristic wavelength of the signal. Our approach to modeling the GPR thin-bed response is a broadband, frequency-dependent computation that utilizes an analytical solution to the three-interface reflectivity and is easy to implement for either transverse electric (TE) or transverse magnetic (TM) polarizations. The AVO curves for TE and TM modes differ …

Reflectivity Modeling Of A Ground-Penetrating-Radar Profile Of A Saturated Fluvial Formation, William P. Clement, Warren Barrash, Michael D. Knoll

Geosciences Faculty Publications and Presentations

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 …

Investigating The Stratigraphy Of An Alluvial Aquifer Using Crosswell Seismic Traveltime Tomography, Geoff J.M. Moret, Michael D. Knoll, Warren Barrash, William P. Clement

CGISS Publications and Presentations

In this study, we investigate the use of crosswell P-wave seismic tomography to obtain spatially extensive information about subsurface sedimentary architecture and heterogeneity in alluvial aquifers. Our field site was a research wellfield in an unconfined aquifer near Boise, Idaho. The aquifer consists of a ~ 20-m-thick sequence of alluvial cobble- and-sand deposits, which have been subdivided into five stratigraphic units based on neutron porosity logs, grainsize analysis, and radar reflection data. We collected crosswell and borehole-to-surface seismic data in wells 17.1 m apart. We carefully considered the impact of well deviation, data quality control, and the choice of inversion …

Traveltime Inversion Of Vertical Radar Profiles, William P. Clement, Michael D. Knoll

Geosciences Faculty Publications and Presentations

Traveltimes of direct arrivals in vertical radar profiles (VRPs) are tomographically inverted to estimate the earth’s electromagnetic (EM) velocity between a surface transmitter and a downhole receiver. We determine the 1D interval velocity model that best fits the first-arrival traveltimes by using a weighted, damped, least-squares inversion scheme. We assess the accuracy of the velocity model using synthetic traveltimes from a known velocity-distribution model simulating an unconfined aquifer. The inverted velocity profile closely matched the velocity profile of the input model in the synthetic examples. Using vertical radar profile data from an unconfined aquifer near Boise, Idaho, we inverted traveltimes …

Image Interpretation Using Appraisal Analysis, Partha S. Routh, Carlyle R. Miller

Geosciences Faculty Publications and Presentations

In geophysical inversion, a significant effort is invested to obtain images of the Earth from finite data. The first step is to obtain an image i.e. solve the inverse problem. This step alone provides significant challenges that are not addressed inthis paper. The next step is to interpret the image in terms of specific questions. For example, what can we say about the average value of a physical property within a certain region of the model? What scale information can we resolve from the data? These questions are problem dependent and may require that inversion be carried out several times …

Applying Reflection Tomography In The Postmigration Domain To Multifold Ground-Penetrating Radar Data, John H. Bradford

CGISS Publications and Presentations

Acquisition and processing of multifold ground-penetrating radar (GPR) data enable detailed measurements of lateral velocity variability. The velocities constrain interpretation of subsurface materials and lead to significant improvement in image accuracy when coupled with prestack depth migration (PSDM). Reflection tomography in the postmigration domain was introduced in the early 1990s for velocity estimation in seismic reflection. This robust, accurate method is directly applicable in multifold GPR imaging. At a contaminated waste facility within the U.S. Department of Energy’s Hanford site in Washington, the method is used to identify significant lateral and vertical velocity heterogeneity associated with infilled waste pits. Using …

Seeing Through The Noise: Seismic Reflection Profiling In Urban Areas, Lee M. Liberty, Thomas L. Pratt, Scott S. Hess

CGISS Publications and Presentations

Studies for urban hazard or resource assessment often take place in densely populated areas characterized by considerable cultural noise. These site conditions can severely compromise seismic reflection data quality. We have collected vibroseis and hammer (weight drop) seismic reflection data in a range of geologic conditions to image stratigraphy and structures in the upper one km along regional highways, city streets,and power line access roads. In addition to the challenges of safety and outreach, acquisition efforts along busy streets and highways often encounter poor receiver coupling and large-amplitude coherent noise from traffic and power lines. Although higher quality seismic reflection …

Crosshole Radar Tomography In An Alluvial Aquifer Nearboise, Idaho, William P. Clement, Warren Barrash

CGISS Publications and Presentations

To determine the distribution of heterogeneities in an unconfined aquifer in Boise, ID, we compute radar 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 sequence of layers 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 …

Relating Damping To Soil Permeability, Paul Michaels

CGISS Publications and Presentations

Published comparisons of complex moduli in dry and saturated soils have shown that viscous behavior is only evident when a sufficiently massive viscous fluid (like water) is present. That is, the loss tangent is frequency dependent for water saturated specimens, but nearly frequency independent for dry samples. While the Kelvin-Voigt (KV) representation of a soil captures the general viscous behavior using a dashpot, it fails to account for the possibly separate motions of the fluid and frame (there is only a single mass element). An alternative representation which separates the two masses, water and frame, is presented here. This Kelvin-Voigt-Maxwell-Biot …

Comparison Of Viscous Damping In Unsaturated Soils, Compression And Shear, Paul Michaels

CGISS Publications and Presentations

Geophysical down-hole surveys can be used to measure the small strain dynamic properties of soils by the effects these properties have on wave propagation. The relevant effects include amplitude decay (corrected for beam divergence) and velocity dispersion. In this paper, down-hole data collected during the GeoInstitute's Denver 2000 field day are presented and analyzed as a Kelvin-Voigt solid. Findings for these unsaturated soils include viscous damping and stiffness which differ significantly for shear and compressional waves. A strong viscous damping is observed in compression, but weak damping is presented in shear. Lumped parameter constitutive models are discussed which mathematically represent …

Wave Field Migration As A Tool For Estimating Spatially Continuous Radar Velocity And Water Content In Glaciers, John H. Bradford, Joel T. Harper

CGISS Publications and Presentations

Normal-moveout velocity analysis can lead to significant overestimates of the velocity structure of temperate glaciers since most englacial reflectors approximate point scatters and violate the assumption of planar flat lying reflectors. Migration velocity analysis (MVA) is a tool that does not depend on the assumption of flat lying reflectors. MVA can provide laterally and vertically continuous velocity estimates from conventional common-offset radar sections. In a study of temperate Bench Glacier, Alaska, we used MVA coupled with dielectric modeling to estimate the distribution of water content along a cross-section of the glacier. We found the glacier has two layers, an upper …

Effects Of Signal Processing And Antenna Frequency On The Geostatistical Structure Of Ground-Penetrating Radar Data, Greg A. Oldenborger, Michael D. Knoll, Warren Barrash

CGISS Publications and Presentations

Recent research has suggested that the geostatistical structure of ground-penetrating radar data may be representative of the spatial structure of hydraulic properties. However, radar images of the subsurface can change drastically with application of signal processing or by changing the signal frequency. We perform geostatistical analyses of surface radar reflection profiles in order to investigate the effects of data processing and antenna frequency on the semivariogram structure of radar reflection amplitudes. Surface radar reflection data collected at the Boise Hydrogeophysical Research Site illustrate the processing- and antenna-dependence of radar semivariograms for a fluvial, cobble-and-sand aquifer. Compensating for signal attenuation and …

Seismic Expression Of Pleistocene Paleoceanographic Changes In The California Borderland From Digitally Acquired 3.5 Khz Subbottom Profiles And Ocean Drilling Program Leg 167 Drilling, Aleksandra Janik, Mitchell W. Lyle, Lee M. Liberty

CGISS Publications and Presentations

We correlate processed 3.5 kHz seismic profiles with physical properties of cores collected during ODP Leg 167 from the Tanner, East Cortes, and San Nicolas Basins through much of the Pleistocene succession. Results indicate that seismic horizons in the unconsolidated Pleistocene sediments (top 50 m) are mainly controlled by density contrasts. Removing of the compaction trend from the density reveals a very interesting relationship between density and composition - the density closely and inversely correlates with organic carbon indicating that large-scale variations in organic carbon are responsible for seismic reflections through their influence on density. This is a significant discovery …

Vsp Traveltime Inversion: Near-Surface Issues, Geoff J. M. Moret, William P. Clement, Michael D. Knoll, Warren Barrash

Geosciences Faculty Publications and Presentations

P-wave velocity information obtained from vertical seismic profiles (VSPs) can be useful in imaging subsurface structure, either by directly detecting changes in the subsurface or as an aid to the interpretation of seismic reflection data. In the shallow subsurface, P-wave velocity can change by nearly an order of magnitude over a short distance, so curved rays are needed to accurately model VSP traveltimes. We used a curved-ray inversion to estimate the velocity profile and the discrepancy principle to estimate the data noise level and to choose the optimum regularization parameter. The curved-ray routine performed better than a straight-ray inversion for …

Multivariate Analysis Of Cross-Hole Georadar Velocity And Attenuation Tomograms For Aquifer Zonation, Jens Tronicke, Klaus Holliger, Warren Barrash, Michael D. Knoll

CGISS Publications and Presentations

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 …

Mae Sae Basin And Wiang Nong Lom: Radiocarbon Dating And Relation To The Active Strike-Slip Mae Chan Fault, Northern Thailand, Spencer H. Wood, Fongsaward S. Singharajwarapan, Tharaporn Bundarnsin, Eric Rothwell

Geosciences Faculty Publications and Presentations

The Mai Sai basin and the Wiang Nong Lom swamp lie north of, and adjacent to, the northeast-trending Mae Chan fault, which is a left -lateral strike-slip fault in northernmost Thailand. The sediment in the Mae Sai basin is about 600 meters thick. The Mae Sai basin has an anticlinal fold below a 390-meter depth. This fold is buried by broadly downwarped sediment. In the Wiang Nong Lom swamp, sediment is about 170 meters thick. Radiocarbon ages from a 2.7-meter auger hole at Ko Mae Mai in the Wiang Nong Lom swamp show that 0.5 meter of sediment accumulated in …

Design Of Geophysical Surveys In Transportation, Paul Michaels

CGISS Publications and Presentations

Designing geophysical investigations for transportation related projects requires special attention to the constraints imposed by right-of-way, irregular topography, noise from traffic, and the need to avoid the interruption of traffic flow. A geophysical engineer needs to be prepared to consider these design issues that are not addressed in a standard procedure such as ASTM D-5777. The author presents design strategies that address these issues, and illustrates the concepts with case histories taken from bridge and highway projects. Beam steering, broadside shooting, and non-traditional designs that preserve alternative analysis options are presented. Transportation engineers who augment traditional subsurface geotechnical surveys with …

Recent Paleoseismic Investigations In Northern And Western Thailand, Clark H. Fenton, Punya Charusiri, Spencer H. Wood

Geosciences Faculty Publications and Presentations

Recent paleoseismic investigations have identified a number of active faults in Northern and Western Thailand. Northern Thailand is an intraplate basin and range province, comprised of north-south-trending Cenozoic intermontane grabens and half grabens, bounded by north- to northwest-striking normal to normal-oblique faults and northeast-striking left-lateral strike-slip faults. The basin-bounding normal faults are marked by steep, linear range fronts with triangular facets and wineglass canyons and have slip rates of 0.1 to 0.8 mm/yr. Based on limited data, the average vertical displacement-per-event is about 1.0 to 1.5 m. These faults are characterized by recurrence intervals of thousands to tens of thousands …

Lower Miocene To Present Stratigraphy Of The Equatorial Pacific Sediment Bulge And Carbonate Dissolution Anomalies, Neil C. Mitchell, Mitchell W. Lyle, Marie B. Knappenberger, Lee M. Liberty

CGISS Publications and Presentations

The bulge is a 600-m-thick regional deposit of pelagic sediment accumulated around the equator. Its stratigraphy reflects a number of factors: how accumulation rates have varied over time, how accumulation has been spatially focused around the equator, how much carbonate dissolution and reworking or nondeposition by bottom currents have occurred, and how much the deposits have been translated northward by motion of the Pacific tectonic plate on which they have accumulated. In order to fully explore the effects of these processes, a spatially continuous stratigraphic database is desirable, as existing cores provide information at only discrete points, and they tend …

Tracer/Time-Lapse Radar Imaging Test At The Boise Hydrogeophysical Research Site, Warren Barrash, Michael D. Knoll, David W. Hyndman, Tom Clemo, Edward C. Reboulet, Elisabeth M. Hausrath

Geosciences Faculty Publications and Presentations

A combined tracer and time-lapse radar imaging experiment was conducted in the unconfined coarse fluvial aquifer at the Boise Hydrogeophysical Research Site in August, 2001. Two tracers (bromide and uranine) were injected to form a plume over a 4-m interval that spanned the contact between hydrostratigraphic units with contrasting permeability. The tracer plume traveled 6.9 m to well B6, passing through well A1 instrumented with 20 sampling zones over a 5-m interval that spanned the injection interval. Radar tomographic data were collected periodically on cross-sectional and longitudinal planes, two of which passed through well A1 for quantitative calibration of radar …

Borehole Radar Attenuation-Difference Tomography During The Tracer/Time-Lapse Test At The Boise Hydrogeophysical Research Site, Sarah E. Goldstein, Tim C. Johnson, Michael D. Knoll, Warren Barrash, William P. Clement

Geosciences Faculty Publications and Presentations

A tracer test and time-lapse radar imaging experiment was conducted at the Boise Hydrogeophysical Research Site to investigate the utility of crosswell radar in imaging an electrically conductive tracer plume. A multilevel water sampling system down gradient from the tracer injection well and in the radar imaging plane was used to collect detailed, 1-dimensional, fluid electrical conductivity data during the tracer test. We compare the spatial and temporal position and concentration variations of the plume as indicated by the fluid conductivity data to those suggested by radar level run attenuation differences, shot-receiver attenuation difference crossplots, and an attenuation-difference tomogram. We …

Hierarchical Geostatistics And Multifacies Systems: Boise Hydrogeophysical Research Site, Boise, Idaho, Warren Barrash, Tom Clemo

CGISS Publications and Presentations

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 …

Equatorial Ocean Circulation In An Extremely Warm Climate, T. C. Moore, D. K. Rea, Mitchell W. Lyle, Lee M. Liberty

CGISS Publications and Presentations

The warm climates of the early Paleogene and the associated diminished near-surface winds should have resulted in a reduction in near-surface ocean circulation. One check on this deduction is the delineation of biogenic sediments associated with an equatorial current system of the early Eocene Pacific. A latitudinal seismic reflection transect across the tropical Pacific along early Paleogene ocean crust reveals a basal high-amplitude reflection package that we take to be the lower Eocene section. This unit varies in thickness by a factor of about two, with the thickest portion forming a low mound some 3˚–4˚ north of the 56 Ma …

Depth Characterization Of Shallow Aquifers With Seismic Reflection, Part Ii—Prestack Depth Migration And Field Examples, John H. Bradford, D. S. Sawyer

CGISS Publications and Presentations

It is common in shallow seismic studies for the compressional-wave velocity in unconsolidated sediments to increase by a factor of four or more at the transition from dry or partial water saturation to full saturation. Under these conditions, conventional NMO velocity analysis fails and leads to large depth and layer thickness estimates if the Dix equation is assumed valid. Prestack depth migration (PSDM) is a means of improving image accuracy. A comparison of PSDM with conventional NMO processing for three field examples from differing hydrogeologic environments illustrates that PSDM can significantly improve image quality and accuracy.

Geologic And Tectonic History Of The Western Snake River Plain, Idaho And Oregon, Spencer H. Wood, Drew M. Clemens

Geosciences Faculty Publications and Presentations

The western Snake River Plain is a Neogene-aged intracontinental rift basin, about 70 km wide and 300 km long, trending northwest across the southern Idaho batholith. Its southeastern end merges with the northeast-trending eastern plain, a structural downwarp associated with extension along the track of the Yellowstone hot spot. Orientation of the western plain rift is parallel to several regional northwest-trending crustal discontinuities, such as the Olympic-Wallowa lineament and the Brothers fault zone, suggesting that the rift failed along zones of lithospheric weakness, as the lithosphere was softened by the passing hot spot. Crustal refraction data and gravity show that …

Identification Of Subsonic P-Waves, Paul Michaels

CGISS Publications and Presentations

A field trial was conducted to test the existence of subsonic (V_p <331 m>/s) P-waves previously reported in the literature. A 1-m-long reverse profile was acquired with three-component (3C) geophones on a sandy silt (unified classification ML). The silt had a porosity of 54%, a degree of water saturation of 63%, and a plasticity index of 10. No subsonic P-waves were observed, although high frequency (up to 1200 Hz) Rayleigh waves were identified by hodogram analysis. These surface waves were observed with horizontal velocities that varied from 40 to 200 m/s. Hodogram observations and theory suggest that …

Depth Characterization Of Shallow Aquifers With Seismic Reflection, Part I—The Failure Of Nmo Velocity Analysis And Quantitative Error Prediction, John H. Bradford

CGISS Publications and Presentations

As seismic reflection data become more prevalent as input for quantitative environmental and engineering studies, there is a growing need to assess and improve the accuracy of reflection processing methodologies. It is common for compressional-wave velocities to increase by a factor of four or more where shallow, unconsolidated sediments change from a dry or partially watersaturated regime to full saturation. While this degree of velocity contrast is rare in conventional seismology, it is a common scenario in shallow environments and leads to significant problems when trying to record and interpret reflections within about the first 30 m below the water …

Use Of Engineering Geophysics To Investigate A Site For A Bridge Foundation, Paul Michaels

CGISS Publications and Presentations

The Idaho Transportation Department (ITD) commissioned a geophysical study to aid in the design of a replacement for an existing concrete span bridge. Because the river current was too swift to place geophones in the river, the solution was to shoot p-wave profiles in a reciprocal geometry (phones on land, shots in the river). In addition to the refraction work, a down-hole seismic profile was acquired to calibrate the surface data. Geotechnical boreholes drilled on the north and south river banks detected a laterally varying soil profile. The south-bank hole encountered 9.1 m of granular overburden, followed by 8.5 m …

Use Of Principal Component Analysis To Determine Down-Hole Tool Orientation And Enhance Sh-Waves, Paul Michaels

CGISS Publications and Presentations

A common problem in down-hole shear wave surveys is the determination of the tool rotation relative to the seismic source polarization direction. This paper reports an algorithm which has been developed to solve for this angle in the horizontal plane. The method employs Principal Component Analysis (PCA) to determine the angle from the large motion of the SH-wave (not first motion). Once the angle is known, the horizontal component data may be rotated so that one component is aligned with the source polarization, and the other orthogonal to it.

Significant findings include the following. First, subtraction of oppositely polarized source …