Abstract:<p>Inverse modeling of large-scale spatially variable parameters fields at fine resolution can be reduced to estimating the projections on dominant principal components of the underlying parameter fields based on principal component analysis of the spatial covariance. However, for unknown or biased prior structural parameters of spatial covariance models, an iterative procedure consisting of two successive steps is usually implemented, i.e., estimation of spatial covariance followed by estimation of the spatially variable parameter fields conditional on the spatial covariance and observations. In this study, we develop an iterative, computationally efficient method to update dominant principal components for nonlinear inverse problems of large-scale spatial fields that adaptively corrects the bias from the initially defined prior spatial covariance. Our algorithm involves two-layer iterations: the inner iteration is to obtain the best estimates of projections on given retained principal components, and the outer iteration implements an efficient rank-one updating to correct the retained principal components using the posterior covariance associated with the best estimates of the projections. Numerical experiments show that inversion results can be significantly improved for large-scale inverse problems with biased structural parameters for spatial covariance. The experiment results show that the iterative correction is essentially to match the distribution patterns of the spatially correlated parameter field with its most dominant principal components. We also investigate the performance of the developed method under different biased covariance model initialization, including model type bias, variance bias and correlation length bias. The correction cannot fundamentally change the smoothness defined by the covariance model type, but can still describe major distribution patterns including anisotropy. Biased variance can be corrected and yields similar best estimates and variance maps, and biased correlation length can be corrected within an applicable range.</p>

A joint reconstruction and model selection approach for large-scale linear inverse modeling (msHyBR v2)

Computationally efficient methods for large-scale atmospheric inverse modeling

Surrogate-Based Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Structural Error.

Efficient iterative methods for hyperparameter estimation in large-scale linear inverse problems

Computational Methods for Large-Scale Inverse Problems: A Survey on Hybrid Projection Methods

Efficient generalized Golub-Kahan based methods for dynamic inverse problems

Bayesian inverse modeling of large-scale spatial fields on iteratively corrected principal components

Generalized hybrid iterative methods for large-scale Bayesian inverse problems

Inexact Generalized Golub-Kahan Methods for Large-Scale Bayesian Inverse Problems

A greedy sensor selection algorithm for hyperparameterized linear Bayesian inverse problems with correlated noise models

A Quasi‐Newton Reformulated Geostatistical Approach on Reduced Dimensions for Large‐Dimensional Inverse Problems

Solving High-Dimensional Inverse Problems with Auxiliary Uncertainty via Operator Learning with Limited Data

Multiscale model reduction method for Bayesian inverse problems of subsurface flow

A Unified Dynamic Approach to Sparse Model Selection

Improve Unscented Kalman Inversion With Low-Rank Approximation and Reduced-Order Model

Multi-Objective Bayesian Optimisation and Joint Inversion for Active Sensor Fusion

Solving Inverse Problems with Model Mismatch using Untrained Neural Networks within Model-based Architectures

Joint learning of nonlinear representation and projection for fast constrained MRSI reconstruction

A novel, finite-element-based framework for sparse data solution reconstruction and multiple choices

Surrogate-Based Bayesian Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Approximation Error

Generalized sparsity-promoting solvers for Bayesian inverse problems: Versatile sparsifying transforms and unknown noise variances