Data Description

Parameters

This SMAP data product contains daily estimates of global ecosystem productivity, including net ecosystem exchange (NEE), gross primary production (GPP), heterotrophic respiration (Rh), and soil organic carbon (SOC), along with quality control metrics. The NEE of CO₂ with the atmosphere is a fundamental measure of the balance between carbon uptake by vegetation GPP, and carbon losses through autotrophic respiration (Ra) and heterotrophic respiration (Rh). The sum of Ra and Rh defines the total ecosystem respiration rate (Rtot), which encompasses most of the annual terrestrial CO₂ efflux to the atmosphere. All parameters are expressed in units of g C m^-2 day^-1. The CO₂ flux state variable outputs are provided in SPL4CMDL files as eight vegetated land-cover classes called Plant Function Types (PFTs). For example, the CO₂ flux state variable outputs are provided in NEE/nee_pft{1..8}_mean, GPP/gpp_pft{1..8}_mean, and RH/gpp_pft{1..8}_mean. The soil carbon pool state variables are provided in SOC/soc_pft{1..8}_mean. Refer to Table 1 for descriptions of the eight PFTs.

Table 1. Plant Function Type (PFT) Classifier Summary

PFT Class	PFT Class Label	PFT Code	PFT Description	PFT Class used in SPL4CMDL
N/A	Water	0	For all ocean and perennial inland water bodies	No
1	Evergreen needleleaf	1	Evergreen needle-leaf trees (mostly conifers)	Yes
2	Evergreen broadleaf	2	Evergreen broadleaf trees	Yes
3	Deciduous needleleaf	3	Deciduous needle-leaf trees	Yes
4	Deciduous broadleaf	4	Deciduous broad-leaf trees	Yes
5	Shrub	5	Shrub (woody perennial)	Yes
6	Grass	6	Grasses (native Graminoids)	Yes
7	Cereal crop	7	Cereal cropland (domesticated agricultural crops such as wheat, oats, barley, rye )	Yes
8	Broadleaf crop	8	Broadleaf crop (domesticated agricultural)	Yes
N/A	Urban and built-up	9	Urban and built-up areas (cities, towns, highways, etc)	No
N/A	Snow and ice	10	Snow and ice (may or may not be perennial)	No
N/A	Barren (rock) or sparsely vegetated	11	Barren, rock, or very sparsely vegetated land	No
N/A	Unclassified	254	Areas otherwise not classified as per above	No
Within each 9 km grid cell the number of 1 km grid cells belonging to each vegetated land class is provided in QA/qa_count_pft{1..8}. Non-vegetated grid cells (i.e. cells where the algorithm is not applied), are determined by combining specified vegetation PFT classes (see Table 1) and long-term MODIS fPAR (MOD15A2) where available. Vegetated PFT grid cells lacking sufficient fPAR retrievals to produce the fPAR climatology and non-vegetated PFT grid cells with otherwise valid fPAR climatology are excluded from SPL4CMDL simulations and QA counts. QA counts are time-static and are therefore identical across files because the PFT classification does not change over the course of data generation within each SPL4CMDL version. Users may use the QA count information to compute total non-vegetated 1 km grid cell coverage, compute percent coverage for each PFT, and account for non-vegetated regions when computing areal averages from SPL4CMDL state variables. For example, when computing the total GPP within a 9 km grid cell, a user would multiply the mean GPP (i.e. /GPP/gpp_mean in g C m-2 d-1) by the vegetated PFT total QA count (i.e. /QA/qa_count).

Refer to the Data Fields document for details on all parameters.

File Information

Format

Data are in HDF5 format. For software and more information, including an HDF5 tutorial, visit the HDF Group's HDF5 website.

File Contents

As shown in Figure 1, each HDF5 file is organized into the following main groups, which contain additional groups and/or data sets:

Data Fields

Each file contains the main data groups summarized in this section. For a complete list and description of all data fields within these groups, refer to the Data Fields document.

All global data arrays have dimensions of 1624 rows and 3856 columns (6,262,144 pixels per layer). Note: The EASE-Grid 2.0 global 1 km reference grid is defined as 14616 lines by 34704 samples (507,233,664 pixels per layer).

EC

Environmental Constraints Data

GEO

Geolocation data, including latitude/longitude coordinate variables in decimal degree units that enable convenient geo-referenced viewing and analysis.

GPP

Gross Primary Production Data

NEE

Net Ecosystem CO₂ Exchange Data

QA

QA includes quality control flags, quality assessment, and valid grid cell counts.

RH

Heterotrophic Respiration Data

SOC

Soil Organic Carbon Data

Metadata Fields

Includes all metadata that describe the full content of each file. For a description of all metadata fields for this product, refer to the Metadata Fields document.

File Naming Convention

Files are named according to the following convention, which is described in Table 2:

SMAP_L4_C_MDL_yyyymmddThhmmss_VLMmmm_NNN.[ext]

For example:

SMAP_L4_C_mdl_20151007T000000_Vv3040_001.h5

Where:

