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

Assessing the Calibration Differences in the Reflective Solar Bands of Terra MODIS and Landsat-7 Enhanced Thematic Mapper Plus

Long-term data records obtained from Earth observing sensors depend not only onthe calibration accuracy of individual sensors but also on the consistency across instruments andplatforms. Hence, sensor calibration intercomparison plays a vital role for a better understandingof various science products. The Moderate Resolution Imaging Spectroradiometer (MODIS)and enhanced thematic mapper plus (ETM+) on the Terra and Landsat 7 platforms have operatedsuccessfully since their launch, collecting measurements in the reflective solar and infrared partsof the spectrum. Terra MODIS has employed a reflectance-based calibration since beginning itsmission. In the case of ETM+, a radiance-based calibration was employed until recent years,when a reflectance-based calibration was introduced. Being in the AM constellation with lessthan 30 min difference in overpass times, near-simultaneous Earth scene measurements can beeffectively used to assess the calibration differences between the spectrally matching bands ofthese two instruments. The pseudoinvariant calibration sites (PICS) in the North African desertare widely used for on-orbit calibration and validation of satellite sensors. Four PICS from thisregion have been employed to assess the multitemporal reflectance differences. Correction forbidirectional reflectance, spectral response function mismatch, and impacts of atmosphericwater-vapor have been incorporated to provide an assessment of the long-term stability ofeach spectral band and reflectance differences amongst them. Results indicate that the spectralbands of both instruments show a long-term stability to within 2% from 2000 to 2017. Thetop-of-atmosphere reflectances between the two instruments postcorrection agree to within 4%.Also included in this paper is a detailed discussion of various parameters contributing to theuncertainties of this cross-calibration. The techniques presented in this paper can be furtherextended to perform similar intercomparison between Landsat 8 Operational Land Imager, AquaMODIS, and Suomi-NPP VIIRS

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

Absolute Radiometric Calibration of Narrow-Swath Imaging Sensors with Reference to Non-Coincident Wide-Swath Sensors

An inter-calibration method is developed to provide absolute radiometric calibration of narrow-swath imaging sensors with reference to non-coincident wide-swath sensors. The method predicts at-sensor radiance using non-coincident imagery from the reference sensor and knowledge of spectral reflectance of the test site. The imagery of the reference sensor is restricted to acquisitions that provide similar view and solar illumination geometry to reduce uncertainties due to directional reflectance effects. Spectral reflectance of the test site is found with a simple iterative radiative transfer method using radiance values of a well-understood wide-swath sensor and spectral shape information based on historical ground-based measurements. At-sensor radiance is calculated for the narrow-swath sensor using this spectral reflectance and atmospheric parameters that are also based on historical in situ measurements. Results of the inter-calibration method show agreement on the 2 5 percent level in most spectral regions with the vicarious calibration technique relying on coincident ground-based measurements referred to as the reflectance-based approach. While the variability of the inter-calibration method based on non-coincident image pairs is significantly larger, results are consistent with techniques relying on in situ measurements. The method is also insensitive to spectral differences between the sensors by transferring to surface spectral reflectance prior to prediction of at-sensor radiance. The utility of this inter-calibration method is made clear by its flexibility to utilize image pairings with acquisition dates differing in excess of 30 days allowing frequent absolute calibration comparisons between wide- and narrow-swath sensors

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

Overview of Intercalibration of Satellite Instruments

Intercalibration of satellite instruments is critical for detection and quantification of changes in the Earth’s environment, weather forecasting, understanding climate processes, and monitoring climate and land cover change. These applications use data from many satellites; for the data to be interoperable, the instruments must be cross-calibrated. To meet the stringent needs of such applications, instruments must provide reliable, accurate, and consistent measurements over time. Robust techniques are required to ensure that observations from different instruments can be normalized to a common scale that the community agrees on. The long-term reliability of this process needs to be sustained in accordance with established reference standards and best practices. Furthermore, establishing physical meaning to the information through robust Système International d’unités traceable calibration and validation (Cal/Val) is essential to fully understand the parameters under observation. The processes of calibration, correction, stabilitymonitoring, and quality assurance need to be underpinned and evidenced by comparison with “peer instruments” and, ideally, highly calibrated in-orbit reference instruments. Intercalibration between instruments is a central pillar of the Cal/Val strategies of many national and international satellite remote sensing organizations. Intercalibration techniques as outlined in this paper not only provide a practical means of identifying and correcting relative biases in radiometric calibration between instruments but also enable potential data gaps between measurement records in a critical time series to be bridged. Use of a robust set of internationally agreed upon and coordinated intercalibration techniques will lead to significant improvement in the consistency between satellite instruments and facilitate accurate monitoring of the Earth’s climate at uncertainty levels needed to detect and attribute the mechanisms of change. This paper summarizes the state-of-the-art of postlaunch radiometric calibration of remote sensing satellite instruments through intercalibration

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

Potential of VIIRS Time Series Data for Aiding the USDA Forest Service Early Warning System for Forest Health Threats: A Gypsy Moth Defoliation Case Study

