Forty-Year Calibrated Record of Earth-Surface Reflected Radiance from Landsat: A Review

Sensors on Landsat satellites have been collecting images of the Earth's surface for nearly 40 years. These images have been invaluable for characterizing and detecting changes in the land cover and land use of the world. Although initially conceived as primarily picture generating sensors, even the early sensors were radiometrically calibrated and spectrally characterized prior to launch and incorporated some capabilities to monitor their radiometric calibration once on orbit. Recently, as the focus of studies has shifted to monitoring Earth surface parameters over significant periods of time, serious attention has been focused toward bringing the data from all these sensors onto a common radiometric scale over this 40-year period. This effort started with the most recent systems and then was extended back in time. Landsat-7 ETM+, the best-characterized sensor of the series prior to launch and once on orbit, and the most stable system to date, was chosen to serve as the reference. The Landsat-7 project was the first of the series to build an image assessment system into its ground system, allowing systematic characterization of its sensors and data. Second, the Landsat-5 TM (still operating at the time of the Landsat-7 launch and continues to operate) calibration history was reconstructed based on its internal calibrator, vicarious calibrations, pseudo-invariant sites and a tie to Landsat-7 ETM+ at the time of the commissioning of Landsat-7. This process was performed in two iterations: the earlier one relied primarily on the TM internal calibrator. When this was found to have some deficiencies, a revised calibration was based more on pseudo-invariant sites, though the internal calibrator was still used to establish the short-term variations in response due to icing build up on the cold focal plane. As time progressed, a capability to monitor the Landsat-5 TM was added to the image assessment system. The Landsat-4 TM, which operated from 1982-1992, was the third system to which the radiometric scale was extended. The limited and broken use of the Landsat-4 TM made this analysis more difficult. Eight-day separated image pairs from Landsat-5 combined with analysis of pseudo invariant sites established this history. The fourth and most challenging effort was making the Landsat-1 to -5 MSS sensors' data internally radiometrically consistent. This effort was particularly complicated by the age of the MSS data, varying formats and processing levels in the archive, limited datasets, and limited documentation available. Ultimately, pseudo-invariant sites were identified in North America and used for this effort. Note that most of the Landsat-MSS archived data had already been calibrated using the MSS internal calibrators, so this processing was imbedded in the result. The final effort was developing an absolute scale for Landsat MSS similar to what was already established for the "TM" sensors. This was accomplished by using simultaneous data from Landsat-5 MSS and Landsat-5 TM, accounting for spectral differences between the sensors using EO-1 Hyperion data. The recalibrated history of the Landsat data and implications to users are discussed. The key result from this work is a consistently calibrated Landsat data archive that spans nearly 40 years with total uncertainties on the order of 10% or less for most sensors and bands

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a Virtual Constellation was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating stable calibration to within 5%the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Students T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process

Worldwide Optimal PICS Search

Pseudo Invariant Calibration Sites (PICS) have proven to be a dependable calibration source for determining degradation of visible and infrared sensor response due to their temporal stability and spatial uniformity. One limit of PICS is that only a handful have been identified, primarily in desert areas of North Africa, Saudi Arabia, and elsewhere. A large number of PICS would not only facilitate calibration of existing and future sensors, but also provide an alternative to internal on-board calibrator data, resulting in significant cost savings and simplification in sensor design. As a result, the process to efficiently identify additional PICS is highly desirable. A relatively straightforward algorithm and processing flow to identify candidate PICS throughout the world has been developed. One goal of the algorithm is to identify PICS with reflectance levels covering more of the sensor dynamic range. As currently implemented, the algorithm makes use of Google Earth Engine to simplify the required image data pre-processing, analysis, and storage, and implements a filtering technique to enhance contiguity of pixels identified as invariant. Application of the proposed algorithm identified not only existing North Africa and Middle East sites with 2% to 2.5% temporal uncertainty, but also sites on other continents with 5% to 6% uncertainty, which can be improved with application of BRDF correction. In general, the algorithm shows potential in providing a means for automated PICS identification

Evaluation of an Extended PICS (EPICS) for Calibration and Stability Monitoring of Optical Satellite Sensors

