Cross-Comparison of Vegetation Indices Derived from Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) Sensors

Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) are currently operational for routine Earth observation. There are substantial differences between instruments onboard both satellites. The enhancements achieved with Landsat-8 refer to the scanning technology (replacing of whisk-broom scanners with two separate push-broom OLI and TIRS scanners), an extended number of spectral bands (two additional bands provided) and narrower bandwidths. Therefore, cross-comparative analysis is very necessary for the combined use of multi-decadal Landsat imagery. In this study, 3,311 independent sample points of four major land cover types (primary forest, unplanted cropland, swidden cultivation and water body) were used to compare the spectral bands of ETM+ and OLI. Eight sample plots with different land cover types were manually selected for comparison with the Normalized Difference Vegetation Index (NDVI), the Modified Normalized Difference Water Index (MNDWI), the Land Surface Water Index (LSWI) and the Normalized Burn Ratio (NBR). These indices were calculated with six pairs of ETM+ and OLI cloud-free images, which were acquired over the border area of Myanmar, Laos and Thailand just two days apart, when Landsat-8 achieved operational obit. Comparative results showed that: (1) the average surface reflectance of each band differed slightly, but with a high degree of similarities between both sensors. In comparison with ETM+, the OLI had higher values for the near-infrared band for vegetative land cover types, but lower values for non-vegetative types. The new sensor had lower values for the shortwave infrared (2.11–2.29 μm) band for OPEN ACCESS Remote Sens. 2014, 6 311 all land cover types. In addition, it also basically had higher values for the shortwave infrared (1.57–1.65 μm) band for non-water land cover types. (2) The subtle differences of vegetation indices derived from both sensors and their high linear correlation coefficient (R > 0.96) demonstrated that ETM+ and OLI imagery can be used as complementary data. (3) LSWI and NBR performed better than NDVI and MNDWI for cross-comparison analysis of satellite sensors, due to the spectral band difference effects.


