The Global Availability of Landsat 5 TM and Landsat 7 ETM+ Land Surface Observations and Implications for Global 30m Landsat Data Product Generation

With the advent of the free U.S. Landsat data policy it is now feasible to consider the generation of global coverage 30 m Landsat data sets with temporal reporting frequency similar to that provided by the monthly Web Enabled Landsat (WELD) products. A statistical Landsat metadata analysis is reported considering more than 800,000 Landsat 5 TM and Landsat 7 ETM + acquisitions obtained from the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center archive. The global monthly probabilities of acquiring a cloud-free land surface observation for December 1998 to November 2001 (2000 epoch) and from December 2008 to November 2011 (2010 epoch) are reported to assess the availability of the Landsat data in the USGS Landsat archive for global multi-temporal land remote sensing applications. The global probabilities of acquiring a cloud-free land surface observation in each of three different seasons with the highest seasonal probabilities of cloud-free land surface observation are reported, considering one, two and three years of Landsat data, to assess the availability of Landsat data for global land cover mapping. The probabilities are derived considering Landsat 5 TM only, Landsat 7 ETM + only, and both sensors combined, to examine the relative benefits of using one or both Landsat sensors. The results demonstrate the utility of combing both Landsat 5 TM and Landsat 7 ETM + data streams to take advantage of their different acquisition patterns and to mitigate the deleterious impact of the Landsat 7 ETM + 2003 scan line failure. Sensor combination provided a greater global acquisition coverage with a 1.7% to 14.4% higher percentage of land locations acquired monthly compared to considering Landsat 7 ETM + data alone. The mean global monthly probability of a cloud-free land surface observation for the combined sensors was up to nearly 1.4 and 6.7 times greater than for ETM + and TM alone respectively. The probability of acquiring a cloud-free Landsat land surface observation in different seasons was greater when more years of data were considered and when both Landsat sensor data were combined. Considering combined sensors and 36 months of data, 86.4% and 84.2% of the global land locations had probabilities ≥ 0.95 for the 2000 and 2010 epochs respectively, with a global mean probability of 0.92 (σ 0.24) for the 2000 epoch and 0.90 (σ 0.28) for the 2010 epoch. These results indicate that 36 months of combined Landsat sensor data will provide sufficient land surface observations for 30 m global land cover mapping using a multi-temporal supervised classification scheme

Alternative Methods to Predict Actual Evapotranspiration Illustrate the Importance of Accounting for Phenology: The Event Driven Phenology Model Part II

Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET – 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET

A New Concept for Simulation of Vegetated Land Surface Dynamics: The Event Driven Phenology Model Part I

Phenologies of the vegetated land surface are being used increasingly for diagnosis and prognosis of climate change consequences. Current prospective and retrospective phenological models stand far apart in their approaches to the subject. We report on an exploratory attempt to implement a phenological model based on a new event driven concept which has both diagnostic and prognostic capabilities in the same modeling framework. This Event Driven Phenological Model (EDPM) is shown to simulate land surface phenologies and phenophase transition dates in agricultural landscapes based on assimilation of weather data and land surface observations from spaceborne sensors. The model enables growing season phenologies to develop in response to changing environmental conditions and disturbance events. It also has the ability to ingest remotely sensed data to adjust its output to improve representation of the modeled variable. We describe the model and report results of initial testing of the EDPM using Level 2 flux tower records from the Ameri- flux sites at Mead, Nebraska, USA, and at Bondville, Illinois, USA. Simulating the dynamics of normalized difference vegetation index based on flux tower data, the predictions by the EDPM show good agreement (RMSE0.8) for maize and soybean during several growing seasons at different locations. This study presents the EDPM used in the companion paper (Kovalskyy and Henebry, 2011) in a coupling scheme to estimate daily actual evapotranspiration over multiple growing seasons

Change and Persistence in Land Surface Phenologies of the Don and Dnieper River Basins

The formal collapse of the Soviet Union at the end of 1991 produced major socio-economic and institutional dislocations across the agricultural sector. The picture of broad scale patterns produced by these transformations continues to be discovered. We examine here the patterns of land surface phenology (LSP) within two key river basins—Don and Dnieper—using AVHRR (Advanced Very High Resolution Radiometer) data from 1982 to 2000 and MODIS (Moderate Resolution Imaging Spectroradiometer) data from 2001 to 2007. We report on the temporal persistence and change of LSPs as summarized by seasonal integration of NDVI (normalized difference vegetation index) time series using accumulated growing degree-days (GDDI NDVI). Three land cover super-classes—forest lands, agricultural lands, and shrub lands—constitute 96% of the land area within the basins. All three in both basins exhibit unidirectional increases in AVHRR GDDI NDVI between the Soviet and post-Soviet epochs. During the MODIS era (2001–2007), different socio-economic trajectories in Ukraine and Russia appear to have led to divergences in the LSPs of the agricultural lands in the two basins. Interannual variation in the shrub lands of the Don river basin has increased since 2000. This is due in part to the better signal-to-noise ratio of the MODIS sensor, but may also be due to a regional drought affecting the Don basin more than the Dnieper basin

