The goal of the MODIS/VIIRS Land Product Subsets project is to provide summaries of selected MODIS/VIIRS Land Products for the community to use for validation of models and remote-sensing products and to characterize field sites.
You may download MODIS Sinusoidal Grid (KMZ) or get more information about the formats and subsets offered by clicking the tabs below.
MODIS Land Subset Products Offered
Detailed information about file naming conventions, QC information, and file formats is provided in the product links. The Land Processes DAAC provides information on all MODIS Land products .
|Product||Acronym||Name||Collection||Frequency (Days)||Resolution (Meter)|
|MCD12Q1||LC||MODIS/Terra+Aqua Land Cover Type Yearly L3||6||annual||500|
|MCD12Q2||LCD||MODIS/Terra+Aqua Land Cover Dynamics Yearly L3||5||annual||500|
|MCD15A2H||LAI/FPAR||MODIS/Terra+Aqua Leaf Area Index/FPAR 8-Day L4||6||8||500|
|MCD15A3H||LAI/FPAR||MODIS/Terra+Aqua Leaf Area Index/FPAR 4-Day L4||6||4||500|
|MCD19A3||MAIAC||MODIS/Terra+Aqua BRDF/Albedo Model Parameters 8-Day L3||6||8||1000|
|MCD43A||BRDF||[Calculated] MODIS/Terra+Aqua BRDF/Albedo 16-Day L3||6||daily||500|
|MCD43A4||NBAR||MODIS/Terra and Aqua Nadir BRDF-Adjusted Reflectance Daily L3||6||daily||500|
|MOD09A1||SREF||MODIS/Terra Surface Reflectance 8-Day L3||6||8||500|
|MOD11A2||TEMP||MODIS/Terra Land Surface Temperature/Emissivity 8-Day L3||6||8||1000|
|MOD13Q1||NDVI/EVI||MODIS/Terra Vegetation Indices 16-Day L3||6||16||250|
|MOD15A2H||LAI/FPAR||MODIS/Terra Leaf Area Index/FPAR 8-Day L4||6||8||500|
|MOD16A2||ET||MODIS/Terra Evapotranspiration 8-Day L4||6||8||500|
|MOD17A2H||GPP||MODIS/Terra Gross Primary Production 8-Day L4||6||8||500|
|MOD17A3H||NPP||MODIS/Terra Net Primary Production Yearly L4||6||annual||500|
|MOD21A2||LSTE||MODIS/Terra Land Surface Temperature/3-Band Emissivity 8-Day L3||6||8||1000|
|MYD09A1||SREF||MODIS/Aqua Surface Reflectance 8-Day L3||6||8||500|
|MYD11A2||TEMP||MODIS/Aqua Land Surface Temperature/Emissivity 8-Day L3||6||8||1000|
|MYD13Q1||NDVI/EVI||MODIS/Aqua Vegetation Indices 16-Day L3||6||16||250|
|MYD15A2H||LAI/FPAR||MODIS/Aqua Leaf Area Index/FPAR 8-Day L4||6||8||500|
|MYD16A2||ET||MODIS/Aqua Evapotranspiration 8-Day L4||6||8||500|
|MYD17A2H||GPP||MODIS/Aqua Gross Primary Production 8-Day L4||6||8||500|
|MYD17A3H||NPP||MODIS/Aqua Net Primary Production Yearly L4||6||annual||500|
|MYD21A2||LSTE||MODIS/Aqua Land Surface Temperature/3-Band Emissivity 8-Day L3||6||8||1000|
Land Surface Temperature and Emissivity (LSTE) products (MOD21A2 and MYD21A2) are temporarily unavailable from the ORNL DAAC’s MODIS Subsetting Tool because of a processing error. Upon resolution of science data quality issues, these data will be reprocessed and made available.
VIIRS Land Subset Products Offered
Detailed information about file naming conventions, QC information, and file formats is provided in the product links. The Land Processes DAAC provides information on all VIIRS Land products .
|Product||Acronym||Name||Collection||Frequency (Days)||Resolution (Meter)|
|VNP09H1||SREF||VIIRS/S-NPP Surface Reflectance 8-Day L3||1||8||500|
|VNP09A1||SREF||VIIRS/S-NPP Surface Reflectance 8-Day L3||1||8||1000|
|VNP13A1||NDVI/EVI||VIIRS/S-NPP Vegetation Indices 16-Day L3||1||16||500|
|VNP15A2H||LAI/FPAR||VIIRS/S-NPP Leaf Area Index/FPAR 8-Day L4||1||8||500|
Daymet Subset Products
Detailed information about file naming conventions, QC information, and file formats is provided in the product links. The ORNL DAAC processes and distributed Daymet Data.
|Daymet||Daymet||Daymet: Daily Surface Weather Data on a 1-km Grid for North America||3||daily||1000|
Tabular Data Format
Format and example data values for [PRODUCT], [BAND], filtered_scaled, and land cover .csv files provided with the subset order.
|Column number||Example values||Description|
|Unique row identifier|
|3||A2000049||Acquisition Date [ A (YYYYDDD) ]|
|4||Lat35.9625Lon-85.35051Samp7Line7||User selected center point coordinates and specified width (Samp) and height (Line) of bounding rectangle in pixels. Width x height denotes number of product values starting in Column 7. (e.g., 7 x 7 =49)|
|5||2015136104550||File Processing Date (YYYYDDDHHMMSS)|
|Product Scientific Data Set (Band): Indicates type of values to follow. Specific values vary by product.|
|7 to n||NDVI Example
[BAND].csv: 4690,4780,4802,4364, … to n
filtered_scaled_[BAND].csv: 0.469,0.478,F,0.4364, … to n
|* Data values of type as specified. Number of data columns as given in Column 4.|
|* Value "F" indicates a filtered pixel due to failed QC criteria or land cover masking|
Time series for all bands of the requested product are provided in one file. These data are the starting point for deriving the other subset data products. The [PRODUCT].csv file is formatted according to Table 1.
Time series for each band are split into a separate file, called [BAND].csv, for each band. The [BAND].csv file is formatted according to Table 1.
Filtered and Scaled Data
Time series for data are provided for science data bands with a scale factor and QC filters applied. Pixels that fail QC criteria are set to "F". Pixels corresponding to land cover different than that of the subset center are set to "F" for files named "*_landcovermask.csv". The filtered_scaled files are formatted according to Table 1.
Land Cover Subsets
Four .csv files containing land cover data from the MCD12Q1 Land Cover Type 1 (IGBP) product are provided with each subset order. Land cover data are formatted following Table 1.
|12q1.csv||Land cover from 2013 for subset area|
|12q1_201x201.csv||Land cover from 2013 for a 201 x 201 km region centered on the user-defined subset center|
|12q1_[1km,500m,250m].csv||Land cover from 12q1.csv resampled to match the spatial resolution of the user-requested product.|
|12q2.csv||MCD12Q2 Land Cover Dynamics timeseries product providing vegetation phenology parameters|
Phenology files are provided with each order. Data are presented, in tabular format, for average day of greenness increase, maximum, decrease, and minimum for the subset area for each detected phenology cycle derived from MCD12Q2. Phenology files may contain up to two cycles per year.
|Phenology.csv||Average day by day of year|
|Phenology_CalendarDates.csv||Average day by calendar date|
Contents of the phenology .csv files are formatted according to Table 2.
|Column header||Description||Example value
|Year||year of row||2001||2001|
|Cycle||full vegetation cycles (maximum of 2)||Cycle 1||Cycle 1|
|Increase||mean day of greenness increase||76||Mar 17|
|Maximum||mean day of greenness maximum||137||May 19|
|Decrease||mean day of greenness decrease||195||Jul 21|
|Minimum||mean day of greenness minimum||308.5||Nov 14|
Statistical Summary of Subset Area
Statistical summaries of the data for each step in the time series for science data bands are provided with the subset order. Pixels corresponding to land cover different than that of the subset center are set to "F" for files named "*_landcovermask.csv". Two .csv files are provided per science data band.
Contents of the statistics .csv files are formatted according to Table 3.
|Column number||Column heading||Description|
|1||date[YYYYDDD]||Acquisition date. Four-digit year and three-digit day of the year|
|2||date||Acquisition calendar date: YYYY-MM-DD|
|4||band||Scientific Data Set. Refers to a specific wavelength range measured by the sensor (Band1, 620-670 nm) or to a specific product derived from several measured bands (MOD15A2: Lai_500m and Fpar_500m; MOD09A1: sur_refl_b01, sur_refl_b02, etc.) or a data variable.|
|5||min||Minimum value of pixels within the subset area|
|6||max||Maximum value of pixels within the subset area|
|7||sum||Sum of pixels within the subset area|
|8||range||Range of pixels within the subset area|
|9||mean||Mean value of pixels within the subset area|
|10||variance||Variance of pixels within the subset area|
|11||standard_deviation||Standard deviation of pixels within the subset area|
|Columns 12-14 for statistics_[BAND].csv|
|12||tot_pixels||Total number of pixels in the requested area|
|13||pixels_pass_qa||Number of pixels in the requested area passing quality filters|
|14||per_cent_pixels_pass_qa||(Number of pixels with acceptable quality in the requested area / total number of pixels) * 100|
|Columns 12-14 for statistics_[BAND]_landcovermask.csv|
|12||tot_pixels||Total number of pixels in the requested area having the same land cover class as the center pixel|
|13||pixels_pass_qa||Number of pixels in the requested area having the same land cover class as the center pixel and passing quality filters|
|14||per_cent_pixels_pass_qa||(Number of pixels with the same land cover as the center pixel and acceptable quality in the requested area / total number of pixels having the same land cover class as the center pixel) * 100|
R Code and Graphics
An R script is provided for each science data band for generating the plots (also provided as .png images) listed in Table 4. R plots are concatenated in a single .pdf file available on the .csv Data page and can be previewed by choosing the "R Code and Graphics" tab at the bottom of the order visualization page.
|R_script[BAND].R||R scripts for generating plots for each science band|
|R_MODIS_Plots_[BAND]_Fig1_timeseries.png||Plot of subset mean values and variation over time|
|R_MODIS_Plots_[BAND]_Fig2_timeseries.png||Plot of subset mean values summarized by month|
|R_MODIS_Plots_[BAND]_Fig3_timeseries.png||Plot of subset mean values summarized by year|
|R_MODIS_Plots_[BAND]_Fig4_timeseries.png||Frequency histogram plot of all subset area mean values|
|R_MODIS_Plots_[BAND]_Fig5_timeseries.png||Time series plot with years stacked and a loess smoothed fit curve for each year|
|R_MODIS_Plots_[BAND]_Fig6_timeseries.png||Plot subset area mean grouped by year and month|
Product GeoTIFF Data and Landcover GeoTIFF Data download links provide compressed tar archive files (GTiff.tar.gz and GTiff_mcd12q.tar.gz) containing subset geotiffs.
Provides one geotiff per band for each date in the time series generated from the rows in [PRODUCT].csv. Files are named as follows:
|[Source MODIS filename]_[BAND].tif||MOD13Q1.A2000049.h11v05.006.2015136104550_250m_16_days_NDVI.tif|
Provides a geotiff for the MCD12Q1 Land Cover Type 1 (IGBP) product from 2013 named as follows:
|[Source MODIS filename]_Land_Cover_Type_1.tif||MCD12Q1.A2013001.h17v02.005.2011085002532_Land_Cover_Type_1.tif|
Citations are provided to allow this specific subset order to be referenced in a publication. Filenames are citation.bib and citation.txt
Aggregate 500m Land Cover (MCD12Q1) to 1km
For the Global subsetting tool, we provide two summaries of data: one is a summary for all pixels in the selected area and the other is the summary for all pixels with the same land cover class as the center pixel (the pixel that contains the field site or flux tower).
The land cover (MCD12Q1) is available at 500m. To obtain a summary of pixels in the selected area with the same land cover class at 1km resolution, the 500m pixels are aggregated to 1km pixels. Below is an image showing the 500m land cover for a 7 x 7 km subset region. The example subset region has a center pixel of "(1) Evergreen Needleef Forest".
The aggregation from 500m to 1km pixel is done by identifying 1km pixels that have 2 or more 500m pixels having the same land cover class as the 1km center pixel. Below is a 1km landcover image showing the aggregation. The black pixels cells are pixels that did not satisfy the aggregation criteria.
Below is a composite image showing the 500m landcover pixels with the aggregated 1-km pixels. The darker regions are 1km pixels that did not satisfy the aggregation criteria.
Method Used to Calculate Albedo from MCD43A1 (MODIS Albedo Model Parameter)
The ALBEDO data are based on the Terra and Aqua(Combined) MODIS BRDF / Albedo Model Parameter Product (MCD43A1). The primary BRDF model parameters from MCD43A1 are used in the following equations to calculate Black-sky, white-sky, and actual (blue-sky) albedo. We have used a solar zenith angle (szn) equal to local solar noon and an optical depth of 0.2 as default values in calculating these data files.
Black-sky Albedo = Parameters_01 + Parameters_02 * (-0.007574 + (-0.070987 * szn2) + (0.307588 * szn3)) + Parameters_03 * (-1.284909 + (-0.166314 * szn2) + (0.041840 * szn3))
White-sky Albedo = Parameters_01 + Parameters_02 * (0.189184) + Parameters_03 * (-1.377622)
Blue-sky (Actual) Albedo = White-sky Albedo * f(optical depth, solar zenith angle, aerosol type, band) + Black-sky Albedo * (1 - f(optical depth, solar zenith angle, aerosol type, band))
Parameters_01 corresponds to the MCD43A1 band "BRDF_Albedo_Parameters_nn.Num_Parameters_01",
Parameters_02 corresponds to the MCD43A1 band "BRDF_Albedo_Parameters_nn.Num_Parameters_02",
Parameters_03 corresponds to the MCD43A1 band "BRDF_Albedo_Parameters_nn.Num_Parameters_03"
for nn=vis (Band_08 Visible Broadband), nn=nir (Band_09 Near Infrared Broadand) and nn=shortwave (Band_10 Shortwave Broadband)
Additional information can be found on the LP DAAC Web site:
and at the BRDF Science Team Web site:
The Fixed Sites Subsets Tool provides customized subsets and visualizations of several MODIS/VIIRS land products and Daymet daily surface weather data for selected field and flux tower sites around the world.Spatial subsets of the land products are available for predefined areas of about 8x8-km centered on the selected sites. The site-specific data can be downloaded in both comma separated value (*.csv) and JSON formats. The data are in the native projection system.
Access using the Network IDs
The Fixed Sites Tool uses its own unique IDs for its 2000+ sites. Most networks (such as FLUXNET, NEON, etc) use their own unique IDs for their sites. For some networks, the data and visualization from the Fixed Sites Tool can be accessed directly using the network IDs. Here is how:
To list all the available sites in a network, go to the following URL:
To access a fixed site page using the network ID, go to the following URL:
To access a product subset using the network ID, go to the following URL:
To access a product subset data for a network site using the Web Service , please send the GET request to:
Additional detail for the web service access is available at the REST API Documentation.
Currently, this feature is available for the following networks:
A network site becomes inactive when they are decommissioned or stops collecting data for a considerable period of time. We will stop updating the sensor data (MODIS, VIIRS) for the inactive sites, and will archive the data for the period the site was operational. If you think a site should not be inactive, please let us know.