Soil moisture estimates at 1 km resolution making a synergistic use of Sentinel data
Authors: Madelon Remi , Rodríguez-Fernández Nemesio J. , Bazzi Hassan , Baghdadi Nicolas , Albergel Clement , Dorigo Wouter , Zribi Mehrez 
Published: March 2023

Abstract. Very high-resolution (∼10–100 m) surface soil moisture (SM) observations are important for applications in agriculture, among other purposes. This is the original goal of the S2MP (Sentinel-1/Sentinel-2-Derived Soil Moisture Product) algorithm, which was designed to retrieve surface SM at the agricultural plot scale by simultaneously using Sentinel-1 (S1) backscatter coefficients and Sentinel-2 (S2) NDVI (Normalized Difference Vegetation Index) as inputs to a neural network trained with Water Cloud Model simulations. However, for many applications, including hydrology and climate impact assessment at regional level, large maps with a high resolution (HR) of around 1 km are already a significant improvement with respect to most of the publicly available SM datasets, which have resolutions of about 25 km. In this study, the S2MP algorithm was adapted to work at 1 km resolution and extended from croplands to herbaceous vegetation types. A target resolution of 1 km also allows the evaluation of the interest in using NDVI derived from Sentinel-3 (S3) instead of S2. Two sets of SM maps at 1 km resolution were produced with S2MP over six regions of ∼104 km2 in Spain, Tunisia, North America, Australia, and the southwest and southeast regions of France for the whole year of 2019. The first set was derived from the combination of S1 and S2 data (S1 + S2 maps), while the second one was derived from the combination of S1 and S3 (S1 + S3 maps). S1 + S2 and S1 + S3 SM maps were compared to each other, to those of the 1 km resolution Copernicus Global Land Service (CGLS) SM and Soil Water Index (SWI) datasets, and to those of the Soil Moisture Active Passive (SMAP) + S1 product. The S2MP S1 + S2 and S1 + S3 SM maps are in very good agreement in terms of correlation (R≥0.9), bias (≤0.04 m3 m−3), and standard deviation of the difference (SDD≤0.03 m3 m−3) over the six domains investigated in this study. In a second step, the S1 + S3 S2MP maps were compared to the other HR maps. S1 + S3 SM maps are well correlated to the CGLS SM maps (R∼0.7–0.8), but the correlations with respect to the other HR maps (CGLS SWI and SMAP + S1) drop significantly over many areas of the six domains investigated in this study. The highest correlations between the HR maps were found over croplands and when the 1 km pixels have a very homogeneous land cover. The bias among the different maps was found to be significant over some areas of the six domains, reaching values of ±0.1 m3 m−3. The S1 + S3 maps show a lower SDD with respect to CGLS maps (≤0.06 m3 m−3) than with respect to the SMAP + S1 maps (≤0.1 m3 m−3) for all the six domains. Finally, all the HR datasets (S1 + S2, S1 + S3, CGLS, and SMAP + S1) were also compared to in situ measurements from five networks across five countries, along with coarse-resolution (CR) SM products from SMAP, SMOS, and the European Space Agency Climate Change Initiative (CCI). While all the CR and HR products show different bias and SDD, the HR products show lower correlations than the CR ones with respect to in situ measurements. The discrepancies in between the different HR datasets, except for the more simple land cover conditions (homogeneous pixels with croplands) and the lower performances with respect to in situ measurement than coarse-resolution datasets, show the remaining challenges for large-scale HR SM mapping.