Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

A sparse reconstruction method for the estimation of multi-resolution emission fields via atmospheric inversion

Journal Article · · Geoscientific Model Development (Online)
Atmospheric inversions are frequently used to estimate fluxes of atmospheric greenhouse gases (e.g., biospheric CO2 flux fields) at Earth's surface. These inversions typically assume that flux departures from a prior model are spatially smoothly varying, which are then modeled using a multi-variate Gaussian. When the field being estimated is spatially rough, multi-variate Gaussian models are difficult to construct and a wavelet-based field model may be more suitable. Unfortunately, such models are very high dimensional and are most conveniently used when the estimation method can simultaneously perform data-driven model simplification (removal of model parameters that cannot be reliably estimated) and fitting. Such sparse reconstruction methods are typically not used in atmospheric inversions. In this work, we devise a sparse reconstruction method, and illustrate it in an idealized atmospheric inversion problem for the estimation of fossil fuel CO2 (ffCO2) emissions in the lower 48 states of the USA. Our new method is based on stagewise orthogonal matching pursuit (StOMP), a method used to reconstruct compressively sensed images. Our adaptations bestow three properties to the sparse reconstruction procedure which are useful in atmospheric inversions. We have modified StOMP to incorporate prior information on the emission field being estimated and to enforce non-negativity on the estimated field. Finally, though based on wavelets, our method allows for the estimation of fields in non-rectangular geometries, e.g., emission fields inside geographical and political boundaries. Our idealized inversions use a recently developed multi-resolution (i.e., wavelet-based) random field model developed for ffCO2 emissions and synthetic observations of ffCO2 concentrations from a limited set of measurement sites. We find that our method for limiting the estimated field within an irregularly shaped region is about a factor of 10 faster than conventional approaches. It also reduces the overall computational cost by a factor of 2. Further, the sparse reconstruction scheme imposes non-negativity without introducing strong nonlinearities, such as those introduced by employing log-transformed fields, and thus reaps the benefits of simplicity and computational speed that are characteristic of linear inverse problems.
Research Organization:
Sandia National Laboratories (SNL), Albuquerque, NM, and Livermore, CA (United States); Sandia National Laboratories (SNL-CA), Livermore, CA (United States); Sandia National Laboratories, Albuquerque, NM (United States)
Sponsoring Organization:
USDOE; USDOE National Nuclear Security Administration (NNSA)
Grant/Contract Number:
AC04-94AL85000
OSTI ID:
1197811
Alternate ID(s):
OSTI ID: 1214233
OSTI ID: 1182960
Report Number(s):
SAND--2014-15988J
Journal Information:
Geoscientific Model Development (Online), Journal Name: Geoscientific Model Development (Online) Journal Issue: 4 Vol. 8; ISSN 1991-9603
Publisher:
Copernicus Publications, EGUCopyright Statement
Country of Publication:
Germany
Language:
English

References (41)

Noiselets journal January 2001
Sparsity-Promoting Solution of Subsurface Flow Model Calibration Inverse Problems book January 2013
Atmospheric inversions for estimating CO2 fluxes: methods and perspectives book January 2010
A Simple Proof of the Restricted Isometry Property for Random Matrices journal January 2008
A sparse Bayesian framework for conditioning uncertain geologic models to nonlinear flow measurements journal September 2010
Locating and quantifying gas emission sources using remotely obtained concentration data journal August 2013
A time-split nonhydrostatic atmospheric model for weather research and forecasting applications journal March 2008
Breakdown of equivalence between the minimal ℓ 1 -norm solution and the sparsest solution journal March 2006
High Resolution Fossil Fuel Combustion CO 2 Emission Fluxes for the United States journal July 2009
A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model journal January 2003
A new global gridded data set of CO 2 emissions from fossil fuel combustion: Methodology and evaluation journal January 2010
Space-based observations of megacity carbon dioxide: SPACE-BASED OBSERVATIONS OF MEGACITY CO journal September 2012
Terrestrial ecosystem production: A process model based on global satellite and surface data journal December 1993
The artificial night sky brightness mapped from DMSP satellite Operational Linescan System measurements journal November 2000
Assessment of ground-based atmospheric observations for verification of greenhouse gas emissions from an urban region journal May 2012
Recent trends in global greenhouse gas emissions:regional trends 1970–2000 and spatial distributionof key sources in 2000 journal June 2005
Sparsity and incoherence in compressive sampling journal April 2007
Matching pursuits with time-frequency dictionaries journal January 1993
The Fukushima inverse problem conference May 2013
Imaging via Compressive Sampling journal March 2008
An Introduction To Compressive Sampling journal March 2008
Bayesian Compressive Sensing Using Laplace Priors journal January 2010
Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information journal February 2006
Compressed sensing journal April 2006
Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? journal January 2006
Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit journal December 2007
Model-Based Compressive Sensing journal April 2010
Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit journal February 2012
Bayesian Compressive Sensing journal June 2008
Exploiting Structure in Wavelet-Based Bayesian Compressive Sensing journal September 2009
Tomographic inversion using ℓ 1 -norm regularization of wavelet coefficients journal July 2007
Regularization of linear and non-linear geophysical ill-posed problems with joint sparsity constraints journal February 2010
Solving or resolving global tomographic models with spherical wavelets, and the scale and sparsity of seismic heterogeneity: Wavelets, sparsity and seismic tomography journal September 2011
Signal reconstruction using sparse tree representations conference August 2005
Practical compressive sensing with Toeplitz and circulant matrices conference July 2010
A very high-resolution (1 km×1 km) global fossil fuel CO 2 emission inventory derived using a point source database and satellite observations of nighttime lights journal January 2011
Assessment of fossil fuel carbon dioxide and other anthropogenic trace gas emissions from airborne measurements over Sacramento, California in spring 2009 journal January 2011
On observational and modelling strategies targeted at regional carbon exchange over continents journal January 2009
A synthesis of carbon dioxide emissions from fossil-fuel combustion journal January 2012
North American CO 2 exchange: inter-comparison of modeled estimates with results from a fine-scale atmospheric inversion journal January 2012
A multiresolution spatial parameterization for the estimation of fossil-fuel carbon dioxide emissions via atmospheric inversions journal January 2014

Similar Records

sparse-msrf:A package for sparse modeling and estimation of fossil-fuel CO2 emission fields
Software · Mon Oct 06 00:00:00 EDT 2014 · OSTI ID:1231949