# Presentation: 181190

## Current Status of the MODIS LST Product and Plans for the Future

- (wan@icess.ucsb.edu)

- MODIS Land-Surface Temperature (LST) products have been generated from Terra MODIS data since late February 2000  and from Aqua MODIS data since June 2002 consistently in version 4 (also collection 4). Its accuracy has been validated with in-situ measurements in wide ranges of atmospheric and surface conditions (Wan et al., 2002 and 2004), as shown in Fig. 1.

**1**. **Accomplishments**

**2. Validation of the MODIS LST Product**

**Fig. 1**, Validation sites (from left): Lake Titicaca, Bolivia & Peru; Walker Lake, NV; rice field, CA; Bridgeport w/o and with snow, CA; soybean field, MS.

**2.1, conventional temperature-based method**

**2.2, advanced radiance-based method**

- TIR radiometers are used to measure the surface radiometric temperature.  Effects of surface emissivity and reflected atmospheric radiation are corrected to obtain the in-situ measured LST using the emissivity value based on land-cover or sample measurements and atmospheric radiative transfer simulations.  Comparisons between the MODIS LSTs and in-situ LSTs give the accuracy.

- Radiosonde balloons are launched to measure the atmospheric temperature and water vapor profiles around the MODIS overpass.  Based on the profiles, MODIS LST and  spectral emissivity measured in the field or from samples or estimated from land-cover, atmospheric radiative transfer simulations are made to calculate brightness temperature T**b** in MODIS band 31. The difference between the calculated T**b** and MODIS T**b**** **values divided by atmospheric transmission in band 31 gives the MODIS LST accuracy.

- This method is limited by the spatial variation in LSTs, especially during daytime, as shown in MAS data clearly.

- The main advantage of this method is that it works for both daytime and nighttime.

- Radiosonding
- 10:29 – 12:29 PDT 29 June 2003
- (17:29 – 19:29 UTC 29 June 2003)
- at 38o29.036 ́N, 115o41.260 ́W
- 1.45km above sea level

- column water vapor = 0.84 cm

**Fig. 2**, Atmospheric temperature (left) and water vapor (right) profiles measured  by radiosonding over Railroad Valley, NV on 6/29/2003.

**Fig. 3**,  Downwelling radiance T**b** at nadir measured by Bomem  spectro-radiomter at 10:30 PDT 6/29/03 in Railroad Valley and the T**b** calculated by MODTRAN4.0 based on measured atmospheric profiles.

- This method requires accurate atmospheric temperature and water vapor profiles, such as in Fig. 2 in the June 2003 Railroad Valley field campaign, and accurate atmospheric radiative transfer code.  The excellent agreements between the measured spectral sky radiance and the radiance calculated with MODTRAN4.0 (Berk et al., 1999) based on measured profiles provide a solid evidence of the good quality of the Bomem TIR spectroradiometer and MODTRAN4.0, as shown in Fig. 3.

- Four TIR radiometers were deployed in a 500m square in the central portion of the Railroad Valley playa.  A Bomem TIR spectroradiometer (MR100), as shown in Fig. 4,  was deployed at 38.4817oN 115.6905oW, by the center of the rectangular consisting of the four TIR radiometers.  Six sounding balloons that carry Vaisala radiosonde RS90-A were launched in clear
- sky days to measure the atmospheric temperature and water vapor profiles, listed in Table I.

**2.2.1, ground measurements in Railroad Valley, NV in June 2003**

**2.2.2, spectral emissivity of playa measured  with sun-shadow method**

- We measured surface-leaving radiance
- under sunshine and shadow conditions
- with the spectroradiometer MR100.  A
- simplified version of the day/night method
- (Wan and Li, 1997) was used to retrieve
- surface temperatures under sunshine and
- shadow conditions, and emissivity values
- in seven bands.  After determining playa
- Surface temperature, the measured
- spectral data were used to calculate the
- spectral emissivity of the playa surface,
- as shown in Fig.5.

**Fig. 4**, Bomem spectroradiometer deployed in RailroadValley, NV.

**Fig. 5**, Spectral emissivity of playa measured by the  sun-shadow method in Railroad Valley, NV.

**2.2.3, Validation Results**

- The radiance-based method was  used to validate both daytime and nighttime LSTs in Terra cases and Aqua cases, in Tables II and III, respectively.  Comparing the errors in the same night granules with two methods indicates that the errors are smaller when the radiance-based method is used.

**Table I**,  Atmospheric temperature and water vapor profiles measured
-  in Railroad Valley under clear-sky conditions in June 26-29, 2003.

| case
no. | date
(m/d/y) | start time PDT
(hh:mm) | duration
(minutes) | height reached
(km) | atmos.
cwv cm) |  Ts-air 
(K) |
| --- | --- | --- | --- | --- | --- | --- |
| A
B
C
D
E
F | 6/26/03
6/27/03
6/27/03
6/28/03
6/29/03
6/29/03 | 10:22
10:14
22:58
02:13
10:29
22:05 | 112
103
102
120
120
120 | 16.1
19.1
23.2
23.3
24.4
24.3 | 0.71
0.78
1.00
0.90
0.84
0.67 | 298.2
299.7
291.8
284.2
305.2
291.6 |

**Table II**,  Radiance-based validation of the V4 Terra MODIS MOD11A1  LSTs.