Pseudo Invariant Calibration Sites (PICS) have been increasingly used as an independent data source for on-orbit radiometric calibration and stability monitoring of optical satellite sensors. Generally, this would be a small region of land that is extremely stable in time and space, predominantly found in North Africa. Use of these small regions, referred to as traditional PICS, can be limited by: i) the spatial extent of an individual Region of Interest (ROI) and/or site; ii) and the frequency of how often the site can be acquired, based on orbital patterns and cloud cover at the site, both impacting the time required to construct a richly populated temporal dataset. This paper uses a new class of continental scaled PICS clusters (also known as Extended PICS or EPICS), to demonstrate their capability in increasing temporal frequency of the calibration time series which ultimately allows calibration and stability assessment at a much finer scale compared to the traditional PICSbased method while also reducing any single location’s potential impact to the overall assessment. The use of EPICS as a calibration site was evaluated using data from Landsat- 8 Operational Land Imager (OLI), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), and Sentinel-2A&B Multispectral Instrument (MSI) images at their full spatial resolutions. Initial analysis suggests that EPICS, at its full potential and with nominal cloud consideration, can significantly decrease the temporal revisit interval of moderate resolution sensors to as much as of 0.33 day (3 collects/day). A traditional PICS is expected to have a temporal uncertainty (defined as the ratio of temporal standard deviation and temporal mean) of 2-5% for TOA reflectance. Over the same time period EPICS produced a temporal uncertainty of 3%. But the advantage to be leveraged is the ability to detect sensor change quicker due to the denser dataset and reduce the impact of any potential ‘local’ changes. Moreover, this approach can be extended to any on-orbit sensor. An initial attempt to quantify the minimum detectable change (a threshold slope value which must be exceeded by the reflectance trend to be considered statistically significant) suggests that the use of EPICS can decrease the time period up to approximately half of that found using traditional PICS-based approach

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a “Virtual Constellation” was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating “stable” calibration to within 5%—the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)–were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Student’s T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process

Reflectance-based Calibration and Validation of the Landsat Satellite Archive

The primary objective of this project was to consistently calibrate the entire Landsat series to a common reflectance scale by performing cross-calibration corrections from Landsat-8 OLI to Landsat- 1 MSS. A consistent radiance-based calibration was already performed from Landsat-8 OLI through Landsat-1 MSS using bright targets and dark targets. The MSS radiance-based calibration results showed an uncertainty of about ±5%. Typically to convert from radiance to reflectance a solar model is used. Unfortunately, there are numerous solar models, all with various levels of accuracies. It was also seen that there is a data format inconsistency for different types of MSS data that impact the radiometric uncertainty of the products when compared to Landsat-8 OLI data. One of the advances Landsat-8 OLI has over to earlier missions is a solar model independent reflectance calibration. Hence, to reduce these uncertainties and remove the dependency on the solar model, direct reflectance-based calibration was performed for all previous missions using Landsat-8 OLI as the “standard”. A consistent cross-calibration of all Landsat sensors was achieved using coincident/near-coincident scene pairs. The work started from cross-calibration of Landsat-8 OLI to Landsat-7 ETM+ and continued through Landsat-1 MSS. Due to the fact each Landsat sensor measures slightly different parts of the electromagnetic spectrum, a spectral band adjustment factor (SBAF) was computed and used prior to the cross-calibration. To determine the significance of the bias derived from cross-calibration, a t-test was performed with a null hypothesis that the bias equals zero at a confidence interval of 95%. From the final calibration equations, it was found that for band 5 of Landsat-1 bias is significant. The effectiveness of these cross-calibration results is discussed by showing a significant improvement in the observed inconsistencies in the absolute calibration of all Landsat sensors for both bright and dark targets. The results show a significant improvement in reflectance calibration, and an overall uncertainty of less than ±3%

Multitemporal Cross-Calibration of the Terra MODIS and Landsat 7 ETM+ Reflective Solar Bands

In recent years, there has been a significant increase in the use of remotely sensed data to address global issues. With the open data policy, the data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Enhanced Thematic Mapper Plus (ETM+) sensors have become a critical component of numerous applications. These two sensors have been operational for more than a decade, providing a rich archive of multispectral imagery for analysis of mutitemporal remote sensing data. This paper focuses on evaluating the radiometric calibration agreement between MODIS and ETM+ using the near-simultaneous and cloud-free image pairs over an African pseudo-invariant calibration site, Libya 4. To account for the combined uncertainties in the top-of-atmosphere (TOA) reflectance due to surface and atmospheric bidirectional reflectance distribution function (BRDF), a semiempirical BRDF model was adopted to normalize the TOA reflectance to the same illumination and viewing geometry. In addition, the spectra from the Earth Observing-1 (EO-1) Hyperion were used to compute spectral corrections between the corresponding MODIS and ETM+ spectral bands. As EO-1 Hyperion scenes were not available for all MODIS and ETM+ data pairs, MODerate resolution atmospheric TRANsmission (MODTRAN) 5.0 simulations were also used to adjust for differences due to the presence or lack of absorption features in some of the bands. A MODIS split-window algorithm provides the atmospheric water vapor column abundance during the overpasses for the MODTRAN simulations. Additionally, the column atmospheric water vapor content during the overpass was retrieved using the MODIS precipitable water vapor product. After performing these adjustments, the radiometric cross-calibration of the two sensors was consistent to within 7%. Some drifts in the response of the bands are evident, with MODIS band 3 being the largest of about 6% over 10 years, a change that will be corrected in Collection 6 MODIS processing

Classification of North Africa for Use as an Extended Pseudo Invariant Calibration Sites (Epics) for Radiometric Calibration and Stability Monitoring of Optical Satellite Sensors