Table 2. File Naming Conventions

Variable	Description
SMAP	Indicates SMAP mission data
L4_C_MDL	Indicates specific product (L4: Level-4; C: Carbon; MDL: Model)
yyyymmddThhmmss	Date/time in Universal Coordinated Time (UTC) of the first data element that appears in the product, where: yyyymmdd 4-digit year, 2-digit month, 2-digit day T Time (delineates the date from the time, i.e. yyyymmddThhmmss) hhmmss 2-digit hour, 2-digit minute, 2-digit second
VLMmmm	Science Version ID, where: Variable Description V Version L Launch Indicator (V: Validated Data) M 1-Digit Major Version Number mmm 3-Digit Minor Version Number Example: Vv3040 indicates a Validated product with a version of 3.040. Refer to the SMAP Data Versions page for version information. Note: The data product Science Version ID (example: Vv3040) consists of the first six characters of the data product Composite Release ID (CRID). The full CRID includes four additional digits that are to be found in individual granule metadata within the DataIdentifcation/DatasetIdentification/CompositeReleaseID field. These additional digits denote minor processing changes, such as runtime configuration and other minor changes that do not impact the science of the data product.
NNN	Number of times the file was generated under the same version for a particular date/time interval (002: 2nd time)
.[ext]	File extensions include: .h5 HDF5 data file .xml XML Metadata file

File Size

Each file is approximately 133 MB.

File Volume

The daily data volume is approximately 133 MB.

Spatial Information

Coverage

Coverage spans from 180°W to 180°E, and from approximately 85.044°N and 85.044°S.

Resolution

Level-4 carbon model inputs include the following spatial resolutions:

• 500 m resolution MODIS-based global PFT classification (from MCD12Q1 Type 5)
• 500 m Fraction of Photosynthetically Active Radiation (fPAR) data (from MOD15A2)
• 9 km resolution SMAP Level-4 soil moisture data (SPL4SMGP)
• ¼ degree pre-processed global, daily averaged meteorology data from the GEOS-5 Forward Processing (FP) system

Level-4 carbon model processing is conducted at 1 km EASE-Grid 2.0 resolution using spatially aggregated MODIS PFT and fPAR inputs. Level-4 carbon model daily global outputs are gridded using a 9 km EASE-Grid 2.0 projection consistent with the SMAP L4 soil moisture data used as input.

Note that while this product has a 9 km spatial resolution, it also retains sub-grid scale heterogeneity information as determined from the 1 km resolution processing using MODIS PFT and fPAR inputs.

For more details regarding inputs used in the carbon model, refer to the Data Sources section.

EASE-Grid 2.0

These data are provided on the global cylindrical EASE-Grid 2.0 (Brodzik et al. 2012). Each grid cell has a nominal area of approximately 9 x 9 km² regardless of longitude and latitude.

EASE-Grid 2.0 has a flexible formulation. By adjusting a single scaling parameter, a family of multi-resolution grids that nest within one another can be generated. The nesting can be adjusted so that smaller grid cells can be tessellated to form larger grid cells. Figure 2 shows a schematic of the nesting to a resolution of 3 km (4872 rows x 11568 columns on global coverage), 9 km (1624 rows x 3856 columns on global coverage) and 36 km (406 rows x 964 columns on global coverage). Note that the grid used for this product has been adjusted using a scaling parameter in order to accommodate a resolution of 1 km.

This feature of perfect nesting provides SMAP data products with a convenient common projection for both high-resolution radar observations and low-resolution radiometer observations, as well as for their derived geophysical products.

Temporal Information

Coverage

Coverage is continuous and spans from 31 March 2015 to present.

SMAP Satellite and Processing Events

Due to instrument maneuvers, data downlink anomalies, data quality screening, and other factors, small gaps in the SMAP time series will occur. Details of these events are maintained on two master lists:

SMAP On-Orbit Events List for Instrument Data Users
Master List of Bad and Missing Data

However, gaps in the SMAP time series do not affect this product. For the analytical variables, the ancillary MODIS fPAR 8-day climatology provides a fallback input source to help ensure there are no spatio-temporal gaps in the modeled data record. 

Latencies

FAQ: What are the latencies for SMAP radiometer data sets?

Resolution

Each Level-4 file is a daily composite. Calculations for this product are conducted at a daily time step in order to provide the necessary precision for resolving dynamic boreal vegetation phenology and carbon cycles (Kimball et al. 2009, Kim et al. 2012).

Data Acquisition and Processing

This section has been adapted from Kimball et al. (2014), the ATBD for this product.

Background

Current capabilities for regional assessment and monitoring of NEE are limited by mismatches between bottom-up and top-down information sources. Atmospheric transport model inversions of CO₂ concentrations from sparsely distributed measurement stations provide information on seasonal patterns and trends in atmospheric CO₂ but little information on underlying processes; these methods are also too coarse to resolve carbon source-sink activity at scales finer than broad latitudinal and continental domains (Piao et al. 2007, Dargaville et al. 2002). Tower CO₂ flux measurement networks provide detailed information on stand-level NEE and associated biophysical processes, but little information regarding spatial variability in these processes over heterogeneous landscapes (Running et al. 1999). Estimates of NEE and component carbon fluxes from satellite remote sensing provide a means for scaling between relatively intensive stand-level measurement and modeling approaches, and top-down assessments from atmospheric model inversions.

