Crustal constraint through complete model space screening for diverse geophysical datasets facilitated by emulation

Roberts, A. W., Hobbs, R. W., Goldstein, M., Moorkamp, Max, Jegen-Kulcsar, Marion D. and Heincke, Björn (2012) Crustal constraint through complete model space screening for diverse geophysical datasets facilitated by emulation Tectonophysics, 572/573 . pp. 47-63. DOI 10.1016/j.tecto.2012.03.006.

[img] Text
Roberts.pdf - Published Version
Restricted to Registered users only

Download (2248Kb) | Contact

Supplementary data:


Deep crustal constraint is often carried out using deterministic inverse methods, sometimes using seismic refraction, gravity and electromagnetic datasets in a complementary or "joint" scheme. With increasingly powerful parallel computer systems it is now possible to apply joint inversion schemes to derive an optimum model from diverse input data. These methods are highly effective where the uncertainty in the system is small. However, given the complex nature of these schemes it is often difficult to discern the uniqueness of the output model given the noise in the data, and the application of necessary regularization and weighting in the inversion process means that the extent of user prejudice pertaining to the final result may be unclear. We can rigorously address the subject of uncertainty using standard statistical tools but these methods also become less feasible if the prior model space is large or the forward simulations are computationally expensive. We present a simple Monte Carlo scheme to screen model space in a fully joint fashion, in which we replace the forward simulation with a fast and uncertainty-calibrated mathematical function, or emulator. This emulator is used as a proxy to run the very large number of models necessary to fully explore the plausible model space. We develop the method using a simple synthetic dataset then demonstrate its use on a joint data set comprising first-arrival seismic refraction. MT and scalar gravity data over a diapiric salt body. This study demonstrates both the value of a forward Monte Carlo approach (as distinct from a search-based or conventional inverse approach) in incorporating all kinds of uncertainty in the modelling process, exploring the entire model space, and shows the potential value of applying emulator technology throughout geophysics. Though the target here is relatively shallow, the methodology can be readily extended to address the whole crust.

Document Type: Article
Additional Information: WOS:000310402200006
Keywords: Bayesian; Statistical methods; Emulation; Joint inversion; Salt diapir; Crustal imaging; MONTE-CARLO METHODS; JOINT INVERSION; VELOCITY; COMPUTATION; GRAVITY
Research affiliation: OceanRep > GEOMAR > FB4 Dynamics of the Ocean Floor > FB4-GDY Marine Geodynamics
OceanRep > The Future Ocean - Cluster of Excellence
Refereed: Yes
DOI etc.: 10.1016/j.tecto.2012.03.006
ISSN: 0040-1951
Projects: Future Ocean
Date Deposited: 06 Dec 2012 10:14
Last Modified: 15 Jun 2016 12:01

Actions (login required)

View Item View Item

Document Downloads

More statistics for this item...