# Presentation: 756672

## MODIS

- Land Surface Temperature and Emissivity Algorithms and Products
- Zhengming Wan
- University of California, Santa Barbara
- January 17-18, 2007

- MODIS Collection 5 Workshop, University of Maryland

## Basic Considerations in MODIS LST Algorithms 

- Basic Considerations in MODIS LST Algorithms

- 1. LST is retrieved from TIR data only in clear-sky conditions.

-   LST is not mixed with cloud-top temperature in the atmospheric product
-   (TIR signal from surface cannot penetrate clouds to reach satellites).

- 2. LST is defined by the radiation emitted by the land surface
-     observed by MODIS at the instantaneous viewing angle.

- Applications may need LST at different angles (nadir or 50o from nadir).

- 3. Proper resolving of the land-atmosphere coupling is the key
-     in retrieving surface & atmospheric properties.

- Use multi-bands in the atmospheric windows for the LST retrieval.

- The values of atmospheric temperature and water vapor are useful
- to improve the LST retrieval.  However, there may be large errors in
- these values.  Use them as indicates of ranges or initial guesses only.

- Integrated retrieval is possible but it takes a lot of computing time.

- 4. Input data: MOD021KM, MOD03, MOD07, MOD10, MOD12, MOD35 & MOD43.

## Ts = C + ( A1 + A2           + A3         ) 

***T******s****** = C + *****(***** A*****1***** + A*****2*****           ******+ A*****3 *****        *****) **

-  **+ ****(***** B*****1***** *****+***** B*****2             *****   ******+****** ******B*****3*****          *****)**

_***1 – ***__***ε***_
***   ******ε***

_***Δε***_
***ε******2***

_***1 – ***__***ε***_
***   ******ε***

_***Δε***_
***ε******2***

_**T**__**11μ**__** – T**__**12μ**_
**       ****2**

***Δε***

***=***

***ε******11******μ***

***ε******12******μ***

- where

- ־

_**T**__**11μ**__** **__***+***__** T**__**12μ**_
**       ****2**

- The generalized split-window algorithm (Wan and Dozier, 1996) in form

- MODIS LST Algorithms (1)

- - emissivities estimated from land cover types (Snyder et al., 1998; Snyder & Wan, 1998).

- A MODIS pixel may cover several 1km grids with different land cover types.

***ε****** = 0.5 ( ******ε******11μ ***+*** ******ε******12μ****** )    ***and

- Emissivities vary even within a land cover type (crop lands may have different
- soils and crops in variable coverage).

- - coefficients Ai, Bi, and C depend on viewing zenith angle (in range of 0-65o).

- - coefficients also depend on ranges of the air surface temperature and column

- water vapor.

- - only process pixels in clear-sky at different MOD35 confidences over land or

- in lakes.

## MODIS LST Algorithms (2)

- The generalized split-window algorithm (continue)

- MODIS LST Algorithms (2)

- The coefficients were obtained from regression analysis of the MODIS
- simulation data created by atmospheric radiative transfer code MODTRAN4
- in wide ranges of surface and atmospheric conditions.

- As a thumb rule, regression models work well only within the space that is
- covered by the data participated in the regression analysis, their performance
- may degrade near the boundary, and they perform poorly outside the space.

- An accurate split-window LST algorithm may be developed only if the
- simulation data are in good quality and they cover a simulation space wider
- than the real physical space of the surface emissivity/temperature changes,
- and the atmospheric temperature and water vapor changes.

- The above statements also apply to a sub-space.  There should be enough
- overlaps between sub-spaces in order to reduce the sensitivity to uncertainties.

## MODIS LST Algorithms (3) 

- MODIS LST Algorithms (3)

-  The MODIS day/night LST algorithm

- (Wan & Li, 1997)

-  retrieve Ts-day, Ts-night, & band emissivities simultaneously

- - The range of viewing zenith angle is separated into 4, 5, or 16

- - be able to adjust the input atmospheric cwv and Ta values.

- is performed for grids larger than MODIS pixels:

- (bands 20, 22, 23, 29, and 31-33).

- - Option for combined use of Terra and Aqua MODIS data in v5.

- - least square-sum fitting 14 observations to solve 13 variables:

- Ts-day, Ts-night, cwv and Ta values for day and night,
- emissivities in the first six bands (small surface effect in b33)
- and a BRDF factor in the first three bands.

- sub-ranges in v3, v4, and v5, respectively.

- with day & night data in seven bands

- - Terrain slope considered in v5 QA.

