Seismic refraction methods are used in environmental and engineering studies to image the shallow subsurface. We present a blind test of inversion and tomographic refraction analysis methods using a synthetic first-arrival-time dataset that was made available to the community in 2010. The data are realistic in terms of the near-surface velocity model, shot-receiver geometry and the data's frequency and added noise. Fourteen estimated models were determined by ten participants using eight different inversion algorithms, with the true model unknown to the participants until it was revealed at a session at the 2011 SAGEEP meeting. The estimated models are generally consistent in terms of their large-scale features, demonstrating the robustness of refraction data inversion in general, and the eight inversion algorithms in particular. When compared to the true model, all of the estimated models contain a smooth expression of its two main features: a large offset in the bedrock and the top of a steeply dipping low-velocity fault zone. The estimated models do not contain a subtle low-velocity zone and other fine-scale features, in accord with conventional wisdom. Together, the results support confidence in the reliability and robustness of modern refraction inversion and tomographic methods.

No abstract

Dense nonaqueous phase liquids (DNAPLs) are a pervasive and persistent category of groundwater contamination. In an effort to better understand their unique subsurface behavior, a controlled and carefully monitored injection of PCE (perchloroethylene), a typical DNAPL, was performed in conjunction with the University of Waterloo at Canadian Forces Base Borden in 1991. Of the various geophysical methods used to monitor the migration of injected PCE, the U.S. Geological Survey collected 500-MHz ground penetrating radar (GPR) data. These data are used in determining calibration parameters for a multiphase flow simulation. GPR data were acquired over time on a fixed two-dimensional surficial grid as the DNAPL was injected into the subsurface. Emphasis is on the method of determining DNAPL saturation values from this time-lapse GPR data set. Interactive full-waveform GPR modeling of regularized field traces resolves relative dielectric permittivity versus depth profiles for pre-injection and later-time data. Modeled values are end members in recursive calculations of the Bruggeman-Hanai-Sen (BHS) mixing formula, yielding interpreted pre-injection porosity and post-injection DNAPL saturation values. The resulting interpreted physical properties of porosity and DNAPL saturation of the Borden test cell, defined on a grid spacing of 50 cm with 1-cm depth resolution, are used as observations for calibration of a 3-D multiphase flow simulation. Calculated values of DNAPL saturation in the subsurface at 14 and 22 hours after the start of injection, from both the GPR and the multiphase flow modeling, are interpolated volumetrically and presented for visual comparison.

Hardware: This program will run on most standard IBM-compatible personal computers with an 80386/80387, 80486, or Pentium central processing unit. It requires 1) a display and graphics adapter that will support the VESA (1991) standard 1024x768x256 SVGA mode, 2) at least 4 Mbyte of memory available to the executable program, and 3) storage capacity (eg. hard disk) of at least 3 Mbyte to hold the required software. Software:The minimum software to run this program includes MS-DOS 3.0 or higher operating system, the executable file, and at least one transmission wavelet file. The format of the transmission wavelet file is described in detail in this report. Other files containing radar field data can be used by the program, and are also described in this report. To change the source code and re-compile the program requires the source code, a program editor, a 32-bit C compiler, Olhoeft's (1994a) graphics library and include file, a 32-bit version of Pinson's (1991) window library and include file, standard C include files, and a C linker, librarian and binder. The program is only known to properly compile and link with the Intel 386/486 C Code Builder Kit, version 1.1 a. Other compiler kits may or may not work.

a b s t r a c t a r t i c l e i n f oA recent evaluation of Martis Creek Dam highlighted the potential for dam failure due to either seepage or an earthquake on nearby faults. In 1972, the U.S. Army Corps of Engineers constructed this earthen dam, located within the Truckee Basin to the north of Lake Tahoe, CA for water storage and flood control. Past attempts to raise the level of the Martis Creek Reservoir to its design level have been aborted due to seepage at locations downstream, along the west dam abutment, and at the base of the spillway. In response to these concerns, the U.S. Geological Survey has undertaken a comprehensive suite of geophysical investigations aimed at understanding the interplay between geologic structure, seepage patterns, and reservoir and groundwater levels. This paper concerns the geologic structure surrounding Martis Creek Dam and emphasizes the importance of a regional-scale understanding to the interpretation of engineering-scale geophysical data. Our studies reveal a thick package of sedimentary deposits interbedded with Plio-Pleistocene volcanic flows; both the deposits and the flows are covered by glacial outwash. Magnetic field data, seismic tomography models, and seismic reflections are used to determine the distribution and chronology of the volcanic flows. Previous estimates of depth to basement (or the thickness of the interbedded deposits) was 100 m. Magnetotelluric soundings suggest that electrically resistive bedrock may be up to 2500 m deep. Both the Polaris Fault, identified outside of the study area using airborne LiDAR, and the previously unnamed Martis Creek Fault, have been mapped through the dam area using ground and airborne geophysics. Finally, as determined by direct-current resistivity imaging, time-domain electromagnetic sounding, and seismic refraction, the paleotopography of the interface between the sedimentary deposits and the overlying glacial outwash plays a principal role both in controlling groundwater flow and in the distribution of the observed seepage.Published by Elsevier B.V.