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.

GRID-SITES: Gridded Solar Iterative Temperature Emission Solver for Fast DEM Inversion

SITES: Solar Iterative Temperature Emission Solver for Differential Emission Measure Inversion of EUV Observations

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

Fast Differential Emission Measure Inversion of Solar Coronal Data

Differential Emission Measures from the Regularized Inversion of Hinode and SDO data

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

A method for global inversion of multi-resolution solar data

Benchmark Test of Differential Emission Measure Codes and Multi-Thermal Energies in Solar Active Regions

Optimizing Grid Computing Configuration and Scheduling for Geospatial Analysis: an Example with Interpolating DEM

Solar multiobject multiframe blind deconvolution with a spatially variant convolution neural emulator

Fast Grid Emissions Sensitivities using Parallel Decentralized Implicit Differentiation

Solar multi-object multi-frame blind deconvolution with a spatially variant convolution neural emulator

Probabilistic Post-Processing of Gridded Atmospheric Variables and Its Application to Site Adaptation of Shortwave Solar Radiation

Solar potential analysis over Indian cities using high-resolution satellite imagery and DEM

Computationally efficient methods for large-scale atmospheric inverse modeling

A SART-Based Iterative Inversion Methodology to Infer the Solar Rotation Rate from Global Helioseismic Data

A high throughput geocomputing system for remote sensing quantitative retrieval and a case study

Grid-enabled High-Performance Quantitative Aerosol Retrieval from Remotely Sensed Data

A fast 3-D inversion for airborne EM data using pre-conditioned stochastic gradient descent

Application of Supervised Descent Method to Transient Electromagnetic Data Inversion