Introduction
Facing an ever-increasing number of observing satellite systems/sensors for the Earth, the joint application of remotely sensed data acquired from various sensors is very important and effective for monitoring global environment changes [1][2][3]. In such a case, the inter-satellite cross-comparison among multiple sensors becomes indispensable for the complementary use of imagery [4,5]. Dinguirard and Slater categorized the approaches of sensor radiometric calibration into three groups, namely preflight, onboard and vicarious calibration [6]. This study focused on the method of vicarious calibration, which refers to comparative analysis with multiple scenes collected during the flight, as widely discussed in other research [5].
By contrast, spectral indices of vegetation, such as the Normalized Difference Vegetation Index (NDVI) [1,7,12,14], the Normalized Burn Ratio (NBR) [10,17] and the Enhanced Vegetation Index (EVI) [17] are widely used as numeric indicators for cross-comparative analysis [3]. Among the broad range of indices, NDVI is the most commonly used for cross-comparative analysis [1,2]. The advantages of using vegetation indices as proxy variables for inter-calibration among existing sensors are the low sensitivity to the uncertainties in atmospheric correction and the variation in the satellite viewing angle [2]. However, visible spectral channels, like the blue, green and red bands, generally lead to apparent spectral band difference effects (SBDEs) on vegetation indices constructed with visible bands [3]. Teillet and Ren pointed out that the SBDE is mainly caused by the interaction between the spectral response profiles of sensors and the spectral dependence of atmospheric gas transmittance [3]. The SBDE of some vegetation indices calculated with the green and blue bands, e.g., the Atmospherically Resistant Vegetation Index and the Modified Triangular Vegetation Index, are especially obvious [3]; Similarly, NDVI, which is derived from the red and near-infrared bands, also results in some SBDE [1]. However, the SBDE on other vegetation indices, such as the Modified Normalized Difference Water Index (MNDWI), the Land Surface Water Index (LSWI) and NBR, was not reported yet.
Currently, two operational Landsat family satellites are available for consistently observing changes on the Earth's surface: Landsat-7 and Landsat-8. Landsat-7 Enhanced Thematic Mapper Plus (ETM+), launched on 15 April 1999, has exceeded its mission lifetime by nearly ten years [18]. However, the Scan Line Corrector (SLC) of the ETM+ sensor has been undergoing a failure since 31 May 2003. Despite the fact that each ETM+ SLC-off scene misses about 22% of the data, these images are some of the most important remotely sensed data in the World [19]. Landsat-8, a significantly improved and upgraded satellite, previously officially called the Landsat Data Continuity Mission (LDCM), launched on 11 February 2013, consists of two distinct sensors: the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS). Landsat-8 will be of great importance in continuing the multi-decadal Landsat data record [18].
As the latest member of the Landsat family of remote sensing satellites, Landsat-8 provides great improvement in the following aspects. First, the Landsat-8 OLI and TIRS instruments use a push-broom design, while Landsat-7 ETM+ and its predecessors all utilized a whisk-broom sensor [18]. Compared with the whisk-broom scanner, the push-broom design can receive stronger signals and improve signal-to-noise performance, due to the usage of long and linear arrays of detectors [18]. However, it also faces some difficulty in co-registering the data of multiple spectral bands, because of the varying sensitivities of the detectors [18]. Second, the Landsat-8 OLI sensor not only keeps the traditional settings of band designations in geometric resolution (30 m), scene size (170 km × 183 km) and revisit cycle (16-day), but also adds to the heritage bands (see Table 1) with new bands, such as coastal and aerosol studies (Band 1: ultra-blue), cirrus cloud detection (Band 9) and a quality assessment band, and adjusts the wavelength of each band [20]. It is worthwhile to note that the spectral range of Landsat-8 OLI is narrower than its predecessor, Landsat-7 ETM+ (Figure 1). The refinement of OLI bands can avoid the atmospheric absorption feature [18]. Table 1 also shows the major adjustments for TIRS in the band number, bandwidth and spatial resolution. According to the US Geological Survey (USGS), both sensors provide geometrically and radiometrically accurate images for repeated monitoring changes on the Earth's surface from 16 to eight days [19,21]. To enhance the utility of the Landsat imagery, it therefore is necessary to cross-compare Landsat-8 and Landsat-7. Relative spectral response (RSR) profiles showing the spectral band difference between Landsat-8 Operational Land Imager (OLI) (solid curve) and Landsat-7 Enhanced Thematic Mapper Plus (ETM+) (short dot curve). Note that these response metadata are available from the web page of the Spectral Characteristics Viewer (http://landsat.usgs.gov/tools_spectralViewer.php) provided by the US Geological Survey (USGS). SWIR, short-wave infrared.
According to the evident differences of both sensors, the pressing questions to be answered are whether and how the substantial changes, especially in the width of visible, near-infrared and shortwave-infrared spectral bands and the entire detector system (shifting from whisk-broom to push-broom design) [18], will affect the ability of consistently detecting and monitoring land surface changes. Taking advantage of the opportunity for providing near-coincident data from both sensors, cross-comparative analyses were conducted with various land cover types in Mainland Southeast Asia in 2013. The objective of this case study was threefold: (1) to compare the differences of spectral bands of Landsat-7 ETM+ and Landsat-8 OLI with randomly selected sample points of four major land cover types; (2) to compare the differences and correlation of four vegetation indices calculated with eight sample plots manually selected within six pairs of ETM+ and OLI images just two days apart; and (3) to compare the performance of different vegetation indices (LSWI, NBR, NDVI and MNDWI) in cross-comparison analysis between ETM+ and OLI. We believe that such an assessment will contribute to better monitoring and mapping of land surface changes by combining images with the eight-day repeat coverage of two sensors.

