# Presentation: 703996

## MODIS C4/C5 LST Status Report

- MODIS C4/C5 LST Status Report
- and Proposed C6 Processing Plan
- Zhengming Wan
- University of California, Santa Barbara

- MODIS/VIIRS ST Meeting, BWI, Maryland, May 14-16, 2008

## 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.

## Basic Considerations in MODIS LST Algorithms (II) 

- Basic Considerations in MODIS LST Algorithms (II)

- 6. The split-window LST algorithm is the primary algorithm used to generate the MODIS LST products including the level-2 product (MOD11_L2 and MYD11_L2) and the 1km level-3 product (MOD11A1 and MYD11A1) because the surface emissivities vary within narrow ranges in the spectral ranges of MODIS bands 31 and 32 for all land-cover types but at different widths.

- 5. The standard MODIS cloudmask product (MOD35 and MYD35) is used         because it is one of the best cloudmask products available and it is good
- in most cases (> 98%). However, mismasking (cloudy pixels as clear-sky pixels at 99% confidence, or clear pixels at 99% confidence as uncertainly clear) exist in all coudmask products.  Therefore, it is important to properly use this cloudmask product and remove LSTs under effects of cloud contaminations after the initial LST generation in the production system.

- 7. Surface emissivities do not significantly change with time because            thermal infrared radiation almost does not penetrate a thin vegetation leaf and its reflectance does not change with the water content in the leaf in the spectral range above 7μm and the skin of sands and soil lands always stay in dry condition normally in clear-sky days.

- Transmittance (in plot A), reflectance and emissivity (both in plot B) spectra of a single ivy leaf (about 0.18 mm thick) in fresh and extremely dry conditions from measurements with a spectroradiometer and integral sphere system by the UCSB MODIS LST group

## 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 99% confidence over snow/ice
 at confidence of >= 95% 
   over non-snow/ice land  <= 2000m
 at confidence of >= 66%
   over non-snow/ice 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.

- This scheme may not work well in wet regions because
- few clear-sky days with valid LSTs exist there.

## Improvements of the C5 LST Products over C4

- 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

## 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
- C4 & C5 MODIS cloudmask products.

- - 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).

## 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 (2008) and Wan and Li (2008)**

