Get Data
DOWNLOADING DATA VIA HTTPS
To learn more about Earthdata Login and register for an account, please see How to Register with Earthdata Login. 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.
SMAP L3 Radiometer Northern Hemisphere Daily 36 km EASE-Grid Freeze/Thaw State, Version 1
This Level-3 (L3) product provides a daily composite of Northern Hemisphere landscape freeze/thaw conditions retrieved by the Soil Moisture Active Passive (SMAP) radiometer from 6:00 a.m. descending and 6:00 p.m. ascending half-orbit passes. SMAP L-band brightness temperatures are used to derive freeze/thaw state and transition data, which are then resampled to an Earth-fixed, Northern Hemisphere azimuthal 36 km Equal-Area Scalable Earth Grid, Version 2.0 (EASE-Grid 2.0).
First public data release
Geographic Coverage
Spatial Coverage: |
|
---|---|
Spatial Resolution: |
|
Temporal Coverage: |
|
Temporal Resolution: | 1 day |
Parameter(s): |
|
Platform(s) | SMAP Observatory |
Sensor(s): | SMAP L-Band Radiometer |
Data Format(s): |
|
Version: | V1 |
Data Contributor(s): | Xu, X., R. S. Dunbar, C. Derksen, A. Colliander, Y. Kim, and J. S. Kimball. |
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.
Xu, X., R. S. Dunbar, C. Derksen, A. Colliander, Y. Kim, and J. S. Kimball. 2016. SMAP L3 Radiometer Northern Hemisphere Daily 36 km EASE-Grid Freeze/Thaw State, Version 1. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: http://dx.doi.org/10.5067/RDEJQEETCNWV. [Date Accessed].Detailed Data Description
Freeze/thaw state and the direction of diurnal freeze/thaw transitions (frozen in the morning to thawed in the afternoon and vice versa) derived from brightness temperatures are output on a fixed Northern Hemisphere azimuthal EASE-Grid 2.0 at 36 km. Freeze/thaw state, the occurrence of freeze/thaw transitions, and the direction of transitions are expressed in boolean values (0 or 1). For freeze/thaw state, 0 indicates thawed conditions and 1 indicates frozen. For freeze/thaw transition state, 0 indicates the a.m. and p.m. FT states are the same (thawed/thawed or frozen/frozen). The transition direction flag is only meaningful if there is a transition (transition state = 1), and is set to 0 for a.m. frozen/p.m. thawed and 1 for a.m. thawed/p.m. frozen.
Also included are brightness temperatures (TBs) in kelvin for a 36 km EASE-Grid 2.0 cell.
Refer to the Data Fields document for details on all parameters.
Data are in HDF5 format. For software and more information, including an HDF5 tutorial, visit the HDF Group's HDF5 Web site.
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-3 radiometer freeze/thaw 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.
Data element arrays are three dimensional, with the exception of transition_direction and transition_state_flag arrays, which are two dimensional. All arrays have dimensions of 500 rows and 500 columns in each a.m. and p.m. layer. For the a.m./p.m. index of the array, the a.m. layer is assigned to the index value 0 and the p.m. layer is assigned to index value 1.
Ancillary Data
Includes all ancillary data, such as landcover classification and open water body fraction.
Freeze/Thaw Retrieval Data
Includes freeze/thaw data, latitude and longitude arrays, and associated quality assessment flags.
Radiometer Data
Includes all radiometer data and associated 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.
Files are named according to the following convention, which is described in Table 1:
SMAP_L3_FT_P_yyyymmdd_RLVvvv_NNN.[ext]
For example:
SMAP_L3_FT_P_20170117_R14010_001.h5
Where:
Variable | Description | ||||||||
---|---|---|---|---|---|---|---|---|---|
SMAP |
Indicates SMAP mission data | ||||||||
L3_FT_P |
Indicates specific product (L3 : Level-3; FT : Freeze/Thaw; P : Passive) |
||||||||
yyyymmdd |
4-digit year, 2-digit month, 2-digit day; date in Universal Coordinated Time (UTC) of the first data element that appears in the product. | ||||||||
RLVvvv |
Composite Release ID (CRID), where:
|
||||||||
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:
|
Each file is approximately 6.1 MB.
The daily data volume is approximately 6.1 MB.
Coverage for this data set spans the Northern Hemisphere for all land regions north of 45°N latitude, and from 180°W to 180°E. The gap in coverage at the North Pole, called a pole hole, has a radius of approximately 400 km. The swath width is 1000 km, enabling nearly complete coverage of the Northern Hemisphere every three days.
Spatial Coverage Map
Figure 2 shows the spatial coverage of this data set.