Ground Reference Data and Testing Sites
Field surveys of land use and land cover types were conducted in the Mekong River Basin, including Laos, Cambodia, Thailand and Vietnam, from 25 February 2013 to 17 March 2013. Using a Casio Exilim EX-H20G (with Hybrid GPS) camera and a Canon GP-E2 Camera GPS receiver, over 10,000 geo-referenced field photos along the field trip route in the four countries were gathered for selecting eight sample plots and sample points across the border area of Myanmar, Laos and Thailand. The thousands of photos were simultaneously taken in a bi-directional manner by a moving vehicle every 2-5 min. The field route traveled by the field team was collected with a handheld GPS receiver (Trimble Juno SB) in a counter-clockwise direction (Figure 2), from Laos, Cambodia to Vietnam and then back to Cambodia, Thailand and Laos. These field photos recorded specific information about land use and land cover types, such as primary or secondary forests, tree plantations (e.g., rubber, eucalyptus, teak, mangrove and fruit orchards), shrubland, grassland, swidden agriculture, croplands (e.g., paddy rice and cassava), village and towns, and so on.
The testing sites chosen for cross-comparison studies between ETM+ and OLI sensors are very important. Three basic rules for selecting sample plots were followed in this study: first, the time difference between both sensors observing the same place should be as short as possible. Usually, the side-by-side overlap between adjacent paths of Landsat-7 and Landsat-8 enables the area to be covered just one day apart. However, the gaps of ETM+ SLC-off images are noticeable along the edge of each scene and gradually reduce towards the scene center. For this reason, it is not appropriate to use the overlapping part of two scenes between adjacent paths. Second, the vegetation target types of sample plots should not encompass pronounced phenological features. For example, deciduous rubber trees start new leaf emergence, and many primary or secondary forests are slashed and burned for swidden cultivation in late March and April. This creates a rapid change in vegetation indices and is not suitable. Third, all the sample plots should be cloud free. Therefore, the border area of Myanmar, Laos and Thailand in Mainland Southeast Asia was chosen as the study area ( Figure 2), because both sensors acquired cloud-free images during the period when Landsat-8 achieved operational orbit. Both sensors almost provide near-coincident data during this period. This region has a tropical monsoon climate. There are three distinct seasons in Mainland Southeast Asia, i.e., the hot season from March to April, the cool and dry season from November to February and the rainy season from May to October [22]. During the dry season, natural landscape types are diverse and the differences among various landscapes are very evident. Land cover types include nine broad categories, such as primary forest, secondary forest (e.g., rubber/teak/eucalyptus plantations), shrubland, grassland, croplands, swidden cultivation (or slashed and burned plots), bare-land, built-up land and water bodies [23,24]. It is worthwhile to mention that swidden cultivation is a dominant land use type in Southeast Asia [25,26]. Swidden cultivation employs a natural or improved fallow phase, which is longer that the cultivation phase of annual crops, sufficiently long to be dominated by woody vegetation, and cleared by means of fire [27]. Additionally, there is infrequent or light cloudy weather during the hot and cool and dry seasons. The weather makes the acquisitions of cloud-free imagery possible for comparative analysis. For example, statistics from the USGS GloVis showed that there are 28 ETM+ or OLI scenes (Worldwide Reference System (WRS)-2 path/row = 130/046) with cloud coverage less than or equal to 10% during 2009-2013. Among them, 25 scenes were acquired from November to April, while only 3 scenes were available between May and October.
In this study, eight sample plots (polygons) covering diverse land cover types were manually selected within six pairs of scenes, which are only two days apart ( Figure 2, Table 2). These polygons were longitudinally located and spanned over 500 km from north to south and about 150 km from west to east. Sample plot 1 is located in the border region of Laos and Thailand. Sample plot 2, 3, 4, 6 and 7 were all located in Northern Thailand, and sample plot 5 and 8 were located in Myanmar. Similarly, the selection of sampling areas abided by three rules: cloud free, located in the overlapping regions between the center part of the ETM+ SLC image and OLI image and diverse land cover types. According to the Global Land Cover Map 2009 products (also called GlobCover 2009) provided by the European Space Agency (ESA), five major land cover types of all sample plots were statistically evaluated ( Figure 3). This shows that all the sample plots were normally closed to open shrubland, unplanted croplands and mosaic croplands/vegetation. Considering the coarse pixel resolution (300 m) of GlobCover 2009, these sample plots actually encompassed even more diverse land cover types. These sampling areas represent the highly different land cover types.   Figure 4, right panel). The sample points of these land cover types were used to compare each spectral band of ETM+ and OLI. All these sample points were not affected by the SLC failure and cloud.