To address these limitations, the primary objectives of the SPL4CMDL product are to:

• Determine NEE regional patterns and temporal behavior (daily, seasonal, and annual) to within the accuracy range of in situ tower measurement estimates of these processes;
• Link NEE estimates with component carbon fluxes (GPP and Rtot) and the primary environmental constraints to ecosystem productivity and respiration.

The SPL4CMDL algorithm supports carbon cycle science objectives by enabling detailed mapping and monitoring of spatial patterns and temporal dynamics of land-atmosphere CO₂ exchange, and the underlying carbon fluxes and environmental drivers of these processes. The SPL4CMDL product also links SMAP land parameter measurements to global terrestrial  CO₂ exchange, including boreal ecosystems, reducing uncertainties about the 'missing sink' on land for atmospheric CO₂ .

Atmospheric transport model inversions of CO₂ concentrations indicate that the Northern Hemisphere terrestrial biosphere is responsible for much of the recent terrestrial sink strength for atmospheric carbon (Dargaville et al. 2002). Variability in land-atmosphere CO₂ exchange is strongly controlled by climatic fluctuations and disturbance, while uncertainty regarding the magnitude and stability of the sink are constrained by a lack of detailed knowledge on the response of underlying processes at regional scales (Denman et al. 2007, Houghton 2003).

The SPL4CMDL product enables quantification and mechanistic understanding of spatial and temporal variations in NEE over a global domain. NEE represents the primary measure of carbon (CO₂) exchange between the land and atmosphere, and the SPL4CMDL product is directly relevant to a range of applications including regional mapping and monitoring of terrestrial carbon stocks and fluxes, climate and drought related impacts on vegetation productivity, and atmospheric transport model inversions of terrestrial source-sink activity for atmospheric CO₂.

For more background information, refer to Section 2.3: Historical Perspective in the ATBD for this product (Kimball et al. 2014).

Acquisition

The following data sources are used as input to calculating this Level-4 carbon product:

• SMAP L4 9 km EASE-Grid Surface and Root Zone Soil Moisture Geophysical Data, Version 4 (SPL4SMGP)
• GMAO GEOS-5 Forward Processing (FP) Model Data: Daily surface meteorology from observation-corrected global atmospheric model analysis
• NASA EOS Terra MODIS fPAR 8-day Data, Version 6 (MOD15A2): Canopy fPAR and land cover classification; if MOD15A2 data are unavailable, the following back-up sources are used to calculate fPAR:
  • SMAP L4 Carbon Model ancillary MODIS fPAR 8-day Climatology: Primary back-up source
  • NASA EOS Aqua MODIS fPAR 8-day Data, Version 6 (MYD15A2): Canopy fPAR and land cover classification: Secondary back-up source
  • VIIRs NDVI Data (VVI3P): Secondary back-up source

In addition, ancillary data sources used as input to calculating this Level-4 carbon product are listed in Table 3.

Table 3. Primary Ancillary Inputs to the SPL4CMDL Algorithm

Parameter	Units	Type	Spatial Resolution	Source
fPAR	%	Dynamic (8-day)	1 km IV	MODIS (MOD15A2 I)
Rsw	MJ m - 2 d - 1	Dynamic (daily)	9 km II	GEOS-5
Tmn	°C	Dynamic (daily)	9 km II	GEOS-5
VPD	Pa	Dynamic (daily)	9 km II	GEOS-5
SM	% Sat.	Dynamic (daily)	9 km	SPL4SMGP
SMrz	% Sat.	Dynamic (daily)	9 km	SPL4SMGP
Ts	°C	Dynamic (daily)	9 km	SPL4SMGP
F/T	Discrete class	Dynamic (daily)	9 km II	GEOS-5 III
Land Cover Class	Discrete class	Static	1 km IV	MODIS (MOD12Q1)
fPAR Climatology	%	Static (8-day)	1 km IV	MODIS (MOD15A2)
Additional Inputs for Algorithm Options
VI (NDVI)	Dimensionless	Dynamic (8-day)	1 km IV	MODIS (MOD13A2, MYD13A2), VIIRS (VVI3P)
Recovery Status	Years	Static	1 km IV	MODIS (MOD13A2, MYD13A2)
I MOD indicates data acquired by the MODIS instrument on the Terra satellite; MYD indicates data acquired by the MODIS instrument on the Aqua satellite. II The native resolution of GEOS-5 FP fields is ¼ degree (latitude) by 3/8 degree (longitude); SPL4CMDL processing internally resamples these fields to 9 km. III Due to the loss of the SMAP radar instrument and operational freeze/thaw (F/T) classification product, SPL4CMDL uses the GMAO GEOS-5-modeled Tsurf parameter to define F/T conditions in the carbon model. IV Derived from finer scale (500 m resolution) MODIS data records and spatially aggregated to 1 km resolution for carbon model processing.

