We are excited to invite you to the 15th seminar of the 2025 Geology and Geophysics Seminar Series, featuring Fabrizio Magrini, who is a Research Fellow/Lecturer at RSES, Australian National University. Fabrizio will be presenting on “Bayesian Inference with Transport Maps“. In this engaging talk, he will show how advances in computational statistics and optimal transport improve MCMC efficiency in Bayesian geophysical inversion, sharing fresh insights and innovative methods in geoscience.
Date: October 22, 2025
Time: 11:00 a.m. – 12:00 p.m., Sydney Time
Location: Room 449 (Conference Room), Madsen Building (F09), School of Geosciences
or Online (Join via zoom)
We look forward to seeing you there in person or joining us online!
https://uni-sydney.zoom.us/j/84796471781?from=addon
Bayesian Inference with Transport Maps
Abstract
Efficient Markov chain Monte Carlo (MCMC) sampling from posterior distributions remains a central challenge in Bayesian geophysical inversion. Recent developments in computational statistics and optimal transport suggest that MCMC efficiency can be improved by reparameterising the sampling problem — specifically, by learning an invertible mapping that recasts the target distribution onto a simpler reference distribution. Here, we introduce a Metropolis–Hastings framework that leverages transport maps parameterised by invertible neural networks. These maps are trained on preliminary MCMC samples from the target distribution and used to propose new samples in a fixed reference space, where proposal design is independent of the target’s structure. The proposed samples are transformed back to the target space via the inverse map, and accepted or rejected according to a modified Metropolis–Hastings criterion. As sampling proceeds, the transport maps are updated, yielding proposals increasingly well adapted to the shape of the target distribution. Across a suite of numerical tests — including a 2-D Rosenbrock distribution, a 3-D earthquake location problem, and Gaussian mixtures up to 16 dimensions — transport-map-driven samplers consistently outperform standard MCMC, reducing integrated autocorrelation times by factors of 2.5 to over 6 (or equivalently, yielding sample sets 2.5–6 times larger for the same number of forward evaluations). This improvement comes at the non-negligible cost of training one or more transport maps, which we quantify systematically. We also provide a quantitative criterion for weighing training cost against sampling speed-up. This shows that transport-map MCMC is advantageous whenever the forward problem is nontrivial, making it a promising approach for Bayesian sampling in geophysics and beyond.
Graphical Abstract
