DOI: 10.5281/zenodo.3241467
title: Geostatistical inverse modeling with large atmospheric data: data files for a case study from OCO-2

description: The files in this data repository provide the inputs required to run an inverse modeling case study. This case study will estimate CO<sub>2</sub> fluxes across North America for July 2015 using synthetic observations that have been created to resemble observations from NASA's Orbiting Carbon Observatory 2 (OCO-2) satellite.

This data repository is specifically linked to a GitHub code repository (http://doi.org/10.5281/zenodo.3241524 or https://github.com/greenhousegaslab/geostatistical_inverse_modeling). That GitHub repository provides scripts for constructing a geostatistical inverse model that will estimate greenhouse gas fluxes or air pollution emissions using atmospheric observations. The GitHub repository includes a case study that can be run out-of-the-box; the case study provides users an opportunity to test out and explore the inverse modeling code. All of the input data files for that case study are provided for download here.

Here is a brief explanation of the different files included in this data repository, but refer to the linked GitHub repository for greater details. All of these files are in a ".mat" file that can be read into Matlab using the <em>load</em> function or can be read into R using the <em>R.matlab</em> package.


	<strong>H_1.mat through H_328.mat</strong>: These files contain the <strong>H</strong> matrix or sensitivity matrix required by the inverse model. These inputs were generated using the Stochastic Time-Inverted Lagrangian Transport (STILT) model as part of NOAA's CarbonTracker-Lagrange program (https://www.esrl.noaa.gov/gmd/ccgg/carbontracker-lagrange/). The <strong>H</strong> matrix is too large to store in a single file. We have therefore split up the matrix into 328 different files. Each file contains a vertical strip of the <strong>H</strong> matrix that corresponds to a different time period of fluxes to be estimated as part of the inverse model.
	<strong>Z.mat</strong>: This file contains the synthetic OCO-2 observations used in the case study. 
	<strong>areas_us.mat</strong>: This file lists the area of each model grid box used in the case study in units of meters<sup>2</sup>. This file only includes grid box area for model grid boxes that fall within the continental United States. We estimate CO<sub>2</sub> fluxes across terrestrial North America on a 1 degree latitude by 1 degree longitude grid as part of the case study. Each of these model grid boxes will have a different area, depending upon the latitude of that model grid box. 
	<strong>distmat.mat</strong>: This file contains a matrix that lists the distance (in kilometers) between the center of each model grid box used in the case study. 
	<strong>land_mask.mat</strong>: We only estimate CO<sub>2</sub> fluxes for terrestrial regions of North America as part of the case study. This land mask is used to convert the fluxes estimated by the inverse model to a latitude-longitude grid that can then be plotted.
description: This work is funded by NASA ROSES grant no. 80NSSC18K0976 and partially supported by NSF grant no. DMS-1720398.