Baseline Alogrithm

The baseline SPL4CMDL algorithm uses daily inputs from the SMAP Level-4 soil moisture stream to define soil moisture and frozen temperature constraints to vegetation productivity, ecosystem respiration, and NEE. The algorithm provides estimates of NEE (g C m^-2 day^-1) and component carbon fluxes for global vegetated land areas at mean daily intervals; the product defines sub-grid scale mean and variability in carbon fluxes for dominant and sub-dominant vegetation classes within each grid cell as determined from finer scale ancillary land cover classification and fPAR inputs. The targeted accuracy for the SPL4CMDL product is a mean annual unbiased RMSE (ubRMSE) accuracy for NEE within 30 g C m^-2 yr^-1 or 1.6 g C m^-2 day^-1, commensurate with the estimated accuracy of in situ tower measurements (Baldocchi et al. 2008, Richardson 2005, Richardson 2008). The baseline 1 km SPL4CMDL spatial resolution is similar to the sampling footprint of CO₂ flux measurements from the global tower network (Running et al. 1999, Baldocchi et al. 2008). Secondary products of scientific value produced during SPL4CMDL processing include surface (<10 cm depth) Soil Organic Carbon (SOC) stocks (g C m^-2), vegetation Gross Primary Production (GPP), heterotrophic soil and litter respiration (Rh), dimensionless (0-100 percent) environmental constraint indices for GPP and Rh, and detailed data Quality Assessment (QA) metrics for NEE.

The SPL4CMDL algorithm consists of Light Use Efficiency (LUE) and terrestrial carbon flux model components used to estimate GPP, respiration, residual NEE carbon fluxes, and underlying SOC pools on a daily basis. The baseline SPL4CMDL algorithm is summarized in Figures 3a and 3b for respective LUE and carbon flux model components. The approach has structural elements similar to the Century (Parton et al. 1987, Ise and Moorcroft 2006) and CASA (Potter et al. 1993) soil decomposition models and the operational MOD17 GPP algorithm (Zhao et al. 2005, Zhao 2008, and Running 2010), but is adapted for use with daily biophysical inputs derived from both global satellite and model analysis data (Kimball et al. 2009, Yi et al. 2013). 

Dynamic daily inputs to the SPL4CMDL algorithms include satellite optical infrared (IR) remote sensing MODIS-based fPAR, GEOS-5 surface meteorology (Rsw, Tmn, VPD) and associated SPL4SMGP soil moisture (SMrz) which provide primary inputs to a LUE algorithm to determine GPP, where Rsw is incoming shortwave solar radiation (MJ m^-2 d^-1); Tmn is minimum daily 2 m air temperature (°C), VPD is atmosphere vapor pressure deficit (Pa), and SMrz is the integrated surface to root zone (0-1 m depth) soil moisture (% Sat.). The SPL4CMDL dynamic inputs also include GEOS-5 surface temperature (Ts, °C), defined frozen temperature (F/T), constraints to GPP, and autotrophic respiration calculations. SMAP Level-4 surface soil moisture (≤ 5 cm depth) and soil temperature are used as primary drivers of the soil decomposition and Rh calculations. Static inputs to the SPL4CMDL algorithms include a global land cover classification, which is used to define the major plant functional types and associated biome-specific Biome Properties Look-Up Table (BPLUT) response characteristics for each vegetated grid cell within the product domain. The BPLUT parameters are defined for up to eight global vegetation (PFT) classes; the model parameters for each global PFT class were calibrated by optimizing carbon model NEE calculations against tower eddy covariance measurement-based daily NEE observations from global FLUXNET sites* representing the major PFT classes (Baldocchi 2008). The land cover classification used for SPL4CMDL processing is consistent with the one used in the production of the fPAR inputs. All model inputs are available as satellite remote sensing derived products or from model (GEOS-5) analysis.

The resulting SPL4CMDL parameters enable characterization of spatial patterns and daily temporal fidelity in NEE, underlying carbon fluxes and SOC pools, and their primary environmental drivers. The resulting fine scale (1 km resolution) SPL4CMDL outputs are spatially aggregated to the coarser 9 km resolution final product grid by weighted linear averaging of outputs according to the fractional cover of individual PFT classes represented within each 9 km grid cell and defined by the underlying 1 km resolution MODIS PFT map. The sub-grid scale means from individual PFT classes are preserved for each 9 km grid cell, while proportional vegetation cover information is included in the qa_counts_pft data field, allowing the coarse resolution data to be decomposed into the relative contributions from individual PFT classes within each cell. These outputs are designed to facilitate improved algorithm and product accuracy over heterogeneous land cover areas, and product outputs that are more consistent with the mean sampling footprint of most tower CO₂ flux measurement sites (Baldocchi 2008, Chen et al. 2012). 

