Towards Routine Mapping of Shallow Bathymetry in Environments with Variable Turbidity: Contribution of Sentinel-2A/B Satellites Mission

Spatiotemporal variability of lake surface water temperature and water quality parameters and its interrelationship with water hyacinth biomass in Lake Tana, Ethiopia

Urbanization, agriculture, and climate change affect water quality and water hyacinth growth in lakes. This study examines the spatiotemporal variability of lake surface water temperature, turbidity, and chlorophyll-a (Chl-a) and their association with water hyacinth biomass in Lake Tana. MODIS Land/ Lake surface water temperature (LSWT), Sentinel 2 MSI Imagery, and in-situ water quality data were used. Validation results revealed strong positive correlations between MODIS LSWT and on-site measured water temperature (R = 0.90), in-situ turbidity and normalized difference turbidity index (NDTI) (R = 0.92), and in-situ Chl-a and normalized difference chlorophyll index (NDCI) (R = 0.84). LSWT trends varied across the lake, with increasing trends in the northeastern, northwestern, and southwestern regions and decreasing trends in the western, southern, and central areas (2001-2022). The spatial average LSWT trend decreased significantly in pre-rainy (0.01 ℃/year), rainy (0.02 ℃/year), and post-rainy seasons (0.01℃/year) but increased non-significantly in the dry season (0.00 ℃/year) (2001-2022, P < 0.05). Spatial average turbidity decreased significantly in all seasons, except in the pre-rainy season (2016-2022). Likewise, spatial average Chl-a decreased significantly in pre-rainy and rainy seasons, whereas it showed a non-significant increasing trend in the dry and post-rainy seasons (2016-2022). Water hyacinth biomass was positively correlated with LSWT (R = 0.18) but negatively with turbidity (R = -0.33) and Chl-a (R = -0.35). High spatiotemporal variability was observed in LSWT, turbidity, and Chl-a, along with overall decreasing trends. The findings suggest integrated management strategies to balance water hyacinth eradication and its role in water purification. The results will be vital in decision support systems and preparing strategic plans for sustainable water resource management, environmental protection, and pollution prevention.

Confronting turbidity, the major challenge for satellite-derived coastal bathymetry

Monitoring the complex seafloor morphology that drives the functioning of shallow coastal ecosystems is vital for assessing marine activities. Satellite-derived bathymetry (SDB) can provide a crucial dataset for creating the bathymetry maps needed to understand hazards and impacts produced by climate change in vulnerable coastal zones. SDB is effective in clear water, but still has limitations in application to areas with some turbidity. Here, using the twin satellites Sentinel-2A/B, we integrate water quality information from the satellite with a multi-temporal compositing method to demonstrate a potential for comprehensively operational bathymetric mapping over a range of environments. The automated compositing method diminishes the turbidity impact in addition to inferring the maximum detectable depth and removing optically deep-water areas. Examining a wide range of conditions along the Caribbean and eastern coast of the U.S. shows detailed bathymetry as deep as 30 m at 10 m spatial resolution with median errors <1 m when compared to high-resolution lidar surveys. These results demonstrate that the model adopted can provide useful bathymetry in areas that do not have consistently clear water and can be extended across multiple geographic regions and optical conditions at local, regional, and national scales.

Performance and Uncertainty of Satellite-Derived Bathymetry Empirical Approaches in an Energetic Coastal Environment

Objectives of this study are to evaluate the performance of different satellite-derived bathymetry (SDB) empirical models developed for multispectral satellite mission applications and to propose an uncertainty model based on inferential statistics. The study site is the Arcachon Bay inlet (France). A dataset composed of 450,837 echosounder data points and 89 Sentinel-2 A/B and Landsat-8 images acquired from 2013 to 2020, is generated to test and validate SDB and uncertainty models for various contrasting optical conditions. Results show that water column optical properties are characterized by a high spatio-temporal variability controlled by hydrodynamics and seasonal conditions. The best performance and highest robustness are found for the cluster-based approach using a green band log-linear regression model. A total of 80 satellite images can be exploited to calibrate SDB models, providing average values of root mean square error and maximum bathymetry of 0.53 m and 7.3 m, respectively. The uncertainty model, developed to extrapolate information beyond the calibration dataset, is based on a multi-scene approach. The sensitivity of the model to the optical variability not explained by the calibration dataset is demonstrated but represents a risk of error of less than 5%. Finally, the uncertainty model applied to a diachronic analysis definitively demonstrates the interest in SDB maps for a better understanding of morphodynamic evolutions of large-scale and complex coastal systems.

