Abstract:Abstract. Atmospheric inverse modeling describes the process of estimating greenhouse gas fluxes or air pollution emissions at the Earth's surface using observations of these gases collected in the atmosphere. The launch of new satellites, the expansion of surface observation networks, and a desire for more detailed maps of surface fluxes has yielded numerous computational and statistical challenges for standard inverse modeling frameworks that were often originally designed with much smaller data sets in mind. In this article, we discuss computationally efficient methods for large-scale atmospheric inverse modeling and focus on addressing some of the main computational and practical challenges. We develop generalized hybrid projection methods, which are iterative methods for solving large-scale inverse problems, and specifically we focus on the case of estimating surface fluxes. These algorithms confer several advantages. They are efficient, in part because they converge quickly, they exploit efficient matrix-vector multiplications, and do not require inverting any matrices. These methods are also robust because they can accurately reconstruct surface fluxes, they are automatic since regularization or covariance matrix parameters and stopping criteria can be determined as part of the iterative algorithm, and they are flexible because they can be paired with many different types of atmospheric models. We demonstrate the benefits of generalized hybrid methods with a case study from NASA's Orbiting Carbon Observatory 2 (OCO-2) satellite. We then address the more challenging problem of solving the inverse model when the mean of the surface fluxes is not known a priori; we do so by reformulating the problem, thereby extending the applicability of hybrid projection methods to include hierarchical priors. We further show that by exploiting mathematical relations provided by the generalized hybrid method, we can efficiently calculate an approximate posterior variance, thereby providing uncertainty information.

A Generalized Variable Projection Algorithm for Least Squares Problems in Atmospheric Remote Sensing

Sparse Representation Based on Projection Method in Online Least Squares Support Vector Machines

An Iterative Implementation of Variable Projection for Separable Nonlinear Optimization Problems

An efficient variable projection formulation for separable nonlinear least squares problems.

Separable nonlinear least squares: the variable projection method and its applications

A Regularized Variable Projection Algorithm for Separable Nonlinear Least-Squares Problems.

Variable projection algorithms with sparse constraint for separable nonlinear models

An Efficient Algorithm for Ill-Conditioned Separable Nonlinear Least Squares

Variable Projection Algorithms: Theoretical Insights and A Novel Approach for Problems with Large Residual

Computationally efficient methods for large-scale atmospheric inverse modeling

Variable projection methods for separable nonlinear inverse problems with general-form Tikhonov regularization

A New Least Squares Method Based Line Generalization in GIS.

Robust Variable Projection Algorithm for the Identification of Separable Nonlinear Models

Modified Gram–Schmidt Method-Based Variable Projection Algorithm for Separable Nonlinear Models

Extensions and uses of the variable projection algorithm for solving nonlinear least squares problems

Convergence Analysis of a Variable Projection Method for Regularized Separable Nonlinear Inverse Problems

Efficient Parameters Estimation Method for the Separable Nonlinear Least Squares Problem

Airborne gravimetry data sparse reconstruction via L1-norm convex quadratic programming

An Iteratively Approximated Gradient Projection Algorithm for Sparse Signal Reconstruction

Parallel Inexact Levenberg-Marquardt Method for Nearly-Separable Nonlinear Least Squares

Accelerating Variable Projection Algorithm by Aitken Technique.