* Note that an augmented FLUXNET global tower site record which includes more calibration sites (335 sites compared to 228 sites for Version 3 SPL4CMDL), expanded tower data records extending to at least 2015, and the addition of new tower sites representing more land cover types has been implemented for this version. 

Algorithm Options

The SPL4CMDL baseline product contains various processing options that are implemented in the algorithm preprocessing stage for handling of the daily model inputs. These processing options are distinct from other options that are more internal to the model algorithms (Kimball et al. 2014). Two major preprocessing options are used in the SPL4CMDL product, including use of estimated clear-sky fPAR inputs for missing or lower quality MODIS fPAR inputs, and use of GEOS-5 surface temperature fields to estimate frozen temperature constraints to the GPP calculations instead of SMAP radar F/T-defined constraints. The use of these preprocessing options are noted in the SPL4CMDL product bit flags as defined in Table 7 of this document and on the Data Fields page.

For more information regarding algorithm options, refer to the ATBD, for this product (Kimball et al. 2014).

Ancillary Data

Ancillary data required as input for the algorithms are summarized in Table 3. For in-depth information on ancillary data, refer to the ATBD, Section 3.2: Ancillary Data Requirements (Kimball et al. 2014).

For more information regarding the algorithm, refer to the ATBD for this product (Kimball et al. 2014).

Processing

Written by the University of Montana's Numerical Terradynamic Simulation Group (NTSG), the SPL4CMDL science code was transferred from NTSG to the NASA Global Modelling and Assimilation Office (GMAO) for translation and implementation as operational code in conjunction with SMAP Level-4 soil moisture production within the GMAO Level-4 SMAP Science Data Processing System (SDS).

To generate the SPL4CMDL product, the processing software: 

1. Ingests SPL4SMGP daily files, MODIS-derived 8-day fPAR files, and GEOS-5 daily surface meteorology data.
2. The ingested data are then inspected for retrievability criteria according to input data quality, ancillary data availability, and land cover conditions.
3. Two pre-processor codes, one for fPAR data and one for global meteorology data, are then executed each day to temporally aggregate and resample these respective inputs for use by the baseline algorithm software. When retrievability criteria are met, the production software invokes the baseline retrieval algorithm to generate the daily carbon model outputs.

SPL4CMDL calculations are conducted at 1 km resolution, benefiting from finer scale (500 m) MODIS fPAR and land cover inputs. The simulations have also been conducted in a consistent global EASE-Grid 2.0 projection format. Model simulations for each 1 km grid cell are conducted using the corresponding (nearest-neighbor) 9 km resolution SMAP Level-4 Soil Moisture and GEOS-5 inputs. The MODIS (MOD/MYD15) fPAR product is produced at 500 m resolution and 8-day temporal fidelity from both NASA EOS Terra and Aqua sensor records.

MODIS fPAR operational products are obtained in a tile-based sinusoidal projection. Preprocessing of these data prior to the SPL4CMDL ingestion involves reprojecting from sinusoidal to 1 km resolution global cylindrical EASE-Grid projection formats, followed by trailing nearest-neighbor temporal interpolation of MOD15A2 good Quality Control (QC; relatively cloud-free with favorable surface conditions) 8-day fPAR series to each daily time step. Missing or low QC 8-day fPAR data are gap filled on a grid cell-wise basis using an ancillary fPAR mean 8-day climatology constructed from the long-term (10+ year) MODIS record. The resulting fPAR data are combined with daily biophysical inputs from GEOS-5 and SPL4SMGP data to estimate NEE, component carbon fluxes (GPP and Rh) and surface SOC pools. SPL4CMDL computes daily Environmental Constraint (EC) indices which influence the GPP and NEE flux calculations, including the estimated bulk environmental reduction to PAR conversion efficiency (ε_mult), low soil moisture and temperature constraints (W_mult, T_mult) to soil decomposition and Rh calculations, and freeze/thaw (F/T) status within each 9 km grid cell. These environmental constraint indices are provided in SPL4CMDL files as the EC/emult_mean, EC/wmult_mean, EC/tmult_mean and EC/frozen_area respective data fields.

Quality, Errors, and Limitations

Error Sources

Many sources of error contribute to the uncertainty in the SPL4CMDL product. The key sources of error or uncertainty in the SPL4CMDL algorithm are:

1. Errors in the ancillary 8-day fPAR inputs
2. Errors in the SPL4SM soil moisture and temperature inputs
3. Errors in the GEOS-5 daily surface meteorology inputs
4. Uncertainty in the internal model parameterization, initialization, and calibration parameters

For more information about error sources refer to the ATBD for this product.

Quality Assessment 

For in-depth details regarding the quality of these Version 3 data, refer to the following reports: 
Validated Assessment Report 
Beta Assessment Report

Quality Overview

SMAP products provide multiple means to assess quality. Each product contains bit flags, uncertainty measures, and file-level metadata that provide quality information. For information regarding the specific bit flags, uncertainty measures, and file-level metadata contained in this product, refer to the Data Fields and Metadata Fields documents.

