As of Thursday, 30 January 2020, NASA has implemented an upgrade to the GEOS Forward Processing (FP) system, which provides ancillary data used in generating SMAP data products. Please see the data news announcement for related impacts.
On Wednesday, February 19th between 9:00 am and 11:00 am MDT, the following data collections may not be available due to planned system maintenance: AMSR-E, Aquarius, ASO, High Mountain Asia, IceBridge, ICESat/GLAS, ICESat-2, MEaSUREs, MODIS, NISE, SMAP, SnowEx, and VIIRS. 
Data Set ID:

SMAP L4 Global Daily 9 km EASE-Grid Carbon Net Ecosystem Exchange, Version 4

The Level-4 (L4) carbon product (SPL4CMDL) provides global gridded daily estimates of net ecosystem carbon (CO2) exchange derived using a satellite data based terrestrial carbon flux model informed by the following: Soil Moisture Active Passive (SMAP) L-band microwave observations, land cover and vegetation inputs from the Moderate Resolution Imaging Spectroradiometer (MODIS), Visible Infrared Imaging Radiometer Suite (VIIRS), and the Goddard Earth Observing System Model, Version 5 (GEOS-5) land model assimilation system. Parameters are computed using an Earth-fixed, global cylindrical 9 km Equal-Area Scalable Earth Grid, Version 2.0 (EASE-Grid 2.0) projection.

This is the most recent version of these data.

Version Summary:

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.

COMPREHENSIVE Level of Service

Data: Data integrity and usability verified; data customization services available for select data

Documentation: Key metadata and comprehensive user guide available

User Support: Assistance with data access and usage; guidance on use of data in tools and data customization services

See All Level of Service Details

Data Format(s):
  • HDF5
Spatial Coverage:
N: 85.044, 
S: -85.044, 
E: 180, 
W: -180
Platform(s):Aqua, SMAP, SUOMI-NPP, Terra
Spatial Resolution:
  • 9 km x 9 km
Temporal Coverage:
  • 31 March 2015
Temporal Resolution1 dayMetadata XML:View Metadata Record
Data Contributor(s):Kimball, J. S., L. A. Jones, T. Kundig, and R. Reichle.

Geographic Coverage

Other Access Options

Other Access Options


As a condition of using these data, you must cite the use of this data set using the following citation. For more information, see our Use and Copyright Web page.

Kimball, J. S., L. A. Jones, T. Kundig, and R. Reichle. 2018. SMAP L4 Global Daily 9 km EASE-Grid Carbon Net Ecosystem Exchange, Version 4. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: [Date Accessed].
7 January 2019
Last modified: 
3 December 2019

Data Description


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 COwith 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 CO2 efflux to the atmosphere. All parameters are expressed in units of g C m-2 day-1. The CO2 flux state variable outputs are provided in SPL4CMDL files as eight vegetated land-cover classes called Plant Function Types (PFTs). For example, the CO2 flux state variable outputs are provided in NEE/nee_pft{1..8}_meanGPP/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 Appendix of this document for details on all parameters.

File Information


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:

File Structure Image SPL4CMDL
Figure 1. Subset of File Contents
For a complete list of file contents for the SMAP Level-4 carbon product, refer to the Appendix of this document. 

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 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).


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 CO2 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.

File Naming Convention

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


For example:



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

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 spans from 180°W to 180°E, and from approximately 85.044°N and 85.044°S.


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

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 km2 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.

Figure 2. Perfect Nesting in EASE-Grid 2.0

Temporal Information


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


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


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.


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 CO2 concentrations from sparsely distributed measurement stations provide information on seasonal patterns and trends in atmospheric CO2 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 CO2 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 CO2 exchange, and the underlying carbon fluxes and environmental drivers of these processes. The SPL4CMDL product also links SMAP land parameter measurements to global terrestrial  CO2 exchange, including boreal ecosystems, reducing uncertainties about the 'missing sink' on land for atmospheric CO2 .

Atmospheric transport model inversions of CO2 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 CO2 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 (CO2) 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 CO2.

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


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 (MOD15A2I)
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 (MOD13A2MYD13A2), VIIRS (VVI3P)
Recovery Status Years Static 1 km IV MODIS (MOD13A2MYD13A2)

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 CO2 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). 

Figure 3a. Baseline Light Use Efficiency (LUE) Carbon Model Structure for Estimating GPP 
(Click image for high-resolution version)
Arrows denote the primary pathways of data flow, while boxes denote the major process calculations. Primary inputs include daily root zone soil moisture (SMrz) and landscape freeze/thaw (FT) status from SMAP Level-4 soil moisture products (in red), and other dynamic ancillary inputs (in green) such as MODIS (MOD15) fPAR and GMAO GEOS-5 daily surface meteorology, as well as vapor pressure deficit (VPD), minimum air temperature (Tmn), and incident solar shortwave radiation (Rsw). Model calculations are performed at 1 km spatial resolution using dominant vegetation class and Biome Properties Look-Up Table (BPLUT) response characteristics for each grid cell defined from a global land cover classification. The resulting GPP calculation is a primary input to the Level-4 carbon terrestrial carbon flux model below (Figure 4b).

