Recent advances in sensor technologies, field methodologies, numerical modeling, and inversion approaches have contributed to unprecedented imaging of hydrogeological properties and detailed predictions at multiple temporal and spatial scales. Nevertheless, imaging results and predictions will always remain imprecise, which calls for appropriate uncertainty quantification (UQ). In this paper, we outline selected methodological developments together with pioneering UQ applications in hydrogeology and hydrogeophysics. The applied mathematics and statistics literature is not easy to penetrate and this review aims at helping hydrogeologists and hydrogeophysicists to identify suitable approaches for UQ that can be applied and further developed to their specific needs. To bypass the tremendous computational costs associated with forward UQ based on full-physics simulations, we discuss proxy-modeling strategies and multi-resolution (Multi-level Monte Carlo) methods. We consider Bayesian inversion for non-linear and non-Gaussian state-space problems and discuss how Sequential Monte Carlo may become a practical alternative. We also describe

[1] Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.Citation: Irving, J., and K. Singha (2010), Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities, Water Resour. Res., 46, W11514,

No abstract

Wavelet dispersion caused by frequency-dependent attenuation is a common occurrence in groundpenetrating 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.

No abstract

No abstract

[1] Over the past decade, significant interest has been expressed in relating the spatial statistics of surface-based reflection ground-penetrating radar (GPR) data to those of the imaged subsurface volume. A primary motivation for this work is that changes in the radar wave velocity, which largely control the character of the observed data, are expected to be related to corresponding changes in subsurface water content. Although previous work has indeed indicated that the spatial statistics of GPR images are linked to those of the water content distribution of the probed region, a viable method for quantitatively analyzing the GPR data and solving the corresponding inverse problem has not yet been presented. Here we address this issue by first deriving a relationship between the 2-D autocorrelation of a water content distribution and that of the corresponding GPR reflection image. We then show how a Bayesian inversion strategy based on Markov chain Monte Carlo sampling can be used to estimate the posterior distribution of subsurface correlation model parameters that are consistent with the GPR data. Our results indicate that if the underlying assumptions are valid and we possess adequate prior knowledge regarding the water content distribution, in particular its vertical variability, this methodology allows not only for the reliable recovery of lateral correlation model parameters but also for estimates of parameter uncertainties. In the case where prior knowledge regarding the vertical variability of water content is not available, the results show that the methodology still reliably recovers the aspect ratio of the heterogeneity.Citation: Irving, J., R. Knight, and K. Holliger (2009), Estimation of the lateral correlation structure of subsurface water content from surface-based ground-penetrating radar reflection images, Water Resour. Res., 45, W12404,

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.