Each HDF5 file contains metadata with Quality Assessment (QA) metadata flags that are set by the GMAO prior to delivery to National Snow and Ice Data Center Distributed Active Archive Center (NSIDC DAAC). A separate metadata file with an .xml file extension is also delivered to NSIDC DAAC with the HDF5 file; it contains the same information as the HDF5 file-level metadata.

Quality Flags

Quality Assessment (QA) fields are also provided with metadata from MODIS fPAR and SPL4SM inputs to the SPL4CMDL algorithms. The QA output incorporates expected model uncertainty propagating from input driver uncertainty including SPL4SMGP, GEOS-5 FP, and MODIS fPAR. QA input error information was assigned by comparing unbiased Root Mean Square Errors (ubRMSE) relative to global historical flux tower benchmark data during SPL4CMDL pre-launch calibration. Input errors are propagated during SPL4CMDL 1 km model calculations using standard error propagation procedures employing the SPL4CMDL model Jacobian and simplifying independence assumptions. Resulting 1 km NEE ubRMSE fields are quadratically averaged to 9 km output fields for each PFT class as defined from 1 km MOD12Q1 land cover and then posted as the NEE QA ubRMSE geophysical variable (g C m^-2 d^-1). The resulting QA information has been evaluated and refined through post-launch SPL4CMDL Cal/Val activities using concurrent eddy covariance CO₂ flux measurements from global tower measurement networks (Baldocchi 2008), comparisons with other similar global carbon products, and algorithm sensitivity studies over the observed range of environmental variability. The above-described QA fields are provided in SPL4CMDL files as the QA/nee_rmse_mean and QA/nee_rmse_pft{1..8}_mean fields. Refer to the SPL4CMDL Product Specification Document Version 2.0 (on the Technical References Tab) for additional details.

Quality control bit flags are provided in SPL4CMDL files to identify retrieval conditions including use of alternative ancillary data sets and exceedance of expected output field value ranges. Alternative ancillary conditions indicated in the QC bit flags include the use of alternative fPAR sources in place of baseline MODIS (MOD15) fPAR inputs, potential gaps in the GEOS-5 input stream, and instances where the ancillary fPAR 8-day climatology is used in place of the dynamic best QC MODIS fPAR input stream to estimate GPP. Expected PFT class-specific range thresholds for each state variable (NEE, GPP, Rh, and SOC) have been established from dynamic algorithm simulations using long-term (10+ year) daily data input records from pre-launch data sources similar to those used for post-launch SPL4CMDL production, including MODIS (MOD15) fPAR, freeze-thaw status (Kim et al. 2012), and MERRA surface meteorology (Yi et al. 2011). These post-launch diagnostics are provided in SPL4CMDL files in the QA/carbon_model_bitflag data field for additional user evaluation. Table 4 indicates the bit-field positions for the above-described flags. A copy of Table 4 is also provided within each file as metadata for quick reference; refer to the QA/carbon_model_bitflag data field.

Table 4. QC Bit Flag Fields, Names, Positions, and Description Metadata

Bit Flag Name	Bit Positions {Start, End}	Number of Bits	Value Range	Description
NEE bit	00 – 00	1	{0|1}	0 = NEE within valid range; 1 = out of valid range
GPP bit	01 – 01	1	{0|1}	0 = GPP within valid range; 1 = out of valid range
Rh bit	02 – 02	1	{0|1}	0 = Rh within valid range; 1 = out of valid range
SOC bit	03 – 03	1	{0|1}	0 = SOC within valid range; 1 = out of valid range
PFT dominant	04 – 07	4	{1..8}	Most frequently occurring (dominant) vegetated PFT class as defined from qa_count
QA score	08 – 11	4	{0,1,2,3}	Relative nee_mean error as ranked by nee_rmse_mean: 0 = (RMSE<1 g C m-2 d-1); 1 = (1<=RMSE<2 g C m-2 d-1); 2= (2=<RMSE<3 g C m-2 d-1); 3 = (RMSE> = 3 g C m-2 d-1)
GPP method	12 – 12	1	{0|1}	0 = derived GPP using 8-day fPAR or NDVI input, 1 = derived GPP via fPAR or NDVI climatology
NDVI method	13 – 13	1	{0|1}	0 = derived GPP using fPAR; 1 = derived GPP using NDVI
F/T method	14 – 14	1	{0|1}	0 = used SPL3SMA F/T; 1 = used GEOS-5 surface temperature
IsFill*	15 – 15	1	{0|1}	0 = is NOT fill value (simulation performed for one or more 1 km grid cells within 9 km grid cell), 1 = is fill value (no 1 km simulation performed within 9 km grid cell). Fill values occur for non-land, non-vegetated, and/or grid cells otherwise lacking valid fPAR data record.
* When IsFill = 1, then all other bit fields will have value 1 and the entire Uint16 integer will evaluate to 65534. Users should therefore check the value of IsFill prior to referencing other bit fields.

