Scientists
Research Grants
Scientific Publications
Informatics Research
NSIDC Scientific Expeditions
Scientists blog from Antarctica and provide a glimpse of what it's like to do research in the field. 
What's hot in the news around climate and sea ice and what are scientists talking about now? 
This Level-2 (L2) soil moisture product provides estimates of global land surface conditions retrieved by the Soil Moisture Active Passive (SMAP) active radar during 6:00 a.m. descending half-orbit passes, as well as ancillary data such as surface temperature and vegetation water content. Input backscatter data used to derive soil moisture are resampled to an Earth-fixed, global, cylindrical 3 km Equal-Area Scalable Earth Grid, Version 2.0 (EASE-Grid 2.0).
Data Set ID:
SPL2SMA
SMAP L2 Radar Half-Orbit 3 km EASE-Grid Soil Moisture, Version 3
Version Summary:
Changes to this version include:
- Transitioned to Validated-Stage 1
- Uses updated SPL1CS0 V3 as input
- Updates to some of the data cubes
- Update of crop database
- Modified iterative bias correction
- Estimation of VWC
Get Data
DOWNLOADING DATA VIA HTTPS
An Earthdata Login account is required to access these data. Please visit the Earthdata Login registration page to register for an account. Once you have logged in, data can be downloaded via a Web browser, command line, or client. For help with downloading data, please see Options Available for Bulk Downloading Data from HTTPS with Earthdata Login.
Earthdata Search: NASA's newest search and order tool for subsetting, reprojecting, and reformatting data.
Geographic Coverage
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
| Parameter(s): |
|
|---|---|
| Spatial Coverage: |
|
| Spatial Resolution: |
|
| Temporal Coverage: |
|
| Temporal Resolution: | 49 minute |
| Data Format(s): |
|
| Platform(s) | SMAP Observatory |
| Sensor(s): | SMAP L-Band Radar |
| Version: | V3 |
| Data Contributor(s): | Kim, S., J. Van Zyl, R. S. Dunbar, E. G. Njoku, J. T. Johnson, M. Moghaddam, and L. Tsang. |
| Metadata XML: | View Metadata Record |
Data Citation
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.
Kim, S., J. Van Zyl, R. S. Dunbar, E. G. Njoku, J. T. Johnson, M. Moghaddam, and L. Tsang. 2016. SMAP L2 Radar Half-Orbit 3 km EASE-Grid Soil Moisture, Version 3. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: https://doi.org/10.5067/J8SGO1E0Y9XZ. [Date Accessed].Collapse All / Open All
Detailed Data Description
Parameter Description
Surface soil moisture (0-5 cm) in cm3/cm3 derived from sigma nought measurements is output on a fixed 3 km EASE-Grid 2.0.
Sigma nought (sigma0), or the backscatter coefficient, is a measure of the strength of radar signals reflected back to the instrument from a target, and is defined as per unit area on the ground. Usually expressed in dB, it is a normalized dimensionless number, comparing the strength observed to that expected from a defined area. The SMAP L-band Radar measures sigma0 using VV, HH, and HV transmit-receive polarizations, and uses separate transmit frequencies for the H (1.26 GHz) and V (1.29 GHz) polarizations. Sigma0 measurements are derived using Synthetic-Aperture Radar (SAR) processing.
Refer to the Data Fields document for details on all parameters.
Background color on
Format
Data are in HDF5 format. For software and more information, including an HDF5 tutorial, visit the HDF Group's HDF5 Web site.
Background color on
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:
For a complete list of file contents for the SMAP Level-2 radar soil moisture product, refer to the Data Fields page.
Data Fields
Each file contains the main data groups summarized in this section. For a complete list and description of all data fields within these groups, refer to the Data Fields document.
All data element arrays are one-dimensional with a size N, where N is the number of valid cells from the radar swath that appear on the grid.
Ancillary Data
Includes all ancillary data, such as surface temperature and vegetation water content.
Radar Data
Includes all radar data, such as cross-polarized sigma0 (sigma nought; also σ0) data.
Soil Moisture Retrieval Data
Includes soil moisture data and quality assessment flags.
Metadata Fields
Includes all metadata that describe the full content of each file. For a description of all metadata fields for this product, refer to the Metadata Fields document.
Background color on
File Naming Convention
Files are named according to the following convention, which is described in Table 1:
SMAP_SPL2SMA_[Orbit#]_D_yyyymmddThhmmss_RLVvvv_NNN.[ext]
For example:
SMAP_SPL2SMA_00934_D_20151225T074951_R13171_002.h5
Where:
| Variable | Description | ||||||||
|---|---|---|---|---|---|---|---|---|---|
SMAP |
Indicates SMAP mission data | ||||||||
SPL2SMA |
Indicates specific product (L2: Level-2; SM: Soil Moisture; A: Active) |
||||||||
[Orbit#] |
5-digit sequential number of the orbit flown by the SMAP spacecraft when data were acquired. Orbit 00000 began at launch. Orbit numbers increment each time the spacecraft flies over the southernmost point in the orbit path. | ||||||||
D |
Descending half-orbit pass of the satellite (where satellite moves from North to South, and 6:00 a.m. is the local solar equator crossing time) | ||||||||
yyyymmddThhmmss |
Date/time in Universal Coordinated Time (UTC) of the first data element that appears in the product, where:
|
||||||||
RLVvvv |
Composite Release ID, where:
Example: R13171 indicates a standard data product with a version of 3.171. Refer to the SMAP Data Versions page for version information. |
||||||||
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:
|
Background color on
File Size
Each half-orbit file is approximately 59 MB.
Background color on
Volume
The daily data volume is approximately 1.5 GB.
Background color on
Spatial Coverage
Coverage spans from 180°W to 180°E, and from approximately 85.044°N and 85.044°S. The gap in coverage at both the North and South Pole, called a pole hole, has a radius of approximately 400 km. The swath width is 1000 km, enabling nearly global coverage every three days.
Spatial Coverage Map
Figure 2 shows the spatial coverage of the SMAP L-Band Radar for one descending half orbit, which comprises one file of this data set.
Background color on
Spatial Resolution
The native spatial resolution of the radar footprint is 1 km. Data are then gridded using the 3 km EASE-Grid 2.0 projection.
Background color on
Projection and Grid Description
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 3 x 3 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 3 shows a schematic of the nesting.
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. For more on EASE-Grid 2.0, refer to the EASE-Grid 2.0 Format Description.
Background color on
Temporal Coverage
Coverage spans from 13 April 2015 to 07 July 2015.
Note: Temporal coverage for this data set is limited due to the premature failure of the SMAP L-Band Radar. On 07 July 2015, the radar stopped transmitting due to an anomaly involving the instrument's high-power amplifier (HPA). For details, refer to the SMAP News Release issued 02 September 2015 by the Jet Propulsion Laboratory (JPL).
Temporal Coverage Gaps
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
Background color on
Temporal Resolution
Each Level-2 half-orbit file spans approximately 49 minutes.
Background color on
Software and Tools
For tools that work with SMAP data, refer to the Tools Web page.
Background color on
Data Acquisition and Processing
Sensor or Instrument Description
For a detailed description of the SMAP instrument, visit the SMAP Instrument page at the JPL SMAP Web site.
Background color on
Data Source
SMAP Level-2 radar soil moisture data (SPL2SMA) are derived from SMAP High-Resolution Radar Sigma Nought, Version 3 (SPL1CS0) data.
Background color on
Theory of Measurements
Retrieval of soil moisture from measured backscatter data typically implies an inversion of the radar forward scattering process. Bare rough surfaces can be characterized in terms of their Root Mean Square (RMS) roughness height, correlation length, and moisture content (a surrogate for dielectric constant). The use of time-series data makes the retrieval a well-constrained estimation problem, under the assumption of a time invariant surface roughness. By taking a co-polarized ratio the soil moisture retrieval becomes insensitive to the correlation length except for very rough surfaces, which enables an accurate retrieval of soil moisture without correlation length information. This approach has been extended to the vegetated surface by introducing a vegetation axis to a lookup table (Kim et al. 2014). A one-axis representation of the vegetation effect is clearly a simplification, considering that different sets of vegetation parameters result in different backscattering coefficients. However, with SMAP’s three measurement channels—HH, VV, and HV—at most three independent parameters can be uniquely estimated, and therefore simplified forward models must be represented in terms of at most three dominant parameters. The simplification results in some errors in soil moisture retrieval, especially in heavily vegetated areas such as forests. Allometric relationships reduce the number of unknowns and may improve the retrievals. The three parameters used to simplify the scattering model are then the dielectric constant of soil, ε, soil surface roughness, s, and VWC.
The SMAP radar HV-channel measurements are reserved for possible use in correcting vegetation effects. The remaining two co-polarized (co-pol) measurements (HH and VV) are not always sufficient to determine s and ε (Kim et al. 2014). One of the main causes is the ambiguity in bare surface scattering: a wet and smooth surface may have the same backscatter as a dry and moderately rough surface. Very often the timescale of the change in s is longer than that of ε (Jackson et al. 1997). Then s may be constrained to be a constant in time, thus resolving the ambiguity (Kim et al. 2014).
The SMAP baseline approach (Kim et al. 2012; Kim et al. 2014) is a multichannel retrieval algorithm that searches for a soil moisture solution such that the difference between modeled and observed backscatter is minimized in the least squares sense. A look-up table representation, or data cube, of a complicated forward model has been demonstrated to be an accurate and fast tool for retrieval (Kim et al. 2012; Kim et al. 2014). The algorithm estimates s first and then retrieves εr using the estimated s. Vegetation effects are included by selecting the sigma0 of the forward model at the VWC level given by an ancillary source. Note that the VWC provided by ancillary information is allowed to vary throughout the time series.
For in-depth information regarding the physics involved in deriving soil moisture from backscatter, refer to the ATBD for this product, Section 2: Physics of the Scattering Problem.
Background color on
Derivation Techniques and Algorithms
SMAP Level-2 active soil moisture data are derived from SMAP High-Resolution Radar Sigma Nought, Version 3 (SPL1CS0) data. Utilizing the soil moisture algorithms discussed below, SMAP sigma0 measurements are converted into an estimate of the surface soil moisture of the 0-5 cm surface soil moisture in units of cm3/cm3 .
The following information has been adapted from Kim et al. 2014.
Soil Moisture Algorithms
The baseline retrieval algorithm inverts a forward scattering model. The equation for forward scattering (Eq. 2-1 in the ATBD) presents four unknowns: RMS height (s), dielectric constant (ε), correlation length (l), and vegetation water content (VWC). It will be shown that mv retrievals are insensitive to errors in knowledge of the correlation length over a wide range of m v, roughness, and correlation lengths (Appendix 5.2) (Kim et al. 2012b). An exponential function is known to describe empirical measurements well (Mattia et al. 1997; Shi et al. 1997), which was used in the forward model and was not considered as an unknown during the retrieval. As an introduction to the discussion of the retrieval algorithm, and to place the discussion of each algorithm component in the proper context, the overall algorithm processing flow is described in the following section. Each portion of the algorithm flow, and the details of the retrieval methods and options, will then be described in detail in the subsequent sections.
Algorithm Flow
The algorithm flow of the SPL2SMA is presented in Figure 11 of the ATBD. The portion of the flowchart through the land surface classification is the initial processing done as a precursor to the actual retrievals. The SPL2SMA processor reads in 1 km resolution sigma0 from the SMAP Level-1C sigma nought product. The 1 km data in natural unit are aggregated onto the 3 km EASE-Grid, during which various quality flags are applied. Three quality flags are derived using the 3 km sigma0 freeze-thaw state (refer to the SPL3FTA ATBD), radar vegetation index (Section 3.4.2 of the SPL3FTA ATBD), and transient water body (Section 3.4.3 of the SPL3FTA ATBD). Static and dynamic ancillary data are sampled for each pixel of the SPL2SMA product. Sigma0 values from the past time are sampled and used by the time-series algorithm. For each 3 km pixel, land cover class information is obtained from the mostly annual ancillary data. The land cover information allows to choose an appropriate data cube for each pixel. The retrieval over the different land cover classes is spatially assembled to create a half-orbit output, followed by the conversion from the dielectric constant to soil moisture.
Baseline Algorithm: Time-series Data-cube Approach
Time-series Retrieval
The SMAP radar provides three independent channels (HH, VV, and HV). HV-channel measurements are reserved for possible use in correcting vegetation effects. The remaining two co-pol measurements (HH and VV) are not always sufficient to determine s and ε (see Section 3.2.4 in the ATBD) (Kim et al. 2012b). One of the main causes is the ambiguity in bare surface scattering: a wet and smooth surface may have the same backscatter as a dry and moderately rough surface has. Very often the time-scale of the change in s is longer than that of ε (Jackson et al. 1997). Then s may be constrained to be a constant in time, thus resolving the ambiguity (Kim et al. 2012b). The concept of ε time-invariant s has also been utilized in other studies (Joseph et al. 2008; Mattia et al. 2009; Verhoest et al. 2007). The SMAP baseline algorithm differs from these studies in that no ancillary or ground measurements or statistical assumptions are required to constrain s.
The SMAP baseline approach (Kim et al. 2012a ; Kim et al. 2012b) is a multichannel retrieval algorithm that searches for a soil moisture solution such that the difference between computed and observed backscatter is minimized in the least squares sense. The algorithm estimates s first and then retrieves εr using the estimated s. Vegetation effects are included by selecting the forward model's sigma0 at the VWC level given by an ancillary source or the SMAP HV measurements. The algorithm retrieves s and real part of dielectric constant (εr) using a time series of N co-pol backscatter measurements: σ0 HH ( t1), σ0 VV (t1), σ0 HH (t2), σ0 VV (t2),..., σ0 HH (t N), and σ0 VV (tN).
There are thus 2 N independent input observations and N+1 unknowns consisting of N εr values and one s value. Note that the VWC provided by ancillary information is allowed to be varying throughout the time series.
Radar backscattering coefficients before conversion to decibels can be modeled as Gaussian random variables (Ulaby et al. 1986) to account for speckle and thermal noise effects. Assuming sufficient integration following power detection, the backscattering coefficient after conversion to decibels can also be modeled as a Gaussian random variable. Because SMAP will observe HH and VV returns at slightly different center frequencies, the effects of speckle and thermal noise on these measurements are statistically independent. Statistical independence of speckle in measurements at differing time steps is also expected. These facts and a maximum likelihood formulation motivate least-squares retrieval approaches based on the average of individual error terms. It is noted that calibration, radio frequency interference, and other error sources may produce correlated error terms. The systematic and correlated components from these sources will be corrected. Any residual s may impact overall algorithm performance and are modeled as uncorrelated Gaussian noise in the ATBD for this product (although they may still contain correlated noise).
Additional details are provided in the ATBD.
Snapshot Retrieval
When the number of time-series input is set to 1, the time-series method becomes a conventional snapshot approach. εr and s are adjusted simultaneously at each time instance so as to minimize the distance between the data cube prediction and the observation. The snapshot retrieval algorithm was also tested with the two sets of field campaign data used in Section 3.2.3 of the ATBD. Figure 14a in the ATBD shows that the bare surface retrieval has an RMSE of 0.055 cm3/cm3, which is larger than the error of the time-series retrieval. Closer examination however indicates that the retrieval e rror increases for larger soil moisture. When the performance of the Monte Carlo simulation of the Kp effect is compared between time-series (Figure 12a of the ATBD) and snapshot (Figure 12c of the ATBD) retrievals, the performance difference is much greater than the in situ validation shows. The radar measurement error of the in situ data (~0.3 dB) is smaller than simulated by the Monte Carlo analysis, which may explain why the Monte Carlo analysis of the snapshot retrieval show the poorer performance.
The snapshot retrieval over the pasture surface is poor. Figure 14b in the ATBD shows a retrieval error of 0.105 cm3/cm3. It is possible that the structure of the grass was important in characterizing the backscatter to be very different from that of the bare surface. The possibility is supported also by the large retrieval error (not shown) when the retrieval was performed using the bare surface data cube (NMM3D) after treating the grass surface as if it is a bare surface. The evaluation results suggest the limited applicability of the snapshot retrieval, and highlights the need for time-series retrievals.
Optional Algorithm
Change Detection by Kim and van Zyl (KVZ)
The full details of this algorithm are described by Kim and van Zyl (2009). A brief summary is presented here. An expression will be derived to relate the backscattering cross section to soil moisture as:
mv = f(σhh, σvv)
where the normalized radar cross sections are expressed in decibels. Since this expression depends on the biomass level, the cross-polarization will be used to compensate the biomass variation over time: how this will be implemented is being studied . One potential choice for f is given by:
(Equation 1)
These two coefficients for each pixel will be determined using the expected minimum and maximum values for soil moisture and the time-series backscattering cross section data. Accommodating different vegetation phenology during the implementation of the above relationship remains as a task. The obvious assumption contained in Equation 1 is that there is a linear relationship between the radar cross-sections and the soil moisture. However, simulated data suggests that this relationship may be nonlinear. Ground radar and airborne data will be used to derive the final expression for the function f(σhh, σvv).
The slowly varying backscattering cross section component due to the biomass variation will be estimated and compensated using the cross-polarization backscattering cross section. Ground radar measurements will help derive the expres sion before the SMAP instrument is launched. However, since the g round measurements are limited, with simulated data various vegetation environments will be considered. When the two extreme conditions (completely dry and wet) of the soil moisture and radar backscatter are known exactly, a test with a corn field experiment shows that the soil moisture may be retrieved with an error of 0.026 cm3/cm3 (Table 3 of Kim and van Zyl). When the errors in the two extreme values are simulated with the random noise, the retrieval error reaches 0.05 cm3/cm3. Alternatively, Kim and van Zyl proposed that C1 at each radar pixel (e.g., at 3 km resolution) is inverted using the mv retrieved at a coarse scale (such as at 36 km resolution) using the radiometer records. The radar-radiometer method is the subject of the ATBD for the SPL2SMAP and SPL3SMAP products.
Additional details are provided in the ATBD.
Note: Algorithms that are currently not being used for these data include Snapshot, Dubois/van Zyl (DVZ), and Shi. Thus, the following data fields contain a fill value of -9999.0:
- soil_moisture_snapshot
- soil_moisture_snapshot_DVZ
- and soil_moisture_snapshot_shi
Algorithm Inputs and Outputs
Table 3 lists the required input data sets to the radar-only soil moisture inversion. These include the primary choices that meet the ancillary data requirements. Ancillary data sets are defined as external data sets that are required as inputs to SMAP retrieval algorithms in the generation of the standard Level-2 and Level-3 active products. Ancillary data needed by the SMAP mission fall into two categories: static ancillary data, which do not change during the mission, and dynamic ancillary data, which require periodic updates in timeframes ranging from seasonally to daily. Static data include parameters such as permanent masks (land/water/forest/urban/mountain), the grid cell average elevation and slope derived from a Digital Elevation Model (DEM), permanent open water fraction, and soils information (primarily sand and clay fraction). All static ancillary data are resampled to the 3 km EASE-Grid 2.0 as the output products. The dynamic ancillary data include land cover, precipitation, and effective soil temperatures. The choice of which ancillary data set to use for a particular SMAP product is based on a number of factors, including its availability and ease of use, and its inherent error and resulting impact on the overall soil moisture or freeze/thaw retrieval accuracy. Latency, spatial resolution, temporal resolution, and global coverage are also important. The choice of a primary source for each of the ancillary data parameters is fully documented in individual SMAP Ancillary Data Reports (refer to Section 6.1 in the ATBD). Brief descriptions are provided below for the input radar backscatter, and for the ancillary parameters generated by the SPL2SMA processor.
| Data Type | Data Source |
|---|---|
| Radar Backscatter (HH/VV/HV) | SPL1CS0 |
| Land Cover Class | MODIS-International Geosphere Biosphere Programme (IGBP) |
| Crop Type | Cropland Data Layer (US) Environment Canada (Canada) ECOCLIMAP (Europe) Global Crop Map (Rest of the world) |
| Vegetation Water Content | SPL2SMA MODIS Normalized Difference Vegetation Index (NDVI) |
| Data Cube | SMAP Science Team |
| Soil Texture | Harmonized World Soil Database + Regional + Food Agricultural Organization (FAO) |
| Soil Surface Temperature | Global Modeling and Assimilation Office (GMAO) or European Center for Medium Range Weather Forecasting (ECMWF) |
| Mountain Area | Combination of Shuttle Radar Topography Mission (SRTM), Alaska United States Geological Survey (USGS) DEM, Canada Geobase DEM, and GTOPO30 (a global DEM with a horizontal grid spacing of 30 arc seconds or approximately 1 km) |
| Freeze/Thaw Flag | GMAO |
| Permanent Ice | MODIS IGBP |
| Snow | Snow & Ice Mapping System (IMS) - NOAA |
| Static Water Body | MODIS MOD44W |
| Transient Water Body | SPL2SMA |
| Precipitation | ECMWF Total Precipitation Forecasts [or Global Precipitation Measurement (GPM) |
| Urban Area | Global Rural-Urban Mapping Project (GRUMP) |
| RFI Contamination | SPL1CS0 |
Radar Backscatter
As previously stated, the primary input to the SPL2SMA algorithm is the Level-1C gridded radar product. This product consists of the calibrated high-resolution (SAR-mode) multiple polarization radar backscatter measurements and associated parameters, Earth-located and mapped to a swath-based, along-track/cross-track, 1 km EASE-Grid 2.0. It is assumed that the backscatter measurements have been fully calibrated and corrected for atmospheric attenuation and ionospheric effects such as Faraday rotation. Information provided in the Level-1C data set or in ancillary calibration files should enable estimates of the sigma nought (σ0) accuracy (Kp) for each measurement based on instrument characteristics.
Vegetation Water Content derived from Radar Vegetation Index
The SPL2SMA algorithm requires information on the VWC. One potential source of vegetation information is a polarimetric parameter derived from SMAP data such as the Radar Vegetation Index (RVI) (Kim and van Zyl 2000).
Alternatively, the vegetation amount may be derived based on low-order polynomial of the NDVI and optical remote-sensing images (refer to the SPL2SMP ATBD). In principle the RVI or any other parameter based on the radar backscatter should correlate with the vegetation amount more strongly than the optical NDVI, because fundamentally optical radiation cannot penetrate vegetation canopies while the L-band microwave signal can. Furthermore the RVI is concurrent with the radar observation. Studies are ongoing with regard to whether the radar backscatter may improve estimates of vegetation effects. The Data Fields Ancillary Data group contains VWC derived from both RVI and NDVI.
Radar Water Bodies
An algorithm for discriminating transient water bodies from land surfaces is also part of the SPL2SMA processor. There are many combinations of water surfaces adjoined by land surfaces, (depending on both land and water surface roughness and presence of vegetation) that affect this discrimination capability. At present, the case of a calm water surface next to lightly vegetated land was studied. It was found that a threshold (3 dB) applied to the ratio of HH to VV effectively separated water from land at SMAP's 40 deg incidence angle. In this study, Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) images were analyzed over fairly flat terrain. Despite the apparent differences between UAVSAR and SMAP, the measurement quality for both instruments is determined by the same statistical processes inherent in SAR imaging. Calibration error and the number of single looks of one UAVSAR pixel are similar to those of one pixel of the SMAP SAR.The algorithm was tested over numerous locations with a detection error of smaller than 10% (Kim et al. 2014).
Refer to the Data Fields document for a description of all data fields.
Background color on
Processing Steps
This data set is generated by the SMAP Science Data Processing System (SDS) at the Jet Propulsion Laboratory (JPL) in Pasadena, California USA. To generate this data set, the processing software ingests only the 6:00 a.m. descending half-orbit files of the SMAP High-Resolution Radar Sigma Nought, Version 3 (SPL1CS0) data set. The SPL1CS0 data set contains swath-gridded radar backscatter observations acquired by the radar and processed using synthetic-aperture processing providing measurements at 1 km resolution. After the 1 km sigma0 data are gridded on the 3 km EASE-Grid 2.0, the gridded data are then tested for retrievability criteria according to input data quality, ancillary data availability, and land cover conditions. These criteria are tested by applying various quality flags.
Two quality flags are derived using the 3 km sigma0: radar vegetation index and radar water body. Static and dynamic ancillary data are sampled for each pixel of the SPL2SMA product. Sigma0 values from the past are sampled and used by the time-series algorithm. For each 3 km pixel, land cover class information is obtained from the mostly annual ancillary data. The land cover information allows one to choose an appropriate data cube for each pixel. The retrieval over the different land cover classes is spatially assembled to create a half-orbit output, followed by the conversion from the dielectric constant to soil moisture.
When retrievability criteria are met, the software invokes the baseline retrieval algorithm to generate soil moisture retrieval. Only cells that are covered by the actual swath for a given projection are written in the product.
For more information on each portion of the algorithm processing flow, and the details of the retrieval methods and options, refer to Section 3: Retrieval Algorithms of the ATBD for this product.
Background color on
Error Sources
Anthropogenic Radio Frequency Interference (RFI), principally from ground-based surveillance radars, can contaminate both radar and radiometer measurements at L-band. Early measurements and results from ESA's Soil Moisture and Ocean Salinity (SMOS) mission indicate that in some regions RFI is present and detectable. The SMAP radar and radiometer electronics and algorithms include design features to mitigate the effects of RFI. The SMAP radar utilizes selective filters and an adjustable carrier frequency to tune to predetermined RFI-free portions of the spectrum while on orbit.
More information about error sources is provided in Section 3.5: Error Budget of the ATBD.
Background color on
Quality Assessment
For in-depth details regarding the quality of these Version 3 Validated data, refer to the following reports:
Validated Assessment Report
Beta Assessment Report
Quality Overview
SMAP products provide multiple means to assess quality. Each product contains bit flags, uncertainty measures, and file-level metadata that provide quality information. For information regarding the specific bit flags, uncertainty measures, and file-level metadata contained in this product, refer to the Data Fields and Metadata Fields documents.
Each HDF5 file contains metadata with Quality Assessment (QA) metadata flags that are set by the Science Data Processing System (SDS) at the JPL prior to delivery to the 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.
A separate QA file with a .qa file extension is also associated with each data file. QA files are ASCII text files that contain statistical information in order to help users better assess the quality of the associated data file.
If a product does not fail QA, it is ready to be used for higher-level processing, browse generation, active science QA, archive, and distribution. If a product fails QA, it is never delivered to NSIDC DAAC.
Data Flags
Note: As SPL2SMA is a research product, the accuracy requirements are tentative and are set to 0.06 cm3/cm3 for soil moisture and up to 5 kg/m2 for vegetation water content (VWC).
Quality Control (QC) is an integral part of the SPL2SMA processing. The QC steps of the SPL2SMA processing are based on the flags that are provided with the input data stream (SPL1CS0), different types of masks, flags, and fractional coverage of other variables provided by ancillary data. The soil moisture retrieval flags are listed in Table 4.
| Sigma0 Flags |
|---|
| Sigma0 Quality (HH) |
| Sigma0 Quality (VV) |
| Sigma0 Quality (HV) |
| Soil Moisture Retrieval Flag: Key Bits |
| Retrieval Success |
| Freeze/Thaw |
| Radar Vegetation Index |
| Radar/Permanent Water Body |
| Surface Quality Flag: Key Bits |
| Dense Vegetation (VWC) |
| Nadir Region |
| Urban |
| Precipitation |
| Snow, Ice, Frozen Surface |
| Mountain |
The overall sigma0 quality flags are deemed good if the following flags are good: range, null data, RFI cleaned, RFI detected, and Kp. Faraday flag is ignored while evaluating the overall sigma0 flag. Generally, 0 is assigned to the favorable condition for soil moisture retrieval. The SPL2SMA processes all data that have favorable conditions for soil moisture retrieval, including no snow cover, no frozen ground, no permenant water, no wetland, no urban, and no snow. As long as sigma0 value in natural unit is greater than 0, soil moisture is retrieved. However, for soil moisture retrieval conducted for regions with high VWC and nadir regions, appropriate flags are added to these data points indicating their susceptibility to high errors. Therefore, there are three flagging states associated with soil moisture retrieval:
- Retrieve and assign good flag (retrieval recommend)
- Retrieve and assign bad flag (retrieval success)
- Do not retrieve and assign bad flag (retrieval attempt)
Radar Backscatter and RFI Flags
The SPL2SMA product also uses the radar backscatter and RFI flag from the SPL1CS0 product.
Freeze/Thaw Flag
The freeze/thaw flag for SPL2SMA is generated from modeled surface temperature information called TSURF from the Global Modeling and Assimilation Office (GMAO). The GMAO-modeled TSURF parameter is used due to the maturity of the SMAP radar freeze/thaw algorithm and availability of radar measurements.
Water Body Flag
The SPL2SMA water body flag is generated from a priori information on permanent open freshwater from the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD44W database. The MOD44W database is used due to the maturity of the SMAP radar open-water algorithm and availability of radar measurements.
Nadir Region Flag
The SPL2SMA product also has a nadir region degradation where the effective resolution in the cross-track direction gradually becomes 36 km. The flag, called the 3 km nadir region flag, is assigned to identify this region.
For more information regarding data flags, refer to the Data Fields document.
Background color on
References and Related Publications
Contacts and Acknowledgments
Investigators
Seungbum Kim, Jakob van Zyl, Roy S. Dunbar, Eni Njoku
Jet Propulsion Laboratory
California Institute of Technology
Pasadena, CA 91109 USA
Joel Johnson
Ohio State University
Columbus, OH 43210 USA
Mahta Moghaddam
University of Michigan
Ann Arbor, MI 48109 USA
Leung Tsang
University of Washington
Seattle, WA 98195 USA
FAQ
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.
Short name
Title
Subsetting... read more
How To
How do I programmatically request data services such as subsetting, reformatting, and reprojection using an API?
The subsetting, reformatting, and reprojection services provided by NSIDC through NASA Earthdata Search can also be accessed programmatically as a synchronous REST interface. This programmatic access is provided via an HTTPS URL containing a series of... read more
How to import and geolocate SMAP Level-1C and Level-2 data in ENVI
The following are instructions on how to import and geolocate SMAP Level-1C HDF5 data in ENVI.
Testing notes
Software: ENVI
Software version: 5.3 and above. If using version 5.3, service pack 5.3.1 is needed.
Platform: Windows 7
Data set: SMAP L1C... read more
How do I search, order, and customize SMAP data using Earthdata Search?
This video tutorial and attached PDF document provide step-by-step instructions on how to search, order, and customize SMAP data using Earthdata Search (https://search.earthdata.nasa.... read more