| case
no. | granule
ID | date & time PDT
(m/d hh:mm) | viewing
zenith
angle (o) | MOD07
cwv, Ts-air
(cm, K) | atmos.
profile
no. (MS) | MODIS
LST
(K) | ΔTb
in b29
(K) |  trans
in b31 |  ΔTb
in b31
(K) |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 1
2
3
4
5
6 | A2003177.1800
A2003178.1840
A2003179.0545
A2003180.0630
A2003180.1830
A2003181.0535 | 6/26 11:02
6/27 11:45
6/27 22:49
6/28 23:32
6/29 10:33
6/29 22:37 | 53.7
11.5
4.6
60.2
12.0
18.4 | 0.67, 299.3
1.90, 305.2
1.30, 291.8
1.30, 293.4
0.96, 301.4
0.67, 291.6 | A (8)
B (8)
C (9)
C (6)
E (9)
F (9) | 320.18
326.78
288.43
287.82
327.62
289.86 | -0.85
-0.18
+0.14
-0.84
-0.58
+0.44 | 0.88
0.92
0.91
0.84
0.92
0.93 | -1.10
-0.59
-0.21
-0.95
-0.81
+0.01 |

**Table III**,  Radiance-based validation of the V4 Aqua MODIS MYD11A1  LSTs.

| case
no. | granule
ID | date & time PDT
(m/d hh:mm) | viewing
zenith
angle (o) | MYD07
cwv, Ts-air
(cm, K) | atmos.
profile no. (MS) | MODIS
LST
(K) | ΔTb
in b29
(K) |  trans
in b31 |  ΔTb
in b31
(K) |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 1
2
3
4
5
6
7 | A2003177.1010
A2003178.0915
A2003178:2020
A2003179.0955
A2003179.2100
A2003180.2005
A2003181.0945 | 6/26 03:11
6/27 02:16
6/27 13:20
6/28 02:59
6/28 14:03
6/29 13:08
6/30 02:47 | 46.9
40.9
44.4
31.9
26.9
55.7
11.4 | 0.50, 284.3
0.48, 281.8
0.54, 335.3
0.83, 283.2
0.72, 338.3
0.94, 330.4
0.41, 284.3 | A (7)
D (6)
B (8)
D (8)
E (6)
E (8)
F (8) | 280.90
281.09
331.00
280.16
331.19
326.47
282.31 | +0.31
+0.15
-0.69
-0.56
-0.63
-0.39
+0.13 | 0.90
0.90
0.90
0.91
0.92
0.87
0.87 | +0.10
+0.20
-0.69
-0.98
-0.70
-0.75
-0.11 |

- Zhengming Wan, Yulin Zhang, and Qincheng Zhang, ICESS, University of California, Santa Barbara, CA

**2.3 Summary of the MODIS LST Validation**

- The MODIS LST products have been validated, most within 1K, with two methods in more than 40 clear-sky cases in wide ranges of surface and atmospheric conditions in 2000-2003, as shown in Fig. 6.

**Fig. 6**, Comparisons between V2 and V4 MODIS LSTs and in-situ values in field campaigns in 2000-03.

- Berk, A., G.P. Anderson, et al., “MODTRAN4 radiative transfer modeling for atmospheric correction”, SPIE, vol. 3756, 348-353, 1999.
- Wan, Z. and Z.-L. Li, “A physics-based algorithm for retrieving land-surface  emissivity and temperature from EOS/MODIS data”, IEEE Trans. Geosci.  Remote Sens., vol. 35, 980-996, 1997.
- Wan, Z., Y. Zhang, Q. Zhang, and Z.-L. Li, “Validation of the land-surface  temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data”, Remote Sens. Environ., vol. 83, 163-180, 2002.
- Wan, Z. Y. Zhang, Q. Zhang, and Z.-L. Li, “Quality assessment and  validation of the MODIS global land surface temperature”, Int. J. Remote Sens., vol. 25, 261-274, 2004.

- This work was supported by NASA’s EOS Program contract NAS5-31370. Thanks to the staffs at Airborne Sensor Facility, NASA Ames Research Center, and at the ER-2 Operations Office, NASA Dryden Flight Research Center, for providing the MAS data,
- and to many others for their supports in the various field campaigns.

**3. Removing Cloud-contaminated LSTs**

**4. Plans for the Future**

- In the MODIS LST PGE, LST is retrieved for pixels in clear-sky conditions at the 99% confidence defined by cloudmask MOD35 or MYD35 over land and in lakes in V3, 99% confidence over land or 66% and higher in lakes in V4.  It is found that cloud-contaminated LSTs exist in both V3 and V4 LST products, e.g., LSTs below freezing point in Lake Tahoe that is open around the year. Cloud-screen schemes were developed with landcover-dependent constraints on temporal-spatial variations in LSTs, as shown in Fig. 7.

**Fig. 7, ** LSTs near the center of Lake Tahoe in the year 2002 MOD11A1 data sets before (left) and after (right) applying the cloud-screen scheme.

**4.1 Validation**

**4.2 Refinements of the MOD11 PGEs in Collection 5 (V5) to Produce CDR Quality LST Products**

- Field campaigns will be conducted in large grasslands in TX and NM to validate the LST and surface emissivity in the MODIS LST products.  Ground truths and station data will be used to evaluate the cloud-screen scheme at various sites of different land-cover types.

- 1) Add an option for fully combined use of Terra and Aqua data in the day/night LST algorithm.

- 2) Change the grid size in MOD11B1 to 6x6 1km grids from 5x5 so the same latitude interval of 0.05o be used in the         MOD11B1 and MOD11C products.

- 3) Separate the whole viewing angle region into six sub-ranges instead of five in V4.

- 4) Use average slope in the east-west direction at 6km grids. Azimuth-dependent MODIS data will be stored for grids with slope larger than 6o, day/night pairs of observations in the same azimuth direction will be used in the day/night algorithm.

- 5) Update the LUT of emissivities in bands 31 and 32 for the split-window algorithm with the values retrieved by the day/night algorithm from Terra MODIS data in 2000-2003.

- 6) Remove cloud-contaminated LSTs in MOD11A1 and MOD11B1.