Note: Although the SPL4CMDL product is global in extent, product accuracy requirements and validation activities were primarily specified for northern (≥45°N) land areas consistent with NRC objectives for better understanding of terrestrial carbon source/sink activity in boreal regions (NRC 2007, Jackson et al. 2011).

For more information, such as algorithm testing procedures, refer to the ATBD (Kimball et al. 2014). For more information regarding data flags, refer to the Data Fields document.

Instrumentation

Description

For a detailed description of the SMAP instrument, visit the SMAP Instrument page at the JPL SMAP website.

Software and Tools

For tools that work with SMAP data, refer to the Tools web page.

Version History

Document Creation Date

October 2015

Document Revision Date

June 2018

Related Data Sets

SMAP Data at NSIDC | Overview

SMAP Radar Data at the ASF DAAC

Related Websites

SMAP at NASA JPL

Contacts and Acknowledgments

Investigators

John Kimball, Lucas Jones, Tobias Kundig
Numerical Terradynamic Simulation Group (NTSG)
College of Forestry & Conservation
The University of Montana
Missoula, MT 59812-1049 USA

Rolf Reichle
NASA Goddard Space Flight Center
Global Modeling and Assimilation Office
Mail Code 610.1
8800 Greenbelt Rd
Greenbelt, MD 20771 USA

References

Baldocchi, D. 2008. Breathing of the Terrestrial Biosphere: Lessons Learned from a Global Network of Carbon Dioxide Flux Measurement Systems. Austr. J. Bot. 56: 1-26.

Brodzik, M. J., B. Billingsley, T. Haran, B. Raup, and M. H. Savoie. 2012. EASE-Grid 2.0: Incremental but Significant Improvements for Earth-Gridded Data Sets. ISPRS Int. J. Geo-Inf. 1(1):32-45. https://doi.org/10.3390/ijgi1010032.

Brodzik, M. J., B. Billingsley, T. Haran, B. Raup, and M. H. Savoie. 2014. Correction: Brodzik, M. J. et al. EASE-Grid 2.0: Incremental but Significant Improvements for Earth-Gridded Data Sets. ISPRS Int. J. Geo-Inf 2012. 1(1):32-45 ISPRS Int. J. Geo-Inf. 3(3):1154-1156. https://doi.org/10.3390/ijgi3031154.

Chen, B., N. C. Coops, D. Fu et al. 2012. Characterizing Spatial Representativeness of Flux Tower Eddy-Covariance Measurements across the Canadian Carbon Prgram Network Using Remote Sensing and Footprint Analysis. Rem. Sens. Environ. 124:742-755.

Dargaville, R. A. D. McGuire, and P. Rayner. 2002. Estimates of Large-Scale Fluxes in High Latitudes from Terrestrial Biosphere Models and an Inversion of Atmospheric CO₂ Measurements. Climatic Change 55:273–285.

