A methodology for elaborating multi-temporal Sentinel-1 and Landsat 8 satellite images for estimating topsoil Soil Moisture Content (SMC) to support hydrological simulation studies is proposed. After pre-processing the remote sensing data, backscattering coefficient, Normalized Difference Vegetation Index (NDVI), thermal infrared temperature and incidence angle parameters are assessed for their potential to infer ground measurements of SMC, collected at the top 5 cm. A non-linear approach using Artificial Neural Networks (ANNs) is tested. The methodology is applied in Western Crete, Greece, where a SMC gauge network was deployed during 2015. The performance of the proposed algorithm is evaluated using leave-one-out cross validation and sensitivity analysis. ANNs prove to be the most efficient in SMC estimation yielding R² values between 0.7 and 0.9. The proposed methodology is used to support a hydrological simulation with the HEC-HMS model, applied at the Keramianos basin which is ungauged for SMC. Results and model sensitivity highlight the contribution of combining Sentinel-1 SAR and Landsat 8 images for improving SMC estimates and supporting hydrological studies.

Soil erosion is one of the main causes of soil degradation among others (salinization, compaction, reduction of organic matter, and non-point source pollution) and is a serious threat in the Mediterranean region. A number of soil properties, such as soil organic matter (SOM), soil structure, particle size, permeability, and Calcium Carbonate equivalent (CaCO₃), can be the key properties for the evaluation of soil erosion. In this work, several innovative methods (satellite remote sensing, field spectroscopy, soil chemical analysis, and GIS) were investigated for their potential in monitoring SOM, CaCO_3, and soil erodibility (K-factor) of the Akrotiri cape in Crete, Greece. Laboratory analysis and soil spectral reflectance in the VIS-NIR (using either Landsat 8, Sentinel-2, or field spectroscopy data) range combined with machine learning and geostatistics permitted the spatial mapping of SOM, CaCO₃, and K-factor. Synergistic use of geospatial modeling based on the aforementioned soil properties and the Revised Universal Soil Loss Equation (RUSLE) erosion assessment model enabled the estimation of soil loss risk. Finally, ordinary least square regression (OLSR) and geographical weighted regression (GWR) methodologies were employed in order to assess the potential contribution of different approaches in estimating soil erosion rates. The derived maps captured successfully the SOM, the CaCO₃, and the K-factor spatial distribution in the GIS environment. The results may contribute to the design of erosion best management measures and wise land use planning in the study region.