Evaluation of Landsat-9 interoperability with Sentinel-2 and Landsat-8 over Europe and local comparison with field surveys

The recent launch of Landsat-9 satellite enriches the opportunities to work with dense time series of multispectral medium-resolution images. The integration of Landsat-9 in a multi-constellation series with Landsat-8 and Sentinel-2 requires a harmonization of the surface reflectance values that can be obtained from the official Level-2 products. This paper proposes the coefficients of the optimal linear transformations for the European continent, which allow to integrate Landsat-9 with the similar operating missions. These coefficients are based on a regression over 30 independent random extractions of 240,000 samples from images of the same areas but acquired by different sensors within two days. The coefficients were validated on an independent dataset. Furthermore, the effects of the proposed harmonization were tested on four popular vegetation indices, by evaluating the distributions of the differences in values obtained from each sensor pair. Finally, a test on a local scale was carried out with a spectroradiometer survey on 16 locations to collect some reference spectra to be compared with the reflectance values provided by the images. The results demonstrate the interoperability of Landsat and Sentinel-2 missions, since reflectance differences are in most cases within the accuracy specifications of the sensors. However, some discrepancies are observed in the blue and SWIR bands, probably due to inconsistencies in the atmospheric correction processes.