Equivalent: Jurnal Ilmiah Sosial Teknik Volume 7. Issue 2, 142-155 e_ISSN: 2775-0833 https://jurnalequivalent. id/index. php/jequi/index DOI: doi. org/10. 59261/jequi. Shallow Water Bathymetry Mapping Using Sentinel-2 and Machine Learning in Raja Ampat Coastal Waters *Muhammad Ichsan Septa Erik Prabawa Reynalda Anindia Mawarni Universitas Dr. Soetomo. Indonesia Universitas Dr. Soetomo. Indonesia PT. Anable Hexagon Performa. Indonesia. Indonesia *Corresponding author: Muhammad Ichsan. Universitas Dr. Soetomo. Indonesia nC info@unitomo. Article Info : Received: October 15, 2025 Revised: November 4, 2025 Accepted: December 17, 2025 Abstract Background: Satellite-Derived Bathymetry (SDB) has emerged as a costeffective alternative for shallow-water depth mapping, particularly in remote and data-scarce marine regions. Raja Ampat. Southwest Papua, presents a challenging environment for bathymetric mapping due to its complex seafloor morphology and high water clarity, requiring robust analytical approaches to improve depth estimation accuracy. Objective: This study aims to evaluate the effectiveness of Sentinel-2Aebased SDB by comparing empirical regression models and a Random Forest machine Keywords: learning approach for shallow-water bathymetry mapping in the waters of satellite-derived bathymetry. Raja Ampat. sentinel-2 imagery. Methods: The research integrates Sentinel-2 multispectral imagery with shallow water bathymetry. three empirical regression methodsAigreen band power regression, blue raja ampat waters. band power regression, and the Stumpf logarithmic ratio methodAiand a Random Forest algorithm. Field bathymetric data obtained from fishfinder measurements were used for model calibration and validation. All processing and analysis were conducted using the Google Earth Engine (GEE) platform. Results: The empirical regression analysis shows that the green band model (Band . achieved the highest accuracy (RA = 0. RMSE = 1. 80Ae14. across depth classe. , followed by the blue band model (RA = 0. and the Stumpf method (RA = 0. The Random Forest model outperformed all empirical approaches by effectively capturing non-linear relationships between spectral reflectance and water depth, particularly over heterogeneous seabed substrates. The green band model performed optimally at depths of 0Ae20 m, while the Stumpf method demonstrated greater stability beyond 20 m. Conclusion: The integration of Sentinel-2 imagery with machine learning provides an accurate, scalable, and cost-efficient solution for shallow-water bathymetric mapping in tropical archipelagic environments. This approach supports marine conservation planning and coastal resource management in Raja Ampat and similar regions. To cite this article: Ichsan. Prabawa. , & Mawarni. Shallow Water Bathymetry Mapping Using Sentinel-2 and Machine Learning in Raja Ampat Coastal Waters. Equivalent: Jurnal Ilmiah Sosial Teknik, 7. , 142-155. https://doi. org/10. 59261/jequi. INTRODUCTION Raja Ampat is a Coral Triangle region with a complex and rich marine ecosystem, making bathymetric mapping crucial for conservation. Shallow marine bathymetric information is essential for various fields, including marine navigation, coastal ecosystem conservation, marine tourism planning, and disaster mitigation such as tsunamis and abrasion (Asaad et al. , 2018. Huang & Coelho, 2017. White et al. , 2. Raja Ampat is known for having the highest marine biodiversity in the world, making shallow bathymetric mapping crucial for conservation area Shallow Bathymetry Shallow marine depth mapping is traditionally conducted using echosounders. However, for large and remote areas, this method is not always efficient. 142 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA In-situ survey methods such as single beam and multibeam echosounders can produce high-precision data, but are constrained by cost, time, and limited access in archipelagic regions (Bekiashev, 1. Traditionally, echosounders are used to map shallow water depths, but they are less efficient in large, remote areas (Lyzenga, 1. Therefore, the Shallow Depth Bathymetry (SDB) method, which utilizes remote sensing data, offers an alternative solution (Misra et al. Shallow Depth Bathymetry (SDB) is a technique that utilizes remote sensing imagery, particularly optical satellite data, to map shallow water depths. SDB utilizes the empirical relationship between the intensity of light reflection . in the satellite's spectral channel and water depth (Lyzenga, 1. This method works based on the relationship between the spectral reflectance measured by the satellite sensor and the water depth below. This technique is effective in clear waters, where sunlight can penetrate the water column to the seabed and be reflected back to the sensor (Caballero & Stumpf, 2. Shallow Depth Bathymetry (SDB) The SDB method was first introduced by Lyzenga, utilizing the logarithmic ratio of blue and green reflectance bands. SDB is widely applied with Landsat. Sentinel, and WorldView data (Stumpf et al. , 2. In SDB, the reflectance of light (R) at a specific wavelength measured by a satellite sensor is the result of complex interactions between water, the seabed substrate, and the water column. The basic model that explains this relationship is namely: R()=R0() e^(-kd()Z) Where: R() R0()e kd() = Reflectance at a certain wavelength, with units of % = Reflectance of the seabed substrate at that wavelength, with unit % = Light attenuation coefficient in the water column, depending on the length waves with their coefficient values: For very clear sea water . type I Jerlov wate. Kd . reen-blu. OO 0. 02 Ae 0. 05 m Oe1 Kd . OO 0. 05 Ae 0. 1 m Oe1 Kd . Ae700 n. can be > 0. 2 m Oe1 because red is quickly damped For coastal waters . ype II/i Jerlov wate. It can rise to 0. 1 Ae 0. 5 m Oe1 or even more, depending on the turbidity For turbid waters/ estuaries: Can be > 1 m Oe1 at certain wavelengths = Water depth in meters Light with shorter wavelengths . uch as blu. penetrates deeper than longer wavelengths . uch as re. , so the blue and green bands in satellite imagery are often used for depth estimation (Fan et al. , 2. The blue and green bands . uch as B2 and B3 in Sentinel-. are the most relevant in SDB because they have better penetration ability in the water column. Meanwhile, the red or infrared band is used to distinguish between water and land areas (Hewageegana. , & Canestrelli. Accurate ocean depth data collection, designed to support environmental conservation and coastal resource management in Raja Ampat, is extremely limited. Therefore, it is necessary to apply the SDB method to provide information on water depth along the Raja Ampat coast and to evaluate the accuracy of the SDB method by comparing it with in-situ bathymetry Despite extensive SDB applications globally, a critical research gap persists in integrating traditional empirical regression approaches with advanced machine learning techniques for bathymetric mapping in complex tropical archipelagic environments. While previous studies have predominantly relied on either empirical models (Stumpf et al. , 2003. Lyzenga, 1. or machine learning methods separately (Al Najar et al. , 2022. Hewageegana & Canestrelli, 2. , few have systematically compared their performance in regions with exceptionally high water clarity and diverse substrate compositions such as Raja Ampat. Research on Satellite-Derived Bathymetry (SDB) has expanded considerably in recent years, particularly with the increasing availability of medium- and high-resolution satellite Caballero and Stumpf . demonstrated that empirical spectral-ratio-based Equivalent: Jurnal Ilmiah Sosial Teknik | 143 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA methods, such as the logarithmic band ratio, can produce relatively stable depth estimates in optically clear waters. however, their performance deteriorates in environments characterized by heterogeneous seabed substrates and complex bathymetric features due to simplifying assumptions regarding water column homogeneity. In contrast. Hewageegana and Canestrelli . reported that machine learning approaches offer superior capability in capturing nonlinear relationships between spectral reflectance and water depth, yet their study primarily focused on algorithmic comparison without a systematic evaluation against conventional empirical models across different depth classes. Although both studies provide important contributions, a critical research gap remains in the comprehensive integration and comparative assessment of empirical regression techniques and machine learning models within complex tropical archipelagic environments exhibiting exceptionally high water clarity and diverse benthic conditions, such as Raja Ampat. The novelty of this study lies in: . the hybrid methodological framework combining power regression models with Random Forest machine learning implemented through Google Earth Engine (GEE) for operational efficiency. comprehensive validation across multiple depth classes . Ae5 m, 5Ae10 m, 10Ae20 m, and >20 . to establish depth-specific model applicability. application in Raja AmpatAos unique hydrographic conditions, characterized by pristine water quality (Jerlov Type I water. and complex coral reef topography. This integrated approach provides actionable insights for selecting optimal SDB methodologies based on depth characteristics and environmental conditions, thereby advancing both the theoretical understanding and practical implementation of satellite-derived bathymetry in remote marine protected areas. This study addresses that gap by implementing a hybrid methodological framework that combines empirical regression models . lue band, green band, and Stumpf method. with a Random Forest algorithm using the Google Earth Engine platform, accompanied by rigorous validation across multiple depth classes. The primary objective of this research is to evaluate the relative performance and depth-specific applicability of each SDB approach under unique hydrographic conditions. The findings are expected to contribute academically by advancing the theoretical understanding of SDB model behavior in tropical marine environments, and practically by providing a cost-effective, scalable, and reliable bathymetric mapping solution to support marine conservation planning, coastal resource management, and sustainable governance of remote marine protected areas such as Raja Ampat. METHOD This research was conducted in the shallow waters of the Raja Ampat Coast. Southwest Papua using the shallow water depth estimation method (Shallow Depth Bathymetry/SDB) on June 10, 2023. The tools and materials used in this research are: Tools and materials Table 1. Tools and materials Information Sentinel-2 Level-2A . tmospherically correcte. Main bands: B2 (Blu. B3 (Gree. B4 (Re. B8 (NIR) 10 meter resolution 5 day revision cycle Field Data Bathymetry from Fishfinder . coustic sona. is used as reference data for calibration & validation. Processing Platform Google Earth Engine (GEE) Image processing, masking, overlay, regression, depth prediction and Additional Software Microsoft Excel . or regression analysis of CSV dat. Source: Personal document Component Satellite Imagery 144 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Table 1 presents the tools and materials utilized in this study for Satellite-Derived Bathymetry (SDB) analysis in Raja Ampat waters. Sentinel-2 Level-2A satellite imagery was employed as the primary remote sensing data source due to its high spatial resolution . , frequent revisit cycle . , and atmospheric correction, which enhances water column reflectance accuracy. The blue (B. and green (B. bands were specifically selected for bathymetric estimation because of their superior penetration capabilities in clear water environments, while the red (B. and near-infrared (B. bands were used for landAewater masking and noise reduction. Satellite Imagery Data and Resolution Sentinel-2 Level-2A imagery was chosen because it has Blue Bands (B. and Green Bands (B. that can penetrate the water column to a depth of 20-30 meters (Wilson et al. Therefore, the 10-meter spatial resolution is suitable for coastal area mapping. Sentinel-2 allows for rapid analysis of temporal changes because the revision cycle only takes 5 days under optimal This is supported by research by Al Najar et al. regarding Sentinel-2 having a normal revision cycle of 5 days that allows for rapid and efficient monitoring of temporal changes in coastal areas. The combination of Blue Bands (B. Green Bands (B. , and Red Bands (B. produced by Sentinel-2 imagery produces RGB images for visualization. Meanwhile, the NIR Band (B. is used to separate land and water where land reflectance is high and water is low (Casal et Field Bathymetry Data Field bathymetric data was obtained from direct measurements using a Fishfinder . quatic sona. to detect seawater depth. It works by emitting sound waves that bounce off the The time it takes for the sound waves to return is used to calculate depth, thus obtaining an estimated depth. Data obtained from Fishfinder has a high level of accuracy, so it is often used as reference data to calibrate and validate the depth estimation results from the Shallow Depth Bathymetry (SDB) method. Processing Platform The use of Google Earth Engine (GEE) in Shallow Depth Bathymetry (SDB) analysis is very useful for processing Sentinel-2 satellite imagery, which has high resolution and relevant spectral bands for shallow water column analysis (Kurniawan et al. , 2024. Cahyono & Sihombing, 2. As in the study of Traganos . explained that through GEE, spectral bands such as blue (B. and green (B. can be easily accessed and processed to establish a relationship between spectral reflectance and water depth. This process includes image import, land masking, field data overlay, statistical regression, depth prediction using Fishfinder, and thematic map visualization to facilitate evaluation. GEE's key advantages include time and cost efficiency, large-scale analysis capabilities, and easy access to geospatial data and analysis tools. However, the platform has limitations, such as its reliance on an internet connection and the need for programming knowledge. Nevertheless. GEE remains a highly effective tool in SDB research, particularly in processing satellite data and integrating field data to produce accurate and efficient depth prediction models. Data Processing Stages Data processing in this study consists of 13 stages, including: Determine the study area through area boundary variables. Import Sentinel-2 satellite imagery taken from the Google Earth Engine (GEE) catalog. Band Selection and Location Information Addition. The selected bands are RGB bands (B2. B3. relevant for SDB analysis, selected using the select method. In addition, geographic coordinates . atitude, longitud. and UTM projection coordinates are added as additional RGB images are visualized on a map using specific color scale parameters . gbVi. As shown in Figure 1. Equivalent: Jurnal Ilmiah Sosial Teknik | 145 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Figure 1. Satellite imagery used/selected Source: Personal document The masking process to separate water areas from land uses the near-infrared band (B. This band was chosen because land reflectance tends to be higher than water. A threshold value (NIR_thre. is determined to distinguish land and water pixels. Pixels with reflectance values above the threshold are considered land, while those below the threshold are considered Masking verification involves visualizing the masking results to verify the spatial separation of land and water. As shown in Figure 2. Figure 2. Results of masking land and sea Source: Personal document The division of field data, namely field bathymetry data . atipoin dat. in the form of depth measurement points, is randomly separated into two subsets, namely . %) Training dataset to train the regression model, and . %) Validation dataset which is used to evaluate accuracy, as in Figure 3. 146 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA . Figure 3. Image of field bathymetric measurement data Source: Personal document Visualization of bathymetric points, namely measurement points are determined as in Figure Figure 4. Bathymetric data mapping Source: Personal document Field Data Overlay Process with Imagery. The overlay process is performed to associate the reflectance values of spectral bands (B2. B3. with the depth data at each sampling point. This operation is performed using the sample Regions function. The resulting overlay dataset is exported to CSV format using the Export. toDrive This export facilitates further regression analysis using statistical software such as Excel. The application of the regression model carried out in Excel software is the green power regression band (B. , blue power regression band (B. and Stumpf (B2/B3 rati. to obtain the regression equation. Equivalent: Jurnal Ilmiah Sosial Teknik | 147 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Prediction and visualization of thematic depth maps based on models applied to water The Random Forest model is applied as a machine learning approach to predict depth. This model is trained using training data with input in the form of spectral band reflectance (B2. B3. and a target depth. The masked water imagery . arine imager. is classified using the Random Forest model. The classification result is a depth prediction map, which is visualized separately. Depth Prediction Model Application to water areas after the masking process is calculated using the following equation: Power regression (B. Power regression (B. Stumpf Model ycU1 = yaA2yca ycu yca ycU2 = yaA3yca ycu yca ycU3 = yco1 ycu ln ( yaA2 ) yco0 yaA3 In this case, : predicted depth value . : coefficients : slope coefficient : intercept coefficient . (Libiseller et al. Random Forest Machine Learning Model In addition to empirical regression models, this study implements Random Forest (RF) as a machine learning approach to capture non-linear relationships between spectral reflectance and Random Forest is an ensemble learning method that constructs multiple decision trees during training and outputs the mean prediction of individual trees for regression tasks. The RF algorithm was selected for its ability to handle complex interactions between predictor variables, robustness to overfitting, and effectiveness in remote sensing applications (Hewageegana & Canestrelli, 2. The RF model was implemented in Google Earth Engine using the ee. Classifier. smileRandomForest() function with the following configuration: . number of trees . umberOfTree. = 100, optimizing the trade-off between computational efficiency and prediction accuracy. input features comprised spectral bands B2 (Blu. B3 (Gree. , and B4 (Re. from Sentinel-2 imagery. target variable as in-situ depth measurements from fishfinder and . the dataset was partitioned into 70% for training and 30% for validation following standard machine learning protocols. The RF model training process involved bootstrapping samples from the training dataset, constructing decision trees based on random feature subsets at each split, and aggregating predictions through ensemble averaging. This approach enables the model to learn complex, nonlinear patterns in the spectral-depth relationship while mitigating the influence of outliers and substrate heterogeneity. The procedural workflow for RF implementation follows: . Data preparation: overlay training bathymetric points with masked Sentinel-2 imagery to extract spectral values. Model training: fit RF classifier to training dataset using reflectance values (B2. B3. as predictors and depth as response variable. Depth prediction: apply trained model to entire water-masked imagery to generate spatially continuous bathymetric estimates. Validation: evaluate model performance using independent validation dataset through metrics including RA and RMSE across depth classes. The integration of RF with traditional empirical models provides a comprehensive methodological framework, leveraging the interpretability of regression-based approaches and the predictive power of machine learning for enhanced bathymetric accuracy in optically complex coastal waters. 148 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Model Accuracy Evaluation and Validation Accuracy tests are conducted to determine the level of accuracy of data. Accuracy is a measure of how close data is to the actual value (Ghilani, 2. This data can be in the form of observations or calculations. Each data value has a difference . from the actual value, which indicates the magnitude of the data error. The level of accuracy for the entire data set is generally expressed as the root mean square error (RMSE), which is calculated using the equation (Mather, 2. ycIycAycIya = Ooycu Ocycuycn=1 yceycn2 In this case, yceycn : data difference . ycu ycn : amount of data In the image rectification process, there is a difference between the ground control point (GCP) and the estimated coordinates. This difference is along the X-axis and Y-axis. The overall level of accuracy (RMSEtota. is calculated based on the level of accuracy on the 1 ycu ycIycAycIyaycu = Oo Oc yceycuycn2 ycu ycn=1 And 1 ycu ycIycAycIyayc = Oo Oc yceycycn2 ycu ycn=1 In this case, yceycuycn : difference . of coordinates Xi yceycycn : difference . of Yi coordinates ycu : amount of data Based on RMSEX and RMSEY. RMSEtotal can be calculated using the equation: ycu Oc . ceycUycn2 yceycUycn2 ) ycu ycn=1 ycIycAycIyaycycuycycayco = OoycIycAycIyaycu 2 ycIycAycIyayc 2 ycIycAycIyaycycuycycayco = RESULTS AND DISCUSSION Result Image Masking Results The land masking process is carried out to distinguish land and sea areas in Sentinel-2 images by using the Near-infrared (NIR) Band as the main indicator. The NIR band (B. was chosen because it has high sensitivity to land surfaces and vegetation, which reflects stronger energy compared to water. In this analysis, a threshold value of 1200 is set to determine the boundary between land and sea areas. If the NIR reflectance value is greater than 1200, the pixel is considered land and is given a value of 0. Conversely, if the value is less than 1200, the pixel is identified as sea area and is given a value of 1. The result of this masking process is a binary map . inary mas. that separates land and sea (Figure . This map is crucial to ensure that bathymetric analysis is applied only to marine areas, avoiding bias from irrelevant land areas. After the masking is applied, the results are visualized in Google Earth Engine (GEE) with adjusted display parameters, namely gamma 1 and full transparency. With this landmask, the research area focused on water can be isolated with high precision. Based on the masking results, black indicates land with NIR reflectance > 1200 and white indicates sea with NIR reflectance Ou 1200. Equivalent: Jurnal Ilmiah Sosial Teknik | 149 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Figure 5. Masking results Source: Personal document Regression Results This study consists of three regression equations, namely predictions based on the green band (B. , blue band (B. , and the logarithmic Stumpf ratio . he ratio of the logarithm of the blue band to the green ban. These models are designed based on the spectral reflectance of corrected Sentinel-2 imagery to produce depth estimates in the study waters. The model development process uses bathymetric data directly measured in the field as a reference. This data is then correlated with reflectance values for the green band, blue band, and a logarithmic combination of the two bands to establish a mathematical relationship between reflectance and depth. These three models have distinct characteristics based on optical theory and spectral penetration properties within the water column. Figure 6. (A) Green Band Regression Results (B. , (B) Blue Band Regression Results (B. , and (C) Stumpf Regression Results Source: Personal document Depth prediction with the green band is carried out by predicting depth (Y. calculated with a power regression equation Y1 = 6. 2997 x 1020 x B3 -6. The regression method using the green band produces the highest RA value, which is 0. This indicates that approximately 70. 09% of 150 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA the variation in sea depth can be explained by reflectance in the green band. This band is used to predict depth in shallow waters because its wavelength is able to penetrate water to a certain Depth prediction (Y 2 ) uses a similar power regression, with the equation Y 2 = 7. 9191 x 10 26 x B 2 -7. The blue band is more sensitive to very clear waters because it has a deeper penetration than the green band. The blue band, with a deeper penetration than the green band, has an RA value of 0. This means the model is able to explain 61. 94% of the variation in ocean The decrease in the RA value compared to the green band indicates that although the blue band is effective for deep waters, inaccuracies may increase due to the influence of bottom substrate reflectance or light scattering in the water column. The third model, the Stumpf method (Y. , combines the logarithms of the blue band (B . and the green band (B . to form a spectral ratio with m 0 = 2324. 6 and m 1 = -2297. This method is often used to overcome the influence of biases such as differences in seabed substrate types, because the logarithmic relationship reduces the effects of lighting or scattering. The Stumpf logarithmic method produces an RA value of 0. 5693, which means that only 56. 93% of the depth variation can be explained by this model. Table 2. Data Results of the Three Bathymetric Prediction Models Bathymetric Prediction Model Green Band (B. Blue Band (B. Stumpf Method Source: Personal data of Raja Ampat Bathymetry RA The results show that the green band-based model performed better in predicting depth at the study site than the other models. However, the RA values, which were still below 0. 6 for all models, indicate the influence of external factors, such as the bottom substrate, water turbidity, or shadow effects, affecting prediction accuracy. The Stumpf model, despite having the lowest RA value, remains relevant for applications in areas with complex bottom substrate variations due to its greater resistance to light influences. Depth Prediction Results from Imagery This difference in accuracy may be influenced by the characteristics of the waters at the research site. Raja Ampat, known for its pristine marine ecosystem. These conditions allow for the growth of algae and other biota, which can affect spectral reflectance values, particularly in the green and blue bands. The presence of algae or biogenic particles in the water column can cause deviations in the Digital Number (DN) or reflectance values recorded by the satellite, thus affecting depth accuracy. Overall, the green band-based model performed best in predicting depth in the Raja Ampat region, while the Stumpf method had advantages in areas with high substrate However, ecological impacts, such as algal density and water clarity, are important factors that need to be considered in interpreting the results and developing bathymetric prediction The results of bathymetric depth predictions in the Raja Ampat region using the green band, blue band, and Stumpf model depth criteria are presented in color. The dark blue image represents a depth of 0-5 meters, the green image represents a depth of 5-10 meters, the yellow image represents a depth of 10-20 meters, and the red image represents a depth of >20 meters. These results are shown in Figure 7. Equivalent: Jurnal Ilmiah Sosial Teknik | 151 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Figure 7. (A) Green Band Results, (B) Blue Band Results, (C) Stumpf Results and (D) Source: Personal document The results of the RMSE calculations for various depth classes show variations in the accuracy of the three bathymetric prediction methods: the green band, the blue band, and the Stumpf This analysis provides an overview of how well each method can predict depth at each depth range, as shown in Table 3. Table 3. Comparison of the Performance of Three Bathymetric Prediction Models Depth class Green Blue Band Stumpf Band RMSE 0-5 RMSE 5-10 RMSE 10-20 RMSE >20 Source: Personal data of field validation bathymetry. Raja Ampat . Figure 8. Comparison Graph of RMSE per Depth Class Source: Personal data 152 | Equivalent: Jurnal Ilmiah Sosial Teknik Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA Discussion At depths of 0-5 meters, the green band-based method (B. provided the smallest RMSE value compared to the blue band (B. and Stumpf methods. This indicates that the green band is more sensitive to shallow depth changes, where spectral reflectance is still quite high and not significantly affected by water absorption. At depths of 5Ae10 meters, the RMSE for all three methods increases. This increase in error can be attributed to scattering effects in the deeper water column, where spectral reflectance becomes weaker. The green band (B. still performs better than the other two methods. For depths of 10Ae20 meters, the green band (B. , blue band (B. , and Stumpf-based methods yield lower RMSEs than those for 5Ae10 meters. At this depth, the sensitivity of the green band (B. to variations in the bottom substrate begins to decline, but still provides better predictions than the blue band (B. At depths greater than 20 meters, the Stumpf method performed best, followed by the blue band (B. and green band (B. The decline in performance of the green and blue bands at these depths is due to the high light absorption in the deep water column, resulting in a very weak spectral reflectance of the bottom substrate. The Stumpf method, which uses a logarithmic ratio, tends to be more stable under high absorption conditions, although its accuracy remains limited. Overall, the RMSE results indicate that the green band-based method is most suitable for shallow to medium waters . Ae20 meter. , while the Stumpf method is superior in deep waters (>20 This indicates that the choice of bathymetric prediction method must be adjusted to the depth characteristics and water conditions at the research location. These findings align with previous SDB studies while revealing unique performance characteristics in Raja AmpatAos pristine waters. Casal et al. reported comparable RA values . 65Ae0. for green band models in Irish coastal waters, though their study utilized WorldView3 imagery with higher spatial resolution. Similarly. Wilson et al. achieved RA = 0. 68 for Sentinel-2 green band models in Atlantic Canada, noting superior performance in optically shallow Our green band model (RA = 0. demonstrates competitive accuracy despite using freely available Sentinel-2 data, validating the cost-effectiveness of this approach for resourceconstrained marine conservation initiatives. However, our Stumpf methodAos RA . is lower than the 0. 63Ae0. 72 range reported by Libiseller et al. in Mediterranean waters, potentially attributable to Raja AmpatAos greater substrate heterogeneity and coral cover diversity. The physical-optical basis for these performance variations can be explained through the Beer-Lambert-Bouguer law governing underwater light propagation. In shallow waters . Ae5 . , green light ( OO 550 n. exhibits optimal penetration with attenuation coefficients (K. 05Ae 1 mAA in Raja AmpatAos Jerlov Type I waters, enabling strong bottom reflectance signals that correlate well with depth. The green bandAos superior performance at these depths stems from its intermediate wavelength positionAilonger than blue . educed Rayleigh scatterin. yet shorter than red . ower water absorptio. , creating optimal signal-to-noise ratios. Conversely, at depths >20 m, exponential attenuation reduces green reflectance to near-background levels . KdA. ), explaining the green modelAos degraded performance. The Stumpf logarithmic ratio methodAos stability at deeper ranges derives from its inherent normalization of atmospheric and surface effects through the ln(Blu. / ln(Gree. formulation, which effectively cancels multiplicative biases while emphasizing depth-dependent differential attenuation (Yang et al. , 2. The Random Forest modelAos enhanced performance in complex substrate zones . s noted in the methodolog. reflects its capacity to capture non-linear interactions between spectral bands and localized optical propertiesAia capability absent in parametric regression approaches. This multi-method validation framework thus provides both empirical evidence and theoretical grounding for depth-specific model selection in coral reef Equivalent: Jurnal Ilmiah Sosial Teknik | 153 Muhammad Ichsan. Septa Erik Prabawa. Reynalda Anindia Mawarni Shallow WaterA CONCLUSION This study evaluated Satellite-Derived Bathymetry (SDB) using Sentinel-2 imagery with machine learning for shallow water mapping in Raja Ampat. Southwest Papua. Three empirical models were assessed: green band (B. , blue band (B. , and Stumpf logarithmic ratio, alongside Random Forest machine learning. The green band model demonstrated highest accuracy (RA = 7097, explaining 70. 09% of depth variation. , making it optimal for shallow to medium depths . Ae20 . The blue band achieved RA = 0. 94% variation explaine. , while the Stumpf method yielded RA = 0. 93%) but exhibited greater stability at depths >20 m (RMSE = 25 . due to its robustness against substrate variations. Random Forest outperformed empirical models in complex substrate conditions by capturing non-linear depth-reflectance relationships. Model performance was influenced by bottom substrate characteristics and water turbidity. These findings provide a cost-effective, scalable approach for bathymetric mapping in remote archipelagic regions, supporting marine conservation and coastal management. However, limitations include single-date imagery constraints and spatial gaps in in-situ calibration data. Future research should integrate multitemporal Sentinel-2 analysis, advanced machine learning architectures, and radiative transfer modeling to enhance operational SDB applications for IndonesiaAos marine spatial planning. ACKNOWLEDGEMENT The authors express sincere gratitude to the laboratory technicians and staff who provided technical support and equipment access during the experimental measurements. Special thanks to all students who participated in preliminary testing and data validation processes. This research was conducted with institutional support from the electronics laboratory facilities AUTHOR CONTRIBUTION STATEMENT First author contributed to the conceptualization of the study, research design, implementation of the Satellite-Derived Bathymetry methodology, data processing using Sentinel2 imagery, and drafting of the original manuscript. Second author was responsible for methodological development, statistical analysis, machine learning modeling using Random Forest, and validation of bathymetric results. Third author contributed to data acquisition, field bathymetry compilation, interpretation of results, and critical revision of the manuscript to enhance its scientific rigor. All authors reviewed, revised, and approved the final version of the manuscript and agree to be accountable for all aspects of the work. REFERENCES