Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification

Collocated magnetotelluric (MT) and seismic profiling is emerging as a necessary combined approach for deep and near-surface imaging but the resulting experimental data are typically interpreted separately since no production programs exist for multidimensional joint inversion of MT and seismic data. We present a joint 2-D inversion approach for imaging collocated MT and seismic refraction data with cross-gradient structural constraints. We describe the main features of the algorithm and first apply it to synthetic data generated for a hypothetical complex geological model. For the synthetic data, we find that the scheme leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and seismic velocity. We apply the scheme to near-surface field data to test the consistency of a previously suggested resistivity-velocity interrelationship and its potential use for subsurface lithofacies discrimination or structural classification. The MT-seismic relationship is found to be in excellent accord with that derived previously for DC resistivity and seismic data set at the test site. Our results suggest that joint MT-seismic cross-gradient imaging leads to improved characterization of heterogeneous geological targets at near-surface to mantle depths.