Premium
Inverse modeling of unsaturated flow parameters using dynamic geological structure conditioned by GPR tomography
Author(s) -
Farmani M. Bagher,
Kitterød NilsOtto,
Keers Henk
Publication year - 2008
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2007wr006251
Subject(s) - vadose zone , geology , ground penetrating radar , inversion (geology) , tomography , inverse problem , infiltration (hvac) , geophysics , soil science , mineralogy , geotechnical engineering , geomorphology , groundwater , radar , meteorology , mathematics , computer science , telecommunications , mathematical analysis , physics , optics , structural basin
A method is presented to estimate flow parameters and geological structure in the vadose zone by combining time‐lapse Ground Penetrating Radar (GPR) traveltime tomography and inverse flow modeling. The traveltime tomography is used to determine the spatial electromagnetic velocity distribution in the vadose zone. These time‐lapse velocity images are converted to time‐lapse volumetric soil water content images using petrophysical relationships. Subsequently, the water content images are used as constraints in the flow inversion. The influence of the tomographic artifacts on the flow inversion is minimized by assigning weights that are proportional to the ray coverage. Our flow inversion algorithm estimates the flow parameters and calibrates the geological structure. The geological structure is defined using a set of control points, the positions of which can be modified during the inversion. After the inversion, the final geological and flow model are used to compute GPR traveltimes to check the consistency between these computed traveltimes and the observed traveltimes.The method is first tested on two synthetic models (a steady state and a transient flow models). Subsequently, the method is applied to characterize a real vadose zone at Oslo Airport Gardermoen, Norway, during the snowmelt in 2005. The flow inversion method is applied to locate and quantify the main geological layers at the site. In particular the inversion method identifies and estimates the location and properties of thin dipping layers with relatively low‐permeability. The flow model is cross validated using an independent infiltration event.