Landsat-7/8 Imagery and Pre-Processing
In order to enable cross-comparative analysis between Landsat-7 ETM+ and Landsat-8 OLI sensors, six pairs of scenes (12 in total) of the Landsat-7 ETM+ SLC-off data and Landsat-8 OLI/TIRS Pre Worldwide Reference System (WRS) 2 data (Table 3) were gathered. Landsat-7 ETM+ SLC-off data refers to the products that have data gaps collected after May 2003 [19]. Landsat-8 OLI/TIRS Pre-WRS-2 data refers to the products acquired before 10 April 2013 [21]. All these data meet the quality standards of the same radiometric and geometric correction, but the geographic extents of each scene differ [19,21]. However, there is an overlap part of over 60 percent of a scene between adjacent paths of two sensors. This overlap area can be covered by scenes from each sensor, which are only two days apart. Considering the SLC failure problem of Landsat-7 ETM+, this research only compares the overlapping areas between the middle of ETM+ scene and OLI images. All twelve Landsat images were available from the Global Visualization Viewer (GloVis, http://glovis.usgs.gov/) developed by the Earth Resource Observation and Science Center (EROS) of the USGS. They are standard level-one terrain-corrected (L1T) products, which have had systematic radiometric and geometric corrections and ensured the overall geometric accuracy with ground control points and a digital elevation model [28].
As earlier studies show, surface reflectance data serves as the prerequisite for comparative analysis of vegetation indices derived from different sensor systems [2]. In this study, each Landsat image has undergone atmospheric correction to convert digital number values into surface reflectance for vegetation indices calculation in the software of Environment for Visualizing Images (ENVI, Version 5.0 SP3, Exelis Visual Information Solutions, Boulder, CO, USA). ENVI has an atmospheric correction module for retrieving spectral reflectance, i.e., the Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) [29]. FLAASH incorporates the Moderate Resolution Transmission (MODTRAN) 4 radiative transfer code [29]. Model parameters describing the tropical atmosphere and rural aerosols together with 2-band (K-T) aerosol retrieval were utilized in ENVI FLAASH. Although the sample areas were comprised of multiple land cover types (Table 2) in the tropical humid region, they could be divided into three groups: natural vegetation, exposed soil (i.e., unplanted cropland and swidden field) and water body. The specific target types can be detected by spectral bands. For instance, the green band works for the identification of water bodies when combined with NIR or short-wave infrared (SWIR) 1 to form the Normalized Difference Water Index (NDWI) [30] and MNDWI [31]. The red and NIR bands are particularly useful for vegetation monitoring [32]. The SWIR 1/2 bands are sensitive to the moisture content of vegetation and soil [33]. Landsat family sensors (since Landsat-5) have all these bands. Therefore, another three commonly used indices, the MNDWI [31], LSWI [34] and NBR [35], as well as NDVI [32] were introduced to compare the difference and correlation of both sensors. The corresponding formulas are given as follows: (1) where ρ green , ρ red , ρ NIR , ρ SWIR1 and ρ SWIR2 refer to the surface reflectance values of Band 2, Band 3, Band 4, Band 5 and Band 7 in the Landsat-7 ETM+ sensor and Band 3, Band 4, Band 5, Band 6 and Band 7 in the Landsat-8 OLI sensor (Table 1), respectively. As opposed to LSWI and MNDWI, NDVI and NBR are more widely used for cross-comparison analysis. LSWI is similar to NBR, as they were both computed with NIR and SWIR bands. It is appropriate to use LSWI for such an analysis. In contrast, MNDWI is less reported and may lead to SBDE.

Statistical Analysis
Statistical analysis methods are widely used in cross-comparison between various satellite sensors [1,2,4,10,15]. Statistical analysis mainly consists of three parts in this study. To begin with, based on the 3,311 sample points, the difference of each band of both sensors was cross-compared.

