Loading Events

« All Events

  • This event has passed.

Markov Chain Monte Carlo for High-dimensional, Nonlinear Problems in Earth Science | HDSI Seminar Series

October 17 @ 1:00 pm - 2:00 pm

Abstract: Earth science nearly always requires estimating models, or model parameters, from data. This could mean to infer the state of the southern ocean from ARGO floats, to compute the state of our atmosphere based on atmospheric observations of the past six hours, or to construct a resistivity model of the Earth’s subsurface from electromagnetic data. All these problems have in common that the number of unknowns is large (millions to hundreds of millions) and that the underlying processes are nonlinear. The problems also all have in common that they can be formulated as the problem of drawing samples from a high-dimensional Bayesian posterior distribution.
Due to the nonlinearity, Markov chain Monte Carlo (MCMC) is a good candidate for the numerical solution of such “inverse problems.” But MCMC is known to be slow when the number of unknowns is large. In this talk, I will argue that an unbiased solution of nonlinear, high-dimensional problems is and remains infeasible, but one can construct efficient and accurate biased estimators that are feasible to apply to high-dimensional problems. I will show examples of biased estimators in action and invert electromagnetic data using an approximate MCMC sampling algorithm called the RTO-TKO (randomize-then-optimize — technical-kock-out).


October 17
1:00 pm - 2:00 pm
Event Categories:




SDSC, the Auditorium
San Diego Supercomputer Center, 10100 Hopkins Dr
La Jolla, 92093 United States
+ Google Map