Alternative methods to predict actual evapotranspiration illustrate the importance of accounting for phenology – Part 2: The event driven phenology model

Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET – 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET

Evaluation of the Event Driven Phenology Model Coupled with the VegET Evapotranspiration Model Through Comparisons with Reference Datasets in a Spatially Explicit Manner

A new model coupling scheme with remote sensing data assimilation was developed for estimation of daily actual evapotranspiration (ET). The scheme represents a mix of the VegET, a physically based model to estimate ET from a water balance, and an event driven phenology model (EDPM), where the EDPM is an empirically derived crop specific model capable of producing seasonal trajectories of canopy attributes. In this experiment, the scheme was deployed in a spatially explicit manner within the croplands of the Northern Great Plains. The evaluation was carried out using 2007-2009 land surface forcing data from the North American Land Data Assimilation System (NLDAS) and crop maps derived from remotely sensed data of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compared the canopy parameters produced by the phenology model with normalized difference vegetation index (NDVI) data derived from the MODIS nadir bi-directional reflectance distribution function (BRDF) adjusted reflectance (NBAR) product. The expectations of the EDPM performance in prognostic mode were met, producing determination coefficient (r2) of 0.8 +/-.0.15. Model estimates of NDVI yielded root mean square error (RMSE) of 0.1 +/-.0.035 for the entire study area. Retrospective correction of canopy dynamics with MODIS NDVI brought the errors down to just below 10% of observed data range. The ET estimates produced by the coupled scheme were compared with ones from the MODIS land product suite. The expected r2=0.7 +/-.15 and RMSE = 11.2 +/-.4 mm per 8 days were met and even exceeded by the coupling scheme0 functioning in both prognostic and retrospective modes. Minor setbacks of the EDPM and VegET performance (r2 about 0.5 and additional 30 % of RMSR) were found on the peripheries of the study area and attributed to the insufficient EDPM training and to spatially varying accuracy of crop maps. Overall the experiment provided sufficient evidence of soundness and robustness of the EDPM and VegET coupling scheme, assuring its potential for spatially explicit applications

Monitoring conterminous United States (CONUS) land cover change with Web-Enabled Landsat Data (WELD)

