landsat 5 thematic mapper reflectance and ndvi 27-year time
TRANSCRIPT
![Page 1: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/1.jpg)
Landsat 5 reflectance and NDVI
27-year time series inconsistencies due to satellite orbit change
Hankui K. Zhang and David P. Roy
Geospatial Sciences Center of Excellence, South Dakota State University,
Brookings, SD 57007, USA
July 28 2016
![Page 2: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/2.jpg)
Lots of Remote Sensing Systems launched since 1972
Landsat-1
Landsat-5
Belward and Skøien, 2014
![Page 3: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/3.jpg)
Landsat 5 local overpass time
• 9:45 a.m Mean Local Time of Descending Node (MLTDN)
- initially orbit maintained by periodic station keeping maneuvers to maintain
ground track and orbit phase with Landsat 4
- initially MLTDN required to not vary by more than ±15 minutes of 9:45 a.m
![Page 4: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/4.jpg)
Landsat 5 local overpass time example
Local time
difference
~1 hour
over 27
years …
(path 27 row 42 somewhere in Texas)
• 9:45 a.m Mean Local Time of Descending Node (MLTDN)
- initially orbit maintained by periodic station keeping maneuvers to maintain
ground track and orbit phase with Landsat 4
- initially MLTDN required to not vary by more than ±15 minutes of 9:45 a.m
![Page 5: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/5.jpg)
Deering, D. W., Eck, T. F., & Banerjee, B. (1999). Characterization of the reflectance anisotropy of three boreal forest canopies in spring – summer. Remote Sensing of Environment, 67, 205 – 229.
Nadir view Reflectance varies with solar zenith
SZ=65°, more shadow, lower ρred SZ=40°, less shadow, higher ρred
![Page 6: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/6.jpg)
NDVI vs Solar Zenith (Pinter 1993, RSE)
![Page 7: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/7.jpg)
Can we see Landsat 5 orbit drift effects in TM images ?
Considers • 3 sites (Vogelmann et al.
2016 RSE), Crater Lake National Park, Oregon
At each site:
• 9 adjacent 30m pixels
• 12 L1T Landsat 5 cloud-free images spanning 1991-2010
• Atmospherically corrected
• summer anniversary date images (no more than ±5 days of Sep. 1)
![Page 8: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/8.jpg)
Blue: sparsely
vegetated pumice
desert
Green: coniferous
forest with gradual
NDVI increase
Orange: coniferous
forest with gradual
NDVI decrease
Open circles
- 9 individual pixel values
Solid colored lines
- median of 9 pixel
values
3 SITE TIME SERIES
![Page 9: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/9.jpg)
Sola
r Z
enit
h
35°
40°
45°
50°
Sola
r Z
enit
h
35°
4
0°
45°
50°
3 SITE TIME SERIES
![Page 10: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/10.jpg)
Landsat 5 TM local overpass times
for 2 sites same path but different rows
Minnesota path 27 row 26 Texas path 27 row 42
2011 reference overpass time (red dots)
![Page 11: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/11.jpg)
𝑡 𝑙𝑜𝑐𝑎𝑙 = 1.36292 × 10−9𝛼5 − 3.15403 × 10−8𝛼4 − 3.15819614 × 10−6𝛼3 +
0.0000652685643𝛼2 + 0.0120604786763𝛼 + 10.06
2011 reference overpass time and solar zenith based on 12 months of non-Antarctic Landsat metadata values
22.14⁰ <= θs <= 89.71⁰ θs mean = 43.23⁰
Zhang, H.K., Roy, D.P., & Kovalskyy, V. (2016). Optimal solar geometry definition for global long-term Landsat time-series bidirectional reflectance normalization. IEEE TGRS, 54 (3), 1410-1418.
where 𝛼 is latitude
q2011 reference = f (tlocal , date, latitude)
![Page 12: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/12.jpg)
Minnesota θobserved – θ2011 reference
Texas θobserved – θ2011 reference
Green: summer difference values
Max = 5.9°
Max = 11.2°
![Page 13: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/13.jpg)
Modeling the reflectance impacts
of overpass time and so solar zenith change over 27 years
forparametersBRDFMODISfNBAR vreferencesreference , , 0 2011 ,
forparametersBRDFMODISfNBAR vobssobs , , 0,
![Page 14: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/14.jpg)
consider
different mean CONUS land cover MODIS BRDF parameters
Roy, D.P., Zhang, H. K., Ju, J., Gomez-Dans, J. L., Lewis, P.E., Schaaf C.B., Sun, Q., Li, J., Huang, H., & Kovalskyy, V., 2016b, A general method to normalize Landsat reflectance data to nadir BRDF adjusted reflectance, Remote Sensing of Environment. 176, 255-271.
Solar Zenith
Nadir view NDVI
![Page 15: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/15.jpg)
NB
AR
red
dif
fere
nce
-0.0
08
-0
.00
4
0
.00
0
NB
AR
NIR
dif
fere
nce
-0
.01
5 -
0.0
10
-
0.0
05
0
.00
0
NB
AR
ND
VI
dif
fere
nce
-0
.01
0
.00
0
.01
0
.02
Minnesota path/row modeled NBAR
difference Red
NIR
NDVI
NBAR obs, l -NBAR 2011 reference, l
closed shrubland class
![Page 16: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/16.jpg)
NB
AR
red
dif
fere
nce
-0.0
08
-0
.00
4
0
.00
0
NB
AR
NIR
dif
fere
nce
-0
.01
5 -
0.0
10
-
0.0
05
0
.00
0
NB
AR
ND
VI
dif
fere
nce
-0
.01
0
.00
0
.01
0
.02
Red
NIR
NDVI Green:
modeled NBAR difference
summer images only
Texas path/row modeled NBAR
difference
NBAR obs, l -NBAR 2011 reference, l
closed shrubland class
![Page 17: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/17.jpg)
NB
AR
red
dif
fere
nce
-0.0
08
-0
.00
4
0
.00
0
NB
AR
NIR
dif
fere
nce
-0
.01
5 -
0.0
10
-
0.0
05
0
.00
0
NB
AR
ND
VI
dif
fere
nce
-0
.01
0
.00
0
.01
0
.02
Minnesota path/row modeled NBAR
difference Red
NIR
NDVI
NBAR obs, l -NBAR 2011 reference, l
closed shrubland class
![Page 18: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/18.jpg)
NB
AR
red
dif
fere
nce
-0.0
08
-0
.00
4
0
.00
0
NB
AR
NIR
dif
fere
nce
-0
.01
5 -
0.0
10
-
0.0
05
0
.00
0
NB
AR
ND
VI
dif
fere
nce
-0
.01
0
.00
0
.01
0
.02
Red
NIR
NDVI Green:
modeled NBAR difference
summer images only
Texas path/row modeled NBAR
difference
NBAR obs, l -NBAR 2011 reference, l
closed shrubland class
![Page 19: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/19.jpg)
NB
AR
red
dif
fere
nce
-0.0
08
-0
.00
4
0
.00
0
NB
AR
NIR
dif
fere
nce
-0
.01
5 -
0.0
10
-
0.0
05
0
.00
0
NB
AR
ND
VI
dif
fere
nce
-0
.01
0
.00
0
.01
0
.02
Red
NIR
NDVI
Texas path/row modeled NBAR
difference
NBAR obs, l -NBAR 2011 reference, l
closed shrubland class
![Page 20: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/20.jpg)
Conclusions • Landsat 5 orbit changed considerably over 27 year life
– Overpass time changed by up to ~1 hour (>> ±15 mins of 09:45 AM MLTDN)
– Solar zenith changed by >10°
• Can see orbit shift in actual Landsat 5 time series
• BRDF modeling findings
– 27-year NDVI change trend 0.0006 NDVI/year, small but not insignificant
– Comparing certain years, i.e., 1995 and 2007 may not be a good idea • NDVI1995 0.11 > NDVI 2007 for anisotropic land cover types
• NDVI 1995 0.05 > NDVI 2007 for average CONUS land cover types
• Further research to develop a Landsat BRDF normalization approach for Landsat 5
orbit drift changes is recommended
• Zhang, H.K and Roy, D.P., Landsat 5 Thematic Mapper reflectance and NDVI 27-year
time series inconsistencies due to satellite orbit change, RSE, In review
![Page 21: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/21.jpg)
Global WELD NEX Version 3.0 September 2009 30m product from 15,058 L1T scenes (7,328 Landsat 5 & 7,730 Landsat 7) Sinusoidal Equal Area Projection
WELD Landsat 5 & 7 observed Solar Zenith used to derive WELD NBAR
![Page 22: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/22.jpg)
WELD Landsat 5 & 7 modeled Solar Zenith used to derive WELD NBAR
Global WELD NEX Version 3.0 September 2009 30m product from 15,058 L1T scenes (7,328 Landsat 5 & 7,730 Landsat 7) Sinusoidal Equal Area Projection
![Page 23: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/23.jpg)
Minnesota (path/row
27/26) modeled NBAR
reflectance and NDVI using
the fixed 12 month mean
CONUS (blue) and closed
shrubland class (red)
spectral BRDF model
parameters and setting the
solar zenith to θobserved.
𝜌𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑𝑟𝑒𝑑
𝜌𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑𝑁𝐼𝑅
𝑁𝐷𝑉𝐼𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑
![Page 24: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/24.jpg)
Same as above but for
Texas (path/row 27/42)
𝜌𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑𝑟𝑒𝑑
𝜌𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑𝑁𝐼𝑅
𝑁𝐷𝑉𝐼𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑
![Page 25: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/25.jpg)
Satellite orbit can change significantly over its life time
![Page 26: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/26.jpg)
Black dots: θreference
plotted every 16 days
for the 2011 Landsat
acquisition dates.
Minnesota θobserved Texas θobserved
Minnesota θobserved - θreference Texas θobserved - θreference
![Page 27: Landsat 5 Thematic Mapper reflectance and NDVI 27-year time](https://reader031.vdocument.in/reader031/viewer/2022030323/589ad81b1a28aba5468ba949/html5/thumbnails/27.jpg)
Minnesota θobserved Texas θobserved
Minnesota θobserved - θreference Texas θobserved - θreference
Black dots: θreference
plotted every 16 days
for the 2011 Landsat
acquisition dates.