Comparison Landsat-7 ETM+ and Landsat-8 OLI Spectral Bands
Considering the chronological order difference of Landsat-7 and Landsat-8 imagery acquisition, two pairs of dates (or images) were utilized for band comparison, i.e., 25 March 2013(ETM+) and 27 March 2013 (OLI) and 1 April 2013(OLI) and 3 April 2013(ETM+). These time-interleaved Landsat scenes help to reduce the impacts of phenological features or weather fluctuation (e.g., showers) on comparative analysis between spectral bands. Then, surface reflectance values were extracted from each spectral band based on the 3,311 sample points of four land cover types. Figures 5 and 6, respectively, display the column plots of the surface reflectance of the spectral bands of both sensors extracted based on the 1,528 and 1,783 sample points. Comparative analysis between spectral bands of the ETM+ and OLI images acquired only two days apart showed that the average surface reflectance values of each band differed slightly. However, the relationship of the surface reflectance of related bands for different land cover types displayed high similarity between both sensors, as they are in analogous spectral regions. For example, the NIR band is very sensitive to growing vegetation (primary forest); SWIR bands are sensitive to soil-exposed land surface, e.g., burned area and unplanted cropland; and a combination of green and SWIR bands perform better in detecting water bodies.
Additionally, some important variations between ETM+ and OLI were obtained in this study. First, among the six spectral bands, only the average surface reflectance values of the SWIR 2 band of the OLI sensor for the four land cover types consistently decreased. Second, the average surface reflectance values of the NIR band of OLI for vegetative types (i.e., primary forest and unplanted cropland) increased slightly, while those for non-vegetative types (e.g., burned area and water body) decreased, on the contrary. It therefore could be concluded that the OLI sensor improved its performance in discriminating vegetative land cover from non-vegetative types. Third, the average surface reflectance values of the SWIR 1 band of OLI for non-water land cover types basically increased, while those for water bodies decreased. Fourth, the spectral bands of burned areas slightly increased over time, while the NIR and SWIR 2 bands further decreased. This highlights the potential for distinguishing burned areas from unburned areas with these bands. Fifth, the NIR and SWIR bands of the OLI sensor for water bodies decreased in comparison with those of the ETM+ sensor, while the visible bands showed completely inverse variations.

Comparison between the Values of Vegetation Indices Derived from ETM+ and OLI Sensors
Vegetation index values of eight sample plots derived from Landsat-8 OLI images acquired on 27 March 2013 and 1 April 2013 were used to compare with those derived from Landsat-7 ETM+ images acquired on 25 March 2013 and 3 April 2013. Pair comparative analysis of four vegetation indices (NDVI, LSWI, MNDWI and NBR) for each sample plot basically demonstrated that there existed subtle differences between the Landsat-7 ETM+ and Landsat-8 OLI sensors (Figure 7). This was clearly illustrated, especially by LSWI, NBR and NDVI, because their average difference values were less than ±0.05, and they had a smaller standard deviation. LSWI might be the optimum parameter (approximated to zero), which demonstrated a very high similarity of both sensors, followed by NBR and NDVI. However, MNDWI performed relatively poorly, because the average difference values reached up to −0.15. As shown in Figure 7, it was very important to note that vegetation indices calculated with non-visible spectral bands, mainly LSWI and NBR, were more stable for comparative analysis between the ETM+ and OLI sensors. The fluctuations of mean values of differenced vegetation indices (LSWI and NBR) among the eight sample plots were very small, less than 0.05. However, the average values of NDVI and MNDWI, calculated with visible and near-infrared or short-wave near-infrared spectral bands, either fluctuated obviously or were non-zero. Our analysis results also validated the possible spectral band difference effects (SBDEs) when using visible spectral bands for calculating vegetation indices [1,3]. Particularly, the SBDEs caused by the red band were far less than the green band. This is evidenced by the large standard deviation of difference values of MNDWI. It therefore could be concluded that the notable fluctuation of difference values of MNDWI among sample plots might be induced by the SBDE.
Another factor that led to the fluctuation of the difference values of vegetation indices might be the different land cover types among eight sample plots. As displayed in Figure 7, the difference values of vegetation indices for unplanted croplands dominated the landscape (i.e., sample plot 1, 4 and 8) and were basically less than ±0.02. However, the difference values of NBR (nearly ±0.05) were greater than those of the other three indices on the whole. This anomaly might have a close connection with cropland burning during the hot season [36,37]. Actually, we indeed noticed that the vast area of croplands were uncultivated because of a shortage of irrigation water during our field trip in northern Thailand. The weeds and crop residuals in the croplands turned yellow. Some fields had been burned, while other more were waiting for burning in the next few days.