Forest cover loss and bare ground gain from 2006 to 2010 for the conterminous United States (CONUS) were quantified at a 30 m spatial resolution using Web-Enabled Landsat Data available from the USGS Center for Earth Resources Observation and Science (EROS) (http://landsat.usgs.gov/WELD.php). The approach related multi-temporal WELD metrics and expert-derived training data for forest cover loss and bare ground gain through a decision tree classification algorithm. Forest cover loss was reported at state and ecoregional scales, and the identification of core forests\u27 absent of change was made and verified using LiDAR data from the GLAS (Geoscience Laser Altimetry System) instrument. Bare ground gain correlated with population change for large metropolitan statistical areas (MSAs) outside of desert or semi-desert environments. Google Earth™ time series images were used to validate the products. Mapped forest cover loss totaled 53,084 km2 and was found to be depicted conservatively, with a user\u27s accuracy of 78% and a producer\u27s accuracy of 68%. Excluding errors of adjacency, user\u27s and producer\u27s accuracies rose to 93% and 89%, respectively. Mapped bare ground gain equaled 5974 km2 and nearly matched the estimated area from the reference (Google Earth™) classification; however, user\u27s (42%) and producer\u27s (49%) accuracies were much less than those of the forest cover loss product. Excluding errors of adjacency, user\u27s and producer\u27s accuracies rose to 62% and 75%, respectively. Compared to recent 2001–2006 USGS National Land Cover Database validation data for forest loss (82% and 30% for respective user\u27s and producer\u27s accuracies) and urban gain (72% and 18% for respective user\u27s and producer\u27s accuracies), results using a single CONUS-scale model with WELD data are promising and point to the potential for national scale operational mapping of key land cover transitions. However, validation results highlighted limitations, some of which can be addressed by improving training data, creating a more robust image feature space, adding contemporaneous Landsat 5 data to the inputs, and modifying definition sets to account for differences in temporal and spatial observational scales. The presented land cover extent and change data are available via the official WELD website (ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/ftp://weldftp.cr.usgs. gov/CONUS_5Y_LandCover/)

Characterization of Landsat-7 to Landsat-8 Reflective Wavelength and Normalized Difference Vegetation Index Continuity

At over 40 years, the Landsat satellites provide the longest temporal record of space-based land surface observations, and the successful 2013 launch of the Landsat-8 is continuing this legacy. Ideally, the Landsat data record should be consistent over the Landsat sensor series. The Landsat-8 Operational Land Imager (OLI) has improved calibration, signal to noise characteristics, higher 12-bit radiometric resolution, and spectrally narrower wavebands than the previous Landsat-7 Enhanced Thematic Mapper (ETM +). Reflective wavelength differences between the two Landsat sensors depend also on the surface reflectance and atmospheric state which are difficult to model comprehensively. The orbit and sensing geometries of the Landsat-8 OLI and Landsat-7 ETM + provide swath edge overlapping paths sensed only one day apart. The overlap regions are sensed in alternating backscatter and forward scattering orientations so Landsat bi-directional reflectance effects are evident but approximately balanced between the two sensors when large amounts of time series data are considered. Taking advantage of this configuration a total of 59 million 30 m corresponding sensor observations extracted from 6317 Landsat-7 ETM + and Landsat-8 OLI images acquired over three winter and three summer months for all the conterminous United States (CONUS) are compared. Results considering different stages of cloud and saturation filtering, and filtering to reduce one day surface state differences, demonstrate the importance of appropriate per-pixel data screening. Top of atmosphere (TOA) and atmospherically corrected surface reflectance for the spectrally corresponding visible, near infrared and shortwave infrared bands, and derived normalized difference vegetation index (NDVI), are compared and their differences quantified. On average the OLI TOA reflectance is greater than the ETM + TOA reflectance for all bands, with greatest differences in the near-infrared (NIR) and the shortwave infrared bands due to the quite different spectral response functions between the sensors. The atmospheric correction reduces the mean difference in the NIR and shortwave infrared but increases the mean difference in the visible bands. Regardless of whether TOA or surface reflectance are used to generate NDVI, on average, for vegetated soil and vegetation surfaces (0 ≤ NDVI ≤ 1), the OLI NDVI is greater than the ETM + NDVI. Statistical functions to transform between the comparable sensor bands and sensor NDVI values are presented so that the user community may apply them in their own research to improve temporal continuity between the Landsat-7 ETM + and Landsat-8 OLI sensor data. The transformation functions were developed using ordinary least squares (OLS) regression and were fit quite reliably (r2 values \u3e 0.7 for the reflectance data and \u3e 0.9 for the NDVI data, p-values \u3c 0.0001)

A General Method to Normalize Landsat Reflectance Data to Nadir BRDF Adjusted Reflectance

The Landsat satellites have been providing spectacular imagery of the Earth\u27s surface for over 40 years. However, they acquire images at view angles ±7.5° from nadir that cause small directional effects in the surface reflectance. There are also variations with solar zenith angle over the year that can cause apparent change in reflectance even if the surface properties remain constant. When Landsat data from adjoining paths, or from long time series are used, a model of the surface anisotropy is required to adjust all Landsat observations to a uniform nadir view (primarily for visual consistency, vegetation monitoring, or detection of subtle surface changes). Here a generalized approach is developed to provide consistent view angle corrections across the Landsat archive. While this approach is not applicable for generation of Landsat surface albedo, which requires a full characterization of the surface bidirectional reflectance distribution function (BRDF), or for correction to a constant solar illumination angle across a wide range of sun angles, it provides Landsat nadir BRDF-adjusted reflectance (NBAR) for a range of terrestrial monitoring applications. The Landsat NBAR is derived as the product of the observed Landsat reflectance and the ratio of the reflectances modeled using MODIS BRDF spectral model parameters for the observed Landsat and for a nadir view and fixed solar zenith geometry. In this study, a total of 567 conterminous United States (CONUS) January and July 2010 Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM+) images that have swath edge overlapping paths sensed in alternating backscatter and forward scattering orientations were used. The average difference between Landsat 5 TM and Landsat 7 ETM+ surface reflectance in the forward and backward scatter directions at the overlapping Landsat scan edges was quantified. The CONUS July view zenith BRDF effects were about 0.02 in the Landsat visible bands, and about 0.03, 0.05 and 0.06, in the 2.1 μm, 1.6 μm and nearinfrared bands respectively. Comparisons of Landsat 5 TM and Landsat 7 ETM+ NBAR derived using MODIS BRDF spectral model parameters defined with respect to different spatial and temporal scales, and defined with respect to different land cover types, were undertaken. The results suggest that, because the BRDF shapes of different terrestrial surfaces are sufficiently similar over the narrow 15° Landsat field of view, a fixed set of MODIS BRDF spectral model parameters may be adequate for Landsat NBAR derivation with little sensitivity to the land cover type, condition, or surface disturbance. A fixed set of BRDF spectral model parameters, derived from a global year of highest quality snow-free MODIS BRDF product values, are provided so users may implement the described Landsat NBAR generation method