- - **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 ****±28****o****(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:

## Radiance-based approach to valid LSTs 

**2.   This approach has been validated by comparisons to the conventional temperature-based approach in large homogeneous fields (lakes, grassland, and rice fields) by the MODIS LST group and César Coll’s group.**

**It uses measured or estimated** **surface emissivity spectra, and atmospheric water vapor** **and temperature profiles in atmospheric radiative transfer (RT) simulations to invert the band-31 brightness temperature (Tb31) of the MODIS observation to a LST value, and compare this value with the LST value in the MOD11 product.**

**3.   It is better than the conventional temperature-based approach because of the large spatial variation in LSTs especially in the daytime and the small horizontal variations in atmospheric temperature and water vapor profiles in clear-sky and flat areas.**

**4.   It is important to use the atmospheric profiles appropriate to MODIS observations by constraints of time (within 2 hours) and distance (≤ 100km).**

**It is important to perform the R-based approach on a series of days and check the quality with ****(Tb31 – Tb32) ****values**: **in successful RT simulations not only Tb31 but also Tb32 of MODIS observations should be well matched simultaneously.**

**7.   The R-based approach has larger uncertainties in wet conditions and it does not work well if the atmospheric profile used in the atmospheric RT simulation is significantly different from the real atmospheric state along the path of MODIS observation. This is the same problem faced by all the T/E methods in the atmospheric correction for extracting emissivity spectra (Li et at., 1999).**

**6.  The NCEP atmospheric profiles may be used In the R-based validation ****with care****.**

## Changes of column water vapor (cwv) over In-salah, Algeria (27.22oN, 2.5oE) during Jan. – Dec. in 2005 and 2007

- (3) The cwv values calculated from NCEP modeled profiles are larger than the values calculated from radio sounding profiles (12Z) by 23-42% in the first 7 cases of 2005. But the difference is smaller in most cases of 2007. Therefore, we should be very careful about the large errors in NCEP profiles in wet conditions.

- Changes of column water vapor (cwv) over In-salah, Algeria (27.22oN, 2.5oE) during Jan. – Dec. in 2005 and 2007

- (1) There is a larger seasonal variation in cwv over this desert site in 2005.

- (2) The cwv values provided by two methods (integrated and direct) in MOD07_L2 may be quite different,
- especially in the wet year 2005.

- (4) If there are large temporal variations
- in the surface emissivities retrieved from
- satellite TIR data over the desert site,
- the major source of errors may be in the
- process of atmospheric correction,
- especially when NCEP modeled profiles
- are used in relatively wet conditions.

## High quality of the Bomem TIR spectroradiometer (MR100) and MODTRAN simulations

**Above, The Bomem TIR spectroradiometer (MR100) was deployed in Railroad Valley, NV in June 2003. It is similar to the instrument used to validate the**
**MODIS SST product. Right, Brightness temperature (Tb) of downwelling radiance at nadir measured by the TIR spectroradiomter at 10:30 PDT 6/29/03 in Railroad Valley and the Tb calculated by MODTRAN4 based on measured atmospheric profiles.**

**High quality of the Bomem TIR spectroradiometer (MR100) and MODTRAN simulations**

## Comparison between temperature-based (in-situ LST) and radiance-based (TOA L31 inverted LST) validations for the Terra MODIS LST product at the rice field in Valencia, Spain, in six ideal clear-sky cases in 2003-2004 (Coll et al., 2005). The atmospheric temperature and water vapor profiles used in the radiance-based validation were measured by radio sounding in Murcia, Spain (38.0o N, 1.17o W), about 160km south of Valencia. The mean and standard deviation of the differences between the in-situ LST and TOA L31 inverted LST values in these ideal clear-sky cases are 0.03K and 0.37K, respectively. The radiance-based approach is also validated by César Coll with more in-situ LSTs and atmospheric profiles measured by radiosode balloons launched at the site in multiple years.

| 
caseno. |   Granule ID    | date &
 time
(m/d hh:mm) | viewing zenithangle (o) | radiosondecwv, Ts-air(in MOD07)(cm, K) | MODISLST
 (δT)(K) |  in-situ LST (K) |  | 
TOA L31
inverted
LST(K) | 
profile-
basedTb31-Tb32(K) |      MODIS      Tb31 -
  Tb32     (K)
 |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 1 | 2003192.1040  |  7/11   10:42 | 27.7 |  1.98, 309.5
(2.16,302.4) |   302.7
   (0.3) | 302.1 |  | 302.7 | 0.23 | 0.44 |
| 2 | 2003224.1040 |  8/12    10:42 | 28.1 |  1.40, 306.8
(1.46,297.5)  |   304.3
   (0.1) | 304.4 |  | 304.1 | 0.05 | 0.49 |
| 3 | 2004190.1020 |  7/08   10:24 | 50.3 |  1.15, 302.8
(1.56,294.4) |   298.6 
   (0.5) | 298.5 |  | 298.9 | 0.14 | 0.52 |
| 4 | 2004209.1050 |  7/27   10:54 | 6.0 |  2.09, 307.1
(1.84,296.9) |   301.5
   (0.2) | 301.1 |  | 301.2 | 0.17 | 0.56 |
| 5 | 2004216.1100 |  8/03   11:00 | 6.0 |  2.18, 308.0
(2.34,298.5) |   303.1
   (0.2) | 303.2 |  | 302.8 | 0.34 | 0.78 |
| 6 | 2004225.1050 |  8/12   10:54 | 5.7 |  2.22, 306.9
(1.94,298.5) |   301.9
   (0.1) | 301.9 |  | 301.7 | 0.26 | 0.58 |

- Comparison between temperature-based (in-situ LST) and radiance-based (TOA L31 inverted LST) validations for the Terra MODIS LST product at the rice field in Valencia, Spain, in six ideal clear-sky cases in 2003-2004 (Coll et al., 2005). The atmospheric temperature and water vapor profiles used in the radiance-based validation were measured by radio sounding in Murcia, Spain (38.0**o** N, 1.17**o** W), about 160km south of Valencia. The mean and standard deviation of the differences between the in-situ LST and TOA L31 inverted LST values in these ideal clear-sky cases are 0.03K and 0.37K, respectively. The radiance-based approach is also validated by César Coll with more in-situ LSTs and atmospheric profiles measured by radiosode balloons launched at the site in multiple years.

## Radiance-based validation of V5 Terra  MODIS LSTs at the CarboEurope site in Hainich, Germany (51.0792°N, 10.4522°E). The atmospheric temperature and water vapour profiles used in the radiance-based validation were measured by radio sounding in Meiningen, Germany (50.56°N, 10.38°E), about 60km south of Hainich.  In column 6, (δT)n indicates the standard deviation of MODIS LSTs in n pixels with valid LSTs in the four pixels surrounding the site. The mean and RMS of errors (in column 8) are      -0.3K and 0.73K. This set of validation results clearly indicates that the LST cannot be validated by longwave radiation data measured at 44m above surface as Wang et al did for the same site in a recent paper published in Remote Sensing of Environment, v112, 623-635, 2008 (bias -2.2K and RMSE 2.5K) .

| 
case
no. | 

Granule ID | 
local solar
date & time
(d/m h:m) | viewing
zenith
angle
(o) | Radiosonde cwv, Ts-air

and values in (M*D07)
(cm, K) | 
MODIS
LST (δT)n
(K) | 
MOD L31 
inverted
LST
(K) | LSTmodis–
MOD L31
inverted
LST
(K) | 
profile-based Tb31-
Tb32
(K) | MODIS
Tb31-Tb32
(K) |  
column 10  - column 9
(K) |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24 | A2004021.2135
A2004023.1015
A2004046.2130
A2004047.1100
A2004077.0935
A2004077.2045
A2004107.0945
A2004107.2100
A2004135.2125
A2004137.1000
A2004166.2040
A2004167.1010
A2004198.2040
A2004199.1010
A2004228.1120
A2004228.2050
A2004258.0955
A2004262.2040
A2004292.0940
A2004292.2050
A2004303.1100
A2004303.2035
A2004350.1015
A2004350.2130 | 21/1 22:22
23/1 10:57
15/2 22:12
16/2 11:42
17/3 10:17
17/3 21:27
16/4 10:27
16/4 21:42
14/5 22:07
16/5 10:42
14/6 21:22
15/6 10:52
16/7 21:22
17/7 10:52
15/8 12:02
15/8 21:32
14/9 10:37
18/9 21:22
18/10 10:22
18/10 21:32
29/10 11:42
29/10 21:17
15/12 10:57
15/12 22:12 | 41.0
20.8
34.3
43.6
56.9
26.3
48.3
9.8
26.5
36.1
33.5
21.3
33.4
21.2
56.5
18.2
42.5
33.3
52.3
18.0
43.4
39.6
12.8
34.1 | 0.37,268.3 (0.46,265.3)
0.70,268.1 (0.65,269.1)
0.90,274.4 (0.92,274.0)
0.72,276.8 (0.96,276.4)
0.74,286.2 (0.94,284.3)
1.22,281.2 (1.34,282.9)
1.02,288.0 (1.00,284.9)
0.89,282.2 (1.06,279.8)
1.08,281.2 (1.09,278.7)
1.03,286.8 (1.60,286.0)
2.28,288.0 (3.25,286.3)
2.00,294.2 (3.39,291.0)
2.20,288.6 (3.09,288.3)
2.62,298.0 (2.94,294.9)
2.00,293.2 (2.59,291.9)
1.67,287.8 (2.33,285.7)
2.23,291.0 (2.35,288.7)
1.22,283.0 (2.05,287.6)
1.10,281.2 (0.94,280.0) 0.78,278.6 (0.88,271.1)
1.41,287.2 (1.61,285.5) 1.72,281.8 (1.30,280.1)
0.48,268.1 (0.59,275.9)
0.85,268.5 (0.67,273.8)  | 264.2  (1.69)3
269.1  (0.62)2
273.0   (---)0
276.6  (0.51)4
285.3  (0.88)4
282.7  (0.45)4
292.5  (0.55)4
280.7  (0.64)4
279.8  (0.22)4
288.5  (0.37)4
286.4  (0.29)4
293.7  (0.76)3
288.5  (0.53)4
296.2 (1.44)4©
290.8  (1.51)3
287.1  (0.09)4
289.1  (0.83)3
286.3  (0.27)4
280.5  (0.60)3
278.7  (0.41)4
287.0  (0.51)4
279.7  (0.43)4
271.8  (0.89)4
271.0  (0.44)4 | 265.0
269.5
273.3
277.2
285.5
283.6
292.7
280.9
279.7
289.0
286.1
293.0
287.8
297.9
291.5
285.5
289.7
286.9
282.0
278.8
287.5
279.6
272.2
271.0 | -0.8
-0.4
-0.3
-0.6
-0.2
-0.9
-0.2
-0.2
0.1
-0.5
0.3
0.7
0.7
-1.7
-0.7
1.6
-0.6
-0.6
-1.5
-0.1
-0.5
0.1
-0.4
0.0 | -0.23
-0.14
-0.13
-0.10
-0.21
-0.10
0.40
-0.17
-0.12
0.22
0.30
0.43
0.07
0.98
0.42
-0.20
0.80
-0.08
0.23
-0.15
0.22
-0.06
-0.26
-0.36 | -0.14
0.07
-0.04
-0.10
0.16
-0.35
0.69
-0.22
-0.01
0.36
0.60
1.00
0.34
0.86
0.83
0.51
0.85
0.03
-0.08
-0.05
0.32
0.18
-0.26
-0.28 | 0.09
0.21
0.09
0.0
0.37
-0.25
0.29
-0.05
0.11
0.14
0.30
0.57
0.27
-0.12
0.41
0.71
0.05
0.11
-0.31
0.10
0.10
0.24
  0.0      0.08 |

- Note: (---)0 indicating that there is no valid LST in the four pixels surrounding the site; © near clouds.

**Radiance-based validation of V5 Terra  MODIS LSTs at the CarboEurope site in Hainich, Germany (51.0792°N, 10.4522°E). The atmospheric temperature and water vapour profiles used in the radiance-based validation were measured by radio sounding in Meiningen, Germany (50.56°N, 10.38°E), about 60km south of Hainich.  In column 6, ****(δT)n ****indicates the standard deviation of MODIS LSTs in n pixels with valid LSTs in the four pixels surrounding the site. The mean and RMS of errors (in column 8) are      -0.3K and 0.73K. ****This set of validation results clearly indicates that the LST cannot be validated by longwave radiation data measured at 44m above surface as Wang et al did for the same site in a recent paper published in Remote Sensing of Environment, v112, 623-635, 2008 ****(bias -2.2K and RMSE 2.5K)** .

## with radio

- sounding
- profiles

## with radio

- sounding
- profiles

## Radiance-based validation of V5 Terra (labelled as T in column 2) and Aqua (labelled as A in column 2) MODIS LSTs at the validation sites in Railroad Valley, Nevada (38.4817oN, 115.6905oW) in cases 1-10 and grassland, Texas (36.2992oN, 102.5706oW) in cases11-14. M*D07 in column 5 represents the atmospheric temperature and water vapour product MOD07 from Terra MODIS data or MYD07 from Aqua MODIS data. 

|  
case
no. | 
Granule ID   (T/A) | 
local solar date & time
(m/d h:m) | 
viewing
zenith
angle
(o) | 
radiosonde
cwv, Ts-air
(M*D07)
(cm, K) | 
MODIS
LST
(δT)
(K) | 
MOD
 L31
inverted
LST (K) | 
LSTmodis–
MOD L31
inverted
LST
(K) | 
profile-based
Tb31-
Tb32
(K) | 
 MODIS
Tb31-Tb32 (K) |  |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
|  |  |  |  |  |  |  |  |  |  |  |
| 1

2

3

4

5

6

7

8

9

10

11

12

13

14 | A2003177.1800 (T)
A2003178.1840 (T)
A2003179.0545 (T)
A2003180.1830 (T)
A2003181.0535 (T)
A2003177.1010 (A)
A2003178.2020 (A)
A2003179.0955 (A)
A2003180.2005 (A)
A2003181.0945 (A)
A2005111.1755 (T)
A2005112.0455 (T)
A2005111.1930 (A)
A2005112.0910 (A) |  6/26 10:19

 6/27 11:02

 6/27 22:06

 6/29 10:50

 6/29 21:54

 6/26 02:28

 6/27 12:37

 6/28 02:16

 6/29 12:25

 6/30 01:54

 4/21 11:05

 4/21 22:09

 4/21 12:40

 4/22 02:20 | 53.7

11.5

4.6

12.0

18.4

46.9

44.4

31.9

55.7

11.4

22.0

5.84

44.9

40.5 | 0.71,298.2
(0.56,299.5)
0.78,299.7
(0.93,304.2)
1.00,289.2
(1.04,287.4)
0.84,305.2
(1.35,305.2)
0.89,295.2
(0.55,290.4)
0.71,288.2
(0.46, 283.8)
0.78,302.7
 (0.41, 293.6)
0.90,284.2
 (0.76, 286.2)
0.84,305.2
 (0.84, 297.7)
0.74,295.2
 (0.40, 288.0)
0.35, 289.5
(0.46,289.4)
0.84, 282.2
(0.62,279.4)
0.58, 293.2
(0.50,288.4)
0.84, 280.2
(0.65,278.9) | 320.4
(1.4)
326.8
(1.7)
288.6
(0.6)
327.5
(1.4)
289.8
(0.6)
281.0
(0.7)
330.3
(1.5)
281.6
(0.4)
326.5
(1.8)
282.5
(0.5)
304.1
(0.4)
278.3
(0.2)
306.5
(0.1)
278.4
(0.1) | 321.5

327.6

288.7

328.5

289.8

281.1

331.1

282.0

327.3

282.4

304.6

278.5

307.1

278.6 | -1.1

-0.8

-0.1

-1.0

0.0

-0.1

-0.8

-0.4

-0.8

0.1

-0.5

-0.2

-0.6

-0.2 | 0.15

0.03

-0.48

0.08

-0.59

-0.78

0.36

-0.66

0.38

-0.69

-0.07

-0.26

0.19

-0.28 | 0.23

0.21

-0.47

0.11

-0.37

-0.70

0.55

-0.48

0.69

-0.53

0.01

-0.40

0.33

-0.27 |  |

## Emissivity spectra of the grassland in TX on 21 April 2005 (left)  and a bare soil site near Salton Sea, CA on 21 June 2006 (right) measured by the sun-shadow method with the Bomem TIR spectroradiometer and comparisons to band emissivities in the V5 MOD11B1 product 

**Emissivity spectra of the grassland in TX on 21 April 2005 (left)  and a bare soil site near Salton Sea, CA on 21 June 2006 (right) measured by the sun-shadow method with the Bomem TIR spectroradiometer and comparisons to band emissivities in the V5 MOD11B1 product**

## Radiance-based validation of V5 Terra and Aqua MODIS LSTs at the validation sites in Salton Sea (33.2oN, 115.75oW) in cases 1-6 and nearby bare soil land (33.25oN, 115.95oW) in cases 7-12. The δT value in column 6 is the standard deviation of MODIS LSTs in the four pixels surrounding the validation site.

| 
case
no. | 
Granule ID   (T/A) | local solar
date & time
(m/d h:m) | viewing
zenith
angle
(o) | radiosonde
cwv, Ts-air
(M*D07)
(cm, K) | MODIS
LST
(δT)
(K) | MOD L31 
inverted
LST
(K) | LSTmodis–
MOD L31
inverted
LST
(K) | profile-based Tb31-
Tb32
(K) | MODIS
 Tb31-Tb32
(K) |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 1

2

3

4

5

6

7

8

9

10

11

12 | A2006171.0545 (T)
A2006171.1000 (A)
A2006171.2100 (A)
A2006172.1830 (T)
A2006173.2050 (A)
A2006174.1820 (T)
A2006171.0545 (T)
A2006171.1000 (A)
A2006171.2100 (A)
A2006172.1830 (T)
A2006173.2050 (A)
A2006174.1820 (T) | 6/19 22:04

6/20 02:19

6/20 13:20

6/21 10:51

6/22 13:08

6/23 10:38

6/19 22:03

6/20 02:19

6/20 13:20

6/21 10:50

6/22 13:08

6/23 10:37 | 8.2

41.4

16.9

1.3

8.0

25.2

9.6

40.4

15.5

3.0

9.5

26.6 | 1.95,300.3
(1.67,303.1)
1.95,300.3
(1.47,301.5)
1.61,308.2
(1.13,303.4)
1.12,309.9
(0.79,302.9)
1.63,312.2
(0.83,303.6)
2.56,308.2
(3.65,304.8)
1.95,305.2
(1.33,304.4)
1.95,300.3
(1.26,300.3)
1.61,315.2
(1.00,297.6)
1.12,309.9
(0.60,309.0)
1.63,312.2
(1.05,298.5)
2.56,308.2
(3.11,313.1) | 300.4
(0.04)
300.3
(0.07)
304.5
(0.1)
303.9
(0.06)
305.4
(0.13)
303.6
(0.23)
302.6
(0.15)
299.1
(0.29)
332.9
(0.25)
327.8
(0.25)
332.7
(0.2)
328.2
(0.15) | 300.4

299.9

304.6

304.1

305.8

303.2

303.6

299.5

335.2

330.1

334.5

331.6 | 0.0

0.4

-0.1

-0.2

-0.4

-0.4

-1.0

-0.4

-2.3

-2.3

-1.8

-3.4 | 0.37

0.40

0.45

0.32

0.55

1.14

-0.28

-0.39

0.92

0.15

1.02

2.13 | 0.31

0.45

0.49

0.33

0.46

1.39

-0.55

-0.51

0.81

-0.06

1.14

2.08 |

## Error Analysis of the LSTs Retrieved by the Split-window Method

- The LST errors versus LST values shown in the last two tables in two groups, one for soil sites in Railroad Valley and near Salton Sea, and another for the grassland and lake sites, in day and night cases, respectively (A). The LST errors versus viewing zenith angle for the cases in the soil group (B), and the LST errors versus column water vapor from measured atmospheric profiles for the cases in the soil day group (C).

- (A)

- (B)

- (C)

**The reason for larger errors in high LST cases with cwv > 1.5cm is that the range of Ts-air ****±**** 16K for the LST values used in the development of the split-window algorithm is not wide enough for the daytime bare soil cases.**

## with radio

- sounding
- profiles

## with NCEP modeled profiles

**with NCEP modeled profiles**

## with radio

- sounding
- profiles

## How does the V5 split-window algorithm perform in deserts?

**How does the V5 split-window algorithm perform in deserts?**

- There is no significant difference in the LSTs at In-salah in the V4 and V5 MOD11A1 products (because the same split-window algorithm has been used in V4 and V5).

- The sensitivity of the split-window algorithm to uncertainties of surface emissivities in bands 31 and 32: if ε**31** reduces by 0.008 and ε**32** increases by 0.008 the retrieved LST value would increase by 2.7K in case of A2007198.1025.

- R-based validation results of LSTs in the V5 MOD11_L2 product at In-salah, Algeria
- (27.22oN, 2.5oE) using radio sounding profiles and NCEP modeled profiles.

- δT = 2.8 ± 0.4K
- when radio sounding profiles
- are used, or

- δT = 3.1 ± 0.5K
- when NCEP modeled profiles
- are used.

- The mean and standard
- deviation of LST errors:

## Reflectance spectra measured from sand and soil samples by the MODIS LST group. There are significant variations in the ranges of

**Reflectance spectra measured from sand and soil samples by the MODIS LST group. There are significant variations in the ranges of**
**3-5 and 8-10****μ****m (left) but the change in the reflectance difference in bands 31 and 32 (at 11 and 12****μ****m shown in right) may cause large errors in LSTs retrieved by the split-window algorithm**

## Comparisons of LST values at In-salah (27.22°N, 2.5°E) retrieved by the day/night and split-window algorithms in the V4 and V5 MOD11B1 products in 2004. In V4, the LSTs retrieved by the day/night method is 1-2K larger at

**Comparisons of LST values at In-salah ****(27.22°N, 2.5°E)** **retrieved by the day/night and split-window algorithms** **in the V4 and V5 MOD11B1 products in 2004. In V4, the LSTs retrieved by the day/night method is 1-2K larger at**
**night or 2-3K larger and dependent on LST values in daytime.**

**However, the LSTs retrieved by the day/night method are only slightly larger in V5 due to the tight bounding with the split-window method as shown in right.**

**But the V4 day/night  algorithm also overestimated LSTs in most regions including lakes and dense vegetation areas. This is not good.**

## Comparisons of retrieved surface emissivity values in bands 29 and 31 at In-salah (27.22°N, 2.5°E) in the V4 and V5 MOD11B1 products in 2004. The

**Comparisons of retrieved surface emissivity values in bands 29 and 31 at In-salah ****(27.22°N, 2.5°E)** **in the V4 and V5 MOD11B1 products in 2004. The**
**low fluctuated band-31 values in V4 is bad. The V5 made improvement in**
**band 31 but failed in band 29 due to the large emis29 values above 0.9 in**
**about half cases.**

## What is wrong with the performance of the V5 day/night algorithm in desert regions? 

| product | C5 (V5) | C4 (V4) |
| --- | --- | --- |
| MOD11A1   day LST | 297 | 189 |
| MOD11A1 night LST | 283 | 262 |
| MOD11B1   day LST | 301 | 195 |
| MOD11B1 night LST | 287 | 274 |

**Numbers of days with valid LST values at In-salah in C5 and C4 MOD11A1 and MOD11B1 products in 2004.**

**2. The definition of clear-sky LSTs in V5 by MOD35 at a confidence of >= 95%  clear over  land  <= 2000m increased the number of daytime LSTs at In-salah by about 50% both in MOD11A1 and MOD11B1 mainly in spring and winter. This C5 change  gives good seasonality in LSTs, as shown in the plot below**
**for day LSTs at In-salah in 2005.**

**1. ****Because the day/night algorithm is tightly bounded with the split-window algorithm in V5, the large errors in LSTs retrieved by the split-window method in desert regions by any reasons also affect the day/night algorithm. So we need to tune the tightness.**** **

## Summary of C5 LST Products

- 1.

- Daytime and nighttime LSTs in the C5 level-2 Terra & Aqua MODIS
-      LST products (M*D11_L2)retrieved by the split-window algorithm have been validated by the temperature-based approach at large homogeneous sites in lakes, rice field and dense vegetation areas. They have been also validated by the radiance-based approach at various sites in arid regions and cold regions including the South Pole. The LST errors are within ±1K in most cases.
-      The LSTs at 6km grids retrieved by the V5 day/night algorithm may be validated indirectly by comparisons to the LSTs from the split-window algorithm aggregated at the 6km grids.

- 2. The land surface emissivities retrieved by the V5 day/night algorithm
-     have been validated only at a few sites. The large fluctuations in the
-     retrieved emissivities in the desert regions may be due to cloud
-     contaminations and large errors in retrieved LSTs.

- 3. Cloud contaminations and the effects of aerosols above average
-     loadings are the major source of errors in the MODIS LST products.
-     The affected LSTs in the C5 level-2 MODIS LST products have not
-     been removed although they were removed from the level-3 products.

##  Proposed Processing Plan for C6 LST Products

- (1), remove cloud-contaminated LSTs not only from level-3 LST products but also from level-2 LST products (MOD11_L2 and MYD11_L2).
- (2), update the coefficient LUT (lst_coef.h) for the split-window algorithm with comprehensive regression analysis of MODIS simulation data in bands 31 and 32 over wide ranges of surface and atmospheric conditions, especially extending the upper boundary for (LST – Ts-air) in arid and semi-arid regions and increasing the overlapping between various sub-ranges in order to reduce the sensitivity of the algorithm to the uncertainties in the input data (i.e., column water vapor and air surface temperature from MOD07 and MYD07).
- (3), make minor adjustments in the classification-based surface emissivity values (band_emis.h), especially for land-cover type of bare soil and rocks.
- (4), tune the day/night algorithm by adjusting weights to improve its performance in desert regions where the incorporated split-window algorithm may not work well.
- (5), generate new LST products for 8-day and monthly at 6km grids (in response to user community requests).

##  Research Issues for C6 LST Products

**(2), develop methods to analyze and correct the effects of thin cirrus clouds and aerosols above the average loading (the MODIS aerosol product was improved by the deep blue algorithm recently), TBD.**

**(1), add flexibility to input options in the daily PGE16 if the need for near real-time processing capability can be better met without scarifying the LST quality, TBD.**

## 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. *Int. 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., (2008). New refinements and validation of the MODIS land-surface temperature/emissivity products. *Remote Sensing of Environment*, 112, 59-74.

- Wan, Z., & Li, Z.-L., (2008).  Radiance-based validation for the V5 MODIS land-surface temperature product. *International Journal of Remote Sensing*, to be published in July.

- Wang, K., Wan, Z., Wang, P., Sparrow M., Liu, J., & Haginoya, S., (2007).  Evaluation and improvement of the MODIS land surface temperature/emissivity products using ground-based measurements at a semi-desert site on the western Tibetan Plateau. *Int. Journal of Remote Sensing*, 28, 2549-2565

- Li, Z.-L., Becker, F., Stoll, M.P., & Wan, Z (1999). Evaluation of six methods for extracting relative emissivity spectra from thermal infrared images. *Remote Sensing of Environment*, 69, 194-214.