- Soil Science Society of America
We address a fundamental problem inherent in least squares based ground penetrating radar tomography problems, and linear inverse Gaussian problems in general: how should the a priori covariance model be chosen? The choice of such a prior covariance model is most often a very subjective task that has major implications on the result of the inversion. We present a method that allows quantification of the likelihood that a given choice of prior covariance model is consistent with the observed tomography data. This is done by comparing statistical properties of samples of the prior and posterior probability density function of the tomographic inverse problem. In essence, if samples of the posterior are unlikely samples of the prior, then such a choice of a priori covariance model is deemed unlikely. This enables one to quantify the consistency of a number of equally probable prior covariance models to data observations. A synthetic data set was used to describe and validate the approach. We determined how a known covariance model could be inferred from a synthetic tomography problem. The methodology was then applied to a nonlinear ground penetrating radar tomography case study. The covariance model deemed most likely was consistent with nearby ground penetrating radar reflection profiles. The method provides useful results even if just a subset as small as 10% of the available data is considered.
- GPR, ground penetrating radar
- HPD, highest probability density
- McMC, Markov chain Monte Carlo
- pdf, probability density function
- Received September 29, 2006.