Abstract. In recent decades, research has been developed to estimate near-shore bathymetry depth values using satellite imagery. Visible and infrared bands are used to derive elevation profile estimates, so to obtain bathymetric in rapid way without mobilisation of persons or equipment and saving the costs. For consequence, Satellite Derived Bathymetry (SDB) is seen as a valid approach for shallow waters survey: strongly supported by the activity of scholars and researchers, multiple methods are available in the literature. This article aims to investigate and compare different SDB methods for sea depth extraction from Sentinel-2 satellite multispectral images, with particular attention to the accuracy of the results. The experiments are conducted on imagery including Blue, Green, Red and Near Infrared bands, with 10 m resolution, concerning the Bay of Pozzuoli (Italy). After removing the glint, the effects caused by the reflection of sunlight through single scattering from sea surface, three methods are applied: Band Ratio method (BRM), 3rd-degree polynomial regression line method (3DPM), and principal component analysis method (PCAM). 3DPM can be seen as a variant of the BRM where the linear law that interprets the correlation between the band ratio values and the depth values is replaced by the third order function. Models are trained using depth data extracted from an Electronic Navigational Chart (ENC) at 1:7,500 scale, which is also used to verify result accuracy. The experiments demonstrate that the 3DPM is better able to obtain a more precise bathymetric model, confirming the greater adaptability of the 3rd order function to interpretate the variability of the interaction of light with water along the water column.