The native spatial resolution of the radiometer footprint is 36 km. Data are then gridded using the 36 km Northern Hemisphere azimuthal EASE-Grid 2.0 projection.
EASE-Grid 2.0
These data are provided on the Northern Hemisphere azimuthal EASE-Grid 2.0 (Brodzik et al. 2012). Each grid cell has a nominal area of approximately 36 x 36 km2 regardless of longitude and latitude. Using this projection, all data arrays have dimensions of 500 rows and 500 columns.
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.

Coverage spans from 31 March 2015 to the present.
Temporal Coverage Gaps
Satellite and Processing Events
Due to instrument maneuvers, data downlink anomalies, data quality screening, and other factors, small gaps in the time series will occur. Refer to the SMAP On-Orbit Events List for Instrument Data Users page for details regarding these gaps.
Latencies
Each Level-3 file is a daily composite of half-orbit files/swaths.
To ensure complete coverage of the freeze/thaw domain in each daily file, a.m. and p.m. data for the current day are combined with a.m. and p.m. data from previous days. A maximum of three days of past data is used, and is necessary only near the southern margin of the freeze/thaw domain.
Software and Tools
For tools that work with SMAP data, refer to the Tools Web page.
Data Acquisition and Processing
For a detailed description of the SMAP instrument, visit the SMAP Instrument page at the JPL SMAP Web site.
SMAP Level-3 radiometer freeze/thaw data (SPL3FTP) are derived from SMAP L1C Radiometer Half-Orbit 36 km EASE-Grid Brightness Temperatures, Version 3 (SPL1CTB).
The SPL3FTP product is derived using a temporal change detection approach that has been previously developed and successfully applied using time-series satellite remote sensing radar backscatter and radiometric brightness temperature data from a variety of sensors and spectral wavelengths. The approach is to identify the landscape freeze/thaw (F/T) state via the temporal response of the normalized polarization ratio (NPR) of the brightness temperature, which is sensitive to changes in the dielectric constant of the landscape that occur as the water within the components transitions between frozen and non-frozen conditions.
This approach assumes that the large changes in dielectric constant occurring between frozen and non-frozen conditions dominates the corresponding NPR temporal dynamics across the seasons, rather than other potential sources of temporal variability such as changes in canopy structure and biomass, large precipitation events, or changes in soil moisture.
This section has been adapted from Dunbar et al. (2016), the Algorithm Theoretical Basis Document (ATBD) for this data set.
The SMAP Level-3 radiometer freeze/thaw data set is derived from SMAP L1C Radiometer Half-Orbit 36 km EASE-Grid Brightness Temperatures, Version 3 (SPL1CTB). The derivation of freeze/thaw from SMAP brightness temperature measurements occurs during an intermediate Level-2 processing step of the input Level-1 brightness temperature data. During the Level-2 processing step, the freeze/thaw algorithm utilizes a seasonal threshold approach to convert SMAP brightness temperature measurements to freeze/thaw state. For an overview of the steps involved in processing this data product, refer to Figure 5 in the Processing Steps section.
Baseline Algorithm
The SPL3FTP freeze/thaw baseline algorithm examines the time series progression of the brightness temperature signature relative to signatures acquired during seasonal reference frozen and thawed states. The algorithm is applied to the normalized polarization ratio (NPR) of SMAP radiometer measurements
NPR=(TBV-TBH)/(TBV+TBH) (Equation 1)
A seasonal scale factor D(t) is defined for an observation acquired at time t as
Dt=(NPR(t)-NPR(fr))/(NPR(th)-NPR(fr)) (Equation 2)
where NPR(t) is the normalized polarization ratio calculated at time t, for which a freeze/thaw classification is sought, and NPR(th) and NPR(fr) are normalized polarization ratios corresponding to the frozen and thawed reference states, respectively. The twenty highest (lowest) NPR values from SMAP radiometer measurements during July and August 2015 (thaw) and January and February 2016 (freeze) for the northern (≥45°N) domain were retained and averaged to create the thaw (freeze) reference. Data were separated by ascending and descending orbit. The methodological approach to NPR freeze and thaw references will be refined in future product releases. In addition, the reference values will be updated following each transition season. Reprocessing of the SMAP data record incorporating annual variations in the SMAP freeze/thaw reference states should improve product accuracy over the use of static reference conditions. The initial SMAP freeze and thaw NPR references are shown in Figure 4 below (and Section 4.2.2 of the ATBD).

