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 Functional 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 Functional 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 Appendix of this 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 Appendix of this document.

All global data fields have dimensions of 1624 rows and 3856 columns (6,262,144 pixels per array).

Environmental Constraints Data

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

Gross Primary Production Data

Net Ecosystem CO₂ Exchange Data

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

Heterotrophic Respiration Data

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 Product Specification Document (Glassy et al., 2016).

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

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 Acquisition and Processing section of this document.

Geolocation

These data are provided on the 9-km global cylindrical EASE-Grid 2.0 equal-area projection. The following tables provide information for geolocating this data set. For more on EASE-Grid 2.0, refer to the EASE Grids website.

Table 3. Geolocation Details for the Global EASE-Grid

Geographic coordinate system	WGS 84
Projected coordinate system	EASE-Grid 2.0 Global
Longitude of true origin	0
Standard Parallel	30
Scale factor at longitude of true origin	N/A
Datum	WGS 84
Ellipsoid/spheroid	WGS 84
Units	meter
False easting	0
False northing	0
EPSG code	6933
PROJ4 string	+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Reference	http://epsg.io/6933

Table 4. Grid Details for the Global EASE-Grid

Grid cell size (x, y pixel dimensions)	9,024.31 m (x) 9,024.31 m (y)
Number of columns	3856
Number of rows	1624
Geolocated lower left point in grid	85.044° S, 180.000 ° W
Nominal gridded resolution	36 km by 36 km
Grid rotation	N/A
ulxmap – x-axis map coordinate of the outer edge of the upper-left pixel	-17367530.45 m
ulymap – y-axis map coordinate of the outer edge of the upper-left pixel	7314540.83 m

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. 

For the period between 19 June and 23 July 2019, an extended gap occurred in the L1 - L3 SMAP products. During this period, the L4 Carbon data set was not informed by SMAP data. For more information on this SMAP outage, users should refer to the SMAP Post-Recovery Notice. 

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 of the ATBD for this product (Kimball et al., 2014).

Instrumentation

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

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 5 (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

Ancillary data sources used as input to calculating this Level-4 carbon product are listed in Table 5.

Table 5. Primary Ancillary Inputs to the SPL4CMDL Algorithm

Parameter	Description	Units	Type	Spatial Resolution	Source
fPAR	fraction of photosynthetically active radiation	%	Dynamic (8-day)	1 km IV	MODIS (MOD15A2I)
Rsw	incident solar shortwave radiation	MJm-2d-1	Dynamic (daily)	9 km II	GEOS-5
Tmn	minimum air temperature	°C	Dynamic (daily)	9 km II	GEOS-5
VPD	vapor pressure deficit	Pa	Dynamic (daily)	9 km II	GEOS-5
SM	surface soil moisture	% Sat.	Dynamic (daily)	9 km	SPL4SMGP
SMrz	root zone soil moisture	% Sat.	Dynamic (daily)	9 km	SPL4SMGP
Ts	surface temperature	°C	Dynamic (daily)	9 km	SPL4SMGP
F/T	freeze/thaw status	Discrete class	Dynamic (daily)	9 km II	GEOS-5 III
Land Cover Class	land cover classification	Discrete class	Static	1 km IV	MODIS (MOD12Q1)
fPAR Climatology	long-term (10+ yrs) mean fPAR	%	Static (8-day)	1 km IV	MODIS (MOD15A2)
Additional Inputs for Algorithm Options
VI (NDVI)	normalized difference vegetation index	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 (~0-5 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). 

* The more recent SPL4CMDL product versions (V4 and later) use an augmented FLUXNET global tower site record which includes more calibration sites (329 sites compared to 228 sites for V3 and earlier), expanded tower data records extending to at least 2015, and additional tower sites representative of more land cover types. 

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 Appendix of this document.

For more information regarding the baseline algorithm and its 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 Section 3.2 of the ATBD.

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 are also conducted in a consistent global EASE-Grid 2.0 projection format. Model simulations for each 1 km grid cell are run 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 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 the SPL4CMDL product, including the latest product release and earlier versions, refer to Jones et al. (2017) and the forthcoming 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 Quality Flags section below, the Appendix of this document, or the Product Specification Document (Glassy et al., 2016).

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. 2017), 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 6 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 6. 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 [SPL3SMA is no longer operational, so bitflag = 1]
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 Appendix of this document.

Software and Tools

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

Version History

Table 7. Summary of Version Changes