This report details one of three experiments performed during FY 2007 for the NASA RPC (Rapid Prototyping Capability) at Stennis Space Center. This RPC experiment assesses the potential of VIIRS (Visible/Infrared Imager/Radiometer Suite) and MODIS (Moderate Resolution Imaging Spectroradiometer) data for detecting and monitoring forest defoliation from the non-native Eurasian gypsy moth (Lymantria dispar). The intent of the RPC experiment was to assess the degree to which VIIRS data can provide forest disturbance monitoring information as an input to a forest threat EWS (Early Warning System) as compared to the level of information that can be obtained from MODIS data. The USDA Forest Service (USFS) plans to use MODIS products for generating broad-scaled, regional monitoring products as input to an EWS for forest health threat assessment. NASA SSC is helping the USFS to evaluate and integrate currently available satellite remote sensing technologies and data products for the EWS, including the use of MODIS products for regional monitoring of forest disturbance. Gypsy moth defoliation of the mid-Appalachian highland region was selected as a case study. Gypsy moth is one of eight major forest insect threats listed in the Healthy Forest Restoration Act (HFRA) of 2003; the gypsy moth threatens eastern U.S. hardwood forests, which are also a concern highlighted in the HFRA of 2003. This region was selected for the project because extensive gypsy moth defoliation occurred there over multiple years during the MODIS operational period. This RPC experiment is relevant to several nationally important mapping applications, including agricultural efficiency, coastal management, ecological forecasting, disaster management, and carbon management. In this experiment, MODIS data and VIIRS data simulated from MODIS were assessed for their ability to contribute broad, regional geospatial information on gypsy moth defoliation. Landsat and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data were used to assess the quality of gypsy moth defoliation mapping products derived from MODIS data and from simulated VIIRS data. The project focused on use of data from MODIS Terra as opposed to MODIS Aqua mainly because only MODIS Terra data was collected during 2000 and 2001-years with comparatively high amounts of gypsy moth defoliation within the study area. The project assessed the quality of VIIRS data simulation products. Hyperion data was employed to assess the quality of MODIS-based VIIRS simulation datasets using image correlation analysis techniques. The ART (Application Research Toolbox) software was used for data simulation. Correlation analysis between MODIS-simulated VIIRS data and Hyperion-simulated VIIRS data for red, NIR (near-infrared), and NDVI (Normalized Difference Vegetation Index) image data products collectively indicate that useful, effective VIIRS simulations can be produced using Hyperion and MODIS data sources. The r(exp 2) for red, NIR, and NDVI products were 0.56, 0.63, and 0.62, respectively, indicating a moderately high correlation between the 2 data sources. Temporal decorrelation from different data acquisition times and image misregistration may have lowered correlation results. The RPC experiment also generated MODIS-based time series data products using the TSPT (Time Series Product Tool) software. Time series of simulated VIIRS NDVI products were produced at approximately 400-meter resolution GSD (Ground Sampling Distance) at nadir for comparison to MODIS NDVI products at either 250- or 500-meter GSD. The project also computed MODIS (MOD02) NDMI (Normalized Difference Moisture Index) products at 500-meter GSD for comparison to NDVI-based products. For each year during 2000-2006, MODIS and VIIRS (simulated from MOD02) time series were computed during the peak gypsy moth defoliation time frame in the study area (approximately June 10 through July 27). Gypsy moth defoliation mapping products from simated VIIRS and MOD02 time series were produced using multiple methods, including image classification and change detection via image differencing. The latter enabled an automated defoliation detection product computed using percent change in maximum NDVI for a peak defoliation period during 2001 compared to maximum NDVI across the entire 2000-2006 time frame. Final gypsy moth defoliation mapping products were assessed for accuracy using randomly sampled locations found on available geospatial reference data (Landsat and ASTER data in conjunction with defoliation map data from the USFS). Extensive gypsy moth defoliation patches were evident on screen displays of multitemporal color composites derived from MODIS data and from simulated VIIRS vegetation index data. Such defoliation was particularly evident for 2001, although widespread denuded forests were also seen for 2000 and 2003. These visualizations were validated using aforementioned reference data. Defoliation patches were visible on displays of MODIS-based NDVI and NDMI data. The viewing of apparent defoliation patches on all of these products necessitated adoption of a specialized temporal data processing method (e.g., maximum NDVI during the peak defoliation time frame). The frequency of cloud cover necessitated this approach. Multitemporal simulated VIIRS and MODIS Terra data both produced effective general classifications of defoliated forest versus other land cover. For 2001, the MOD02-simulated VIIRS 400-meter NDVI classification produced a similar yet slightly lower overall accuracy (87.28 percent with 0.72 Kappa) than the MOD02 250-meter NDVI classification (88.44 percent with 0.75 Kappa). The MOD13 250-meter NDVI classification had a lower overall accuracy (79.13 percent) and a much lower Kappa (0.46). The report discusses accuracy assessment results in much more detail, comparing overall classification and individual class accuracy statistics for simulated VIIRS 400-meter NDVI, MOD02 250-meter NDVI, MOD02-500 meter NDVI, MOD13 250-meter NDVI, and MOD02 500-meter NDMI classifications. Automated defoliation detection products from simulated VIIRS and MOD02 data for 2001 also yielded similar, relatively high overall classification accuracy (85.55 percent for the VIIRS 400-meter NDVI versus 87.28 percent for the MOD02 250-meter NDVI). In contrast, the USFS aerial sketch map of gypsy moth defoliation showed a lower overall classification accuracy at 73.64 percent. The overall classification Kappa values were also similar for the VIIRS (approximately 0.67 Kappa) versus the MOD02 (approximately 0.72 Kappa) automated defoliation detection product, which were much higher than the values exhibited by the USFS sketch map product (overall Kappa of approximately 0.47). The report provides additional details on the accuracy of automated gypsy moth defoliation detection products compared with USFS sketch maps. The results suggest that VIIRS data can be effectively simulated from MODIS data and that VIIRS data will produce gypsy moth defoliation mapping products that are similar to MODIS-based products. The results of the RPC experiment indicate that VIIRS and MODIS data products have good potential for integration into the forest threat EWS. The accuracy assessment was performed only for 2001 because of time constraints and a relative scarcity of cloud-free Landsat and ASTER data for the peak defoliation period of the other years in the 2000-2006 time series. Additional work should be performed to assess the accuracy of gypsy moth defoliation detection products for additional years.The study area (mid-Appalachian highlands) and application (gypsy moth forest defoliation) are not necessarily representative of all forested regions and of all forest threat disturbance agents. Additional work should be performed on other inland and coastal regions as well as for other major forest threats

