This study evaluates the accuracy of bathymetric maps generated from multispectral satellite datasets acquired from different multispectral sensors, namely the Worldview 2, PlanetScope, and the Sentinel 2, in the bay of Elounda in Crete. Image pre-processing steps were implemented before the use of the three empirical methods for estimating bathymetry. A dedicated correction and median filter have been applied to minimize noise from the sun glint and the sea waves. Due to the spectral complexity of the selected study area, statistical correlation with different numbers of bands was applied. The analysis indicated that blue and green bands obtained the best results with higher accuracy. Then, three empirical models, namely the Single Band Linear Algorithm, the Multiband Linear Algorithm, and the Ratio Transform Algorithm, were applied to the three multispectral images. Bathymetric and error distribution maps were created and used for the error assessment of results. The accuracy of the bathymetric maps estimated from different empirical models is compared with on-site Single beam Echo Sounder measurements. The most accurate bathymetric maps were obtained using the WorldView 2 and the empirical model of the Ratio Transform algorithm, with the RMSE reaching 1.01 m.

Shallow bathymetry inversion algorithms have long been applied in various types of remote sensing imagery with relative success. However, this approach requires that imagery with increased radiometric resolution in the visible spectrum be available. The recent developments in drones and camera sensors allow for testing current inversion techniques on new types of datasets with centimeter resolution. This study explores the bathymetric mapping capabilities of fused RGB and multispectral imagery as an alternative to costly hyperspectral sensors for drones. Combining drone-based RGB and multispectral imagery into a single cube dataset provides the necessary radiometric detail for shallow bathymetry inversion applications. This technique is based on commercial and open-source software and does not require the input of reference depth measurements in contrast to other approaches. The robustness of this method was tested on three different coastal sites with contrasting seafloor types with a maximum depth of six meters. The use of suitable end-member spectra, which are representative of the seafloor types of the study area, are important parameters in model tuning. The results of this study are promising, showing good correlation (R² > 0.75 and Lin’s coefficient > 0.80) and less than half a meter average error when they are compared with sonar depth measurements. Consequently, the integration of imagery from various drone-based sensors (visible range) assists in producing detailed bathymetry maps for small-scale shallow areas based on optical modelling.

Geomorphic indices can be used to examine the geomorphological and tectonic processes responsible for the development of the drainage basins. Such indices can be dependent on tectonics, erosional processes and other factors that control the morphology of the landforms. The inter-relationships between geomorphic indices can determine the influence of regional tectonic activity in the shape development of drainage basins. A Multi-Criteria Decision Analysis (MCDA) procedure has been used to perform an integrated cluster analysis that highlights information associated with the dominant regional tectonic activity. Factor Analysis (FA) and Analytical Hierarchy Process (AHP) were considered within that procedure, producing a representation of the distributed regional tectonic activity of the drainage basins studied. The study area is western Crete, located in the outer fore-arc of the Hellenic subduction zone, one of the world's most tectonically active regions. The results indicate that in the landscape evolution of the study area (especially the western basins) tectonic controls dominate over lithological controls.