Version	Date	Version Changes
Version 1	October 2015	First public data release
Version 2	April 2016	Changes to this version include: Transitioned to Validated-Stage 2 Using SPL4SMAU V2 Validated and SPL4SMGP V2 Validated data as input Update to process radiometer data from 2015-03-31 to 2015-04-12 Some data fields renamed from *_av to *_mean Updated to have have continuous RMSE-based "quality" fields instead of the categorical quality flag in V1
Version 3	July 2017	Changes to this version include: Uses dynamic 8-day fPAR inputs obtained from the latest (Collection 6) MODIS fPAR record at 500 m resolution. The preprocessor was updated to handle the finer resolution MODIS Collection 6 inputs, which are interpolated to 1 km resolution prior to model processing. The prior (Version 2) processor used MODIS Collection 5 fPAR inputs, which were derived at 1 km resolution. Updated the ancillary MODIS fPAR 8-day climatology used for fPAR gap-filling as an L4C model preprocessing step to reflect new MODIS Collection 6 fPAR inputs. The fPAR climatology is derived from a longer 14-year (2000-2014) MODIS record relative to the original 12-year (2000-2012) Collection 5 fPAR record used in Version 2 processing. For each grid cell, a sine-curve-based seasonal fPAR climatology curve is now used to identify and screen anomalous 8-day fPAR variations in the preprocessor. This change reduces impacts of anomalous fPAR temporal variations that may not be captured by the MODIS fPAR product quality control (QC) flags, particularly during seasonal transitions at northern latitudes. Updated and recalibrated the ancillary Biome Properties Look-Up Table (BPLUT) and re-initialized the model initial global soil organic carbon (SOC) pools to reflect new MODIS Collection 6 fPAR inputs. The BPLUT calibration was conducted using global historical FLUXNET in situ tower eddy covariance CO2 flux measurement records for representative global land cover types using a similar step-wise calibration procedure employed for the Version 2 product. A minor bug fix to the post-processor was made to ensure that all grid cell no-data fill values are identified with a consistent -9999 notation; the prior Version 2 product erroneously assigned some no-data values as -999900.
V4	June 2018	Changes to this version largely affect model inputs and ancillary files rather than changes to the internal model structure or code. Note that Version 4 is slightly better than Version 3 in RMSE terms, with improvement generally larger for drier sites. Specific changes include: The carbon model biome properties lookup table (BPLUT) has been calibrated using an augmented FLUXNET global tower site record which includes more calibration sites (335 sites compared to 228 sites for V3), expanded tower data records extending to at least 2015, and the addition of new tower sites representing more land cover types. Revised Level-4 carbon global model calibration and SOC initialization using an extended (2000-2017) MODIS fPAR (V006) record and the latest SMAP Nature Run (NRv7.2) climate data records. Implemented minor changes to spatial weighting of calibration tower sites within a model grid cell and reduced the outlier influence on model response curve fitting.
V5	August 2020	Changes to this version primarily represent changes to the model inputs and model Biome Properties Look-up Table (BPLUT) calibration rather than changes to the internal model structure or code. The Version 5 product accuracy and performance is largely consistent with Version 4. Specific changes include: Revised Level-4 carbon global model calibration and soil organic carbon (SOC) initialization using the latest SMAP Nature Run (NRv8.1) soil moisture and soil temperature data, and MERRA-2 surface meteorology. The carbon model BPLUT has been recalibrated using a FLUXNET2015 global tower site record that includes 329 sites representing all major global plant functional type (PFT) classes. Implemented minor changes to spatial weighting of FLUXNET2015 tower sites in the Level-4 carbon model BPLUT calibration. Northern (≥ 45°N) tower sites now have greater impact in defining global PFT response characteristics; whereby, the adjusted weighting is designed to enhance product performance in boreal-arctic biomes. It also incorporated a more detailed set of PFT realistic range constraints for the BPLUT parameters prior to the model calibration and optimization. Enhanced an attribute in the legend_carbon_model_bitflag data field to include descriptions of the bit values.

Related Data Sets

SMAP Data at NSIDC | Overview

SMAP Radar Data at the ASF DAAC

Related Websites

SMAP at NASA JPL

Contacts and Acknowledgments

John Kimball, Lucas Jones, Tobias Kundig
Numerical Terradynamic Simulation Group
The University of Montana
Missoula, MT

Rolf Reichle
NASA Goddard Space Flight Center
Greenbelt, MD

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.

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, J. V. Ardizzone, G.-K. Kim, R. A. Lucchesi, E. B. Smith, and B. H. Weiss. 2016. 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. (see Technical References or PDF).

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.