Denman, K. L., G. Brasseur, A. Chidthaisong, P. Ciais, P. M. Cox, R. E. Dickinson, D. Hauglustaine, C. Heinze, E. Holland, D. Jacob, U. Lohmann, S. Ramachandran, P. L. da Silva Dias, S. C. Wofsy, and X. Zhang. 2007. Couplings Between Changes in the Climate System and Biogeochemistry. In: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [Solomon, S., D. Qin, M. Manning, Z. Chen, M. Marquis, K. B. Averyt, M. Tignor, and H. L. Miller (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA.

Entekhabi, D. et al. 2014. SMAP Handbook–Soil Moisture Active Passive: Mapping Soil Moisture and Freeze/Thaw from Space. Pasadena, CA USA: SMAP Project, JPL CL#14-2285, Jet Propulsion Laboratory.

Glassy, J., J. S. Kimball, L. Jones, R. H., Reichle, R. A. Lucchesi, J. V. Ardizzone, G.-K. Kim, and B. H. Weiss. 2016. Soil Moisture Active Passive (SMAP) Mission Level 4 Carbon (L4_C) Product Specification Document. GMAO Office Note No. 11, Version 2.0. NASA Goddard Space Flight Center, Greenbelt, MD, USA. (SMAP_L4C_PSD_v2_0_20160421.pdf, 1.3 MB)

Houghton, R. A. 2003. Why are Estimates of the Terrestrial Carbon Balance so Different? Global Change Biol. 9:500-9.

Ise, T., and P. R. Moorcroft. 2006. The Global-Scale Temperature and Moisture Dependencies of Soil Organic Carbon Decomposition: An Analysis Using a Mechanistic Decomposition Model. Biogeochemistry 80: 217–231.

Jackson, T. J., J. S. Kimball, R. Reichle, W. Crow, A. Colliander, and E. Njoku. 2011. SMAP Science Calibration and Validation Plan. SMAP Science Document, No. 014, Version 1.2 (Preliminary Release), JPL D-52544:1-93.

Jones, L. A. et al. 2017. The SMAP Level 4 Carbon Product for Monitoring Ecosystem Land–Atmosphere CO₂ Exchange. IEEE Transactions on Geoscience and Remote Sensing, 55(11):6517-6532. https://doi.org/10.1109/TGRS.2017.2729343.

Kim, Y., J. Kimball, K. Zhang, K. McDonald. 2012. Satellite Detection of Increasing Northern Hemisphere Non-Frozen Seasons from 1979 to 2008: Implications for Regional Vegetation Growth. Remote Sens. Environ., 121:472-487.

Kimball, J. S., Jones L. A.,  Zhang K., Heinsch F. A., McDonald K. C., and Oechel W. C. 2009. A Satellite Approach to Estimate Land-Atmosphere CO₂ Exchange for Boreal and Arctic Biomes using MODIS and AMSR-E. IEEE Geosci. Remote Sens. 47:569-87.

Kimball, J. S., L. A. Jones, J. P. Glassy, and R. Reichle. 2014. SMAP Algorithm Theoretical Basis Document, Release A: L4 Carbon Product. SMAP Project, JPL D-66484, Jet Propulsion Laboratory, Pasadena CA. 76 pp. (PDF, 2.7 MB; see Technical References)

Kimball, J. S., L. A. Jones, J. Glassy, E. N. Stavros, N. Madani, R. H. Reichle, T. Jackson, and A. Colliander. 2015. Soil Moisture Active Passive (SMAP) Project Calibration and Validation for the L4_C Beta-Release Data Product. NASA/TM–2015-104606, 42:1-37. (PDF, 13.7 MB; see Technical References)

Kimball, J. S., L. A. Jones, J. Glassy, E. N. Stavros, N. Madani, R. H. Reichle, T. Jackson, and A. Colliander. 2016. Soil Moisture Active Passive (SMAP) Mission Assessment Report for the Version 2 Validated Release L4_C Data Product. GMAO Office Note No. 13 (Version 1.0):1-37. NASA Goddard Space Flight Center, Greenbelt, MD, USA. (PDF, 2.25 MB; see Technical References)

National Research Council (NRC). 2007. Earth Science and Applications from Space: National Imperatives for the Next Decade and Beyond (Executive Summary). National Academy of Sciences, National Academies Press, Washington DC. 1-35. http://www.nap.edu/catalog/11820.html.

Parton, W. J., D. S. Schimel, C. V. Cole, and D. S. Ojima. 1987. Analysis of Factors Controlling Soil Organic Matter Levels in Great Plains Grasslands. Soil Sci. Soc. Am. J. 51: 1173–1179.

Piao, S., P. Ciais, P. Friedlingstein, P. Peylin, M. Reichstein, S. Luyssaert, H. Margolis, J. Fang, A. Barr, A. Chen, A. Grelle, D. Y. Hollinger, T. Laurila, A. Lindroth, A. D. Richardson, and T. Vesala, 2007. Net Carbon Dioxide Losses of Northern Ecosystems in Response to Autumn Warming. Nature 451:49-52, https://doi.org/doi:10.1038/nature06444.

Potter, C. S., J. T. Randerson, C. B. Field, P. A. Matson, P. M. Vitousek, H. A. Mooney, and S. A. Klooster. 1993. Terrestrial Ecosystem Production: A Process Model Based on Global satellite and Surface Data. Global Biogeochemical Cycles 7(4):811-841.

Richardson, A. D., D. Y. Hollinger. 2005. Statistical Modeling of Ecosystem Respiration Using Eddy Covariance Data: Maximum Likelihood Parameter Estimation, and Monte Carlo Simulation of Model and Parameter Uncertainty Applied to Three Simple Models. Ag. For. Meteor. 131:191-208.

Richardson, A. D, M. D. Mahecha, E. Falge, et al., 2008. Statistical Properties of Random CO₂ Flux Measurement Uncertainty Inferred from Model Residuals. Ag. For. Meteor. 148:38-50.

Running, S. W., D. D. Baldocchi, D. P. Turner, S. T. Gower, P. S. Bakwin, and K. A. Hibbard. 1999. A Global Terrestrial Monitoring Network Integrating Tower Fluxes, Flask Sampling, Ecosystem Modeling and EOS Satellite Data. Remote Sensing of Environment 70:108-127.

Yi, Y., J. S. Kimball, L. A. Jones, R. H. Reichle, R. Nemani, and H. A. Margolis. 2013. Recent Climate and Fire Disturbance Impacts on Boreal and Arctic Ecosystem Productivity Estimated Using a Ssatellite-Based Terrestrial Carbon Flux Model. J. Geophys. Res. Biogeosci. 118:1-17.

Zhao, M., F. A. Heinsch, R. R. Nemani, and S. W. Running. 2005. Improvements of the MODIS Terrestrial Gross and Net Primary Production Global Data Set. Remote Sensing of Environment 95(2):164-175.

Zhao, M., and S. W. Running. 2010. Drought-Induced Reduction in Global Terrestrial Net Primary Production from 2000 through 2009. Science 329(5994):940-943.

Technical References

For additional references, such as ATBDs, refer to the Technical References tab at the top of this user guide.