H11M-03
GPR waveguide and full-waveform inversion for the hydrogeophysical characterization of soils and aquifers

Monday, 14 December 2015: 08:30
3016 (Moscone West)
Jan Van Der Kruk1, Adam R Mangel2, Nils Gueting3, Sebastian Busch4, Anja Klotzsche3, Stephen M Moysey5, Johan Alexander Huisman6 and Harry Vereecken7, (1)Agrosphere Institute (IBG-3), Forschungszentrum Jülich, Jülich, Germany, (2)Clemson University, Environmental Engineering and Earth Science, Clemson, SC, United States, (3)Agrosphere Institute (IBG-3) Forschungszentrum Jülich, Deutschland, Germany, (4)Forschungszentrum Jülich, Jülich, Germany, (5)Clemson University, Clemson, SC, United States, (6)Forschungszentrum Jülich, Agrosphere (IBG 3), Jülich, Germany, (7)Forschungszentrum Julich GmbH, Julich, Germany
Abstract:
Ground-penetrating radar (GPR) is a powerful tool for characterizing hydrologic processes. Coupled hydrogeophysical inversion can be used to invert time-lapse GPR data to obtain soil hydraulic properties. However, the inversion will fail if the hydrological model is not appropriately conceptualized, or when the ray-based methods that are often used to limit the computation time are not appropriate. Ray-based approaches cannot be used in the case of wetting fronts in the vadose zone due to precipitation/infiltration or thawing events, where low-velocity layers with high water content can trap the GPR waves and act as a waveguide such that multiple internal reflections cause dispersion. Utilizing the dispersion, manifested as a singled appearance of arrivals on multi-offset data, we can invert for waveguide properties. Single- or two-layer waveguide inversion approaches return average water contents but are incapable of representing the gradational nature of the water content distribution of the shallow subsurface. Recently, a shuffled complex evolution algorithm was implemented that used a piece-wise linear function to closely match the shallow gradational water content profile present for early-time infiltration events for different n-parameter values of the Mualem–van Genuchten equation using HYDRUS-1D. Similar waveguide phenomena can arise in the presence of high porosity layers in saturated aquifers for crosshole GPR. These high porosity layers result in significant late arrival high-amplitude elongated wave trains that were detected in crosshole data from the Widen, Boise and Krauthausen aquifers. Full-waveform inversion of these data that uses an accurate forward model to explain the measured data is able to return decimeter scale resolution of the porosity based on the GPR velocity. With increasing computer power, it is expected that these advanced processing can be soon included in the coupled inversion approaches.