## Major refinements implemented in the V5 daily LST code (PGE16) 

|    Specification / Action |   in V4 |                       in V5	 |
| --- | --- | --- |
| grid size of LST/emissivities in M*D11B1 retrieved from the day/night algorithm | 5km x 5km (exactly 4.63km) | 6km x 6km (exactly 5.56km) |
| number of sub-ranges of zenith view angles | 5 for the whole scan swath | 2x8 for the whole scan swath |
| effect of slope in the M*D11B1 grid | not considered | considered in the QA |
| temporal averaging in the 1km LST product | yes | no |
| option of combined use of Terra and Aqua data in the day/night algorithm | no | yes |
| incorporate the split-window method into the day/night algorithm | partially with landcover-based
em31, em32 and initial Ta, cwv | fully with em31, em32, Ta and cwv as variables in the iterations |
| clear-sky pixels defined by MODIS cloudmask | at 99% confidence over land
  at 66% confidence over lakes |  at confidence of >= 95% 
   over  land  <= 2000m
 at confidence of >= 66%
   over  land > 2000m
 at confidence of >= 66%
   over lakes |
| removing cloud-contaminated LSTs | not implemented |  implemented for M*D11A1
    and M*D11B1 |
| empirical optical-leak correction to band 32 | not implemented |  made for the last four pixels each scan
 line in Terra MODIS L1B granules |

- Major refinements implemented in the V5 daily LST code (PGE16)

## Remove cloud-contaminated LSTs in C5 level-3 LST products with 

- Remove cloud-contaminated LSTs in C5 level-3 LST products with

- constraints (δΤ) on the temporal variations in clear-sky LSTs (in PGE16C)

