Empirical cross sensor comparison of Sentinel-2A and 2B MSI, Landsat-8 OLI, and Landsat-7 ETM+ top of atmosphere spectral characteristics over the conterminous United States

Remote sensing landscape monitoring approaches frequently benefit from a dense time series of observations. To enhance these time series, data from multiple satellite systems need to be integrated. Landsat image data is a valuable 30-meter resolution source of spatial information to assess forest conditions over time. Together both operational Landsat satellites—7 and 8—provide a revisit frequency of 8 days at the equator. This moderate temporal frequency provides essential information to detect annual large area abrupt land cover changes. However, the ability to measure subtle and short lived intraseasonal changes is challenged by gaps in Landsat imagery at key points in time. The first Sentinel-2 satellite mission was launched by the European Space Agency in 2015. This moderate resolution data stream provides an opportunity to supplement the Landsat data record. The objective of this study is to assess the potential for integrating top of atmosphere Landsat and Sentinel 2 image data archived in the Google Earth Engine compute environment. In this paper we assess absolute and proportional differences in near-contemporaneous observations for six bands with comparable spectral response functions and spatial resolution between the Sentinel-2 Multi Spectral Instrument and Landsat Operational Land Imager and Enhanced Thematic Mapper Plus imagery. We assessed differences using absolute difference metrics and major axis linear regression between over 10,000 image pairs across the conterminous United States and present cross sensor transformation models. Major axis linear regression results indicate that Sentinel MSI data are as spectrally comparable to the two types of Landsat image data as the Landsat sensors are with each other. Root-mean-square deviation (RMSD) values ranging from 0.0121 to 0.0398 were obtained between MSI and Landsat spectral values, and RMSD values ranging from 0.0124 and 0.0372 were obtained between OLI and ETM+. Despite differences in their spatial, spectral, and temporal characteristics, integration of these datasets appears to be feasible through the application of bandwise linear regression corrections.