MOD11_User_Guide_V6_2
MOD11_User_Guide_V6_2
Zhengming Wan
ERI, University of California, Santa Barbara
December 2013
Note:
This document is a living document that describes the Collection-6 (C6) MODIS Land
Surface Temperature (LST) products. It is revised as progress is made in the development
and assessment of the LST products. Described is the current state of the MODIS LST
products. The purpose of the document is to give the potential user of LST products an
understanding of the MODIS LST products and the data in those products.
In the early days, the MODIS LST products MOD11_L2, MOD11A1, and MOD11B1
had been validated at stage 1 with in situ measurements in more than 50 clear-sky cases
in the temperature range from -10oC to 58oC and the column water vapor range of 0.4-
4cm, most of them presented in published papers (Wan et al., 2002 and 2004; Coll et al.,
2005; Wan, 2006; Wan, 2008; Wan and Li, 2008). Detailed validation of the C6 MODIS
LST product is given in the most recently published paper (Wan, 2014).
Please use the C6 MODIS LST product in your applications because its quality is much
better than the qualities of the C4.1 LST product (generated by the V4 algorithm from C5
input data) and C5 products (generated by the V5 algorithm from C5 input data).
As shown in the following plots, the viewing angle dependence in the band 29 emissivity
values [red in daytime and green in nighttime data records] at the north Tassili site (latitude
27.0o, longitude 7.65o) in the C6 MYD11B1 LST product in 2007 [right] is much better than in
the C4.1 product [left] because the C6 LST/E product provides more stable results. Wan and Li
(2011) explained with spectral measurement results that the spectral transmittance of a
single ivy leaf (roughly 0.18-mm thickness) are negligibly small (<0.02) in the whole TIR
spectral region in its fresh condition and very small (<0.03) in the 6–14μm region in its
dry condition and that the emissivity of the single ivy leaf varies significantly with water
content only in the 3.5–6μm spectral region, and its variation in the 8–14μm region is
negligible. It is reasonable to expect the most top thin layer (with thickness in a few
tenths of one mm) of land surface in bare soil regions is staying in dry conditions during
the daytime in clear-sky days so that the surface emissivity will not change significantly
during daytime in clear-sky conditions.
Table of Contents
1. INTRODUCTION
1.1. File Format of LST Products
1.2. Sequence of LST Products
The major refinements implemented in the V6 daily LST Product Generation Executive
(PGE) code are given in (Wan, 2014). This V6 PGE code has been used in the
Collection-6 (C6) reprocessing and forward processing of the MODIS LST product.
Brief descriptions of the LST data products are given here to give perspective to the
sequence. Expanded descriptions of the LST products are given in following sections.
The first product, MOD11_L2, is a LST product at 1km spatial resolution for a swath.
This product is the result of the generalized split-window LST algorithm (Wan and
Dozier, 1996). Geolocation data (latitude and longitude) at a coarse resolution (5 lines by
5 samples) is also stored in the product. The second product, MOD11A1, is a tile of daily
LST product at 1km spatial resolution. It is generated by mapping the pixels in the
MOD11_L2 products for a day to the Earth locations on the sinusoidal projection. The
third product, MOD11B1, is a tile of daily LST and emissivities at 6km spatial resolution.
It is generated by the day/night LST algorithm (Wan and Li, 1997). The fourth product,
MOD11B2, is an eight-day LST product by averaging from two to eight days of the
MOD11B1 product. The fifth product, MOD11B3, is a monthly LST product generated
from the daily MOD11B1 product. The sixth product, MOD11A2, is an eight-day LST
product by averaging from two to eight days of the MOD11A1 product. The seventh
product, MOD11C1, is a daily global LST product in a geographic projection. It is
created by assembling the MOD11B1 daily tiles together and resampling the SDSs at
6km grids to the 0.05° spatial resolution of the Climate Modeling Grid (CMG) cells. The
eighth product, MOD11C2, is an eight-day composite of LST at the same resolution as
MOD11C1. The ninth product, MOD11C3, is a monthly composite of LST at the same
resolution as MOD11C2. Similarly, there are nine C6 LST products starting with MYD11
for the Aqua MODIS.
* The number 65535 in uint16 may be shown as -1s in 16-bit integer by some software
toolkits, for example, by ncdump in the HDF toolkit.
The effective calibration formula for the "LST" SDS is
LST = the SDS data in uint16 * 0.02, giving a value in the range of 150-1310.7K.
The four QAFraction PSAs are especially useful to granules in ocean regions where only
a small number of island pixels exist. Because the total number of land and coastal pixels
is highly variable in granules covering both land and ocean, the values of QAPercent and
QAFraction PSAs are calculated on the base of the total number of all pixels in a
granule. Therefore, we can always calculate how many pixels with LST in good quality
and other quality from these PSA values, even for island pixels.
The ECS global attribute ArchiveMetadata.0 contains information relevant to production
of the data product. It also contains an alternate bounding of geographic coverage of the
swath. These data may be useful in determining what version of the algorithm was used
to generate the product. Contents are described in Table 6.
Table 6. Listing of objects in the global attribute ArchiveMetadata.0 in MOD11_L2.
The StructMetadata.0 global attribute is used by the HDF-EOS toolkit to specify the
mapping relationships between the geolocation data and the LST data (SDSs). Mapping
relationships are unique in HDF-EOS and are stored in the product using HDF structures.
Description of the mapping relationships is not given here. Use of HDF-EOS toolkit,
other EOSDIS supplied toolkits may be used to geolocate the data.
Table 7. Listing of objects in the global attribute StructMetadata.0 in MOD11_L2.
Object Definition
DIMENSION_1 along_swath_lines_1km (10*nscans)
DIMENSION_2 Cross_swath_pixels_1km (Max_EV_frames)
DIMENSION_3 Coarse_swath_lines_5km (2*nscans)
DIMENSION_4 Coarse_swath_pixels_5km (Max_EV_frames/5)
GeoDimension= Coarse_swath_lines_5km
DIMENSIONMAP_1 DataDimension= along_swath_lines_1km
Offset=2 Increment=5
GeoDimension= Cross_swath_pixels_1km
DIMENSIONMAP_2 DataDimension= Cross_swath_pixels_1km
Offset=2 Increment=5
GEOFIELD_1 GeoFieldName=Latitude
GEOFIELD_2 GeoFieldName=Longitude
DATAFIELD_1 DataFieldName=LST
DATAFIELD_2 DataFieldName=QC
DATAFIELD_3 DataFieldName=Error_LST
DATAFIELD_4 DataFieldName=Emis_31
DATAFIELD_5 DataFieldName=Emis_32
DATAFIELD_6 DataFieldName=View_angle
DATAFIELD_7 DataFieldName=View_time
2.5. Quality Assurance
Indicators of quality are given in metadata objects in the CoreMetadata.0 global attribute
QA and in a quality control (QC) SDS, generated during production, or in post-product
scientific and quality checks of the data product. QA metadata objects in the
CoreMetadata.0 global attribute are the AutomaticQualityFlag and the
ScienceQualityFlag and their corresponding explanations. The AutomaticQualityFlag is
set according to rules based on data conditions encountered during a run of the LST
algorithm. Setting of this QA flag is fully automated. The rules used to set it are liberal;
nearly all of the data or intermediate calculations would have to be anomalous for it to be
set to "Failed". Typically, it will be set to "Passed". The ScienceQualityFlag is set post
production either after an automated QA program is run on the data product or after the
data product is inspected by a qualified LST investigator. Content and explanation of this
flag are dynamic so it should always be examined if present. A sampling of products will
be inspected. Sampling may be random, in support of field campaigns, or event driven.
The QC SDS in the data product provides additional information on algorithm results for
each pixel. The QC SDS unsigned 16-bit data are stored as bit flags in the SDS. This QC
information can be extracted by reading the bits in the 16-bit unsigned integer. The
purpose of the QC SDS is to give the user information on algorithm results for each pixel
that can be viewed in a spatial context. The QC information tells if algorithm results were
nominal, abnormal, or if other defined conditions were encountered for a pixel. The QC
information should be used to help determine the usefulness of the LST data for a user's
needs. The bit flags in the QC SDS are listed in Table 8.
Table 8. Bit flags defined in the QC SDS in the MOD11_L2 product. Note that bit 0 is the least
significant bit.
bits Long Name Key
00=Pixel produced, good quality, not necessary to examine
more detailed QA
01=Pixel produced, unreliable or unquantifiable quality,
1&0 Mandatory QA flags recommend examination of more detailed QA
10=Pixel not produced due to cloud effects
11=Pixel not produced primarily due to reasons other than cloud
(such as ocean pixel, poor input data)
00=good data quality of L1B in bands 31 and 32
01=missing pixel
3&2 Data quality flag
10=fairly calibrated
11=poorly calibrated, LST processing skipped
00=cloud free pixel
01=pixel only with thin cirrus
5&4 Cloud flag
10=fraction of sub-pixel clouds<= 2/16
11=LST affected by nearby clouds
00=generalized split-window method
01=day/night method
7&6 LST model number
10=high LST w/o atmospheric & emis corrections
11=cirrus effects corrected
00=no multi-method comparison
01=multi-method comparison done
9&8 LST quality flag
10=fair consistency
11=good consistency
00=inferred from land cover type
01=MODIS retrieved
11 & 10 Emissivity flag
10=TBD
11=default value used
00=emis quality not checked
01=emis quality checked with land cover type
13 & 12 Emis quality flag
10=emis quality checked with NDVI
11=emis view-angle dependence checked
00=error in emis_31 emis_32 <= 0.01
01=error in emis_31 emis_32 <= 0.02
15 & 14 Emis error flag
10=error in emis_31 emis_32 <= 0.04
11=error in emis_31 emis_32 > 0.04
It should be noted that fillvalue 0 listed for the SDS QC in Table 9 is valid for the bit
flags only when a fillvalue 0 is present in the SDS LST pixels (so the 00-01 bits in the
QC pxels have a value of 10 or 11). A value of 0 in the QC bit flags means good data
quality, cloud free, or small error in emis_31 and emis_32, and etc, if a pixel has a valid
LST value. We do not discriminate fillvalue 0 from valid value 0 for all bit flags in the
QC in order to minimize the data volume. Users should read SDSs LST and QC at the
same time in order to properly interpret their values in an easy way.
The ECS global attribute ArchiveMetadata.0 contains information relevant to the input
data, an alternate geographic coverage bounds, and information relevant to version of the
algorithm and product. A listing of objects along with sample values is given in Table 11.
Table 11. ArchiveMetadata.0 of the MOD11A1 data product.
Object Name Sample Value Comment
HorizontalTileNumber "08"
VerticalTileNumber "05"
AlgorithmPackageAcceptanceDate "102004"
AlgorithmPackageMaturityCode "Normal" Algorithm
AlgorithmPackageName "MOD_PR11A" Descriptors
AlgorithmPackageVersion "6"
"Moderate-Resolution Imaging
InstrumentName
SpectroRadiometer"
ProcessingDateTime "2006-07-27T23:09:54.000Z"
"MODIS Level-3 1km Land Surface
LongName
Temperature and Emissivity"
ProcessingCenter "MODAPS"
SPSOParameters "2484 and 3323"
only keep the
LocalInputGranuleID "1920,0535,0540,0715,1740,1915"
time IDs
EastBoundingCoordinate -103.941403
WestBoundingCoordinate -130.540731 Extent of the
NorthBoundingCoordinate 39.995833 tile coverage.
SouthBoundingCoordinate 30.004167
The StructMetadata.0 global attribute (Table 12) is used by the HDF-EOS toolkit to
create the mapping relationships between the defined grid and data (SDSs).
Table 12. Listing of objects in the global attribute StructMetadata.0 in MOD11A1.
Object Definition
DIMENSION_1 XDim = 1200
DIMENSION_2 YDim = 1200
DataField_1 DataFieldName=LST_Day_1km
DataField_2 DataFieldName=QC_Day
DataField_3 DataFieldName=Day_view_time
DataField_4 DataFieldName=Day_view_angle
DataField_5 DataFieldName=LST_Night_1km
DataField_6 DataFieldName=QC_Night
DataField_7 DataFieldName=Night_view_time
DataField_8 DataFieldName=Night_view_angle
DataField_9 DataFieldName=Emis_31
DataField_10 DataFieldName=Emis_32
DataField_11 DataFieldName=Clear_day_cov
DataField_12 DataFieldName=Clear_night_cov
Table 13. Bit flags defined for SDSs QC_day and QC_Night in MOD11A1. Note that bit 0 is the
least significant bit.
bits Long Name Key
00=LST produced, good quality, not necessary to examine
more detailed QA
01=LST produced, other quality, recommend examination of
1 & 0 Mandatory QA flags
more detailed QA
10=LST not produced due to cloud effects
11=LST not produced primarily due to reasons other than cloud
00=good data quality
01=other quality data
3 & 2 Data quality flag
10=TBD
11=TBD
00=average emissivity error <= 0.01
01=average emissivity error <= 0.02
5 & 4 Emis Error flag
10=average emissivity error <= 0.04
11=average emissivity error > 0.04
00=average LST error <= 1K
01=average LST error <= 2K
7 & 6 LST LST Error flag
10=average LST error <= 3K
11=average LST error > 3K
4. MOD11A2 Eight-day LST
An eight-day compositing period was chosen because twice of such period is the exact
ground track repeat period of the Terra platform. LST over eight days is the averaged
LSTs of the MOD11A1 product over eight days.
Note: * a negative sign before the zenith view angle indicates that the MODIS views the
Earth surface from east.
** The 03-00 four bits are for view angle flag of the companion observation (which is
night observation if LST_Day_6km valid, or day observation otherwise), with key 0-15
as index of view angle sub-range of the companion observation: 0-7 if view from east (0
at the west end of scan line); 8-15 if view from west (15 at the east end of scan line). The
06-04 three bits are used for a flag of time difference between the day and night MODIS
observations, with key 0-7 as couples of days, 7 representing the range from 7 to 16 (i.e.,
14-32 days). Bit 07 is for DEM slope flag, with key 0=DEM slope not considered or
1=DEM slope considered in the day/night algorithm. There are eight view angle sub-
ranges in the whole range of viewing zenith angle at the surface from 0-65o and they are
0-10, 10-20, 20-30, 30-39, 39-47, 47-54, 54-60, and 60-65 degrees.
5.3. Local Attributes
The local attributes for SDSs LST_Day_6km and LST_Night_6km are similar to those in
Table 5.
Terra/Aqua 0=no
3
Combined-use flag 1=yes
00=average emissivity error <= 0.01
01=average emissivity error <= 0.02
5 & 4 Emis Error flag
10=average emissivity error <= 0.04
11=average emissivity error > 0.04
00=average LST error <= 1K
01=average LST error <= 2K
7 & 6 LST LST Error flag
10=average LST error <= 3K
11=average LST error > 3K
* Each bit in the 8-bit unsigned integer indicates clear-sky (1) or not (0) in the
corresponding day or night. Bit 00 is for the first day or night, and bit 07 is for the last
day or night in the 8-day period.
Z. Wan and J. Dozier, 1996, "A generalized split-window algorithm for retrieving land-
surface temperature from space", IEEE Trans. Geosci. Remote Sens., v34, n4, pp. 892-
905.
Z. Wan and Z.-L. Li, 1997, "A physics-based algorithm for retrieving land-surface
emissivity and temperature from EOS/MODIS data", IEEE Trans. Geosci. Remote Sens.,
v35, n4, pp. 980-996.
W. C. Snyder, Z. Wan, Y. Zhang and Y.-Z. Feng, 1998, "Classification-based emissivity
for land surface temperature measurement from space", Int. J. Remote Sens., v19, n14,
pp. 2753-2774.
W. C. Snyder and Z. Wan, 1998, "BRDF models to predict spectral reflectance and
emissivity in the thermal infrared", IEEE Trans. Geosci. Remote Sens., v36, n1, pp. 214-
225.
Z. Wan, Y. Zhang, X. Ma, M. D. King, J. S. Myers, and X. Li, 1999, "Vicarious
calibration of the Moderate-Resolution Imaging Spectroradiometer Airborn Simulator
thermal infrared channels", Appl. Optics, v38, n20, pp. 6294-6306.
Z. Wan, 2002, "Estimate of noise and systematic error in early thermal infrared data of
the Moderate Resolution Imaging Spectroradiometer (MODIS)", Remote Sens. Environ.,
80, 47-54.
Z. Wan, Y. Zhang, Z.-L. Li, R. Wang, V.V. Salomonson, A. Yves, R. Bosseno, and J. F.
Hanocq, 2002, "Preliminary estimate of calibration of the Moderate Resolution Imaging
Spectroradiometer (MODIS) thermal infrared data using Lake Titicaca", Remote Sens.
Environ., 80, 497-515.
Z. Wan, Y. Zhang, Q. Zhang, and Z.-L. Li, 2002, "Validation of the land-surface
temperature products retrieved from Terra Moderate Resolution Imaging
Spectroradiometer data", Remote Sens. Environ., 83, 163-180.
Z. Wan, Y. Zhang, Q. Zhang, and Z.-L. Li, 2004, "Quality assessment and validation of
the MODIS land surface temperature", Int. J. Remote Sens., 25, 261-274.
C, Coll, V. Caselles, J.M. Galve, E. Valor, R. Niclos, J.M. Sanchez, and R. Rivas, 2005,
"Ground measurements for the validation of land surface temperatures derived from
AATSR and MODIS data", Remote Sens. Environ., 97, 288-300.
Z. Wan, 2008, "New refinements and validation of the MODIS land-surface
temperature/emissivity products ", Remote Sens. Environ., 112, 59-74.
Z. Wan, 2014, "New refinements and validation of the Collection-6 MODIS land-surface
temperature/emissivity products ", Remote Sens. Environ., 140, 36-45.
Z. Wan, & Z.-L. Li, 2008, “Radiance-based validation of the V5 MODIS land-surface
temperature product”, International Journal of Remote Sensing, 29, 5373–5393.
Z. Wan, & Z.-L. Li, 2011, Chapter 25, MODIS land surface temperature and emissivity.
In B. Ramachandran, C. O. Justice, & M. J. Abrams (Eds.), Land remote sensing and
global environmental change, NASA's Earth observing system and the science of ASTER
and MODIS, NASA's Earth observing system and the science of ASTER and MODIS.
New York Dordrecht Heidelberg London: Springer, http://dx.doi.org/10.1007/978-1-
4419-6749-7.