Kim, Y., J. S. Kimball, J. Glassy, J. Du. 2017. An extended global Earth system data record on daily landscape freeze–thaw status determined from satellite passive microwave remote sensing. Earth System Science Data, 9(1): 133–147. https://doi.org/10.5194/essd-9-133-2017

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: L4 Carbon Product, Release A. SMAP Project, JPL D-66484, Jet Propulsion Laboratory, Pasadena CA. 76 pp. (see Technical References or PDF).

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.

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. 

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.

Appendix - Data Fields

This appendix provides a description of all data fields within the SMAP L4 Global Daily 9 km Carbon Net Ecosystem Exchange product. The data are provided in the following main HDF5 groups:

• EC — Environmental Constraints
• GEO — Geolocation Information
• GPP — Gross Primary Production
• NEE — Net Ecosystem Exchange
• QA — Quality Assurance
• RH — Heterotrophic Respiration
• SOC — Soil Organic Carbon
• Metadata

For a description of metadata fields for this product, refer to the Product Specification Document (Glassy et al., 2016).

SMAP Level-4 carbon files include HDF5 groups that contain six primary variables: EC, GEO, GPP, NEE, RH, and SOC. With the exception of EC and GEO, each group includes data fields for eight Plant Functional Type (PFT) classes (refer to Table A8). Note that the QA variables also includes PFT classes. All output variables are summarized in the following tables.

Environmental Constraints (EC)

Table A1. Data Fields for EC Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
emult_mean	Float32	2	percent	0.0	100.0	-9999.0	1624 x 3856
frozen_area	Float32	2	percent	0.0	100.0	-9999.0	1624 x 3856
tmult_mean	Float32	2	percent	0.0	100.0	-9999.0	1624 x 3856
wmult_mean	Float32	2	percent	0.0	100.0	-9999.0	1624 x 3856

Geolocation Information (GEO)

Table A2. Data Fields for GEO Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
latitude	Float32	2	decimal degree	-89.999	89.999	-9999.0	1624 x 3856
longitude	Float32	2	decimal degree	-179.999	179.999	-9999.0	1624 x 3856

Gross Primary Production (GPP)

Table A3. Data Fields for GPP Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
gpp_mean	Float32	2	gCm-2day-1	0.0	30.0	-9999.0	1624 x 3856
gpp_pft[1-8]_mean	Float32	2	gCm-2day-1	0.0	30.0	-9999.0	1624 x 3856
gpp_std_dev	Float32	2	gCm-2day-1	0.0	30.0	-9999.0	1624 x 3856

Net Ecosystem Exchange (NEE) 

Table A4. Data Fields for NEE Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
nee_mean	Float32	2	gCm-2day-1	-30.0	20.0	-9999.0	1624 x 3856
nee_pft[1-8]_mean	Float32	2	gCm-2day-1	-30.0	20.0	-9999.0	1624 x 3856
nee_std_dev	Float32	2	gCm-2day-1	-30.0	20.0	-9999.0	1624 x 3856

Quality Assurance (QA)

Table A5. Data Fields for QA Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
carbon_model_bitflag	Uint16	2	dimensionless	0.0	65534	65534	1624 x 3856
nee_rmse_mean	Float32	2	gCm-2day-1	0.0	20.0	-9999.0	1624 x 3856
nee_rmse_pft[1-8]_mean	Float32	2	gCm-2day-1	0.0	20.0	-9999.0	1624 x 3856
qa_count	Uint8	2	dimensionless	0	81.0	254	1624 x 3856
qa_count_pft[1-8]	Uint8	2	dimensionless	0	81.0	254	1624 x 3856

Heterotrophic Respiration (RH)

Table A6. Data Fields for RH Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
rh_mean	Float32	2	gCm-2day-1	0.0	20.0	-9999.0	1624 x 3856
rh_pft[1-8]_mean	Float32	2	gCm-2day-1	0.0	20.0	-9999.0	1624 x 3856
rh_std_dev	Float32	2	gCm-2day-1	0.0	20.0	-9999.0	1624 x 3856

Soil Organic Carbon (SOC)

Table A7. Data Fields for SOC Group

Data Field Name	Type	Bit	Unit	Valid Min	Valid Max	Fill/Gap Value	Dimensions
soc_mean	Float32	2	gCm-2	0.0	25000.0	-9999.0	1624 x 3856
soc_pft[1-8]_mean	Float32	2	gCm-2	0.0	25000.0	-9999.0	1624 x 3856
soc_std_dev	Float32	2	gCm-2	0.0	25000.0	-9999.0	1624 x 3856

Data Field Definitions