Figure 3b. Terrestrial Carbon Flux Model for Estimating NEE
(Click image for high-resolution version)
Primary algorithm inputs (in red) include daily GPP from the LUE model (Figure 4a), and surface soil moisture (SM) and surface temperature (Ts) from the SMAP Level-4 product. NEE is the primary (validated) output, while GPP, respiration (Rh + Ra), and SOC are secondary (research) outputs. 

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

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).


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 (Wmult, Tmult) 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 Appendix of this document or the Product Specification Document.

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 CO2 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}

of Bits

Value Range


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



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

Table 5. 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.

Related Data Sets

SMAP Data at NSIDC | Overview

SMAP Radar Data at the ASF DAAC

Related Websites


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


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.

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.

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 CO2 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. 

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 CO2 Exchange. IEEE Transactions on Geoscience and Remote Sensing, 55(11):6517-6532.

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 CO2 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. 

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.

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,

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 CO2 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. ResBiogeosci. 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:

  • Root Data - time and space coordinate data sets are not stored in a particular group
  • 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.

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 Function Type (PFT) classes (refer to Table A9). Note that the QA variables also includes PFT classes. All output variables are summarized in the following tables.

Root Data 

Table A1. Data Fields for Root Data
Data Field Name Type Unit Valid Min Valid Max Fill/Gap Value
EASE2_global_projection string N/A N/A N/A N/A
x Float64 m -17367531 17367531 0.0
y Float64 m -7342231 7342231 0.0

Environmental Constraints (EC)

Table A2. 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 A3. 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 A4. 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 A5. 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 A6. 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 A7. 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 A8. 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

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

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.

Mean value of GPP for pixels classified as PFT class [1-8]. PFT class designations are listed in Table A9.
Table A9. 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 A9.


Global daily 9km Net Ecosystem Exchange standard deviation.


Carbon model quality bitflags. Refer to Table A7 in the SPL4CMDL user guide for information on bit flags. A copy of Table A7 is also is 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 A9.


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 A9.


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 A9.


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 A10 lists the acronyms and abbreviations used in this document.

Table A10. 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 Function 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

How To

Programmatic access to data with services such as subsetting, reformatting, and reprojection
This article provides a step-by-step getting started guide to utilizing an Application Programming Interface, or API, for programmatic access to data from the NSIDC Distributed Active Archive Center (DAAC) based on spatial and temporal filters, as well as subsetting, reformatting, and... read more
How to import and geolocate SMAP Level-3 and Level-4 data in ENVI
The following are instructions on how to import and geolocate SMAP Level-3 Radiometer Soil Moisture HDF5 data in ENVI. Testing notes Software: ENVI Software version: 5.3 Platform: Windows 7 Data set: SMAP L3 Radiometer Global Daily 36 km EASE-Grid Soil... read more
How do I visualize SMAP data in Worldview?
This video tutorial provides step-by-step instructions on how to visualize SMAP data in Worldview ( NASA Worldview is a map-based application... read more
How to Import SMAP HDF Data Into ArcGIS
Selected SMAP L4, Version 4 HDF data (SPL4SMAU, SPL4SMGP, & SPL4SMLM) can be added to ArcGIS with a simple drag/drop or using the 'Add Data' function. These data can be imported and visualized but not geolocated. In order to import, project, and scale these data and other SMAP L3 and L4 HDF... read more
How do I access data using OPeNDAP?
Data can be programmatically accessed using NSIDC’s OPeNDAP Hyrax server, allowing you to reformat and subset data based on parameter and array index. For more information on OPeNDAP, including supported data sets and known issues, please see our OPeNDAP documentation: ... read more
How to learn more about SMAP ancillary data
SMAP Ancillary data sets are used to produce SMAP Level-1, -2, -3, and -4 standard data products. Several of these ancillary data sets are produced by external organizations, such as NOAA, the NASA Global Modeling and Assimilation... read more
How to extract point and area data samples using AppEEARS
This step-by-step tutorial demonstrates how to access MODIS and SMAP data using the Application for Extracting and Exploring Analysis Ready Samples (AppEEARS). AppEEARS allows users to access, explore, and download point and area data with spatial, temporal, and parameter subsets. Interactive... read more
Visualize NSIDC data as WMS layers with ArcGIS and Google Earth
NASA's Global Imagery Browse Services (GIBS) provides up to date, full resolution imagery for selected NSIDC DAAC data sets. ... read more
Search, order, and customize NSIDC DAAC data with NASA Earthdata Search
NASA Earthdata Search is a map-based interface where a user can search for Earth science data, filter results based on spatial and temporal constraints, and order data with customizations including re-formatting, re-projecting, and spatial and parameter subsetting. Thousands of Earth science data... read more


What are the latencies for SMAP radiometer data sets?
The following table describes both the required and actual latencies for the different SMAP radiometer data sets. Latency is defined as the time (# days, hh:mm:ss) from data acquisition to product generation. Short name Title Latency Required Actual (mean1) SPL1AP SMAP L1A... read more
What data subsetting, reformatting, and reprojection services are available for SMAP data?
The following table describes the data subsetting, reformatting, and reprojection services that are currently available for SMAP data via the NASA Earthdata Search tool and a Data Subscription... read more