Computationally efficient methods for large-scale atmospheric inverse modeling

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 have 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 they do not require inversion of 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.

Quantifying methane emissions from the global scale down to point sources using satellite observations of atmospheric methane

“…Numerical solution for min(J (x)) using the adjoint of the atmospheric transport model or other variational methods optimizes a state vector of any dimension by avoiding explicit construction of the full Jacobian matrix K and may use various procedures to estimate Ŝ (Bousserez et al, 2015;Cho et al, 2022). Analytical solution provides a closed-form expression for Ŝ but requires the computationally expensive construction of K column-by-column with n perturbation runs of the atmospheric transport model.…”

Assimilating the dynamic spatial gradient of a bottom-up carbon flux estimation as a unique observation in COLA (v2.0)

“…can be inferred from atmospheric CO2 measurements using "top-down" (hereafter quotation marks will be omitted) techniques of the Bayesian synthesis (e.g., Rodenbeck et al, 2003;Zammit-Mangion et al, 2022;Cho et al, 2022) and data assimilation (DA) techniques (e.g., Peters et al, 2007;Feng et al, 2009;Chevallier et al, 2010;J. Liu et al, 2014;.…”