Refinement of the method for using pseudo-invariant sites for long term calibration trending of Landsat reflective bands

The long term calibration history of the Landsat 5 TM instrument has recently been defined using a time series of desert sites in Northern Africa. This correction is based on the assumption that the atmosphere is invariant and the reflectance of each site is approximately constant and Lambertian over time. As a result, the top of the atmosphere reflection is assumed constant when corrected for variations in the solar elevation angle and earth-sun distance. While this is true to first order and is the basis for all current temporal calibration, there are multiple known sources of residual error in the data. A methodology is presented for reducing the variation in pseudo-invariant site trending data based on correction for the BRDF. This work establishes a means to use DIRSIG to model the L5 calibration site. It combines a digital elevation map and desert atmosphere with a surface BRDF to reduce the residual errors in the calibration data. A set of Landsat 7 ETM+ calibration days is utilized to optimize the surface reflectance properties used in DIRSIG. These optimized parameters are then used to model the L5 TM calibration days. The results of the DIRSIG modeling are compared to the solar elevation angle and time of year trends of the original data and analyzed for their effectiveness at describing and reducing the residual errors. A major goal of this effort is to understand the contribution that BRDFs make to the current calibration errors and to develop methods that are robust enough to be applicable to a wider range of sites to enable extension of the methodology to earlier data sets (e.g. Landsat MSS). Additionally, while Landsat has a 30 m reflective resolution, the pseudo-invariant site calibration approach is valid for all spatial resolutions. Depending on another instrument\u27s field of view, the BRDF error reduction technique used by L5 TM could either be used on the same desert calibration site or on a subsection of the area

Remote Sensing of Environment: Current status of Landsat program, science, and applications

Formal planning and development of what became the first Landsat satellite commenced over 50 years ago in 1967. Now, having collected earth observation data for well over four decades since the 1972 launch of Landsat- 1, the Landsat program is increasingly complex and vibrant. Critical programmatic elements are ensuring the continuity of high quality measurements for scientific and operational investigations, including ground systems, acquisition planning, data archiving and management, and provision of analysis ready data products. Free and open access to archival and new imagery has resulted in a myriad of innovative applications and novel scientific insights. The planning of future compatible satellites in the Landsat series, which maintain continuity while incorporating technological advancements, has resulted in an increased operational use of Landsat data. Governments and international agencies, among others, can now build an expectation of Landsat data into a given operational data stream. International programs and conventions (e.g., deforestation monitoring, climate change mitigation) are empowered by access to systematically collected and calibrated data with expected future continuity further contributing to the existing multi-decadal record. The increased breadth and depth of Landsat science and applications have accelerated following the launch of Landsat-8, with significant improvements in data quality. Herein, we describe the programmatic developments and institutional context for the Landsat program and the unique ability of Landsat to meet the needs of national and international programs. We then present the key trends in Landsat science that underpin many of the recent scientific and application developments and followup with more detailed thematically organized summaries. The historical context offered by archival imagery combined with new imagery allows for the development of time series algorithms that can produce information on trends and dynamics. Landsat-8 has figured prominently in these recent developments, as has the improved understanding and calibration of historical data. Following the communication of the state of Landsat science, an outlook for future launches and envisioned programmatic developments are presented. Increased linkages between satellite programs are also made possible through an expectation of future mission continuity, such as developing a virtual constellation with Sentinel-2. Successful science and applications developments create a positive feedback loop—justifying and encouraging current and future programmatic support for Landsat