Next Article in Journal
Data Integration for Climate Vulnerability Mapping in West Africa
Previous Article in Journal
Evaluation of the Consistency of MODIS Land Cover Product (MCD12Q1) Based on Chinese 30 m GlobeLand30 Datasets: A Case Study in Anhui Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Effect of Temporal Interval Length on the Blending of Landsat-MODIS Surface Reflectance for Different Land Cover Types in Southwestern Continental United States

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
Urban Science Department, College of Applied Arts and Science of Beijing Union University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2015, 4(4), 2542-2560; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi4042542
Submission received: 18 September 2015 / Revised: 4 November 2015 / Accepted: 9 November 2015 / Published: 17 November 2015

Abstract

:
Capturing spatial and temporal dynamics is a key issue for many remote-sensing based applications. Consequently, several image-blending algorithms that can simulate the surface reflectance with high spatial-temporal resolution have been developed recently. However, the performance of the algorithm against the effect of temporal interval length between the base and simulation dates has not been reported. In this study, our aim was to evaluate the effect of different temporal interval lengths on the accuracy using the widely used blending algorithm, Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM), based on Landsat, Moderate-resolution Imaging Spectroradiometer (MODIS) images and National Land Cover Database (NLCD). Taking the southwestern continental United States as the study area, a series of experiments was conducted using two schemes, which were the assessment of STARFM with (i) a fixed base date and varied simulation date and (ii) varied base date and specific simulation date, respectively. The result showed that the coefficient of determination (R2), Root Mean Squared Error (RMSE) varied, and overall trend of R2 decreased along with the increasing temporal interval between the base and simulation dates for six land cover types. The mean R2 value of cropland was lowest, whereas shrub had the highest value for two schemes. The result may facilitate selection of an appropriate temporal interval when using STARFM.

1. Introduction

According to the fundamentals of satellite sensor design, a trade-off must be made among the spatial, temporal, and spectral resolutions. A single-sensor is constrained by the specific aim, which constitutes this sensor-specific data framework [1]; therefore, there is no single satellite sensor that can produce multispectral/hyperspectral images with both fine spatial and temporal resolution. To resolve this constraint, several spatial-temporal image fusion models have been developed to produce high spatial and temporal resolution reflectance, which has aroused great interest within the remote-sensing community [2]. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) was initially proposed, and has been documented to have the capacity of blending Landsat Thematic Mapper (TM)/Enhanced TM Plus (ETM+) images (30 m, 16-day) and Moderate-resolution Imaging Spectroradiometer (MODIS) images (500 m, daily) to simulate the daily surface reflectance at the Landsat spatial resolution and MODIS temporal frequency [3]. Subsequently, this model has been applied in many research fields, such as forest disturbance [4], evapotranspiration [5], net ecosystem exchange [6], gross primary production [7], land cover classification [8], phenology [9], leaf area index [10], and public health [11]. An enhanced STARFM, called ESTARFM [12], allowed overcoming some limitations of the original in heterogeneous landscape but was not superior in homogeneous regions with high temporal variance [1]. Meanwhile, several other improved and derived versions based on STARFM have subsequently been developed [4,13,14,15].
However, the spatial-temporal image fusion models have limitations when the simulation dates are distance from the available image pair at base date [16]. Until now, the temporal interval length between the base date and simulation date (hereafter referred to as Tbs) has been determined empirically using the available data when the spatial-temporal image fusion model has been used. No quantitative assessment of spatial-temporal image fusion model has been undertaken when the Tbs is changed. Therefore, if we need to simulate surface reflectance with high spatial resolution, which does not exist on a specific date, the issue arises of how long Tbs is suitable, with a predefined accuracy, using spatial-temporal image fusion. To our knowledge, there is no published literature evaluating the performance of spatial-temporal image fusion against the effect of Tbs.
In this study, we assessed the performance of spatial-temporal image fusion with an increasing Tbs and analyzed the accuracy with different Tbs values for different land cover types. The STARFM was used with one-pair of Landsat-MODIS at a base date for the study. Combined with land cover data, the main objective of this study was to give a quantitative analysis on the performance of STARFM with (i) a fixed base date, varied simulation dates and (ii) varied base dates, a specific simulation date.

2. Data and Methods

2.1. Study Area

The study area was located in the southwest of the continental United States with an area of 36 × 30 km2 (Figure 1), which is for the Central Arizona-Phoenix Long-Term Ecological Research (CAP LTER) project [17]. The study area has a dry (mean annual precipitation of 180 mm) and warm (mean summer temperature of 30.8 °C) climate with two distinct wet seasons (one is in summer, the other is in winter). The native vegetation is dominated by two subdivisions of the Sonoran desert scrub: One is Arizona Upland subdivision with Paloverde-Mixed Cacti series (composed of Cercidium microphyllum, Olneya tesota, Simmondsia chinensis, Larrea tridentata, Encelia farinosa, Fouquieria splendens, Carnegiea gigantea and Opuntia sp.), the other is Lower Colorado River subdivision with Creosotebush-Bursage series (Larrea tridentata and Ambrosia dumosa) [18]. The other native riparian vegetation is characterized by riparian scrubland along minor drainages. The managed vegetation for the study area is mainly cropland [19]. The land cover types in this region include cropland, grassland, shrub, urban, water, and wetland, which were reclassified from the eight broad classes based on National Land Cover Database (NLCD) classification system [20,21,22]. Specifically, cropland type includes cultivated crops; grassland type includes grassland/herbaceous and pasture/hay; shrub type includes dwarf scrub and shrub/scrub; urban type includes (i) developed, open space, (ii) developed, low intensity, (iii) developed, medium intensity, and (iv) developed high intensity; water type includes open water and perennial ice/snow; and wetland type includes woody wetlands and emergent herbaceous wetlands. The other land cover types comprised less than 0.1% of the study area, consequently, they were not taken into account in this study. More than half of the cropland and wetland has two growing seasons, while grassland and shrub has one growing season per year [23]. The peak of growth for shrub, wetland is in early spring, early summer to early autumn, respectively. Cropland has the shortest growth length and fastest green-up rate. Grassland has the longest growth length, which never drops below a certain value. Shrub growth starts from early August or November, and ends in early summer. Wetland, with two growing seasons, green up in mid-February and defoliate in July. The date of beginning and end growth for cropland are closer to wetland than shrub [23].

