Premium
Fast Stokes Flow Simulations for Geophysical‐Geodynamic Inverse Problems and Sensitivity Analyses Based On Reduced Order Modeling
Author(s) -
OrtegaGelabert O.,
Zlotnik S.,
Afonso J. C.,
Díez P.
Publication year - 2020
Publication title -
journal of geophysical research: solid earth
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.983
H-Index - 232
eISSN - 2169-9356
pISSN - 2169-9313
DOI - 10.1029/2019jb018314
Subject(s) - markov chain monte carlo , inverse problem , uncertainty quantification , data assimilation , computer science , inversion (geology) , algorithm , surrogate model , sensitivity (control systems) , monte carlo method , computation , mathematical optimization , mathematics , bayesian probability , geophysics , geology , physics , machine learning , artificial intelligence , mathematical analysis , paleontology , statistics , structural basin , electronic engineering , meteorology , engineering
Markov chain Monte Carlo (MCMC) methods have become standard in Bayesian inference and multi‐observable inversions in almost every discipline of the Earth sciences. In the case of geodynamic and/or coupled geophysical‐geodynamic inverse problems, however, the computational cost associated with the solution of large‐scale 3‐D Stokes forward problems has rendered probabilistic formulations impractical. Here we present a novel and extremely efficient method to produce ultrafast solutions of the 3‐D Stokes problem for MCMC simulations. Our approach combines the individual benefits of Reduced Basis techniques, goal‐oriented error formulations, and MCMC algorithms to produce an accurate and computationally efficient surrogate for the forward problem. Importantly, the surrogate adapts itself during the MCMC simulation according to the history of the chain and the goals of the inversion. This maximizes the efficiency of the forward problem and removes the need for preinversion off‐line computations to build a surrogate. We demonstrate the benefits and limitations of the method with several numerical examples and show that in all cases the computational cost is of the order of <1% compared to a traditional MCMC approach. The method is general enough to be applied to a range of problems, including uncertainty quantification/propagation, adjoint‐based geodynamic inversions, sensitivity analyses in mantle convection problems, and in the creating surrogate models for complex forward problems (e.g., heat transfer, seismic tomography, and magnetotellurics).