Correlation Analysis of Vegetation Indices Derived from ETM+ and OLI
Scatter plots with correlation coefficients of vegetation indices (NDVI, LSWI, MNDWI and NBR) consistently showed that the imagery acquired two days apart has very high correlation between Landsat-7 ETM+ and Landsat-8 OLI (Figures 8-11). In each plot, Landsat-8 OLI-derived vegetation indices are displayed on the Y-axis, while Landsat-7 ETM+-derived vegetation indices are shown on the X-axis. Coefficients of determination (R 2 ) of each plot were greater than 0.96 on the whole. They demonstrated that vegetation indices obtained from Landsat-7 ETM+ and Landsat-8 OLI were highly linearly correlated. Therefore, it can be concluded that the imagery of both sensors can be used as complementary data. For all eight sample plots, the mean values of slope and the Y-intercept were used to form the inter-calibration equations of vegetation indices between Landsat-7 ETM+ and Landsat-8 OLI imagery as follows: However, there still existed some differences among the four vegetation indices. To begin with, the two indices (LSWI and NBR), computed from the near-infrared (NIR) and short-wave infrared (SWIR 1/2) bands, had slightly larger R 2 values (over 0.99) than those of NDVI and MNDWI, which are computed using visible bands (the red and green bands) and the NIR/SWIR bands for each sample plot. The NBR obtained the highest average R 2 value (0.999), while the MNDWI got the lowest (0.983) among eight sample plots. Then, the LSWI and NBR were highly correlated and basically followed the 1:1 line, while the NDVI and MNDWI also had high correlation coefficients, but did not always follow the 1:1 line. The fluctuations of NDVI and MNDWI values extracted from both sensors further proved the potential spectral band difference effects (SBDEs). Based on Equations (5)(6)(7)(8)

Conclusions
Landsat-8 restarted the multi-decadal Landsat data record of global surface changes after Landsat-5 discontinued its Earth observation, and Landsat-7 has continuous problems in SLC and limited fuel [3,4,38]. Currently, Landsat-8 and Landast-7 will assure an eight-day repeat coverage for monitoring the entire Earth, especially for mapping crop growth and forest disturbances. Consequently, comparison analysis between Landsat-8 and its predecessors will be of significant importance for detecting spatial and temporal land surface changes, since the Landsat satellites have collected invaluable global imagery for over 40 years [39].
In this study, practical and repeatable comparative analyses of Landsat-7 ETM+ and Landsat-8 OLI were carried out from two aspects: spectral bands comparison with sample points and vegetation indices comparison with sample plots. We used Landsat imagery (ETM+ SLC-off and OLI Pre-WRS-2 data) that was just two days apart, which helps to reduce the impacts of phenology and weather fluctuations. Spectral bands comparison of both sensors showed that the spectral bands differed slightly, but with high similarities. Specifically, in comparing with ETM+, first, the OLI had higher values of the near-infrared band for vegetative land cover types, but lower values for non-vegetative types. Second, the new sensor had lower values for the shortwave infrared (2.11-2.29 µm) band for all land cover types. Third, it also basically had higher values for the shortwave infrared (1.57-1.65 µm) band for non-water land cover types. The difference of vegetation indices comparison showed that there were subtle differences between both sensors, which demonstrated high similarity. LSWI, NBR and NDVI had average difference values less than ±0.05 and a smaller standard deviation. LSWI might be the optimum parameter (approximated to zero), followed by NBR and NDVI. However, MNDWI performed relatively poorly, because the average difference values reached up to −0.15. Additionally, correlation analysis of the vegetation indices indicated that both sensors had a very high linear correlation coefficient, with R 2 greater than 0.96. The subtle differences and high correlation of vegetation indices demonstrated that ETM+ and OLI imagery can be used as complementary data. Lastly, the NIR and SWIR bands-calculated vegetation indices (i.e., LSWI and NBR) performed better than NDVI and MNDWI for cross-comparison analysis of satellite sensors, due to the spectral band difference effects.
To our knowledge, this study may be the first comparative analysis between Landsat-8 OLI and Landsat-7 ETM+. It preliminarily compares the difference and correlation of spectral bands and vegetation indices from both sensors. However, it is worthwhile to emphasize that more comparison analysis between Landsat-8 OLI/TIRS and other sensors should be carried out in the future. As the band wavelength gets shorter and the entire detector system changes from whisk-broom to push-broom design [18], related comparative studies between Landsat-8 OLI/TIRS and previous Landsat satellites will deserve more attention. Such research work will lay the foundation for integrating historical Landsat imagery over 40 years for monitoring the spatio-temporal change, pattern and process of the Earth ecosystem.