2.2. Data

The datasets used in this study were as follows: (i) Landsat TM/ETM+ surface reflectance (2001–2012), (ii) MODIS reflectance (2001–2012) and (iii) land cover data (2001, 2006 and 2011).
The Landsat TM/ETM+ surface reflectance product (path/row: 037/037) covering the study area was acquired form EarthExplorer of United States Geological Survey (USGS) (http://earthexplorer.usgs.gov/) during the period 2001 to 2012. Only Landsat surface reflectance data with overall cloud cover of less than 20% were selected as candidates, which can cover the most of the Landsat data within the extent of study area without cloud contaminated. Candidate data were then clipped to the extent of the study area (36 × 30 km2), and the clipped data with cloud cover of more than 15% were not used based on the corresponding cfmask band. A total of 299 scenes of Landsat data were finally used (TM: 142, ETM+: 157), and information of frequency and day of year (DOY) distribution is shown in Figure 2. The gap within Landsat 7 ETM+ Scan Line Corrector (SLC)-off data was not taken into account for the performance assessment of spatial-temporal image fusion.
Figure 1. Land cover map (National Land Cover Database (NLCD), year: 2001) of the study area.
Figure 1. Land cover map (National Land Cover Database (NLCD), year: 2001) of the study area.
Ijgi 04 02542 g001
Figure 2. Distribution of the Landsat data used in the study. (a) day of year (DOY), (b) frequency from the period 2001 to 2012. L5 denotes Landsat TM data, L7 denotes Landsat ETM+ data.
Figure 2. Distribution of the Landsat data used in the study. (a) day of year (DOY), (b) frequency from the period 2001 to 2012. L5 denotes Landsat TM data, L7 denotes Landsat ETM+ data.
Ijgi 04 02542 g002
For MODIS data, MOD09A1 product data (8-day reflectance, 500 m, Collection 5, tile: h08v05) for the study area were obtained from the National Aeronautics and Space Administration (NASA) Earthdata portal (https://search.earthdata.nasa.gov/) for the period 2001 to 2012. The MODIS data were re-projected, resampled and clipped to the same spatial extent (36 × 30 km2) and resolution as the Landsat projection using the MODIS Reprojection Tool (MRT) (available at URL https://lpdaac.usgs.gov/tools/modis_reprojection_tool). Any invalid MOD09A1 data were eliminated using the quality assurance (QA) layer included in the product [6].
For land cover data, NLCD with a spatial resolution of 30 m were obtained from the Multi-Resolution Land Characteristics Consortium (MRLC) (http://www.mrlc.gov/index.php) for the years 2001, 2006, and 2011. We assumed that the status of the study area during 2001–2003 could be represented by NLCD 2001; the status of study area during 2004–2008 could be represented by NLCD 2006; and the status of study area during 2009–2012 could be represented by NLCD 2011. Nationally, the overall accuracy of NLCD 2001, and 2006 is 79% and 78%, respectively [24]. For NLCD 2011, the accuracy assessment is currently underway [25].

2.3. STARFM

STARFM, which was developed by Gao, Masek, Schwaller and Hall [3], was used in the study. Only one pair of Landsat and MODIS images acquired on the base date and one MODIS observation on the simulation date were taken into account.
Landsat-like surface reflectance was obtained on the simulation date using the following steps. First, the spectral similar neighboring pixels within a local moving window of Landsat data were identified. Second, a weight Wijk was calculated for each spectral similar neighboring pixel based on: (i) the spectral difference between the Landsat and MODIS data on the base date, (ii) the temporal difference of the MODIS data between the base and simulation date, and (iii) the spatial Euclidean distance between the neighbor and the central pixel within the local moving window. Finally, the surface reflectance of the central pixel was calculated as follows:
L ( x w / 2 , y w / 2 , t s ) = i = 1 w j = 1 w b = 1 n W i j k × ( M ( x i , y j , t s ) + L ( x i , y j , t b ) M ( x i , y j , t b ) )
where L and M indicate the surface reflectance of Landsat and MODIS, respectively; L(xw/2, yw/2, ts) is the surface reflectance of the central pixel (xw/2, yw/2) on the simulation date ts for Landsat; M(xi, yj, tb) is the surface reflectance of pixel (xi, yj) within local moving window on the base date (tb) for MODIS; and w is the size of the local moving window. For more detailed information of STARFM, see Gao, Masek, Schwaller and Hall [3].

2.4. Evaluation with a Fixed Base Date

To evaluate the performance of STARFM with a fixed base date, we changed the value of Tbs. Specifically, one date was selected as the base date for every month of 2001, 2003 and 2008. For each month, only the Landsat data with the lowest cloud cover percentage were selected as the base data, and the corresponding date was defined as base date (Table 1). The number of base dates in 2008 was 20, which is more than 12. When evaluating the performance of STARFM with each fixed base date, the value of Tbs increased until it reached the last simulation date in 2012.
Four pairs of Landsat-MODIS data on the base dates of 9 February, 16 May, 28 August, and 8 November 2001, were selected. Land cover data from the NLCD for 2001, 2006, and 2011 were used to assess the performance of STARFM on a fixed base date. Based on four fixed base dates (9 February, 16 May, 28 August, and 8 November 2001), the performance of STARFM was evaluated according to six land cover types. We defined four seasons and the corresponding months as winter (December, January and February), spring (March, April and May), summer (June, July and August), and autumn (September, October and November). We assumed that each Landsat-MODIS pair on the base date represented the average status of the study areas in the corresponding season.
Table 1. Details of the base date for the spatial-temporal image fusion. Abbreviations: January (Jan.), February (Feb.), March (Mar.), April (Apr.), June (Jun.), July (Jul.), August (Aug.), September (Sep.), October (Oct.), November (Nov.), December (Dec.). For column of Sensor, L5 denotes Landsat TM, and L7 denotes Landsat ETM+. *The base date for February 2003 is 30 January 2003, because no Landsat data is available in February 2003 and the Landsat data on 30 January 2003 is the closest one to February 2003.
Table 1. Details of the base date for the spatial-temporal image fusion. Abbreviations: January (Jan.), February (Feb.), March (Mar.), April (Apr.), June (Jun.), July (Jul.), August (Aug.), September (Sep.), October (Oct.), November (Nov.), December (Dec.). For column of Sensor, L5 denotes Landsat TM, and L7 denotes Landsat ETM+. *The base date for February 2003 is 30 January 2003, because no Landsat data is available in February 2003 and the Landsat data on 30 January 2003 is the closest one to February 2003.
YearMonthDayDOYSensorDay since Start of the Dataset
Fixed Base Date (2001)Specific Simulation Date (2009)
Feb. 9May 16Aug. 28Nov. 8Jan. 14Apr. 12Jul. 9Oct. 21
2001Jan.2424L5----2912300030883192
Feb.940L5----2896298430723176
Mar.1372L532---2864295230403144
Apr.14104L564---2832292030083112
May16136L596---2800288829763080
Jun.17168L512832--2768285629443048
Jul.19200L516064--2736282429123016
Aug.28240L7200104--2696278428722976
Sep.21264L522412824-2672276028482952
Oct.23296L525616056-2640272828162920
Nov.8312L527217672-2624271228002904
Dec.18352L7312216112402584267227602864
2003Jan.1414L57046085044322192228023682472
Feb.30*30L57206245204482176226423522456
Mar.1170L77606645604882136222423122416
Apr.20110L58007046005282096218422722376
May22142L58327366325602064215222402344
Jun.23174L58647686645922032212022082312
Jul.9190L58807846806082016210421922296
Aug.10222L59128167126401984207221602264
Sep.11254L59448487446721952204021282232
Oct.13286L59768807767041920200820962200
Nov.14318L510089128087361888197620642168
Dec.16350L510409448407681856194420322136
2008Jan.2020L72536244023362264360448536640
Feb.536L72552245623522280344432520624
2960L52576248023762304320408496600
Mar.2484L72600250424002328296384472576
Apr.17108L52624252824242352272360448552
25116L72632253624322360264352440544
May3124L52640254424402368256344432536
11132L72648255224482376248336424528
2008Jun.12164L72680258424802408216304392496
20172L52688259224882416208296384488
Jul.22204L52720262425202448176264352456
30212L72728263225282456168256344448
Aug.23236L52752265625522480144232320424
Sep.24268L52784268825842512112200288392
Oct.10284L5280027042600252896184272376
18292L7280827122608253688176264368
Nov.3308L7282427282624255272160248352
11316L5283227362632256064152240344
Dec.5340L7285627602656258440128216320
29364L5288027842680260816104192296

2.5. Evaluation for a Specific Simulation Date

To evaluate the performance of STARFM on a specific simulation date, the value of Tbs was increased until the farthest based date was reached. Specifically, four simulation dates in 2009 were selected for the study: 14 January, 12 April, 9 July, and 21 October 2009. Therefore, for each simulation date, 40 base date Landsat-MODIS pairs were used to produce 40 simulations (2001, 2003, and 2008) (Table 1). The land cover data from NLCD for 2001, 2006, and 2011 were also used to assess the performance of STARFM on specific simulation dates. We assumed that each simulated reflectance at the simulation date represented the average status of the study areas in the corresponding season.

2.6. Accuracy Assessment

We assessed the accuracy of simulations by comparing the simulated Landsat-like data to a Landsat observation on the same date. A linear regression model (observed versus simulated surface reflectance) was used to assess the performance of STARFM, with reference to the statistical parameter, R2 and Root Mean Squared Error (RMSE), which were used to measure the fitness of the linear regression, differences between simulated and observed reflectance, respectively.

3. Results

3.1. Landsat-Like Surface Reflectance with a Fixed Base Date

The R2 and RMSE value varied seasonally, and R2 value was generally decreased as the value of the temporal interval between the base and simulation dates (Tbs) increased. The RMSE value for band 1 is lowest compared to other bands. Generally, the DOYs of simulation dates for best/worst performance (maximum R2 and minimum RMSE/ minimum R2 and maximum RMSE) at each simulation year were near to/far from the DOY of base date (Figure 3 and Figure 4).
Figure 3. Distribution of the day of year (DOY) for maximum R2 (not filled marker) and minimum RMSE (filled marker) value at each year from 2001 to 2012. The corresponding base date (dash line) is (a) 9 February (DOY: 40), (b) 16 May (DOY: 136), (c) 28 August (DOY: 240) and (d) 8 November (DOY: 312) for year 2001, respectively. The markers for six Landsat-like bands are circle, triangle, diamond, square, hexagram and pentagram, respectively. The color for each land cover type is black (WAT: water), red (URB: urban), green (SHR: shrub), blue (GRA, grassland), violet (CRO: cropland) and brown (WET: wetland). The meaning of marker and color for symbols is the same for Figure 4.
Figure 3. Distribution of the day of year (DOY) for maximum R2 (not filled marker) and minimum RMSE (filled marker) value at each year from 2001 to 2012. The corresponding base date (dash line) is (a) 9 February (DOY: 40), (b) 16 May (DOY: 136), (c) 28 August (DOY: 240) and (d) 8 November (DOY: 312) for year 2001, respectively. The markers for six Landsat-like bands are circle, triangle, diamond, square, hexagram and pentagram, respectively. The color for each land cover type is black (WAT: water), red (URB: urban), green (SHR: shrub), blue (GRA, grassland), violet (CRO: cropland) and brown (WET: wetland). The meaning of marker and color for symbols is the same for Figure 4.
Ijgi 04 02542 g003
Specifically, with one pair of Landsat-MODIS images as base data on 9 February 2001 (DOY: 40), the high R2 and low RMSE values often occurred at the beginning (base date ± 39 days) or end (last month) of the year (Figure 3a). The DOY of simulated reflectance with low R2 and high RMSE value at each year was mainly distributed from 140 to 300 (Figure 4a). When the temporal distance was within ±70 days between simulation date and base date, the high R2 and low RMSE value for water and grass can also be obtained. Among each land cover type, the performance of STARFM was worst for cropland, which median and mean R2 values for the six bands were lower than 0.25, the RMSE value for band 4 had reached to 0.08 (Figures S1a, S2a, and S3a). The best performance of STARFM was simulated result for shrub, the median and mean R2 values for each Landsat band was higher than any other land cover type, the simulated reflectance of band 5 was closer (R2) to the observed reflectance than to the other bands, and result of band 1 had lowest bias compared to other bands (Figures S1c, S2c, and S3c).
Figure 4. Distribution of the day of year (DOY) for minimum R2 (not filled marker) and maximum RMSE (filled marker) value at each year from 2001 to 2012. The corresponding base date (dash line) is (a) 9 February (DOY: 40), (b) 16 May (DOY: 136), (c) 28 August (DOY: 240) and (d) 8 November (DOY: 312) for year 2001, respectively.
Figure 4. Distribution of the day of year (DOY) for minimum R2 (not filled marker) and maximum RMSE (filled marker) value at each year from 2001 to 2012. The corresponding base date (dash line) is (a) 9 February (DOY: 40), (b) 16 May (DOY: 136), (c) 28 August (DOY: 240) and (d) 8 November (DOY: 312) for year 2001, respectively.
Ijgi 04 02542 g004
When the base date was moved to 16 May 2001 (DOY: 136), the highest R2 and low RMSE often occurred in the middle (base date ± 30 days) of the year for most land cover types. However, the high R2 and low RMSE value can also been obtained when the simulation DOY was extended to the beginning of the year for water, cropland, wetland and the band 5 and 7 of urban (Figure 3b). The DOY of simulated reflectance with low R2 and high RMSE value at each year was mainly distributed at two extents (150–300 and 340–365) (Figure 4b). The simulated reflectance of shrub still had high R2 and low RMSE value (Figures S4c, S5c, and S6c). For cropland, the situation is reversed (Figures S4a, S5a, and S6a). For urban areas, the R2 value had two significant declines, one from 2003 to 2004, and the other from 2008 to 2009. The upper 25% of data are quite sparse in comparison with the lower 25% for urban land cover type, the RMSE value had a slight increase from year 2003 to 2004 (Figures S4d, S5d, and S6d).
For base date of 28 August 2001 (DOY: 240), the distribution of DOY for simulated reflectance with high R2 and low RMSE value at each year was different compared to previous two base dates. The extent of DOY with high R2 (base date ± 40 days) was different from low RMSE value for water, urban, cropland and wetland (Figure 3c). The high RMSE for water, urban and cropland can be found near the base date (Figure 4c). The R2 value near the base date was higher than other three base dates. Specifically, for cropland, the R2 value was 0.7 or higher, with the exception of band 4. An obvious seasonal variation was also apparent (Figures S7a and S8a). For grassland, the R2 value was 0.9 or even higher for each band. The trend was similar to the results of the other three base date (Figures S7b, S8b, and S9b). For shrub, the simulated reflectance was closer to the observed reflectance, with R2 values higher than 0.9 and almost RMSE values lower than 0.04 (Figures S7c, S8c, and S9c). For other land cover types, the trends of the R2 value were similar to that of the R2 value for the base date of 16 May 2001, but there was a higher R2 value for the base date 28 August 2001 (Figures S7d,e,f, S8d,e,f, and S9d,e,f).
For the base date of 8 November 2001 (DOY: 312), the DOY for high R2 and low RMSE often distributed at the end (300–365) or beginning (1–60) of the year (Figure 3d). The R2 of cropland was satisfactory at the beginning of the simulation date, and was 0.8 or higher for bands 2 and 3. The corresponding RMSE value is lower than 0.04. The simulated reflectance for band 4 still held the lowest median R2 value (<0.1) (Figures S10a, S11a, and S12a). For the other five land cover types, the trend in the R2 and RMSE value was similar to that for the R2 and RMSE value for the base date of 9 February 2001, but with higher R2 value for the base date of 8 November 2001 (Figures S10, S11, and S12).

3.2. Landsat-Like Surface Reflectance on A Specific Simulation Date

Similarly, the DOYs of base dates for best/worst performance (maximum R2 and minimum RMSE/ minimum R2 and maximum RMSE) at each base year were near to/far from the DOY of simulation date (Figure 5 and Figure 6). When the simulation date was 14 January 2009 (DOY: 14), the R2 and RMSE value varied as the Tbs increased, and the overall trend of R2 become more apparent. The high R2 and low RMSE value often occurred at the beginning and end of the year (Figure 5a), and the R2 value at the end was higher than at the beginning of the year 2008 for six land cover types (Figures S13 and S14). The low R2 and high RMSE value mainly distributed at the middle of the year (DOY: 60–210) (Figure 6a). For cropland, the R2 values of the first three and last two bands exceeded 0.70 and 0.60 for the nearest two base dates (5 December and 29 December 2008), and the corresponding Tbs was 40 and 16 days, respectively (Figure S13a). The corresponding magnitude of RMSE for cropland at year 2001, 2003 and 2008 were similar (Figure S14a). The patterns of R2 distributions for grassland and shrub were also similar, but the mean R2 value of shrub was higher than grassland (Figure S13b,c). The R2 value of urban for year 2008 was much higher than for years 2001 and 2003; this is also reflected in the lower RMSE value for year 2008 compared to other two years (Figures S13d and S14d). The trend of R2 for water and wetland was similar, but the trend of RMSE was different (Figures S13e,f and S14e,f).
Figure 5. Distribution of the day of year (DOY) for maximum R2 (not filled marker) and minimum RMSE (filled marker) value for year 2001, 2003 and 2008. The corresponding simulation date (dash line) is (a) 14 January (DOY: 14), (b) 12 April (DOY: 102), (c) 9 July (DOY: 190) and (d) 21 October (DOY: 294) for year 2009, respectively. The size of marker from small to big denotes years 2001, 2003 and 2008, respectively. The markers for six Landsat-like bands are circle, triangle, diamond, square, hexagram and pentagram, respectively. The color for each land cover type is black (WAT: water), red (URB: urban), green (SHR: shrub), blue (GRA, grassland), violet (CRO: cropland) and brown (WET: wetland). The meaning of marker and color for symbols is the same for Figure 6.
Figure 5. Distribution of the day of year (DOY) for maximum R2 (not filled marker) and minimum RMSE (filled marker) value for year 2001, 2003 and 2008. The corresponding simulation date (dash line) is (a) 14 January (DOY: 14), (b) 12 April (DOY: 102), (c) 9 July (DOY: 190) and (d) 21 October (DOY: 294) for year 2009, respectively. The size of marker from small to big denotes years 2001, 2003 and 2008, respectively. The markers for six Landsat-like bands are circle, triangle, diamond, square, hexagram and pentagram, respectively. The color for each land cover type is black (WAT: water), red (URB: urban), green (SHR: shrub), blue (GRA, grassland), violet (CRO: cropland) and brown (WET: wetland). The meaning of marker and color for symbols is the same for Figure 6.
Ijgi 04 02542 g005
When the simulation date was 12 April 2009 (DOY: 102), the R2 and RMSE value varied with the increasing of Tbs, and the overall trend of the accuracy was flat compared to the simulation date of 14 January 2009. The high R2 values often occurred in the middle of the year (Figure 5b). Two peaks of R2 existed for six land cover types except simulated reflectance of band 4 for cropland and wetland. Meanwhile, the magnitude of RMSE was lower at the middle of the year (Figures S15 and S16). For cropland, no R2 value exceeded 0.60 for each band, and RMSE value was slightly higher than the result of first simulation date, 14 January 2009 (Figures S15a and S16a). The trend of R2 for six land cover types was similar compared to the first simulation date.
When the simulation date was 9 July 2009 (DOY: 190), the result pattern (R2, RMSE) was similar to second simulation date, 12 April 2009 (Figure 5c, Figure S17 and Figure S18). Combining the R2 and RMSE, the accuracy of simulated reflectance of cropland was lower than the second simulation date (Figures S17a and S18a). The RMSE of band 4, 5, and 7 for water was still higher than the first three bands (Figure S18e).
Figure 6. Distribution of the day of year (DOY) for minimum R2 (not filled marker) and maximum RMSE (filled marker) value for year 2001, 2003 and 2008. The corresponding simulation date (dash line) is (a) 14 January (DOY: 14), (b) 12 April (DOY: 102), (c) 9 July (DOY: 190) and (d) 21 October (DOY: 294) for year 2009, respectively. The size of marker from small to big denotes year 2001, 2003 and 2008 respectively.
Figure 6. Distribution of the day of year (DOY) for minimum R2 (not filled marker) and maximum RMSE (filled marker) value for year 2001, 2003 and 2008. The corresponding simulation date (dash line) is (a) 14 January (DOY: 14), (b) 12 April (DOY: 102), (c) 9 July (DOY: 190) and (d) 21 October (DOY: 294) for year 2009, respectively. The size of marker from small to big denotes year 2001, 2003 and 2008 respectively.
Ijgi 04 02542 g006
When the simulation date was set as 21 October 2009 (DOY: 294), the R2 at the near side of the simulation date had higher value (Figure 5d). The magnitude of RMSE was also similar to the other simulation date (Figures S19 and S20). However, the shape of R2 value was different to the first simulation date. In the other words, the high value of R2 occurred at the edge of the year for the first simulation date (14 January 2009), but the high R2 value occurred at the near side of the simulation date for the last simulation date (21 October 2009).

4. Discussion

For the six land cover types within the study area, the performance of STARFM was evaluated by two ways, one was a fixed base date and the other was a fixed simulation date.
Four base dates (9 February, 16 May, 28 August and 8 November 2001) were selected from the total 44 base dates, and were assumed to represent the average status of winter, spring, summer and autumn. The selection of different base dates affected the accuracy of the simulated reflectance. This phenomenon was also observed in the results using a specific simulation date. The high R2 and corresponding low RMSE value was often located at the same season. For the R2 value in the same season (base or other year), the longer the temporal distance was, the lower R2 value obtained. Meanwhile, the R2 value at the same season of other year might be higher than at the different season of the same year. The reason for high R2 value is the similar vegetation phenology or land surface status between the base and simulation date [26,27]. For example, high R2 value always occurred near the beginning or end of the simulation year when the base dates of 9 February 2001 (Figure 7) and 8 November 2001 (Figure 3) were selected. A high R2 and low RMSE value always occurred in the middle of the simulation year when the base dates of 16 May 2001 and 28 August 2001 were selected (Figure 3). For the results obtained using a specific simulation date, a similar status was also found (Figure 5). This may indicate that a good result of STARFM can be obtained in the situations when the base and simulation dates are in the same season. However, the temporal distance from base/simulation date is varied for different land cover types.
Figure 7. Visual comparison between the observed and simulated Landsat reflectance (Near-infrared (NIR)-red-green composite) using the base Landsat-MODIS pair data at 9 February 2001. The upper/lower row is the observed/simulated Landsat reflectance.
Figure 7. Visual comparison between the observed and simulated Landsat reflectance (Near-infrared (NIR)-red-green composite) using the base Landsat-MODIS pair data at 9 February 2001. The upper/lower row is the observed/simulated Landsat reflectance.
Ijgi 04 02542 g007
Table 2. Percentage (%) of land cover change for the three periods. The abbreviations CRO, GRA, SHR, URB, WAT, WET and OTH denote cropland, grassland, shrub, urban, water, wetland and others, respectively. The subscripts ‘from’ and ‘to’ denote the original and latest land cover type, respectively.
Table 2. Percentage (%) of land cover change for the three periods. The abbreviations CRO, GRA, SHR, URB, WAT, WET and OTH denote cropland, grassland, shrub, urban, water, wetland and others, respectively. The subscripts ‘from’ and ‘to’ denote the original and latest land cover type, respectively.
CROtoGRAtoSHRtoURBtoWATtoWETtoOTHtoTotalto
2001
➔2006
CROfrom25.79-0.734.080.040.00-30.65
GRAfrom-0.820.020.050.010.00-0.90
SHRfrom0.09-33.092.030.06--35.27
URBfrom--0.0028.97---28.97
WATfrom--0.000.010.310.01-0.33
WETfrom--0.060.050.013.69-3.81
OTHfrom0.00- 0.00-0.07-0.08
Totalfrom25.880.8233.9135.190.433.770.00100.00
2006
➔2011
CROfrom21.77-0.034.000.09--25.88
GRAfrom-0.770.020.010.010.00-0.81
SHRfrom--32.351.530.03--33.91
URBfrom---35.190.00--35.19
WATfrom--0.010.000.42--0.43
WETfrom---0.02-3.75-3.77
OTHfrom-------0.00
Totalfrom21.770.7732.4040.750.553.750.00100.00
2001
➔2011
CROfrom21.68-0.668.170.130.00-30.65
GRAfrom0.000.770.040.070.020.00-0.90
SHRfrom0.08-31.633.480.080.06-35.34
URBfrom--0.0028.970.000.00-28.97
WATfrom--0.000.010.300.01-0.33
WETfrom--0.060.060.013.67-3.81
OTHfrom0.00--0.00--0.000.01
Totalfrom21.770.7732.4040.750.553.750.00100.00
More detailed information was obtained by combining the land cover types within the study area (Table 2, Figure 8). For cropland, the overall R2 values of the six Landsat bands for four base dates were low, especially for band 4 (near-infrared band). The reason for the low accuracy for simulated reflectance could be attributed to the rapid variation of the cropland canopy structure, especially during the growing season. For the base date of 9 February 2001, the accuracy of simulated reflectance is quite low, the maximum R2 value is 0.50 for band 3; possibly, the reflectance of winter is essentially different compared to other seasons, which lead to the lower R2 value. For a specific simulation date, when the base date was similar to the first simulation date (14 January 2009), the R2 value can exceeded 0.6, whereas the R2 value of the six Landsat bands was lower than 0.6 using other three simulation dates. It is possible that the crops had not been planted yet, and most of this area therefore had no vegetation. The surface ground status would then be almost identical between the base and simulation dates. It can also been reflected from the DOY distribution of lowest R2 and highest RMSE value of each year (Figure 4 and Figure 6). The land cover change of the cropland is another reason for the low accuracy of the simulated reflectance, with 4.08%, 4% and 8.17% of the total study area for cropland area converted to urban area from 2001 to 2006, 2006 to 2011, and 2001 to 2011, respectively. The percentage of overall change for the combination of all three periods was lower than the sum of the separate periods because some change occurred twice in the same pixel. In this case, the most recent land cover class was used [20]. Grassland appears invariant compared to cropland and urban areas. However, the area of grassland is small, the representativeness of the result is limited, and needs to be improved. For urban areas, there are two significant decreases in the R2 value at beginning of 2004 and 2009. The reason for this phenomenon is the assumption of representativeness of the NLCD 2001, 2006, and 2011. The spatial distribution of urban area for year 2001–2003, 2004–2008, and 2009–2012 can be represented by NLCD 2001, NLCD 2006 and NLCD 2011, respectively. With the development of the city within the study area, the urban area is expanded. Consequently, the R2 value for the periods of 2004–2008 and 2009–2012 can be improved when the base date is set within these two periods, respectively. For shrub, the R2 value had an obvious seasonal variation, and remained at a good level. Perhaps the canopy structure for shrub is stable, and this leads to good result of simulated surface reflectance compared to other land cover types. For water, because the reflectivity of water is small, the overall R2 values of the six Landsat bands were not good. The R2 value, which was higher than 0.6, could only be obtained when the simulation date is near to the base date. For wetland, the accuracy of simulated reflectance for band 4 was worse than other bands (Figures S2f, S5f, S9f, and S11f), which might be because wetland contains water surface, and the reflectivity was low. A significant decrease in the R2 was apparent at the beginning of 2009; possibly, the land cover type of wetland had changed, whereas it is not reflected using NLCD 2011.
Figure 8. Spatial distribution of land cover changes for the study area for (a) 2001 to 2006, (b) 2006 to 2011, and (c) 2001 to 2011. The base map is band 4 (near-infrared band) of Landsat data on 8 May 2001. The land cover shown on the map is the latest land cover type.
Figure 8. Spatial distribution of land cover changes for the study area for (a) 2001 to 2006, (b) 2006 to 2011, and (c) 2001 to 2011. The base map is band 4 (near-infrared band) of Landsat data on 8 May 2001. The land cover shown on the map is the latest land cover type.
Ijgi 04 02542 g008
There are several limitations and constraints for this study. First, the STARFM cannot accurately simulate the short term, transient changes that are not recorded in either base date or simulation date [3,12]. Second, the use of eight-day surface reflectance from the MODIS product (MOD09A1) may affect the accuracy of simulated surface reflectance. The phenology of vegetation or other associated factors may change within the eight-day period. Meanwhile, Bidirectional Reflectance Distribution Function (BRDF) effect of MODIS eight-day surface reflectance may also have some influence of the simulated result. Third, the performance of STARFM will be affected if the land cover type of the study area had changed from base date to simulation date. Last, the sensor degradation is not considered for the study, which may affect the accuracy of result with long temporal interval. Therefore, further studies should focus on the following issues to improve the accuracy of the simulated results: (i) the usage of improved versions based on STARFM, and (ii) a combination of data assimilation and STARFM.

5. Conclusion

The performance of the spatial-temporal image fusion model (STARFM) was assessed for temporal interval of different sizes using Landsat and MODIS images from the southwest of the continental United States. Two schemes were applied for the evaluation, one using a fixed base date and the other using a specific simulation date. In addition, the land cover data at three periods (2001, 2006 and 2011) were used in the study. The accuracy of simulated reflectance for six Landsat bands was assessed for each of the six land cover types. Although it is expected that the accuracy will reduce with the increased interval from the base date, this study provides a thorough discussion of possible reasons affecting the accuracy of the simulated data and warns readers about limitations of the study.
Generally, there was an obvious seasonal variation for simulated reflectance for each land cover type. To obtain the simulated reflectance at a particular date, the selection of nearer base date during the same season or the same season from a different year will produce results with higher accuracy compared to the rest of the dates. However, the corresponding Tbs value for higher accuracy is different for different land cover types. According to the result of the study area, the land cover with most variation is cropland, which need data with the nearest date (especially at growing season) during the same year to obtain a simulated reflectance with high accuracy. For a land cover type with a stable structure, such as urban, shrub and grassland, the Tbs value can be more than 100 days, especially during the growing season. For water and wetland, the Tbs value is varied, and the performance of STARFM needs to be improved for these land cover types. The results provided guidance on how to select an appropriate temporal interval for specific applications when using the STARFM method.

Supplementary Files

Supplementary File 1

Acknowledgments

This research is supported by the research grant funded by Natural Science Foundation of China (grant No. 41371362), the research grant funded by China Postdoctoral Science Foundation (grant No. 2014M550871), the research grant of the Direct Youth Foundation in Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences (grant No. Y5SJ2000CX), Major Special Project—the China High-Resolution Earth Observation System. We thank the USGS Earth Resources Observation Systems (EROS) data center for providing Landsat data and the Land Processes Distributed Active Archive Center (LP-DAAC) and MODIS science team for providing free MODIS products.

Author Contributions

Dongjie Fu, Lifu Zhang and Hao Chen conceived and designed the experiments; Dongjie Fu, Hao Chen and Juan Wang performed the experiments; Dongjie Fu, Xuejian Sun and Taixia Wu analyzed the data; and Dongjie Fu drafted the paper. All authors read and edited drafts, and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emelyanova, I.V.; McVicar, T.R.; Van Niel, T.G.; Li, L.T.; van Dijk, A.I.J.M. Assessing the accuracy of blending Landsat-MODIS surface reflectances in two landscapes with contrasting spatial and temporal dynamics: A framework for algorithm selection. Remote Sens. Environ. 2013, 133, 193–209. [Google Scholar] [CrossRef]
  2. Gevaert, C.M.; García-Haro, F.J. A comparison of STARFM and an unmixing-based algorithm for Landsat and MODIS data fusion. Remote Sens. Environ. 2015, 156, 34–44. [Google Scholar] [CrossRef]
  3. Gao, F.; Masek, J.; Schwaller, M.; Hall, F. On the blending of the Landsat and MODIS surface reflectance: Predicting daily Landsat surface reflectance. IEEE Trans. Geosci Remote Sens. 2006, 44, 2207–2218. [Google Scholar]
  4. Hilker, T.; Wulder, M.A.; Coops, N.C.; Linke, J.; McDermid, G.; Masek, J.G.; Gao, F.; White, J.C. A new data fusion model for high spatial- and temporal-resolution mapping of forest disturbance based on Landsat and MODIS. Remote Sens. Environ. 2009, 113, 1613–1627. [Google Scholar] [CrossRef]
  5. Cammalleri, C.; Anderson, M.C.; Gao, F.; Hain, C.R.; Kustas, W.P. Mapping daily evapotranspiration at field scales over rainfed and irrigated agricultural areas using remote sensing data fusion. Agric. For. Meteorol. 2014, 186, 1–11. [Google Scholar] [CrossRef]
  6. Fu, D.; Chen, B.; Zhang, H.; Wang, J.; Black, T.A.; Amiro, B.D.; Bohrer, G.; Bolstad, P.; Coulter, R.; Rahman, A.F.; et al. Estimating landscape net ecosystem exchange at high spatial-temporal resolution based on Landsat data, an improved upscaling model framework, and eddy covariance flux measurements. Remote Sens. Environ. 2014, 141, 90–104. [Google Scholar] [CrossRef]
  7. Singh, D. Generation and evaluation of gross primary productivity using Landsat data through blending with MODIS data. Int. J. Appl. Earth Obs. Geoinform. 2011, 13, 59–69. [Google Scholar] [CrossRef]
  8. Watts, J.D.; Powell, S.L.; Lawrence, R.L.; Hilker, T. Improved classification of conservation tillage adoption using high temporal and synthetic satellite imagery. Remote Sens. Environ. 2011, 115, 66–75. [Google Scholar] [CrossRef]
  9. Walker, J.J.; de Beurs, K.M.; Wynne, R.H.; Gao, F. Evaluation of Landsat and MODIS data fusion products for analysis of dryland forest phenology. Remote Sens. Environ. 2012, 117, 381–393. [Google Scholar] [CrossRef]
  10. Wu, M.; Wu, C.; Huang, W.; Niu, Z.; Wang, C. High-resolution leaf area index estimation from synthetic Landsat data generated by a spatial and temporal data fusion model. Comput. Electron. Agric. 2015, 115, 1–11. [Google Scholar] [CrossRef]
  11. Liu, H.; Weng, Q.H. Enhancing temporal resolution of satellite imagery for public health studies: A case study of West Nile Virus outbreak in Los Angeles in 2007. Remote Sens. Environ. 2012, 117, 57–71. [Google Scholar] [CrossRef]
  12. Zhu, X.L.; Chen, J.; Gao, F.; Chen, X.H.; Masek, J.G. An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions. Remote Sens. Environ. 2010, 114, 2610–2623. [Google Scholar] [CrossRef]
  13. Fu, D.; Chen, B.; Wang, J.; Zhu, X.; Hilker, T. An improved image fusion approach based on enhanced spatial and temporal the adaptive reflectance fusion model. Remote Sens. 2013, 5, 6346–6360. [Google Scholar] [CrossRef]
  14. Huang, B.; Wang, J.; Song, H.; Fu, D.; Wong, K. Generating high spatiotemporal resolution land surface temperature for urban heat island monitoring. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1011–1015. [Google Scholar] [CrossRef]
  15. Weng, Q.; Fu, P.; Gao, F. Generating daily land surface temperature at Landsat resolution by fusing Landsat and MODIS data. Remote Sens. Environ. 2014, 145, 55–67. [Google Scholar] [CrossRef]
  16. Gao, F.; Hilker, T.; Zhu, X.; Anderson, M.; Masek, J.; Wang, P.; Yang, Y. Fusing Landsat and MODIS data for vegetation monitoring. IEEE Geosci. Remote Sens. Mag. 2015, 3, 47–60. [Google Scholar] [CrossRef]
  17. Grimm, N.; Redman, C. Approaches to the study of urban ecosystems: The case of Central Arizona—Phoenix. Urban. Ecosyst. 2004, 7, 199–213. [Google Scholar] [CrossRef]
  18. Buyantuyev, A.; Wu, J. Urbanization alters spatiotemporal patterns of ecosystem primary production: A case study of the Phoenix metropolitan region, USA. J. Arid Environ. 2009, 73, 512–520. [Google Scholar] [CrossRef]
  19. Buyantuyev, A.; Wu, J.; Gries, C. Estimating vegetation cover in an urban environment based on landsat ETM+ imagery: A case study in Phoenix, USA. Int. J. Remote Sens. 2007, 28, 269–291. [Google Scholar] [CrossRef]
  20. Jin, S.; Yang, L.; Danielson, P.; Homer, C.; Fry, J.; Xian, G. A comprehensive change detection method for updating the national Land Cover Database to circa 2011. Remote Sens. Environ. 2013, 132, 159–175. [Google Scholar] [CrossRef]
  21. Fry, J.A.; Xian, G.; Jin, S.; Dewitz, J.A.; Homer, C.G.; LIMIN, Y.; Barnes, C.A.; Herold, N.D.; Wickham, J.D. Completion of the 2006 National Land Cover Database for the conterminous United States. Photogramm. Eng. Remote Sens. 2011, 77, 858–864. [Google Scholar]
  22. Homer, C.; Dewitz, J.; Fry, J.; Coan, M.; Hossain, N.; Larson, C.; Herold, N.; McKerrow, A.; VanDriel, J.N.; Wickham, J. Completion of the 2001 National Land Cover Database for the counterminous United States. Photogramm. Eng. Remote Sens. 2007, 73, 337–341. [Google Scholar]
  23. Buyantuyev, A.; Wu, J. Urbanization diversifies land surface phenology in arid environments: Interactions among vegetation, climatic variation, and land use pattern in the Phoenix metropolitan region, USA. Landsc. Urban Plan. 2012, 105, 149–159. [Google Scholar] [CrossRef]
  24. Wickham, J.D.; Stehman, S.V.; Gass, L.; Dewitz, J.; Fry, J.A.; Wade, T.G. Accuracy assessment of NLCD 2006 land cover and impervious surface. Remote Sens. Environ. 2013, 130, 294–304. [Google Scholar] [CrossRef]
  25. Wickham, J.; Homer, C.; Vogelmann, J.; McKerrow, A.; Mueller, R.; Herold, N.; Coulston, J. The multi-resolution land characteristics (MRLC) consortium—20 years of development and integration of USA national land cover data. Remote Sens. 2014, 6, 7424–7441. [Google Scholar] [CrossRef]
  26. Sexton, J.O.; Urban, D.L.; Donohue, M.J.; Song, C. Long-term land cover dynamics by multi-temporal classification across the Landsat-5 record. Remote Sens. Environ. 2013, 128, 246–258. [Google Scholar] [CrossRef]
  27. Bhandari, S.; Phinn, S.; Gill, T. Preparing Landsat image time series (LITS) for monitoring changes in vegetation phenology in Queensland, Australia. Remote Sens. 2012, 4, 1856–1886. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Fu, D.; Zhang, L.; Chen, H.; Wang, J.; Sun, X.; Wu, T. Assessing the Effect of Temporal Interval Length on the Blending of Landsat-MODIS Surface Reflectance for Different Land Cover Types in Southwestern Continental United States. ISPRS Int. J. Geo-Inf. 2015, 4, 2542-2560. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi4042542

AMA Style

Fu D, Zhang L, Chen H, Wang J, Sun X, Wu T. Assessing the Effect of Temporal Interval Length on the Blending of Landsat-MODIS Surface Reflectance for Different Land Cover Types in Southwestern Continental United States. ISPRS International Journal of Geo-Information. 2015; 4(4):2542-2560. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi4042542

Chicago/Turabian Style

Fu, Dongjie, Lifu Zhang, Hao Chen, Juan Wang, Xuejian Sun, and Taixia Wu. 2015. "Assessing the Effect of Temporal Interval Length on the Blending of Landsat-MODIS Surface Reflectance for Different Land Cover Types in Southwestern Continental United States" ISPRS International Journal of Geo-Information 4, no. 4: 2542-2560. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi4042542

Article Metrics

Back to TopTop