SMAP radiometer (a) freeze and (b) thaw references; (c) reference difference between panels (a) and (b). Units are NPR scaled by 100.
A threshold level T is then defined such that
D(t) > T
D(t) <=T (Equation 3)
defines the thawed and frozen landscape states, respectively. This series of equations (1-3) are run on a grid cell-by-cell basis for unmasked portions of the FT domain. The output from Equation (3) is a dimensionless binary state variable designating either frozen or thawed conditions for each unmasked grid cell. The threshold values can be optimized on a grid cell by grid cell basis, but are fixed at 0.5 for this release. Optimization approaches will be evaluated in advance of future product releases.
Following the pixel-wise determination of freeze/thaw state, two additional processing steps are applied to mitigate summer season false freeze and winter season false thaw retrievals. First, if the brightness temperature magnitude at either V or H pol is greater than 273 K, the pixel is set to thaw regardless of the retrieval. Second a temporally fixed 'never frozen' mask calculated from daily AMSR-E freeze/thaw maps (using the approach described in Kim et al., 2012) is applied to remove obviously false summer freeze flags. False freeze retrievals occur in some regions of the F/T domain because of small differences between the reference freeze and thaw values (see Figure 4c). Implementation of daily AMSR-E-derived ‘never frozen’ and ‘never thawed’ masks in future product releases will address this issue.
Ancillary Data
Ancillary data sets are used to:
- Support initialization of the thresholds employed in the algorithm
- Set flags that indicate potential problem regions
- Define masks where no retrievals should be performed
Ancillary data used in SPL3FTP processing includes data sets of inland open water, permanent ice and snow, and urban areas in order to derive masks so that no retrievals occur over these regions. Ancillary data sets of mountainous areas, fractional open water cover, and precipitation are used to derive flags so that a confidence interval can be associated with the retrieval. All ancillary data sets are resampled to a spatial scale and geographic projection that matches the SPL3FTP product in accordance with the guidelines of the SMAP mission.
Ancillary data sets used for SPL3FTP data processing were in place prior to launch, with no need for periodic updates during post-launch operations. A continuous surface map of fractional area of open water was used to represent fractional water coverage within a grid consistent with the resolution and projection of the SPL3FTP product. For the SPL3FTP development, the lake fraction threshold within a grid cell was set to 50%. Determination of a physically-based lake fraction will be finalized for a forthcoming SPL3FTP release. Table 3 lists the ancillary data employed in support of SPL3FTP production. Similar ancillary data were used for production of the SMAP radar freeze/thaw (SPL3FTA) product.
Data Type | Data Source(s) | Frequency | Resolution | Extent | Use |
---|---|---|---|---|---|
Vegetation Type | Moderate Resolution Imaging Spectroradiometer International Geosphere Biosphere Programme (MODIS-IGBP) | Once | 250 m | Global | Sensitivity Analysis |
Precipitation |
European Center for Medium Range Weather Forecasting (ECMWF) Forecasts | Time of Acquisition | 0.25 degrees | Global | Sensitivity Analysis |
Static Water Bodies |
MODIS Land-Water Mask (MODIS44W) | Once | 250 m | Global | Mask/Flag |
Mountainous Areas | NASA Global Digital Elevation Model (DEM) | Once | 30 m | Global | Mask/Flag |
Permanent Ice and Snow | MODIS-IGBP Permanent Ice and Snow Class | Once | 500 m | Global | Mask/Flag |
Seasonal Snow | National Oceanic and Atmospheric Administration/National Ice Center Interactive Multisensor Snow and Ice Mapping System (NOAA IMS) | Daily | 1 km | Northern Hemisphere | Flag |
Never-Thawed/Never-Frozen Masks | Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E) | Annual | 25 km | Northern Hemisphere | False Flag Mitigation |
For more information, refer to the ATBD for this product.
This product is generated by the SMAP Science Data Processing System (SDS) at the Jet Propulsion Laboratory (JPL) in Pasadena, California USA. Figure 5 shows the processing sequence for generation of the SMAP L3 freeze/thaw (F/T) radiometer product (SPL3FTP).
The derivation of freeze/thaw from SMAP brightness temperature measurements occurs during an intermediate Level-2 processing step of the input Level-1 brightness temperature data. During the Level-2 processing step, the freeze/thaw algorithm utilizes a seasonal threshold approach to convert SMAP brightness temperature measurements to freeze/thaw state.
To generate this product, the processing software:
- Ingests one day's worth of Level-1 files and creates individual Northern Hemisphere composites as two-dimensional or three-dimensional arrays for each output parameter defined in the Level-1 data
- Intermediate Level-2 processing step: Converts SMAP brightness temperature measurements to freeze/thaw state. Classifies frozen and thawed landscape states on a grid cell-by-cell basis for unmasked portions of the FT domain by:
- Utilizing the NPR of SMAP radiometer measurements during seasonal reference frozen and thawed states
- Applying a fixed threshold of 0.5 to determine either frozen or thawed conditions relative to the reference states
- Employing ancillary data sets to set flags for potential problem regions, and define masks where no retrievals should be performed
- Mitigating summer season false freeze and winter season false thaw retrievals by:
- Designating pixels as 'thaw' when TB magnitude at V or H pol is greater than 273 K
- Applying a fixed AMSR-E derived 'never frozen' mask
- The processing software then combines a.m. and p.m. data for the current day with a.m. and p.m. data from previous days to ensure complete coverage of the freeze/thaw domain in each daily file. Note that a maximum of three days of past data is used, and is necessary only near the southern margin of the freeze/thaw domain. Wherever data overlap occurs, as is typical at high latitudes, data which were acquired closest to 6:00 a.m. and 6:00 p.m. local solar times are chosen.
For details regarding each of these processing steps, refer to the Derivation Techniques and Algorithms section of this document.

