Nikolaj K Christensen

USGS2021

Model structural uncertainty quantification and hydrogeophysical data integration using airborne electromagnetic data

A irborne electromagnetic (AEM) data are used to estimate large - scale model structural geometry, i.e. the spatial distribution of different lit hological units based on assumed or estimated resistivity - lithology relationships, and the uncertainty in those structures given imperfect measurements. Geophysically derived estimates of model structural uncertainty are then combined with hydrologic obse rvations to assess the impact of model structural error on hydrologic calibration and prediction errors. Using a synthetic numerical model, we describe a sequential hydrogeophysical approach that: (1) uses Bayesian Markov chain Monte Carlo (McMC) methods to produce a robust estimate of uncertainty in electrical resistivity parameter s , (2) combines geophysical parameter uncertainty estimates with borehole observations of lithology to produce probabilistic estimates of model structural uncertainty over the e ntire AEM survey area using geostatistical sequential indicator simulation algorithms, and (3) uses model structural estimates along with hydrologic observations to quantify both hydrologic parameter and prediction uncertainty using a second McMC sampling algorithm. Results of simulations will be presented that illustrate the complete workflow from geophysical parameter uncertainty analysis to the impact of model structural uncertainty on hydrologic parameter estimates.

Conference Paper↗

USGS2017

Generation of 3-D hydrostratigraphic zones from dense airborne electromagnetic data to assess groundwater model prediction error

We present a new methodology to combine spatially dense high-resolution airborne electromagnetic (AEM) data and sparse borehole information to construct multiple plausible geological structures using a stochastic approach. The method developed allows for quantification of the performance of groundwater models built from different geological realizations of structure. Multiple structural realizations are generated using geostatistical Monte Carlo simulations that treat sparse borehole lithological observations as hard data and dense geophysically derived structural probabilities as soft data. Each structural model is used to define 3-D hydrostratigraphical zones of a groundwater model, and the hydraulic parameter values of the zones are estimated by using nonlinear regression to fit hydrological data (hydraulic head and river discharge measurements). Use of the methodology is demonstrated for a synthetic domain having structures of categorical deposits consisting of sand, silt, or clay. It is shown that using dense AEM data with the methodology can significantly improve the estimated accuracy of the sediment distribution as compared to when borehole data are used alone. It is also shown that this use of AEM data can improve the predictive capability of a calibrated groundwater model that uses the geological structures as zones. However, such structural models will always contain errors because even with dense AEM data it is not possible to perfectly resolve the structures of a groundwater system. It is shown that when using such erroneous structures in a groundwater model, they can lead to biased parameter estimates and biased model predictions, therefore impairing the model's predictive capability.

Water Resources Research↗