Assessing the Ability to Quantify Bathymetric Change over Time Using Solely Satellite-Based Measurements

Coastal regions are undergoing rapid change, due to natural and anthropogenic forcings. A current constraint in understanding and modeling these changes is the lack of multi-temporal bathymetric data, or recursive observations. Often, it is difficult to obtain the repeat observations needed to quantify bathymetric change over time or events. However, the recent availability of ICESat-2 bathymetric lidar creates the option to map coastal bathymetry from solely space-based measurements via satellite-derived bathymetry with multispectral imagery (IS-2/SDB). This compositional space-based bathymetric mapping technique can assess temporal change along the coasts without other remote sensing or in situ data. However, questions exist as to the accuracy of the technique relative to both quantitative uncertainties and the ability to resolve the spatial patterns of erosion and deposition in the nearshore environment, indicative of geomorphologic change. This paper addresses the concept using data from the Florida panhandle (Northern Gulf of Mexico) collected by Sentinel-2 and ICESat-2 at two epochs to assess the feasibility of using IS-2/SDB for bathymetric change detection at scientifically relevant scales, spatial resolutions and accuracies. The comparison of the satellite-only result is compared to airborne data collected at similar epochs to reveal both quantitatively and qualitatively the utility of this technique.

Satellite-derived bathymetry using Landsat-8 and Sentinel-2A images: assessment of atmospheric correction algorithms and depth derivation models in shallow waters

Satellite-derived bathymetry (SDB) has an extensive prospect in nearshore bathymetry for its high efficiency and low costs. Atmospheric correction and bathymetric modeling are critical processes in SDB, and examining the performance of related algorithms and models will contribute to the formulation of reliable bathymetry strategies. This study explored the effectiveness of three general atmospheric correction algorithms, namely Second Simulation of a Satellite Signal in the Solar Spectrum (6S), Atmospheric correction for OLI 'lite' (ACOLITE), and QUick Atmospheric Correction (QUAC), in depth retrieval from Landsat-8 and Sentinel-2A images using different SDB models over Ganquan Island and Oahu Island. The bathymetric Light Detection and Ranging (LiDAR) data was used for SDB model training and accuracy verification. The results indicated that the three atmospheric correction algorithms could provide effective corrections for SDB. For the SDB models except log-transformed band ratio model (LBR) and support vector machine (SVM), the impact of different atmospheric corrections on bathymetry was basically the same. Furthermore, we assessed the performance of six different SDB models: Lyzenga's model (LM), generalized additive model (GAM), LBR, SVM, multilayer perceptron (MLP), and random forest (RF). The bathymetric accuracy, consistency of bathymetric maps and generalization ability were considered for the assessment. Given sufficient training data, the accuracy of the machine learning models (SVM, MLP, RF) was generally superior to that of the empirical inversion models (LM, GAM, LBR), with the root mean square error (RMSE) varied between 0.735 m to 1.177 m. MLP achieved the best accuracy and consistency. When the depth was deeper than 15 m, the bathymetry error of all the SDB models increased sharply, and LM, LBR and SVM reached the upper limit of depth retrieval capability at 20-25 m. In addition, LM and LBR were demonstrated to have better adaptability in heterogeneous environment without training data.

Improved Filtering of ICESat-2 Lidar Data for Nearshore Bathymetry Estimation Using Sentinel-2 Imagery

