The Reduced-Order Hybrid Monte Carlo Sampling Smoother

Abstract

Hybrid Monte-Carlo (HMC) sampling smoother is a fully non-Gaussian four-dimensional data assimilation algorithm that works by directly sampling the posterior distribution formulated in the Bayesian framework. The smoother in its original formulation is computationally expensive due to the intrinsic requirement of running the forward and adjoint models repeatedly. Here we present computationally efficient versions of the HMC sampling smoother based on reduced-order approximations of the underlying model dynamics. The schemes developed herein are tested numerically using the shallow-water equations model on Cartesian coordinates. The results reveal that the reduced-order versions of the smoother are capable of accurately capturing the posterior probability density, while being significantly faster than the original full order formulation.

32 pages, 2 figures

Cited by in corpus (8)

• Goal-Oriented Optimal Design of Experiments for Large-Scale Bayesian Linear Inverse Problems
• Structure-Preserving Galerkin POD Reduced-Order Modeling of Hamiltonian Systems
• Stochastic Learning Approach to Binary Optimization for Optimal Design of Experiments
• A Bayesian Approach to Multivariate Adaptive Localization in Ensemble-Based Data Assimilation with Time-Dependent Extensions
• An Optimal Experimental Design Framework for Adaptive Inflation and Covariance Localization for Ensemble Filters
• A Machine Learning Approach to Adaptive Covariance Localization
• Medical Image Retrieval Based On the Parallelization of the Cluster Sampling Algorithm
• Investigation of Nonlinear Model Order Reduction of the Quasigeostrophic Equations through a Physics-Informed Convolutional Autoencoder