Magnetic resonance tomography constrained by ground-penetrating radar for improved hydrogeophysical characterization

Geophysical methods can characterize aquifer systems noninvasively and are particularly helpful to image the complex depositional architecture of the subsurface. Among these, ground-penetrating radar (GPR) is an effective tool for detailed investigations of shallow subsurface geometry, but it provides only limited information on hydraulic properties. Magnetic resonance tomography (MRT) provides parameters such as water content (porosity) and relaxation time/hydraulic conductivity, but it suffers from resolution limits. Furthermore, it requires knowledge of subsurface electrical resistivity, which can be obtained by electrical resistivity tomography (ERT) also suffering from resolution limits. To overcome the limitations in resolution, we have incorporated GPR reflectors as structural information into the ERT and MRT data inversion. We test the methodology on a synthetic example and find improved imaging properties compared to standard inversion, particularly at greater depths, where the resolution is limited. We apply the methodology to a test site that is characterized by a complex depositional architecture. The Quaternary deposits consist of interbedded meltwater deposits (aquifers) and till (aquitards), overlain by aeolian deposits. To image the subsurface depositional architecture in three dimensions, a 200 x 250 m area was surveyed by GPR. The use of GPR constraints clearly improves the resolution and zonation of the subsurface image, which is validated by drill-core analyses. We develop a workflow to combine GPR, MRT, and ERT, leading the way to high-resolution hydro geologic models that can be used for groundwater studies.