The accurate estimation of nearshore bathymetry is necessary for multiple aspects of coastal research and practices. The traditional shipborne single-beam/multi-beam echo sounders and Airborne Lidar bathymetry (ALB) have a high cost, are inefficient, and have sparse coverage. The Satellite-derived bathymetry (SDB) method has been proven to be a promising tool in obtaining bathymetric data in shallow water. However, current empirical SDB methods for multispectral imagery data usually rely on in situ depths as control points, severely limiting their spatial application. This study proposed a satellite-derived bathymetry method without requiring a priori in situ data by merging active and passive remote sensing (SDB-AP). It realizes rapid bathymetric mapping with only satellite remotely sensed data, which greatly extends the spatial coverage and temporal scale. First, seafloor photons were detected from the ICESat-2 raw photons based on an improved adaptive Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm, which could calculate the optimal detection parameters for seafloor photons by adaptive iteration. Then, the bathymetry of the detected seafloor photons was corrected because of the refraction that occurs at the air–water interface. Afterward, the outlier photons were removed by an outlier-removal algorithm to improve the retrieval accuracy. Subsequently, the high spatial resolution (0.7 m) ICESat-2 derived bathymetry data were gridded to match the Sentinel-2 data with a lower spatial resolution (10 m). All of the ICESate-2 gridded data were randomly separated into two parts: 80% were employed to train the empirical bathymetric model, and the remaining 20% were used to quantify the inversion accuracy. Finally, after merging the ICESat-2 data and Sentinel-2 multispectral images, the bathymetric maps over St. Thomas of the United States Virgin Islands, Acklins Island in the Bahamas, and Huaguang Reef in the South China Sea were produced. The ICESat-2-derived results were compared against in situ data over the St. Thomas area. The results showed that the estimated bathymetry reached excellent inversion accuracy and the corresponding RMSE was 0.68 m. In addition, the RMSEs between the SDB-AP estimated depths and the ICESat-2 bathymetry results of St. Thomas, Acklins Island, and Huaguang Reef were 0.96 m, 0.91 m, and 0.94 m, respectively. Overall, the above results indicate that the SDB-AP method is effective and feasible for different shallow water regions. It has great potential for large-scale and long-term nearshore bathymetry in the future.

Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength

In this study, we present a new method called BathySent to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentinel 2 spectral bands, acquired quasi-simultaneously, from a single satellite dataset. Our method was based on the linear dispersion law, which related water depth to wave celerity and wavelength: when the water depth was less than about half the dominant wavelength, the wave celerity and wavelength decreased due to decreasing water depth (h) as the waves propagated towards the coast. Instead of using a best weighted (c,λ) fit with the linear dispersion relation to retrieve h, we proposed solving the linear dispersion relation for each (c,λ) pair to find multiple h-values within the same resolution cell. Then, we calculated the weighted averaged h-value for each resolution cell. To improve the precision of the final bathymetric map, we stacked the bathymetry values from N-different datasets acquired from the same study area on different dates. We first tested the algorithm on a set of images representing simulated ocean waves, then we applied it to a real set of Sentinel 2 data obtained of our study area, Gâvres peninsula (France, 47°,67 lat.; −3°35 lon.). A comparison with in situ bathymetry yielded good results from the synthetic images (r2 = 0.9) and promising results with the Sentinel 2 images (r2 = 0.7) in the 0–16 m depth zone.

ICESat-2 Elevation Retrievals in Support of Satellite-Derived Bathymetry for Global Science Applications

Bathymetry retrievals from 2D, multispectral imagery, referred to as Satellite-Derived Bathymetry (SDB), afford the potential to obtain global, nearshore bathymetric data in optically clear waters. However, accurate SDB depth retrievals are limited in the absence of "seed depths." The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) space-based altimeter has proven capable of accurate bathymetry, but methods of employing ICESat-2 bathymetry for SDB retrievals over broad spatial extents are immature. This research aims to establish and test a baseline methodology for generating bathymetric surface models using SDB with ICESat-2. The workflow is operationally efficient (17-37 min processing time) and capable of producing bathymetry of sufficient vertical accuracy for many coastal science applications, with RMSEs of 0.96 and 1.54 m when using Sentinel-2 and Landsat 8, respectively. The highest priorities for further automation have also been identified, supporting the long-range goal of global coral reef habitat change analysis using ICESat-2-aided SDB.

UAV-Derived Multispectral Bathymetry

