Bayesian joint inversion of controlled source electromagnetic and magnetotelluric data to image freshwater aquifer offshore New Jersey

Abstract

SUMMARY Joint inversion of multiple electromagnetic data sets, such as controlled source electromagnetic and magnetotelluric data, has the potential to significantly reduce uncertainty in the inverted electrical resistivity when the two data sets contain complementary information about the subsurface. However, evaluating quantitatively the model uncertainty reduction is made difficult by the fact that conventional inversion methods—using gradients and model regularization—typically produce just one model, with no associated estimate of model parameter uncertainty. Bayesian inverse methods can provide quantitative estimates of inverted model parameter uncertainty by generating an ensemble of models, sampled proportional to data fit. The resulting posterior distribution represents a combination of a priori assumptions about the model parameters and information contained in field data. Bayesian inversion is therefore able to quantify the impact of jointly inverting multiple data sets by using the statistical information contained in the posterior distribution. We illustrate, for synthetic data generated from a simple 1-D model, the shape of parameter space compatible with controlled source electromagnetic and magnetotelluric data, separately and jointly. We also demonstrate that when data sets contain complementary information about the model, the region of parameter space compatible with the joint data set is less than or equal to the intersection of the regions compatible with the individual data sets. We adapt a trans-dimensional Markov chain Monte Carlo algorithm for jointly inverting multiple electromagnetic data sets for 1-D earth models and apply it to surface-towed controlled source electromagnetic and magnetotelluric data collected offshore New Jersey, USA, to evaluate the extent of a low salinity aquifer within the continental shelf. Our inversion results identify a region of high resistivity of varying depth and thickness in the upper 500 m of the continental shelf, corroborating results from a previous study that used regularized, gradient-based inversion methods. We evaluate the joint model parameter uncertainty in comparison to the uncertainty obtained from the individual data sets and demonstrate quantitatively that joint inversion offers reduced uncertainty. In addition, we show how the Bayesian model ensemble can subsequently be used to derive uncertainty estimates of pore water salinity within the low salinity aquifer.

Seismic Tomography Using Variational Inference Methods

“…Variational inference provides a different way to solve a Bayesian inference problem: Within a predefined family of probability distributions, one seeks an optimal approximation to a target distribution, which in this case is the Bayesian posterior pdf. This is achieved by minimizing the Kullback‐Leibler (KL) divergence (Kullback & Leibler, )—one possible measure of the difference between two given pdfs (Blatter et al, ), in our case the difference between approximate and target pdfs (Bishop, ; Blei et al, ). Since the method casts the inference problem into an optimization problem, it can be computationally more efficient than either MC sampling or neural network methods and provides better scaling to higher‐dimensional problems.…”

Marine electrical imaging reveals novel freshwater transport mechanism in Hawai‘i

“…The substitution of conductive seawater with freshwater will increase the electrical resistivity of any geological formation ( 28 ). Various marine CSEM techniques were proven successful in imaging the electrical structure of continuous offshore freshened groundwater in different coastal sediment environments such as those of New Zealand ( 29 ), the U.S. Atlantic coast ( 2 , 30 ), and nearshore Israel ( 31 , 32 ). In volcanic geology, where seawater-saturated basalts have resistivities of <10 ohm·m ( 33 ), submarine freshwater-saturated basalts will manifest as 600– to 1100–ohm·m resistive anomalies ( 34 ) embedded in a conductive background of seawater-saturated basalts.…”

Development of the Inversion Method for Transient Electromagnetic Data

“…It cannot only provide the estimation of model parameters that are often provided by deterministic inversion, but it also can determine the uncertainty of parameter estimation under the influence of multiple solutions of the inverse problem, observation error, and noise, so as to provide more information for data interpretation. In addition, because the prior distribution of Bayesian inversion makes it easy to introduce multiple model spaces, it is suitable for joint inversion of multiple parameters, such as parallel sampling of resistivity and charge rate parameters in TEM data [36], joint inversion of artificial-source electromagnetic data and MT data [65].…”