Landsat-7 ETM+, Landsat-8 OLI, and Sentinel-2 MSI surface reflectance cross-comparison and harmonization over the Mediterranean Basin area

In the Mediterranean area, vegetation dynamics and phenology analysed over a long time can have an important role in highlighting changes in land use and cover as well as the effect of climate change. Over the last 30 years, remote sensing has played an essential role in bringing about these changes thanks to many types of observations and techniques. Satellite images are to be considered an important tool to grasp these dynamics and evaluate them in an inexpensive and multidisciplinary way thanks to Landsat and Sentinel satellite constellations. The integration of these tools holds a dual potential: on the one hand, allowing us to obtain a longer historical series of reflectance data, while on the other hand making data available with a higher frequency even within a specific timeframe. The study aims to conduct a comprehensive cross-comparison analysis of long-time-series pixel values in the Mediterranean regions. For this scope comparisons between Landsat-7 (ETM+), Landsat-8 (OLI), and Sentinel-2 (MSI) satellite sensors were conducted based on surface reflectance products. We evaluated these differences using Ordinary Least Squares (OLS) and Major Axis linear regression (RMA) analysis on points extracted from over 15,000 images across the Mediterranean basin area from 2017 to 2020. Minor but consistent differences were noted, necessitating the formulation of suitable adjustment equations to better align Sentinel-2 reflectance values with those of Landsat-7 or Landsat-8. The results of the analysis are compared with the most-used harmonization coefficients proposed in the literature, revealing significant differences. The root-mean-square deviation, the mean difference and the orthogonal distance regression (ODR) slope show an improvement of the parameters for both models used (OLS and RMA) in this study. The discrepancies in reflectance values leads to corresponding variations in the estimation of biophysical parameters, such as NDVI, showing an increase in the ODR slope of 0.3. Despite differences in spatial, spectral, and temporal characteristics, we demonstrate that integration of these datasets is feasible through the application of band-wise regression corrections for a sensitive and heterogeneous area like those of the Mediterranean basin area.