Scalable Bayesian transport maps for high-dimensional non-Gaussian spatial fields

by   Matthias Katzfuss, et al.

A multivariate distribution can be described by a triangular transport map from the target distribution to a simple reference distribution. We propose Bayesian nonparametric inference on the transport map by modeling its components using Gaussian processes. This enables regularization and accounting for uncertainty in the map estimation, while still resulting in a closed-form and invertible posterior map. We then focus on inferring the distribution of a nonstationary spatial field from a small number of replicates. We develop specific transport-map priors that are highly flexible and are motivated by the behavior of a large class of stochastic processes. Our approach is scalable to high-dimensional fields due to data-dependent sparsity and parallel computations. We also discuss extensions, including Dirichlet process mixtures for marginal non-Gaussianity. We present numerical results to demonstrate the accuracy, scalability, and usefulness of our methods, including statistical emulation of non-Gaussian climate-model output.


page 2

page 9

page 14

page 18


Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

Multivariate spatially-oriented data sets are prevalent in the environme...

Greedy inference with layers of lazy maps

We propose a framework for the greedy approximation of high-dimensional ...

Transport Monte Carlo

In Bayesian inference, transport map is a promising alternative to the c...

Simple approximate MAP Inference for Dirichlet processes

The Dirichlet process mixture (DPM) is a ubiquitous, flexible Bayesian n...

Bayesian nonstationary and nonparametric covariance estimation for large spatial data

In spatial statistics, it is often assumed that the spatial field of int...

A Distributed Framework for the Construction of Transport Maps

The need to reason about uncertainty in large, complex, and multi-modal ...

A Scalable Partitioned Approach to Model Massive Nonstationary Non-Gaussian Spatial Datasets

Nonstationary non-Gaussian spatial data are common in many disciplines, ...

Code Repositories


Scalable Bayesian transport maps for high-dimensional non-Gaussian spatial fields

view repo