As a result, the output Level-3 radiometer freeze/thaw product distinguishes four levels of freeze/thaw conditions determined from the ascending 6:00 a.m. and descending 6:00 p.m. SPL1CTB data, including:
- Frozen (from both a.m. and p.m. overpass times)
- Non-frozen (a.m. and p.m.)
- Transitional (a.m. frozen; p.m. non-frozen)
- Inverse-transitional (a.m. non-frozen; p.m. frozen)
For more information on the algorithm processing flow, refer to the ATBD for this product, Section 2.2: L3_FT_P Production.
Anthropogenic Radio Frequency Interference (RFI), principally from ground-based surveillance radars, can contaminate both radar and radiometer measurements at L-band. The SMAP radar and radiometer electronics and algorithms include design features to mitigate the effects of RFI. The SMAP radiometer utilizes selective filters and an adjustable carrier frequency to tune to predetermined RFI-free portions of the spectrum while on orbit.
The landscape freeze/thaw state retrieval represented by the SPL3FTP algorithm and products characterizes the predominant frozen or non-frozen state of the land surface within the sensor Field of View (FOV) and does not distinguish freeze/thaw characteristics among different landscape elements, including surface snow, soil, open water or vegetation. The lower frequency L-band retrievals from SMAP are expected to have greater sensitivity to surface soil freeze/thaw conditions under low to moderate vegetation cover. Microwave freeze/thaw sensitivity is strongly constrained by intervening vegetation biomass, soil moisture levels, and snow wetness. Ambiguity in relating changes in the radiometer signal to these specific landscape components is a challenge to validation of the freeze/thaw product (Colliander et al., 2012). In northern boreal and tundra landscapes, L-band penetration depth is greater under frozen conditions when land surface liquid water levels are low, and markedly reduced under thawed conditions due to characteristically moist surface organic layer and soil active layer conditions, even under relatively low tundra vegetation biomass levels (Du et al. 2014).
Note that spatial classification error is expected to be larger in regions with small differences between frozen and thawed NPR references. This includes areas where: freeze/thaw is ephemeral and densely vegetated areas due to vegetation scattering effects on microwave emissivity, which reduces TB V and H-polarization differences and NPR dynamic range. In regions of complex terrain, freeze/thaw heterogeneity is greater which also adversely impacts retrieval performance. In arid regions, the small amount water present in the thawed state makes the soil permittivity close to the frozen state, which can cause false freeze retrieval errors.
Finally, a major assumption of the seasonal threshold based temporal change freeze/thaw classification is that the major temporal shifts in brightness temperature are caused by land surface dielectric changes from temporal freeze/thaw transitions. This assumption generally holds for higher latitudes and elevations where seasonal frozen temperatures are a significant part of the annual cycle and a large constraint to land surface water mobility and ecosystem processes (e.g., Kim et al. 2012). However, freeze/thaw classification accuracy is expected to be reduced where other environmental factors may cause large temporal shifts in brightness temperature, including large rainfall events and surface inundation, and changes in vegetation biomass (e.g. phenology, disturbance and land cover change). While there is a strong NPR response to freeze/thaw transitions, NPR is not stable during summer due to the influence of vegetation, soil moisture, etc. Depolarization of summer season measurements leads to false freeze retrievals that must be mitigated. Winter season false thaw in areas of complex terrain are due to uncertainty in the references due to sub-grid heterogeneity.
For an assessment of algorithm performance and sources of uncertainty using in situ observations, and other satellite data sets, refer to the Validated Assessment Report for this product.
For in-depth details regarding the quality of these Version 1 Validated data, refer to the following report:
Validated Assessment Report
Quality Overview
The SPL3FTP product has sufficient fidelity and accuracy to identify the primary seasonal freeze and thaw transitions, and distinguish diurnal freeze/thaw state changes common during seasonal transitions.
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.
References and Related Publications
Contacts and Acknowledgments
Investigators
Xiaolan Xu, Scott Dunbar, Andreas Colliander
Jet Propulsion Laboratory
California Institute of Technology
Pasadena, CA 91109 USA
Chris Derksen
Climate Research Division
Environment and Climate Change Canada
Toronto, ON M3H 5T4 Canada
John Kimball, Youngwook Kim
Numerical Terradynamic Simulation Group (NTSG)
College of Forestry & Conservation
The University of Montana
Missoula, MT 59812-1049 USA
FAQ
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.
... read moreThe following table describes the data subsetting, reformatting, and reprojection services that are currently available for SMAP data via the NASA Earthdata Search tool.