Wavelet dispersion caused by frequency-dependent attenuation is a common occurrence in ground-penetrating radar (GPR) data, and is displayed in the radar image as a characteristic “blurriness” that increases with depth. Correcting for wavelet dispersion is an important step that should be performed before GPR data are used for either qualitative interpretation or the quantitative determination of subsurface electrical properties. Over the bandwidth of a GPR wavelet, the attenuation of electromagnetic waves in many geological materials is approximately linear with frequency. As a result, the change in shape of a radar pulse as it propagates through these materials can be well described using one parameter, Q∗, related to the slope of the linear region. Assuming that all subsurface materials can be characterized by some Q∗ value, the problem of estimating and correcting for wavelet dispersion becomes one of determining Q∗ in the subsurface and deconvolving its effects using an inverse-Q filter. We present a method for the estimation of subsurface Q∗ from reflection GPR data based on a technique developed for seismic attenuation tomography. Essentially, Q∗ is computed from the downshift in the dominant frequency of the GPR signal with time. Once Q∗ has been obtained, we propose a damped-least-squares inverse-Q filtering scheme based on a causal, linear model for constant-Q wave propagation as a means of removing wavelet dispersion. Tests on synthetic and field data indicate that these steps can be very effective at enhancing the resolution of the GPR image.

To obtain the highest-resolution ray-based tomographic images from crosshole ground-penetrating radar (GPR) data, wide angular ray coverage of the region between the two boreholes is required. Unfortunately, at borehole spacings on the order of a few meters, high-angle traveltime data (i.e., traveltime data corresponding to transmitter-receiver angles greater than approximately 50° from the horizontal) are notoriously difficult to incorporate into crosshole GPR inversions. This is because (1) low signal-to-noise ratios make the accurate picking of first-arrival times at high angles extremely difficult, and (2) significant tomographic artifacts commonly appear when high- and low-angle ray data are inverted together. We address and overcome thesetwo issues for a crosshole GPR data example collected at the Boise Hydrogeophysical Research Site (BHRS). To estimate first-arrival times on noisy, high-angle gathers, we develop a robust and automatic picking strategy based on crosscorrelations, where reference waveforms are determined from the data through the stacking of common-ray-angle gathers. To overcome incompatibility issues between high- and low-angle data, we modify the standard tomographic inversion strategy to estimate, in addition to subsurface velocities, parameters that describe a traveltime ‘correction curve’ as a function of angle. Application of our modified inversion strategy, to both synthetic data and the BHRS data set, shows that it allows the successful incorporation of all available traveltime data to obtain significantly improved subsurface velocity images.

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. This underlines the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raises some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.