|       δΤ (K) | Description of land-cover  (type #) |
| --- | --- |
| 3.0 |  water (0) |
| 7.6 |  evergreen needleleaf forest (1) |
| 7.2 |  evergreen broadleaf forest  (2) |
| 7.2 |  deciduous needleleaf forest (3) |
| 7.0 |  deciduous broadleaf forest (4) |
| 7.0 |  mixed forest (5) |
| 8.0 |  closed shrublands (6) |
| 9.0 |  open shrublands (7) |
| 8.4 |  woody savannas (8) |
| 9.0 |  savannas (9) |
| 9.0 |  grasslands (10) |
| 5.0 |  permanent wetlands (11) |
| 8.0 |  croplands (12) |
| 8.0 |  urban and built-up (13) |
| 8.0 |  cropland and mosaics (14) |
| 4.0 |  snow and ice (15) |
| 11 |  bare soil and rocks (16) |
| 10 |  unclassified (17) |

- - In step 1, remove the worst LSTs that are different from the 32-day maximum by more than 4 times the δΤ value or different from the 16-day maximum by more than 3 times the δΤ value.
- - In step 2, remove the LSTs that are different from the 8-day maximum by more than 2 times the δΤ value, then calculate the 8-day average value of the remaining LSTs.
- - In step 3, remove the LSTs that are different from the 8-day average value by more than the δΤ value.
- - In order to consider the larger natural temporal changes in clear-sky LSTs in the growing and drying seasons, and in cold regions, the δΤ values are adjusted on the basis of statistical values of mean and standard deviation of LSTs in four periods of 8 days for each land-cover type in the tile under processing.

## Improvements of the C5 LST Products over C4 (I)

- due to using cloudmask combined with surface elevation

- C5

- C4

- shown in example of MOD11B1 in tile h25v05 retrieved
- from Terra MODIS data acquired on 21 January 2003.

- daytime LST

- nighttime LST

- emissivity RGB

## Improvements of the C5 LST Products over C4 (II)

- due to applying an empirical correction for optical leak to band 32 in the last four pixels each scan line in the Terra MODIS L1B granules, where the leak cannot be corrected
- by a physical model, shown in nighttime LSTs in MOD11A1.A2003194.h11v04. Note that LST values in the dark stripe (right) are cooler than their neighboring by 2-4K.

- C5

- C4

## Improvements of the C5 LST Products over C4 (III)

- Cloud-contaminated LSTs were removed in level-3 C5 LST products. Daytime LSTs in MOD11B1.A2003203.h17v05 (7/22/03) are shown below. The minimum daytime LST value is 264.68K in the right image (C4) or 284.26K in the left image (C5). Note that cloud-contaminated
- LSTs in the C5 level-2 LST products are not removed yet.

- C5

- C4

## Validation of the C5 LST Products generated in V5 tests

- By comparisons of LST values in the C5 MOD11_L2 and MYD11_L2 products with the in-situ values in Wan et al., 2002; Wan et al., 2004; Coll et al., 2005, and radiance-based validation results over Railroad Valley, NV in June 2003 and a grassland in northern TX in April 2005. LST errors < 1K in most cases.

- See details in Wan (2007)

- - In M*D11_L2, if valid LSTs are available in both C4 &
- C5, their difference is less than 0.2-0.4K in most cases.

- - In M*D11A1 within latitude ±28o(MODIS orbits w/o overlapping), if valid LSTs are available in both C4 & C5, their difference is less than 0.2-0.4K in most cases. Outside the latitude region, if valid LSTs are available in both C4 & C5 and at the same view time (indicating temporal average not applied in C4), their difference is less than 0.2-0.4K in most cases. Users should remove cloud-contaminated LSTs in the C4 product before using them in applications.

- - LSTs severely contaminated by clouds were removed from level-3 C5 products, but not from all C4 products.
- It is very difficult to remove such LSTs from the 8-day C4 M*D11A2 products because the cloud contamination effect may be reduced in the 8-day averaging.

- Notes for applications of C4 & C5 LST products:

## Comparison of the C4 and C5 CMG LST Products (M*D11C)

- The LSTs in M*D11C products are based on the LSTs retrieved by the day/night algorithm and supplemented by the LSTs retrieved by the split-window algorithm.

- LST_day in C4 MOD11A1.A2001089

- LST_day in C5 MOD11A1.A2001089 (3/30)

- Highlights:

- - C4 & C5 daytime LSTs have different spatial
- distributions because V4 & V5 PGEs use the
- MODIS cloudmask differently.

- - The mean and standard deviation of the differences between LSTs retrieved by the two algorithms are less than 0.2K and 0.5K in V5 so the 6km LSTs from the day/night algorithm can be validated in-directly.

- - However, they are larger (about 1.5K and 1.8K) in V4 (the effects of aerosol and cloud contaminations propagate into clear-sky days thru the initialization with the affected lower emis values in the V4 day/night algorithm).

## Possible Future Enhancements of LST Products in C6

-  refine the existing algorithms and V5 PGEs for the LST/E products from Terra and Aqua MODIS data to improve the stability for accuracy and flexibility in multiple options to inputs so that the code will be suitable for the near real-time processing of MODIS and NPP/NPOESS VIIRS data in order to generate consistent long-term ESDR/CDR of LST/E products.
-  remove cloud-contaminated data records from levels 2 and 3 LST/E products.
-  develop methods to analyze and correct the effects of thin cirrus clouds and aerosols above the average loading (the MODIS aerosol product will be improved by the deep blue algorithm in C6).
-  measure surface emissivity spectra in the field with the sun-shadow method, and make more validation of the LST products with the radiance-based approach.

## A short list of references

- Coll, C., Caselles, V., Galve, J.M.., Valor, E., Niclòs, R., Sánchez, J.M., &   Rivas, R. (2005). Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data. *Remote Sensing of Environment*, 97, 288-300.

- Snyder, W.C., Wan, Z., Zhang, Y., & Feng, Y.-Z. (1998). Classification-based emissivity
- for land surface temperature measurement from space. *International Journal of Remote*
*Sensing*, 19, 2753-2574.

- Wan, Z., & Dozier, J. (1996). A generalized split-window algorithm for retrieving land-surface temperature from space. *IEEE Trans. Geoscience and Remote Sensing,* 34, 892–905.

- Wan, Z., & Li,  Z.-L. (1997). A physics-based algorithm for retrieving land-surface
- emissivity and temperature from EOS/MODIS data. *IEEE Trans. Geoscience and*
*Remote Sensing*, 35,  980-996.

- Wan, Z., Zhang, Y., Zhang, Y.Q., & Li, Z.-L. (2002). Validation of the land-surface temperature products retrieved from Moderate Resolution Imaging Spectroradiometer data. *Remote Sensing of Environment*, 83, 163-180.

- Wan, Z., Zhang, Y., Zhang, Y.Q., & Li, Z.-L. (2004). Quality assessment and validation of the global land surface temperature. *International Journal of Remote Sensing*, 25, 261-274.

- http://www.icess.ucsb.edu/modis/LstUsrGuide/usrguide.html

- Wan, Z., (2007). New refinements and validation of the MODIS land-surface temperature/
- emissivity products. *Remote Sensing of Environment*, in press.