Package sparse-msrf
--------------------

This package contains the Matlab codes required to model a
non-stationary fossil-fuel CO2 (ffCO2) emission field using wavelets
and then solve an inverse problem to estimate time-variant ffCO2
emission fields using limited measurements of ffCO2
concentrations. The theory behind these codes are in:

1. J. Ray, V. Yadav, J. Lee, A. M. Michalak, S. Lefantzi, B. van Bloemen
   Waanders and S. A. McKenna, "A multiresolution spatial parametrization
   for the estimation of fossil-fuel carbon dioxide emissions via
   atmospheric inversions", Sandia Technical Report,
   SAND2013-2919. Printed April 2013. Sandia National Laboratories,
   Livermore, CA.

2. J. Ray, V. Yadav, A. M. Michalak, B. van Bloemen Waanders and
   S.A. McKenna, "A multiresolution spatial parameterization for the
   estimation of fossil-fuel carbon dioxide emissions via atmospheric
   inversions", accepted, Geoscientific Model Development,
   2014. SAND2013-8688J. Geoscientific Model Development Discussions
   version: doi:10.5194/gmdd-7-1277-2014.

Contents:
--------

data/edgar2005.1x1 : ffCCO2 emissions, from EDGAR, for 2005, for US,
                     aggregated to 1degree resolution.

data/nlights-0.1x0.1 : radiances of lights at night, along with
                       lattitude & longitude of each grid cell

data/tower.obs: The names of sensors where ffCO2 measurements are made

data/US_coords.dat : Polygon file describe the border of the Lower48
                     states

data/vulcan.8day.1x1: Emissions from the Vulcan inventory, averaged
                      over 8 days and aggregated to 1 degree
                      resolution. For 2002. Raw data
		      at http://vulcan.project.asu.edu/

doc/ : User-manual and description of the code

examples/ex01 : Matlab codes that demonstrate how Vulcan emissions can
	        be decomposed on different types of wavelets and what
	        kind of sparsities they lead to. The codes generate
	        Figs. 2, 3 and 4 in [2].

examples/ex02 : Matlab codes that demonstrate (1) creation of the
	        Multiscale Random Field model using Haars, the setting
	        up of the inverse problem and its solution using
	        sparse reconstruction.

H-matrices/ : Directory containing the transport matrices. Unzip the
	    matrices.zip file in this folder.

msrf-selection/ : Matlab codes that allow a person to modify wavelet
		  coefficients at each level. It also has codes that
		  allow one to take forward and inverse wavelet
		  transforms on a sparsified wavelet basis set.


sparse-inversion/ Codes that solve the ffCO2 inverse problem using
		  sparse reconstruction.

incoherence/	  Codes that compute the incoherene between the MsRF and 
		  (1) transport model / H matrix and (2) standard random 
		  sampling matrices.