An increasing number of Earth-observing satellite sensors are being launched to meet the insatiable demand for timely and accurate data to help the understanding of the Earth’s complex systems and to monitor significant changes to them. The quality of data recorded by these sensors is a primary concern, as it critically depends on accurate radiometric calibration for each sensor. Pseudo Invariant Calibration Sites (PICS) have been extensively used for radiometric calibration and temporal stability monitoring of optical satellite sensors. Due to limited knowledge about the radiometric stability of North Africa, only a limited number of sites in the region are used for this purpose. This work presents an automated approach to classify North Africa for its potential use as an extended PICS (EPICS) covering vast portions of the continent. An unsupervised classification algorithm identified 19 “clusters” representing distinct land surface types; three clusters were identified with spatial uncertainties within approximately 5% in the shorter wavelength bands and 3% in the longer wavelength bands. A key advantage of the cluster approach is that large numbers of pixels are aggregated into contiguous homogeneous regions sufficiently distributed across the continent to allow multiple imaging opportunities per day, as opposed to imaging a typical PICS once during the sensor’s revisit period. In addition, this work proposes a technique to generate a representative hyperspectral profile for these clusters, as the hyperspectral profile of these identified clusters are mandatory in order to utilize them for performing cross-calibration of optical satellite sensors. The technique was used to generate the profile for the cluster containing the largest number of aggregated pixels. The resulting profile was found to have temporal uncertainties within 5% across all the spectral regions. Overall, this technique shows great potential for generation of representative hyperspectral profiles for any North African cluster, which could allow the use of the entire North Africa Saharan region as an extended PICS (EPICS) dataset for sensor cross-calibration. Furthermore, this work investigates the performance of extended pseudo-invariant calibration sites (EPICS) in cross-calibration for one of Shrestha’s clusters, Cluster 13, by comparing its results to those obtained from a traditional PICS-based cross-calibration. The use of EPICS clusters can significantly increase the number of cross-calibration opportunities within a much shorter time period. The cross-calibration gain ratio estimated using a cluster-based approach had a similar accuracy to the cross-calibration gain derived from region of interest (ROI)-based approaches. The cluster-based cross-calibration gain ratio is consistent within approximately 2% of the ROI-based cross-calibration gain ratio for all bands except for the coastal and shortwave-infrared (SWIR) 2 bands. These results show that image data from any region within Cluster 13 can be used for sensor crosscalibration. Eventually, North Africa can be used a continental scale PICS

Normalization of Pseudo-invariant Calibration Sites for Increasing the Temporal Resolution and Long-Term Trending

Given their low level of temporal, spatial, and spectral variability, Pseudo-Invariant Calibration Sites (PICS) have been increasingly desired as data sources for radiometric calibration of Earth imaging satellite sensors. The temporal resolution for PICS data acquired by any sensor is limited by the amount of time required for it to make subsequent passes over the site. Consequently, for any given PICS, it can take many years of imaging to develop a sufficient amount of cloud-free data to perform radiometric calibration; this can be especially problematic for sensors in their early years after launch. This thesis presents techniques to combine Landsat-8; normally acquiring data for every 16 days, image data from multiple PICS into a single dataset with increased temporal resolution and is called “PICS Normalization Process” or PNP. Landsat-8 Operational Land Imager (OLI) data from six Saharan desert sites were normalized to the Libya-4 reference. The normalized data were then merged into a “Super PICS” dataset, and the estimation of calibration drift was derived. The results of the Super PICS dataset show that the temporal resolution of the calibration dataset can be increased by approximately a factor of three to four times. The normalization process was performed on radiometrically and geometrically corrected image data (“L1T” product), and also on the same image data corrected for BRDF effects using a quadratic function of the solar zenith angle and TOA reflectance over a region of interest. An additional uncertainty analysis was performed using the BRDF corrected image data based on the following parameters which are involved in this whole BRDF PICS Normalization Process: Worst-case histogram bin analysis, Temporal Uncertainty of each PICS, BRDF Super PICS uncertainty. The resulting uncertainties are within the currently accepted satellite calibration range, within 3% for all spectral bands. Overall, the process indicates a calibration drift for OLI within 0.15% per year, agreeing quite well with the calibration drift derived from the on-board calibrators

Vicarious Methodologies to Assess and Improve the Quality of the Optical Remote Sensing Images: A Critical Review

Over the past decade, number of optical Earth observing satellites performing remote sensing has increased substantially, dramatically increasing the capability to monitor the Earth. The quantity of remote sensing satellite increase is primarily driven by improved technology, miniaturization of components, reduced manufacturing, and launch cost. These satellites often lack on-board calibrators that a large satellite utilizes to ensure high quality (e.g., radiometric, geometric, spatial quality, etc.) scientific measurement. To address this issue, this work presents “best” vicarious image quality assessment and improvement techniques for those kinds of optical satellites which lacks on-board calibration system. In this article, image quality categories have been explored, and essential quality parameters (e.g., absolute and relative calibration, aliasing, etc.) have been identified. For each of the parameters, appropriate characterization methods are identified along with its specifications or requirements. In cases of multiple methods, recommendation has been made based-on the strengths and weaknesses of each method. Furthermore, processing steps have been presented, including examples. Essentially, this paper provides a comprehensive study of the criteria that needs to be assessed to evaluate remote sensing satellite data quality, and best vicarious methodologies to evaluate identified quality parameters such as coherent noise, ground sample distance, etc