Bathymetry is considered an important component in marine applications as several coastal erosion monitoring and engineering projects are carried out in this field. It is traditionally acquired via shipboard echo sounding, but nowadays, multispectral satellite imagery is also commonly applied using different remote sensing-based algorithms. Satellite-Derived Bathymetry (SDB) relates the surface reflectance of shallow coastal waters to the depth of the water column. The present study shows the results of the application of Stumpf and Lyzenga algorithms to derive the bathymetry for a small area using an Unmanned Aerial Vehicle (UAV), also known as a drone, equipped with a multispectral camera acquiring images in the same WorldView-2 satellite sensor spectral bands. A hydrographic Multibeam Echosounder survey was performed in the same period in order to validate the method’s results and accuracy. The study area was approximately 0.5 km2 and located in Tuscany (Italy). Because of the high percentage of water in the images, a new methodology was also implemented for producing a georeferenced orthophoto mosaic. UAV multispectral images were processed to retrieve bathymetric data for testing different band combinations and evaluating the accuracy as a function of the density and quantity of sea bottom control points. Our results indicate that UAV-Derived Bathymetry (UDB) permits an accuracy of about 20 cm to be obtained in bathymetric mapping in shallow waters, minimizing operative expenses and giving the possibility to program a coastal monitoring surveying activity. The full sea bottom coverage obtained using this methodology permits detailed Digital Elevation Models (DEMs) comparable to a Multibeam Echosounder survey, and can also be applied in very shallow waters, where the traditional hydrographic approach requires hard fieldwork and presents operational limits.

Intertidal Bathymetry Extraction with Multispectral Images: A Logistic Regression Approach

In this study, a methodology to estimate the intertidal bathymetry from multispectral remote sensing images is presented. The technique is based on the temporal variability of the water and the intertidal zone reflectance and their correlation with the tidal height. The water spectral behavior is characterized by high absorption at the infrared (IR) band or radiation with higher wavelengths. Due to tidal cycles, pixels on the intertidal zone have higher temporal variability on the near IR spectral reflectance. The variability of IR reflectivity in time is modeled through a sigmoid function of three parameters, where the inflection parameter corresponds to the pixel elevation. The methodology was tested at the Tagus river estuary in Lisbon, Portugal, and at the Bijagós archipelago, in the West African nation of Guinea-Bissau. Multispectral images from Sentinel-2 satellites were used, after atmospheric corrections from ACOLITE processor and the derived bathymetric model validated with in situ data. The presented method does not require additional depth data for calibration, and the output can generate intertidal digital elevation models at 10 m spatial resolution, without any manual editing by the operator. The results show a standard deviation of 0.34 m at the Tagus tidal zone, with −0.50 m bias, performing better than the Stumpf ratio transform algorithm, also applied to the test areas to derive intertidal bathymetry. This methodology can be used to update intertidal elevation models with clear benefits to monitoring of intertidal dynamics, morphodynamic modeling, and cartographic update.

Atmospheric correction for satellite-derived bathymetry in the Caribbean waters: from a single image to multi-temporal approaches using Sentinel-2A/B

Different atmospheric correction (AC) procedures for Sentinel-2 satellites are evaluated for their effectiveness in retrieving consistent satellite-derived bathymetry (SDB) over two islands in the Caribbean (Buck and Culebra). The log-ratio method for SDB, which allows use of minimal calibration information from lidar surveys (25 points in this study), is applied to several Sentinel-2A/B scenes at 10 m spatial resolution. The overall performance during a one-year study period depends on the image quality and AC. Three AC processors were evaluated: ACOLITE Exponential model (EXP), ACOLITE Dark Spectrum Fitting model (DSF), and C2RCC model. ACOLITE EXP and ACOLITE DSF produce greater consistency and repeatability with accurate results in a scene-by-scene analysis (mean errors ∼1.1 m) for depths up to 23 m (limit of lidar surveys). In contrast, C2RCC produces lower accuracy and noisier results with generally higher (>50%) errors (mean errors ∼2.2 m), but it is able to retrieve depth for scenes in Buck Island that have moderately severe sunglint. Furthermore, we demonstrate that a multi-temporal compositing model for SDB mapping, using ACOLITE for the input scenes, could achieve overall median errors <1 m for depths ranging 0-23 m. The simple and effective compositing model can considerably enhance coastal SDB estimates with high reliability and no missing data, outperforming the traditional single image approaches and thus eliminating the need to evaluate individual scenes. The consistency in the output from the AC correction indicates the potential for automated application of the multi-scene compositing technique, which can apply the open and free Sentinel-2 data set for the benefit of operational and scientific investigations.