Defines the parameters of the cylindrical 9 km Earth-fixed EASE-Grid 2.0 projection and the mapping from latitude/longitude to grid-native coordinates

x
The x coordinate values from the cylindrical 9 km Earth-fixed EASE-Grid 2.0 projection

y
The y coordinate values from the cylindrical 9 km Earth-fixed EASE-Grid 2.0 projection

Environmental-constraint multiplier. PFT class designations are listed in Table A8.

Frozen area percentage (environmental-constraint). Mean frozen_area, per 9 x 9 km pixel, not distinguished by PFT (units = percent constraint; min=0 fully constrained; max = 100.0 not constrained).

Environmental-constraint on temperature.

Environmental-constraint on moisture.

Latitude (north..south) coordinates in decimal degrees.

Longitude (west..east) coordinates in decimal degrees.

Global daily 9km Gross Primary Productivity (GPP) mean.

Table A8. PFT Classes

Class	Designation
1	Evergreen Needleleaf
2	Evergreen Broadleaf
3	Deciduous Needleleaf
4	Deciduous Broadleaf
5	Shrub
6	Grass
7	Cereal Crop
8	Broadleaf Crop

Global daily 9km Gross Primary Productivity (GPP) standard deviation.

Global daily 9km Net Ecosystem Exchange (mean).

Mean value of NEE for pixels classified in PFT classes 1 through 8. The PFT class designations are listed in Table A8.

Global daily 9km Net Ecosystem Exchange standard deviation.

Carbon model quality bitflags. Refer to Table 6 in this User Guide for information on bit flags. A copy of Table 6 is also included as an attribute called legend_carbon_model_bitflag data field for quick reference. 

NEE unbiased RMSE estimate (units: gCm^-2day^-1).

NEE unbiased RMSE estimate (units: gCm^-2day^-1) for PFT class [1-8]. PFT class designations are listed in Table A8.

Carbon model QA count 0 to 81 (overall).

QA count for PFT class [1-8]. PFT class designations are listed in Table A8.

Mean heterotrophic respiration (Rh) value, not distinguished by PFT; in gCm^-2day^-1 units, (valid_min=0.0,valid_max=20.0).

Mean value of Rh for pixels classified as PFT class [1-8]. PFT class designations are listed in Table A8

Global daily 9km heterotrophic respiration standard deviation.

Global daily 9km soil organic carbon mean.

Mean value of SOC for pixels classified as PFT class [1-8]. PFT class designations are listed in Table A8.

Global daily 9km soil organic carbon standard deviation.

Fill/Gap Values 

SMAP data products employ fill and gap values to indicate when no valid data appear in a particular data field. Fill values ensure that data fields retain the correct dimension. Gap values locate portions of a data stream that do not appear in the output data file. Fill values appear in the SPL4CMDL data product over ocean and water surfaces. SMAP data products employ a specific set of data values to connote that a field is fill. The selected values that represent fill are dependent on the data type.

No valid value in the SPL4CMDL data product is equal to the values that represent fill. If any exceptions should exist in the future, the SPL4CMDL content will provide a means for users to discern between fields that contain fill and fields that contain genuine data values. Operationally, software should not attempt to ever test two floating point values for equality, but instead test two given values using a relational inequality operator and a tolerance-based difference in absolute values (e.g. if abs(a-b) <= 1E-14) method.

For discrete categorical variables such as the bit-flags (carbon_bitflag and carbon_qual_flag_pft[1-8]), all values are defined, therefore no missing data values are expected with these variables. For the analytical variables, the fPAR climatology provides a fallback input source to help assure there are no spatio-temporal gaps in the modeled data record. Therefore no gaps are expected to occur in the SPL4CMDL output data stream.

Acronyms and Abbreviations

Table A9. Acronyms and Abbreviations

Abbreviation	Definition
Char	8-bit character
Int8	8-bit (1-byte) signed integer
Int16	16-bit (2-byte) signed integer
Int32	32-bit (4-byte) signed integer
EC	Environmental Constraints
Float32	32-bit (4-byte) floating-point integer
Float64	64-bit (8-byte) floating-point integer
fPAR	Fraction of Photosynthetically Active Radiation
GEO	Geolocation Information
GPP	Gross Primary Production
N/A	Not Applicable
NEE	Net Ecosystem Exchange
PFT	Plant Functional Type
QA	Quality Assurance
RH	Heterotrophic Respiration
SOC	Soil Organic Carbon
Uint8	8-bit (1-byte) unsigned integer
Uint16	16-bit (2-byte) unsigned integer
UTC	Universal Coordinated Time