A simplified and robust surface reflectance estimation method (SREM) for use over diverse land surfaces using multi-sensor data

Surface reflectance (SR) estimation is the most critical preprocessing step for deriving geophysical parameters in multi-sensor remote sensing. Most state-of-the-art SR estimation methods, such as the vector version of the Second Simulation of the Satellite Signal in the Solar Spectrum (6SV) radiative transfer (RT) model, depend on accurate information on aerosol and atmospheric gases. In this study, a Simplified and Robust Surface Reflectance Estimation Method (SREM) based on the equations from 6SV RT model, without integrating information of aerosol particles and atmospheric gasses, is proposed and tested using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper plus (ETM+), and Landsat 8 Operational Land Imager (OLI) data from 2000 to 2018. For evaluation purposes, (i) the SREM SR retrievals are validated against in situ SR measurements collected by Analytical Spectral Devices (ASD) from the South Dakota State University (SDSU) site, USA; (ii) cross-comparison between the SREM and Landsat spectral SR products, i.e., Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Landsat 8 Surface Reflectance Code (LaSRC), are conducted over 11 urban (2013-2018), 13 vegetated (2013-2018), and 11 desert/arid (2000 to 2018) sites located over different climatic zones at a global scale; (iii) the performance of the SREM spectral SR retrievals for low to high aerosol loadings is evaluated; (iv) spatio-temporal cross-comparison is conducted for six Landsat paths/rows located in Asia, Africa, Europe, and the United States of America from 2013 to 2018 to consider a large variety of land surfaces and atmospheric conditions; (v) cross-comparison is also performed for the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), and the Soil Adjusted Vegetation Index (SAVI) calculated from both the SREM and Landsat SR data; (vi) the SREM is also applied to the Sentinel-2A and Moderate Resolution Imaging Spectrometer (MODIS) data to explore its applicability; and (vii) errors in the SR retrievals are reported using the mean bias error (MBE), root mean squared deviation (RMSD), and mean systematic error (MSE). Results depict significant and strong positive Pearson's correlation (r), small MBE, RMSD, and MSE for each spectral band against in situ ASD data and Landsat (LEDAPS and LaSRC) SR products. Consistency in SREM performance against Sentinel-2A (r = 0.994, MBE = -0.009, and RMSD = 0.014) and MODIS (r = 0.925, MBE = 0.007, and RMSD = 0.014) data suggests that SREM can be applied to other multispectral satellites data. Overall, the findings demonstrate the potential and promise of SREM for use over diverse surfaces and under varying atmospheric conditions using multi-sensor data on a global scale.